This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: xyouts,0.025,20,'!mt!x* = 10',charsize=1.8,color=3 xyouts,0.025,2,'!mt!x* = 1',charsize=1.8,color=11 x ;calculate the altitude at the weighting function maximum print,' ' print,' Tau Altitude of Max' for i=0,2 do begin x=wherE(w(i,*) eq max(w(i,*))) print,tau(i),' ',z(x(0)) endfor print,' ' p device,/close set_plot,'win' s end e ;The code will print: ; Tau Altitude of Max ; 1 0.000000 ; 10 18.0000 ; 100 37.0000 ; ;The altitude at maximum W, zmax, which is obtained by using dW/dz = 0, ;is given by: tau* exp(-zmax/H) = 1 (for m=1) (see Petty eq. 7.60). ;Thus: zmax = H ln(tau*) ;0 km (tau*=1), 18.4 km (tau*=10), 36.8 km (tau*=100) ;The IDL values were constrained to integer altitudes since z ;was defined in 1-km increment...
View Full Document
- Fall '10
- Height, Tau Altitude, /nodata,subtitle='Transmission Weighting Functions, weighting function maximum, different optical depths