; IDL Version 4.0.1 (Win32 x86) ; Journal File for user@Microsoft Windows Host ; Working directory: C:\RSI\IDL40 ; Date: Thu Feb 08 15:02:50 1996 t=dindgen(31)*0.1D0 help,t z=dindgen(61)*0.01d0 q=0.03*t dn=0.04d0 d0=0.1d0 l1=0.60d0 d=(dn*(l1-z)+d0*z)/l1 & plot,z,d a=!pi*d*d/4d0 plot,z,a v=q/a plot,t,v help,q,a wdelete,0 hlp,a ; 1 Double array (61). Min = 0.0012566371, Max = 0.0078539819 ainv=1.d0/a & help,ainv ainv=1.d0/a & hlp,ainv ; 1 Double array (61). Min = 127.32395, Max = 795.77469 q=ainv#transpose(q) help,q contour,q hlp,q ; 1 Double array (61, 31). Min = 0.00000000, Max = 71.619721 q=0.03*t hlp,q ; 1 Double array (31). Min = 0.00000000, Max = 0.089999998 v=ainv#transpose(q) contour,v,lev=[0,20,40,60] contour,v,lev=[0,10,20,30,40,50,60,70] contour,v,z,t,lev=[0,10,20,30,40,50,60,70] contour,v,z,t,lev=[0,10,20,30,40,50,60,70],/follow v=-ainv#transpose(q) contour,v,z,t,lev=(-1)*[0,10,20,30,40,50,60,70],/follow ; % CONTOUR: Contour levels must be in increasing order. contour,v,z,t,lev=(-1)*reverse([0,10,20,30,40,50,60,70]),/follow dvdz=(8.d0/!dpi)*transpose(q)*l1*l1*(d0-dn)##(1.d0/(d*d*d)) hlp,dvdz ; 1 Double array (61, 31). Min = 0.00000000, Max = 77.349301 contour,dvdz,z,t,lev=[0,10,20,30,40,50,60,70],/follow contour,dvdz,z,t,lev=[0,2,4,6,8,10,20,30,40,50,60,70],/follow convacc=v*dvdz & hlp,convacc ; 1 Double array (61, 31). Min = -5539.7353, Max = 0.00000000 contour,convacc,z,t,/follow,lev=[100*findgen(10),1000*findgen(6)] ; % CONTOUR: Contour levels must be in increasing order. contour,convacc,z,t,/follow,lev=[100*findgen(10),1000*(findgen(5)+1)] contour,convacc,z,t,/follow,lev=(-1)*reverse([100*findgen(10),1000*(findgen(5)+1)]) print,convacc(30,20) ; -150.00863 print,v(30,20) ; -15.590688 print,dvdz(30,20) ; 9.6216817 wdelete,0 print,(8.d0/!dpi)*transpose(q)*l1*l1*(d0-dn) ; 0.00000000 ; 0.00016501184 ; 0.00033002368 ; 0.00049503552 ; 0.00066004737 ; 0.00082505921 ; 0.00099007105 ; 0.0011550829 ; 0.0013200947 ; 0.0014851066 ; 0.0016501184 ; 0.0018151303 ; 0.0019801421 ; 0.0021451539 ; 0.0023101658 ; 0.0024751776 ; 0.0026401895 ; 0.0028052013 ; 0.0029702131 ; 0.0031352250 ; 0.0033002368 ; 0.0034652487 ; 0.0036302605 ; 0.0037952724 ; 0.0039602842 ; 0.0041252960 ; 0.0042903079 ; 0.0044553197 ; 0.0046203316 ; 0.0047853434 ; 0.0049503552 print,transpose(q) ; 0.00000000 ; 0.0029999999 ; 0.0059999999 ; 0.0089999998 ; 0.012000000 ; 0.015000000 ; 0.018000000 ; 0.021000000 ; 0.023999999 ; 0.026999999 ; 0.029999999 ; 0.032999999 ; 0.035999999 ; 0.038999999 ; 0.041999999 ; 0.044999999 ; 0.047999999 ; 0.050999999 ; 0.053999999 ; 0.056999999 ; 0.059999999 ; 0.062999999 ; 0.065999999 ; 0.068999998 ; 0.071999998 ; 0.074999998 ; 0.077999998 ; 0.080999998 ; 0.083999998 ; 0.086999998 ; 0.089999998 print,d0-dn ; 0.060000000 print,l1 ; 0.60000000 print,((8.d0/!dpi)*transpose(q)*l1*l1*(d0-dn))(20) ; 0.0033002368 print,d(30) ; 0.070000000 dvdz=((8.d0/!dpi)*transpose(q)*(d0-dn)/l1)##(1.d0/(d*d*d)) print,dvdz(30,20) ; 44.544823 convacc=v*dvdz & print,convacc(30,20) ; -694.48441 dvdt=(-4.d0/!dpi)*(transpose(q/t))##(1.d0/(d*d)) ; % Program caused arithmetic error: Floating illegal operand dvdt=(-4.d0/!dpi)*(replicate(0.03d0,1,30))##(1.d0/(d*d)) print,dvdt(30,20) ; -7.7953442 locacc=dvdt totacc=convacc+locacc & print,totacc(30,20) ; -702.27976 rho=1600.d0 g=9.81d0 g=-9.81d0 dpdz=-rho*(totacc+g) hlp,totacc ; 1 Double array (61, 30). Min = -23989.498, Max = -3.8197186 surface,alog10(totacc),z,t ; % SURFACE: X, Y, or Z array dimensions are incompatible ; % Program caused arithmetic error: Floating illegal operand dvdt=(-4.d0/!dpi)*(replicate(0.03d0,1,31))##(1.d0/(d*d)) locacc=dvdt totacc=convacc+locacc & print,totacc(30,20) ; -702.27976 dpdz=-rho*(totacc+g) hlp,dpdz,totacc ; 1 Double array (61, 31). Min = 21807.550, Max = 41088970. ; 2 Double array (61, 31). Min = -25670.796, Max = -3.8197186 surface,alog10(dpdz),z,t dz=shift(z,-1)-z print,z ; 0.00000000 0.010000000 0.020000000 0.030000000 0.040000000 0.050000000 ; 0.060000000 0.070000000 0.080000000 0.090000000 0.10000000 0.11000000 ; 0.12000000 0.13000000 0.14000000 0.15000000 0.16000000 0.17000000 ; 0.18000000 0.19000000 0.20000000 0.21000000 0.22000000 0.23000000 ; 0.24000000 0.25000000 0.26000000 0.27000000 0.28000000 0.29000000 ; 0.30000000 0.31000000 0.32000000 0.33000000 0.34000000 0.35000000 ; 0.36000000 0.37000000 0.38000000 0.39000000 0.40000000 0.41000000 ; 0.42000000 0.43000000 0.44000000 0.45000000 0.46000000 0.47000000 ; 0.48000000 0.49000000 0.50000000 0.51000000 0.52000000 0.53000000 ; 0.54000000 0.55000000 0.56000000 0.57000000 0.58000000 0.59000000 ; 0.60000000 print,dz ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; -0.60000000 dz=(shift(z,-1)-z)>0 print,dz ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 0.010000000 ; 0.00000000 p=dpdz-dpdz & hlp,p ; 1 Double array (61, 31). Min = 0.00000000, Max = 0.00000000 for i=1,60 do p(i,0)=int_tabulated(z(0:i,0),dpdz(0:i,0)) plot,z(*,0),p(*,0) wdelete,0 plot,z(*,0),p(*,0) hlp,z(*,0),p(*,0) ; 1 Double array (61). Min = 0.00000000, Max = 0.60000000 ; 2 Double array (61). Min = 0.00000000, Max = 18584.925 print,g*rho*l1 ; -9417.6000 plot,findgen(10),findgen(10)^2 erase plot,findgen(10),findgen(10)^2 wdelete,0 plot,findgen(10),findgen(10)^2 save,/all,file='\ep324\prob5_15.dat'