Index: /trunk/mesoscale/PLOT/SPEC/LES/compturb.pro =================================================================== --- /trunk/mesoscale/PLOT/SPEC/LES/compturb.pro (revision 112) +++ /trunk/mesoscale/PLOT/SPEC/LES/compturb.pro (revision 113) @@ -10,5 +10,5 @@ ;smoothampl=4. smoothampl=0. -;smoothampl=5. +smoothampl=5. ;;;;;;; ; @@ -41,219 +41,219 @@ ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;; pour article, decommenter ici et commenter l'appel au deuxieme fichier +;;; pour article, decommenter ici et commenter l'appel au deuxieme fichier +;;;com +;;!p.multi = [0,2,1] +; +; path = paths(0) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;; +;;; COEFFICIENT POUR LOIS DE SIMILITUDE +;;; +;;@coeffhf +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = a_wt_bot +;zey = a_h +;set_name = saveplot+'comp_A_HFLUX.ps' +;set_title = "Similarity: vert. eddy heat flux on Mars" +;set_titlex = "Dimensionless vertical heat flux /!Dmax!N" +;set_titley = 'Dimensionless height z/z!Di!N' +;set_subtitle = '' +;set_xrange = [-0.5,1.] +;set_yrange = [0.,1.5] +;set_tickx = 0.25 +;set_ticky = 0.1 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zesim = 1 - 1.2 * a_h +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] +;;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )] +;;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )] +;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin +;for ll = 1, n_elements(localtimes)-1 do begin +; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin +;endfor +;;wherex = set_xrange(0) +;;wherey = set_yrange(1) +;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zecolor = 0. & loadct, 33 +;for i=1, n_elements(paths)-1 do begin +;; zecolor = 50. + zecolor +;zecolor = 0. & loadct, 0 +; path = paths(i) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; zefield = a_wt_bot +; zey = a_h +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; for ll = 0, n_elements(localtimes)-1 do begin +; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor +; endfor +; ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;endfor +;loadct, 0 +;;oplot, [0.,0.], set_yrange, linestyle=0, color=0 +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP +; ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_End, /PNG ;;com -;!p.multi = [0,2,1] - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;; -;; COEFFICIENT POUR LOIS DE SIMILITUDE -;; -;@coeffhf - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = a_wt_bot -zey = a_h -set_name = saveplot+'comp_A_HFLUX.ps' -set_title = "Similarity: vert. eddy heat flux on Mars" -set_titlex = "Dimensionless vertical heat flux /!Dmax!N" -set_titley = 'Dimensionless height z/z!Di!N' -set_subtitle = '' -set_xrange = [-0.5,1.] -set_yrange = [0.,1.5] -set_tickx = 0.25 -set_ticky = 0.1 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zesim = 1 - 1.2 * a_h -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] -;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )] -;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )] -;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin -for ll = 1, n_elements(localtimes)-1 do begin - user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) - if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin -endfor -;wherex = set_xrange(0) -;wherey = set_yrange(1) -;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin -; zecolor = 50. + zecolor -zecolor = 0. & loadct, 0 - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = a_wt_bot - zey = a_h - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - for ll = 0, n_elements(localtimes)-1 do begin - user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) - if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor - endfor - ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -endfor -loadct, 0 -;oplot, [0.,0.], set_yrange, linestyle=0, color=0 - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP - ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG -;com - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;; -;; COEFFICIENT POUR LOIS DE SIMILITUDE -;; -;@coeffw - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = a_vel_bot -zey = a_h -set_name = saveplot+'comp_A_VEL.ps' -set_title = "Similarity: vert. vel. variance on Mars" -set_titlex = "Dimensionless vert. vel. variance /w!D*!N!U2!N" -set_titley = 'Dimensionless height z/z!Di!N' -set_subtitle = '' -set_xrange = [0.,1.2] -set_yrange = [0.,1.5] -set_tickx = 0.25 ;0.1 -set_ticky = 0.1 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] -;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;com -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin -for ll = 1, n_elements(localtimes)-1 do begin - user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) - if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin -endfor -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin -zecolor = 0. & loadct, 0 -; zecolor = 50. + zecolor - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = a_vel_bot - zey = a_h - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - for ll = 0, n_elements(localtimes)-1 do begin - user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) - if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor - endfor - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -endfor -loadct, 0 - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP - ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG - +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; path = paths(0) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;; +;;; COEFFICIENT POUR LOIS DE SIMILITUDE +;;; +;;@coeffw +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = a_vel_bot +;zey = a_h +;set_name = saveplot+'comp_A_VEL.ps' +;set_title = "Similarity: vert. vel. variance on Mars" +;set_titlex = "Dimensionless vert. vel. variance /w!D*!N!U2!N" +;set_titley = 'Dimensionless height z/z!Di!N' +;set_subtitle = '' +;set_xrange = [0.,1.2] +;set_yrange = [0.,1.5] +;set_tickx = 0.25 ;0.1 +;set_ticky = 0.1 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] +;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;com -;stop - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE) -set_name = saveplot+'comp_W_STAR.ps' -set_title = '';"Free convection velocity scale W!D*!N (m.s!U-1!N)" -set_titlex = 'Local Time (h)' -set_titley = 'Free convection velocity scale (m s!U-1!N)' -set_subtitle = '' -set_xrange = [08.,17.] -set_yrange = [00.,06.] -set_tickx = 1.0 -set_ticky = 0.5 -wherex = 8.5;11.0 -wherey = 5.5;1.75 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -altlin=0 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path - plots, wherex + 2.*dddx/1.1, wherey - plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin +;for ll = 1, n_elements(localtimes)-1 do begin +; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin +;endfor +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin - ;zecolor = float(255)/n_elements(paths) + zecolor - zecolor = 0 - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = localtime - zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - oplot, zefield, zey, color=zecolor, linestyle=altlin+i - xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor - plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy - plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -endfor -loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;for i=1, n_elements(paths)-1 do begin +;zecolor = 0. & loadct, 0 +;; zecolor = 50. + zecolor +; path = paths(i) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; zefield = a_vel_bot +; zey = a_h +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; for ll = 0, n_elements(localtimes)-1 do begin +; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor +; endfor +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;endfor +;loadct, 0 +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP +; ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_End, /PNG +; +;;;com +;;stop +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; path = paths(0) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = localtime +;zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE) +;set_name = saveplot+'comp_W_STAR.ps' +;set_title = '';"Free convection velocity scale W!D*!N (m.s!U-1!N)" +;set_titlex = 'Local Time (h)' +;set_titley = 'Free convection velocity scale (m s!U-1!N)' +;set_subtitle = '' +;set_xrange = [08.,17.] +;set_yrange = [00.,06.] +;set_tickx = 1.0 +;set_ticky = 0.5 +;wherex = 8.5;11.0 +;wherey = 5.5;1.75 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;altlin=0 +;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path +; plots, wherex + 2.*dddx/1.1, wherey +; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zecolor = 0. & loadct, 33 +;for i=1, n_elements(paths)-1 do begin +; ;zecolor = float(255)/n_elements(paths) + zecolor +; zecolor = 0 +; path = paths(i) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; zefield = localtime +; zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor +; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;endfor +;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_End, /PNG +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; path = paths(0) & print, path @@ -271,10 +271,10 @@ set_titley = 'Boundary layer height (km)' set_subtitle = '' -;set_xrange = [09.,18.] -set_yrange = [00.,07.] +set_xrange = [08.,18.] +set_yrange = [00.,08.] ;[00.,07.] set_tickx = 1. set_ticky = 1. wherex = 09.;13.25 ;12.75 -wherey = 06.;2.5 ;2.0 +wherey = 06.5;2.5 ;2.0 ; set_xrange = [13.5,18.] & set_yrange = [06.4,08.0] & set_tickx = 0.5 & set_ticky = 0.2 ; wherex=15.30 @@ -321,4 +321,61 @@ if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;; path = paths(0) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;;set_name = saveplot+'comp_dHEIGHT.ps' +;;set_title = "Boundary layer height growth dz!Di!N/dt (km/hour)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Boundary layer height growth (km/hour)' +;;set_subtitle = '' +;;set_xrange = [9.,18.] +;;set_yrange = [-0.4,1.6] +;;set_tickx = 1. +;;set_ticky = 0.2 +;;wherex = 10.30 +;;wherey = 0.1 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;altlin = 0 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path +;; plots, wherex + 2.*dddx/1.1, wherey +;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;zecolor = 0. & loadct, 33 +;;for i=1, n_elements(paths)-1 do begin +;; ;zecolor = float(255)/n_elements(paths) + zecolor +;; zecolor = 0 +;; path = paths(i) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; zefield = localtime +;; zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor +;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;endfor +;;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_End, /PNG +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -333,16 +390,17 @@ ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;zefield = localtime -;zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) -;set_name = saveplot+'comp_dHEIGHT.ps' -;set_title = "Boundary layer height growth dz!Di!N/dt (km/hour)" +;zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) +;zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) +;set_name = saveplot+'comp_WMAX.ps' +;set_title = '';"Maximum and minimum vertical velocity (m.s!U-1!N)" ;set_titlex = 'Local Time (h)' -;set_titley = 'Boundary layer height growth (km/hour)' +;set_titley = 'Maximum and minimum vertical velocity (m s!U-1!N)' ;set_subtitle = '' -;set_xrange = [9.,18.] -;set_yrange = [-0.4,1.6] +;;set_xrange = [09.0,18.0] +;set_yrange = [-12.0,18.0];[-15.0,20.0] ;set_tickx = 1. -;set_ticky = 0.2 -;wherex = 10.30 -;wherey = 0.1 +;set_ticky = 2. +;wherex = 8.5;11.00 +;wherey = 16.5;5.5 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_Start, FILENAME=set_name @@ -350,13 +408,14 @@ ;altlin = 0 ;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path +;oplot, zefield, zey2, color=0, linestyle=altlin +;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) ; plots, wherex + 2.*dddx/1.1, wherey ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;zecolor = 0. & loadct, 33 ;for i=1, n_elements(paths)-1 do begin ; ;zecolor = float(255)/n_elements(paths) + zecolor -; zecolor = 0 +; zecolor = 0. ; path = paths(i) & print, path ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -365,8 +424,10 @@ ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; zefield = localtime -; zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +; zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) +; zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i -; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor +; oplot, zefield, zey2, color=zecolor, linestyle=altlin+i +; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue @@ -377,184 +438,5 @@ ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_End, /PNG - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) -zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) -set_name = saveplot+'comp_WMAX.ps' -set_title = '';"Maximum and minimum vertical velocity (m.s!U-1!N)" -set_titlex = 'Local Time (h)' -set_titley = 'Maximum and minimum vertical velocity (m s!U-1!N)' -set_subtitle = '' -;set_xrange = [09.0,18.0] -set_yrange = [-12.0,18.0];[-15.0,20.0] -set_tickx = 1. -set_ticky = 2. -wherex = 8.5;11.00 -wherey = 16.5;5.5 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -altlin = 0 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -oplot, zefield, zey2, color=0, linestyle=altlin -dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) - plots, wherex + 2.*dddx/1.1, wherey - plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin - ;zecolor = float(255)/n_elements(paths) + zecolor - zecolor = 0. - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = localtime - zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) - zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - oplot, zefield, zey, color=zecolor, linestyle=altlin+i - oplot, zefield, zey2, color=zecolor, linestyle=altlin+i - xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor - plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy - plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -endfor -loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax -set_name = saveplot+'comp_TKE.ps' -set_title = '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[++]!Dmax!N (m!U2!N.s!U-2!N)" -set_titlex = 'Local Time (h)' -set_titley = 'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)' -set_subtitle = '' -;set_xrange = [09.0,18.0] -set_yrange = [00.0,18.0];[00.0,26.0] -set_tickx = 1. -set_ticky = 2. -wherex = 8.5;11.70 -wherey = 16.;5.4 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -altlin = 0 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0) - plots, wherex + 2.*dddx/1.1, wherey - plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin - ;zecolor = float(255)/n_elements(paths) + zecolor - zecolor = 0. - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = localtime - zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - oplot, zefield, zey, color=zecolor, linestyle=altlin+i - xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitles(i), color=zecolor - plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy - plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -endfor -loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax -set_name = saveplot+'comp_HFBOT.ps' -set_title = '';"Maximum heat flux !Dmax!N (K.m.s!U-1!N)" -set_titlex = 'Local Time (h)' -set_titley = 'Maximum vertical eddy heat flux (K m s!U-1!N)' -set_subtitle = '' -;set_xrange = [09.0,18.0] -set_yrange = [00.0,02.0];[00.0,02.8] -set_tickx = 1. -set_ticky = 0.2 -wherex = 8.5;10.30 -wherey = 1.8;0.7 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -altlin = 0 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) - plots, wherex + 2.*dddx/1.1, wherey - plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin - ;zecolor = float(255)/n_elements(paths) + zecolor - zecolor = 0. - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/getturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = localtime - zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - oplot, zefield, zey, color=zecolor, linestyle=altlin+i - xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor - plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy - plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -endfor -loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG - - +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -569,20 +451,77 @@ ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;zefield = localtime -;zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top -;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) -;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) -;set_name = saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps' -;set_title = "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)" +;zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax +;set_name = saveplot+'comp_TKE.ps' +;set_title = '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[++]!Dmax!N (m!U2!N.s!U-2!N)" ;set_titlex = 'Local Time (h)' -;set_titley = 'Potential temperature time derivative (K.hour!U-1!N)' +;set_titley = 'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)' ;set_subtitle = '' -;set_xrange = [ 09., 18.] -;set_yrange = [-01., 06.] +;;set_xrange = [09.0,18.0] +;set_yrange = [00.0,18.0];[00.0,26.0] ;set_tickx = 1. -;set_ticky = 1. -;wherex = 13.25 -;wherey = 5.5 +;set_ticky = 2. +;wherex = 8.5;11.70 +;wherey = 16.;5.4 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;altlin = 0 +;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0) +; plots, wherex + 2.*dddx/1.1, wherey +; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zecolor = 0. & loadct, 33 +;for i=1, n_elements(paths)-1 do begin +; ;zecolor = float(255)/n_elements(paths) + zecolor +; zecolor = 0. +; path = paths(i) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; zefield = localtime +; zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitles(i), color=zecolor +; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;endfor +;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_End, /PNG +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; path = paths(0) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = localtime +;zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax +;set_name = saveplot+'comp_HFBOT.ps' +;set_title = '';"Maximum heat flux !Dmax!N (K.m.s!U-1!N)" +;set_titlex = 'Local Time (h)' +;set_titley = 'Maximum vertical eddy heat flux (K m s!U-1!N)' +;set_subtitle = '' +;;set_xrange = [09.0,18.0] +;set_yrange = [00.0,02.0];[00.0,02.8] +;set_tickx = 1. +;set_ticky = 0.2 +;wherex = 8.5;10.30 +;wherey = 1.8;0.7 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name ;!P.Charsize = 1.2 ;altlin = 0 @@ -591,5 +530,190 @@ ; plots, wherex + 2.*dddx/1.1, wherey ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zecolor = 0. & loadct, 33 +;for i=1, n_elements(paths)-1 do begin +; ;zecolor = float(255)/n_elements(paths) + zecolor +; zecolor = 0. +; path = paths(i) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/getturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; zefield = localtime +; zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor +; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;endfor +;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_End, /PNG +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;; path = paths(0) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top +;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;;set_name = saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps' +;;set_title = "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Potential temperature time derivative (K.hour!U-1!N)' +;;set_subtitle = '' +;;set_xrange = [ 09., 18.] +;;set_yrange = [-01., 06.] +;;set_tickx = 1. +;;set_ticky = 1. +;;wherex = 13.25 +;;wherey = 5.5 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;altlin = 0 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) +;; plots, wherex + 2.*dddx/1.1, wherey +;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;zecolor = 0. & loadct, 33 +;;for i=1, n_elements(paths)-1 do begin +;; ;zecolor = float(255)/n_elements(paths) + zecolor +;; zecolor = 0. +;; path = paths(i) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; zefield = localtime +;; zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top +;; ;zey = DERIV(localtime, t_bot) +;; zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;; ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor +;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;endfor +;;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;oplot, set_xrange, [0,0], linestyle=2 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_End, /PNG +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;; path = paths(0) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top +;;set_name = saveplot+'comp_TTOP.ps' +;;set_title = "Temperature at the top of the boundary layer (K)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Temperature at the top of the boundary layer (K)' +;;set_subtitle = '' +;;set_xrange = [09.,18.] +;;set_yrange = [195.,225.] +;;set_tickx = 1. +;;set_ticky = 5. +;;wherex = 13.00 +;;wherey = 223.5 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;altlin = 0 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path +;; plots, wherex + 2.*dddx/1.1, wherey +;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;zecolor = 0. & loadct, 33 +;;for i=1, n_elements(paths)-1 do begin +;; ;zecolor = float(255)/n_elements(paths) + zecolor +;; zecolor = 0. +;; path = paths(i) & print, path +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; file = path+'/compturb.dat' & restore, filename=file +;; file = path+'/getturb.dat' & restore, filename=file +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; zefield = localtime +;; zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i +;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor +;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy +;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue +;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;endfor +;;loadct, 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_End, /PNG +; +; +;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;; +;lev = 4 +;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;; +; +; path = paths(0) & print, path +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; file = path+'/compturb.dat' & restore, filename=file +; file = path+'/addturb.dat' & restore, filename=file +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = localtime +;zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) +;set_name = saveplot+'comp_VELMAX.ps' +;set_title = '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)" +;set_titlex = 'Local Time (h)' +;set_titley = "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)" +;set_subtitle = '' +;;set_xrange = [09.0,18.0] +;set_yrange = [00.0,14.0] +;set_tickx = 1. +;set_ticky = 2. +;wherex = 8.5;10.30 +;wherey = 12.;0.7 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;print, h(lev) +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;altlin = 0 +;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) +; plots, wherex + 2.*dddx/1.1, wherey +; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;zecolor = 0. & loadct, 33 @@ -600,11 +724,8 @@ ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; file = path+'/compturb.dat' & restore, filename=file -; file = path+'/getturb.dat' & restore, filename=file +; file = path+'/addturb.dat' & restore, filename=file ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; zefield = localtime -; zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top -; ;zey = DERIV(localtime, t_bot) -; zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) -; ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) +; zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i @@ -613,131 +734,10 @@ ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; ;endfor ;loadct, 0 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;oplot, set_xrange, [0,0], linestyle=2 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_End, /PNG - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -; path = paths(0) & print, path -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; file = path+'/compturb.dat' & restore, filename=file -; file = path+'/getturb.dat' & restore, filename=file -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = localtime -;zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top -;set_name = saveplot+'comp_TTOP.ps' -;set_title = "Temperature at the top of the boundary layer (K)" -;set_titlex = 'Local Time (h)' -;set_titley = 'Temperature at the top of the boundary layer (K)' -;set_subtitle = '' -;set_xrange = [09.,18.] -;set_yrange = [195.,225.] -;set_tickx = 1. -;set_ticky = 5. -;wherex = 13.00 -;wherey = 223.5 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;altlin = 0 -;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path -; plots, wherex + 2.*dddx/1.1, wherey -; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;zecolor = 0. & loadct, 33 -;for i=1, n_elements(paths)-1 do begin -; ;zecolor = float(255)/n_elements(paths) + zecolor -; zecolor = 0. -; path = paths(i) & print, path -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; file = path+'/compturb.dat' & restore, filename=file -; file = path+'/getturb.dat' & restore, filename=file -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; zefield = localtime -; zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; oplot, zefield, zey, color=zecolor, linestyle=altlin+i -; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor -; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy -; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue -; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;endfor -;loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_End, /PNG - - -;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;; -lev = 4 -;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;; - - path = paths(0) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/addturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) -set_name = saveplot+'comp_VELMAX.ps' -set_title = '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)" -set_titlex = 'Local Time (h)' -set_titley = "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)" -set_subtitle = '' -;set_xrange = [09.0,18.0] -set_yrange = [00.0,14.0] -set_tickx = 1. -set_ticky = 2. -wherex = 8.5;10.30 -wherey = 12.;0.7 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -print, h(lev) -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -altlin = 0 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) - plots, wherex + 2.*dddx/1.1, wherey - plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zecolor = 0. & loadct, 33 -for i=1, n_elements(paths)-1 do begin - ;zecolor = float(255)/n_elements(paths) + zecolor - zecolor = 0. - path = paths(i) & print, path - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - file = path+'/compturb.dat' & restore, filename=file - file = path+'/addturb.dat' & restore, filename=file - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - zefield = localtime - zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - oplot, zefield, zey, color=zecolor, linestyle=altlin+i - xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor - plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy - plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -endfor -loadct, 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_End, /PNG Index: /trunk/mesoscale/PLOT/SPEC/LES/getturb.pro =================================================================== --- /trunk/mesoscale/PLOT/SPEC/LES/getturb.pro (revision 112) +++ /trunk/mesoscale/PLOT/SPEC/LES/getturb.pro (revision 113) @@ -17,7 +17,7 @@ ;; constantes ;; -;p0=610. & t0=220. & r_cp=1/4.4 & grav=3.72 & R=192. -p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191. -print, 'ATTENTION ATTENTION R/cp !!!!', r_cp +p0=610. & t0=220. & r_cp=1.0/4.4 & grav=3.72 & R=192. +;p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191. +;print, 'ATTENTION ATTENTION R/cp !!!!', r_cp ; @@ -171,5 +171,5 @@ OPENR, 23, 'input_more' & READF, 23, hgtu, tsurfu & CLOSE, 23 -@include_ustar.pro +;@include_ustar.pro ;@include_hfx.pro @@ -344,67 +344,47 @@ ;*******************; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = a_vel_bot -;zey = a_h -;set_name = 'A_VEL.ps' -;set_title = "Scaled vertical velocity variance (m.s!U-1!N)" -;set_titlex = 'Scaled vertical velocity variance (m.s!U-1!N)' -;set_titley = 'Scaled height (km)' -;set_subtitle = '' -;set_xrange = [0.,1.2] -;set_yrange = [0.,1.2] -;set_tickx = 0.1 -;set_ticky = 0.1 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] -;;localtimes = localtime -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;altlin=0 -;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin, /ISOTROPIC -;for ll = 1, n_elements(localtimes)-1 do begin -; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin -;endfor -;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0 -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = localtime -zey = pbl_height -set_name = 'HEIGHT.ps' -set_title = 'Boundary layer height (km)' -set_titlex = 'Local Time (h)' -set_titley = 'Boundary layer height (km)' -set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' -set_xrange = [8.,17.] -set_yrange = [0.,7.] -set_tickx = 1. -set_ticky = 1. -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;oplot, zefield, zey, psym=5 -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = a_vel_bot +;;zey = a_h +;;set_name = 'A_VEL.ps' +;;set_title = "Scaled vertical velocity variance (m.s!U-1!N)" +;;set_titlex = 'Scaled vertical velocity variance (m.s!U-1!N)' +;;set_titley = 'Scaled height (km)' +;;set_subtitle = '' +;;set_xrange = [0.,1.2] +;;set_yrange = [0.,1.2] +;;set_tickx = 0.1 +;;set_ticky = 0.1 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] +;;;localtimes = localtime +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;altlin=0 +;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin, /ISOTROPIC +;;for ll = 1, n_elements(localtimes)-1 do begin +;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin +;;endfor +;;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0 +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;zefield = localtime -;zey = w_star_bot -;set_name = 'W_STAR.ps' -;set_title = "Free convection velocity scale (m.s!U-1!N)" +;zey = pbl_height +;set_name = 'HEIGHT.ps' +;set_title = 'Boundary layer height (km)' ;set_titlex = 'Local Time (h)' -;set_titley = 'Free convection velocity scale (m.s!U-1!N)' -;set_subtitle = '' -;set_xrange = [8.,18.] -;set_yrange = [0.,6.] +;set_titley = 'Boundary layer height (km)' +;set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' +;set_xrange = [11.,17.] ;[8.,17.] +;set_yrange = [0.,8.] ;set_tickx = 1. -;set_ticky = 0.5 +;set_ticky = 1. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_Start, FILENAME=set_name @@ -414,119 +394,103 @@ ;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = localtime -;zey = hfbot -;zey2 = hfsurf -;zey3 = hfbotmax -;set_name = 'HFBOT.ps' -;set_title = "Heat flux at the top of the radiative layer (K.m.s!U-1!N)" -;set_titlex = 'Local Time (h)' -;set_titley = 'Heat flux at the top of the radiative layer (K.m.s!U-1!N)' -;set_subtitle = '' -;set_xrange = [8.,18.] -;set_yrange = [0.,2.5] -;set_tickx = 1. -;set_ticky = 0.5 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;oplot, zefield, zey2, linestyle=1 -;oplot, zefield, zey3, linestyle=2 -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = transpose(tke) -zex = localtime -zey = h -set_name = 'TKE.ps' -set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[++] -set_titlex = 'Local Time (h)' -set_titley = 'Altitude above surface (km)' -set_subtitle = '' ;'Mean over the simulation domain' -set_xrange = [8.,17.] -set_yrange = [0.,7.] ;; [0.,200.] -set_tickx = 1. -set_ticky = 1. ;; 50. -minval = 0. -maxval = 15. -nlev = maxval-minval -pal = 22 -rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -;; 0. levels -lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;; 1. background -loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 - ;; 2. color field - loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) - contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;; 4. choose output -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = 2.*transpose(ztke) -;set_title = "Vertical wind variance (m!U2!N.s!U-2!N)" -;minval = -6. -;maxval = 12. -;pal = 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = transpose(ztke) -zex = localtime -zey = h -set_name = 'zTKE.ps' -set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] -set_titlex = 'Local Time (h)' -set_titley = 'Altitude above surface (km)' -set_subtitle = '' ;'Mean over the simulation domain' -set_xrange = [8.,17.] -set_yrange = [0.,7.] -set_tickx = 1. -set_ticky = 1. -minval = 0. -maxval = 8. -nlev = maxval-minval -pal = 22 -rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -;; 0. levels -lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;; 1. background -loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;; 2. color field -loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) - contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;; 4. choose output -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;restore, filename='addturb.dat' -;modvar = SMOOTH(TEMPORARY(modvar), [0,smoothampl], /EDGE_TRUNCATE) -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = transpose(modvar) -;zex = localtime +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = w_star_bot +;;set_name = 'W_STAR.ps' +;;set_title = "Free convection velocity scale (m.s!U-1!N)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Free convection velocity scale (m.s!U-1!N)' +;;set_subtitle = '' +;;set_xrange = [8.,18.] +;;set_yrange = [0.,6.] +;;set_tickx = 1. +;;set_ticky = 0.5 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;;oplot, zefield, zey, psym=5 +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = hfbot +;;zey2 = hfsurf +;;zey3 = hfbotmax +;;set_name = 'HFBOT.ps' +;;set_title = "Heat flux at the top of the radiative layer (K.m.s!U-1!N)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Heat flux at the top of the radiative layer (K.m.s!U-1!N)' +;;set_subtitle = '' +;;set_xrange = [8.,18.] +;;set_yrange = [0.,2.5] +;;set_tickx = 1. +;;set_ticky = 0.5 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;oplot, zefield, zey2, linestyle=1 +;;oplot, zefield, zey3, linestyle=2 +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = transpose(tke) +;zex = localtime ;zey = h -;set_name = 'modvar.ps' -;set_title = "Horizontal wind speed variance (m!U2!N.s!U-2!N)" +;set_name = 'TKE.ps' +;set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[++] ;set_titlex = 'Local Time (h)' ;set_titley = 'Altitude above surface (km)' ;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [7.,17.] +;set_xrange = [11.,17.] ;[8.,17.] +;set_yrange = [0.,8.] ;; [0.,200.] +;set_tickx = 1. +;set_ticky = 1. ;; 50. +;minval = 0. +;maxval = 20. +;nlev = maxval-minval +;pal = 22 +;rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;!P.Charsize = 1.2 +;;; 0. levels +;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;;; 1. background +;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +; ;; 2. color field +; loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) +; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;;; 3. contour field +;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;;; 4. choose output +;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = 2.*transpose(ztke) +;;set_title = "Vertical wind variance (m!U2!N.s!U-2!N)" +;;minval = -6. +;;maxval = 12. +;;pal = 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = transpose(ztke) +;zex = localtime +;zey = h +;set_name = 'zTKE.ps' +;set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] +;set_titlex = 'Local Time (h)' +;set_titley = 'Altitude above surface (km)' +;set_subtitle = '' ;'Mean over the simulation domain' +;set_xrange = [8.,17.] ;set_yrange = [0.,7.] ;set_tickx = 1. ;set_ticky = 1. ;minval = 0. -;maxval = 10. +;maxval = 8. ;nlev = maxval-minval ;pal = 22 @@ -547,155 +511,225 @@ ;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = transpose(wt) -zex = localtime -zey = h -set_name = 'HF.ps' -set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux " ; -set_titlex = 'Local Time (h)' -set_titley = 'Altitude above surface (km)' -set_subtitle = '' ;'Mean over the simulation domain' -set_xrange = [8.,17.] -set_yrange = [0.,7.] ;; [0.,200.] -set_tickx = 1. -set_ticky = 1. ;; 50. -minval = -2. -maxval = 2. -nlev = floor(maxval-minval)*10 -pal = 33 -rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;minval = -0.8 -;maxval = 1.2 -;pal = 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -;; 0. levels -lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;; 1. background -loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;; 2. color field -loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) - ;;;-------------------------------------------------------------------------------------------------------------------------------- - ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone - ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5 - nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) - TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b - ;;;-------------------------------------------------------------------------------------------------------------------------------- - contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;; 4. choose output -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;restore, filename='tpot_profB' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = t -;zey = ht -;set_name = 'T.ps' -;set_title = "Potential Temperature (K)" -;set_titlex = 'Potential Temperature (K)' +; +;;restore, filename='addturb.dat' +;;modvar = SMOOTH(TEMPORARY(modvar), [0,smoothampl], /EDGE_TRUNCATE) +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = transpose(modvar) +;;zex = localtime +;;zey = h +;;set_name = 'modvar.ps' +;;set_title = "Horizontal wind speed variance (m!U2!N.s!U-2!N)" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Altitude above surface (km)' +;;set_subtitle = '' ;'Mean over the simulation domain' +;;set_xrange = [7.,17.] +;;set_yrange = [0.,7.] +;;set_tickx = 1. +;;set_ticky = 1. +;;minval = 0. +;;maxval = 10. +;;nlev = maxval-minval +;;pal = 22 +;;rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;;; 0. levels +;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;;;; 1. background +;;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;;;; 2. color field +;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) +;; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;;;; 3. contour field +;;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;;;; 4. choose output +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = transpose(wt) +;zex = localtime +;zey = h +;set_name = 'HF.ps' +;set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux " ; +;set_titlex = 'Local Time (h)' ;set_titley = 'Altitude above surface (km)' ;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [190.,230.] ;[min(t),max(t)] -;set_yrange = [0.,7.] -;set_tickx = 5. -;set_ticky = 1. -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;localtimes = [9,10,11,12,13,14,15,16,17,18] -;localtimes = [7,9,11,13,15,17] -;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] -;;localtimes = [15.0] +;set_xrange = [8.,17.] +;set_yrange = [0.,8.] ;; [0.,200.] +;set_tickx = 1. +;set_ticky = 1. ;; 50. +;minval = -1.5 ;-2. +;maxval = 2.5 ;2. +;nlev = floor(maxval-minval)*10 +;pal = 33 +;rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;minval = -0.8 +;;maxval = 1.2 +;;pal = 0 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_Start, FILENAME=set_name ;!P.Charsize = 1.2 +;;; 0. levels +;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;;; 1. background +;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;;; 2. color field +;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) +; ;;;-------------------------------------------------------------------------------------------------------------------------------- +; ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone +; ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5 +; nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) +; TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b +; ;;;-------------------------------------------------------------------------------------------------------------------------------- +; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;;; 3. contour field +;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;;; 4. choose output +;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;restore, filename='tpot_profB' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = t +;;zey = ht +;;set_name = 'T.ps' +;;set_title = "Potential Temperature (K)" +;;set_titlex = 'Potential Temperature (K)' +;;set_titley = 'Altitude above surface (km)' +;;set_subtitle = '' ;'Mean over the simulation domain' +;;set_xrange = [190.,230.] ;[min(t),max(t)] +;;set_yrange = [0.,7.] +;;set_tickx = 5. +;;set_ticky = 1. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;localtimes = [9,10,11,12,13,14,15,16,17,18] +;;localtimes = [7,9,11,13,15,17] +;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] +;;;localtimes = [15.0] +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;;altlin=0 & loadct, 0 +;;altlin=1 & loadct, 0 +;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;;plot, zefield(*,nntt), les_column, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;for ll = 1, n_elements(localtimes)-1 do begin +;; CASE altlin OF +;; 0: altlin=1 +;; 1: altlin=0 +;; ENDCASE +;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin +;;; if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin +;;endfor +;;;oplot, ro_tpot, ro_column, psym=7 +;;;xyouts, 192.0, 0.25, '09:00' +;;;xyouts, 205.5, 0.25, '11:00' +;;;xyouts, 214.5, 0.25, '13:00' +;;;xyouts, 219.5, 0.25, '17:00' +;;;xyouts, 220.5, 2.30, '17:00 [RO]' +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = stst +;;zey = ht +;;set_name = 'STST.ps' +;;set_title = "Static stability (K.m!U-1!N)" +;;set_titlex = 'Static stability (K.m!U-1!N)' +;;set_titley = 'Altitude above surface (km)' +;;set_subtitle = 'Mean over the simulation domain' +;;set_xrange = [-5.,5.] +;;set_yrange = [0.,7.] +;;set_tickx = 1. +;;set_ticky = 1. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;localtimes = [9,11,13,15,17] +;;localtimes = [12,13,14,15,16] +;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] +;;;localtimes = [15.0] +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 ;;altlin=0 & loadct, 0 -;altlin=1 & loadct, 0 -;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;;plot, zefield(*,nntt), les_column, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;for ll = 1, n_elements(localtimes)-1 do begin -; CASE altlin OF -; 0: altlin=1 -; 1: altlin=0 -; ENDCASE -; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin -;; if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin -;endfor -;;oplot, ro_tpot, ro_column, psym=7 -;;xyouts, 192.0, 0.25, '09:00' -;;xyouts, 205.5, 0.25, '11:00' -;;xyouts, 214.5, 0.25, '13:00' -;;xyouts, 219.5, 0.25, '17:00' -;;xyouts, 220.5, 2.30, '17:00 [RO]' -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = stst -;zey = ht -;set_name = 'STST.ps' -;set_title = "Static stability (K.m!U-1!N)" -;set_titlex = 'Static stability (K.m!U-1!N)' -;set_titley = 'Altitude above surface (km)' -;set_subtitle = 'Mean over the simulation domain' -;set_xrange = [-5.,5.] -;set_yrange = [0.,7.] -;set_tickx = 1. -;set_ticky = 1. -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;localtimes = [9,11,13,15,17] -;localtimes = [12,13,14,15,16] -;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] -;;localtimes = [15.0] -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;altlin=0 & loadct, 0 -;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;oplot, zefield(*,nntt), zey(*,nntt), psym=5 -;for ll = 1, n_elements(localtimes)-1 do begin -; CASE altlin OF -; 0: altlin=1 -; 1: altlin=0 -; ENDCASE -; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) -; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin -; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5 -;endfor -;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - - -if (n_elements(xtke) eq 0) then stop - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;zefield = 2.*transpose(xtke) -;;minval = -4 -;;maxval = 8. ;12. -;;pal = 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;oplot, zefield(*,nntt), zey(*,nntt), psym=5 +;;for ll = 1, n_elements(localtimes)-1 do begin +;; CASE altlin OF +;; 0: altlin=1 +;; 1: altlin=0 +;; ENDCASE +;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) +;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin +;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5 +;;endfor +;;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; +; +;if (n_elements(xtke) eq 0) then stop +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;zefield = 2.*transpose(xtke) +;;;minval = -4 +;;;maxval = 8. ;12. +;;;pal = 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zex = localtime +;;zey = h +;;set_name = 'xTKE.ps' +;;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Altitude above surface (km)' +;;set_subtitle = '';'Mean over the simulation domain' +;;set_xrange = [8.,17.] +;;set_yrange = [0.,7.] +;;set_tickx = 1. +;;set_ticky = 1. +;;minval = 0. +;;maxval = 8. +;;nlev = maxval-minval +;;pal = 22 +;;rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;;; 0. levels +;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;;;; 1. background +;;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;;;; 2. color field +;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) +;; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;;;; 3. contour field +;;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;;;; 4. choose output +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = transpose(ytke) + transpose(xtke) ;zex = localtime ;zey = h -;set_name = 'xTKE.ps' -;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] +;set_name = 'hTKE.ps' +;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] ;set_titlex = 'Local Time (h)' ;set_titley = 'Altitude above surface (km)' -;set_subtitle = '';'Mean over the simulation domain' +;set_subtitle = '' ;'Mean over the simulation domain' ;set_xrange = [8.,17.] -;set_yrange = [0.,7.] +;set_yrange = [0.,1.] ;[0.,7.] ;; [0.,200.] ;set_tickx = 1. -;set_ticky = 1. +;set_ticky = 0.1 ;1. ;50. ;minval = 0. -;maxval = 8. +;maxval = 12. ;10. ;nlev = maxval-minval ;pal = 22 @@ -716,40 +750,6 @@ ;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = transpose(ytke) + transpose(xtke) -zex = localtime -zey = h -set_name = 'hTKE.ps' -set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] -set_titlex = 'Local Time (h)' -set_titley = 'Altitude above surface (km)' -set_subtitle = '' ;'Mean over the simulation domain' -set_xrange = [8.,17.] -set_yrange = [0.,1.] ;[0.,7.] ;; [0.,200.] -set_tickx = 1. -set_ticky = 0.1 ;1. ;50. -minval = 0. -maxval = 12. ;10. -nlev = maxval-minval -pal = 22 -rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -;; 0. levels -lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;; 1. background -loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;; 2. color field -loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) - contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;; 4. choose output -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - -if (n_elements(wmax) eq 0) then stop +; +;if (n_elements(wmax) eq 0) then stop ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -762,10 +762,10 @@ set_titley = 'Altitude above surface (km)' set_subtitle = '' -set_xrange = [8.,18.] -set_yrange = [0.,7.] +set_xrange = [11.,17.] ;[8.,18.] +set_yrange = [0.,8.] ;[0.,8.] set_tickx = 1. set_ticky = 1. minval = 0. -maxval = 16. ;15.;14.;13. +maxval = 18. ;16. ;15.;14.;13. nlev = maxval-minval pal = 22 @@ -776,4 +776,5 @@ ;; 0. levels lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)] ;;; 1. background loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 @@ -782,5 +783,5 @@ contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill ;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +loadct, 0 & contour, zefield, zex, zey, levels=levhalf, c_labels=findgen(n_elements(levhalf))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (levhalf LT 0.0) ;; 4. choose output if (saveps eq 'true') then PS_End, /PNG @@ -796,10 +797,10 @@ set_titley = 'Altitude above surface (km)' set_subtitle = '' -set_xrange = [8.,18.] -set_yrange = [0.,7.] +set_xrange = [11.,17.] ;[8.,18.] +set_yrange = [0.,8.] set_tickx = 1. set_ticky = 1. minval = 0. -maxval = 9.;12. +maxval = 10.;12. nlev = maxval-minval pal = 22 @@ -810,4 +811,5 @@ ;; 0. levels lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)] ;;; 1. background loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 @@ -816,67 +818,67 @@ contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill ;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +loadct, 0 & contour, zefield, zex, zey, levels=levhalf, c_labels=findgen(n_elements(levhalf))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (levhalf LT 0.0) ;; 4. choose output if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - -restore, filename='addturb.dat' -velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -zefield = transpose(velmax) -zex = localtime -zey = h -set_name = 'velmax.ps' -set_title = "Maximum horizontal wind speed (m.s!U-1!N)" -set_titlex = 'Local Time (h)' -set_titley = 'Altitude above surface (km)' -set_subtitle = '' ;'Mean over the simulation domain' -set_xrange = [8.,17.] -set_yrange = [0.,4.] -set_tickx = 1. -set_ticky = 1. -minval = 0. -maxval = 30. ;12. -nlev = maxval-minval -pal = 22 -rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -if (saveps eq 'true') then PS_Start, FILENAME=set_name -!P.Charsize = 1.2 -;; 0. levels -lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;; 1. background -loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;;; 2. color field -loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) - contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;; 3. contour field -loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;; 4. choose output -if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; - - -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = localtime -;zey = depressions -;set_name = 'DEPR.ps' -;set_title = "Percentage of depressions in the domain" +; +; +;restore, filename='addturb.dat' +;velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;zefield = transpose(velmax) +;zex = localtime +;zey = h +;set_name = 'velmax.ps' +;set_title = "Maximum horizontal wind speed (m.s!U-1!N)" ;set_titlex = 'Local Time (h)' -;set_titley = 'Percentage' -;set_subtitle = '' -;set_xrange = [8.,18.] -;set_yrange = [0.,3.] +;set_titley = 'Altitude above surface (km)' +;set_subtitle = '' ;'Mean over the simulation domain' +;set_xrange = [11.,17.] ;[8.,17.] +;set_yrange = [0.,4.] ;set_tickx = 1. ;set_ticky = 1. +;minval = 0. +;maxval = 30. ;12. +;nlev = maxval-minval +;pal = 22 +;rrr = 'no' ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;if (saveps eq 'true') then PS_Start, FILENAME=set_name ;!P.Charsize = 1.2 -;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;oplot, zefield, zey, psym=5 -; oplot, localtime, abs(psmin), linestyle=1 +;;; 0. levels +;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;;; 1. background +;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;;;; 2. color field +;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) +; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;;; 3. contour field +;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;;; 4. choose output ;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +; +; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;zefield = localtime +;;zey = depressions +;;set_name = 'DEPR.ps' +;;set_title = "Percentage of depressions in the domain" +;;set_titlex = 'Local Time (h)' +;;set_titley = 'Percentage' +;;set_subtitle = '' +;;set_xrange = [8.,18.] +;;set_yrange = [0.,3.] +;;set_tickx = 1. +;;set_ticky = 1. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;if (saveps eq 'true') then PS_Start, FILENAME=set_name +;;!P.Charsize = 1.2 +;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;;oplot, zefield, zey, psym=5 +;; oplot, localtime, abs(psmin), linestyle=1 +;;if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; stop Index: /trunk/mesoscale/PLOT/SPEC/LES/report_tasi.inc =================================================================== --- /trunk/mesoscale/PLOT/SPEC/LES/report_tasi.inc (revision 112) +++ /trunk/mesoscale/PLOT/SPEC/LES/report_tasi.inc (revision 113) @@ -1,10 +1,19 @@ + saveplot = './' -cases = ['MERIDIANI_tau050', 'MERIDIANI_tau100', 'MERIDIANI_tau200', 'MERIDIANI_tau500'] +cases = ['baseline', 'hill', 'little_hill', 'realistic_extreme'] paths = cases -tttitle = ['tau = 0.5' $ - ,'tau = 1.0' $ - ,'tau = 2.0' $ - ,'tau = 5.0' $ - ] +tttitle = cases tttitles = tttitle shorttt = tttitle + + +;saveplot = './' +;cases = ['MERIDIANI_tau050', 'MERIDIANI_tau100', 'MERIDIANI_tau200', 'MERIDIANI_tau500'] +;paths = cases +;tttitle = ['tau = 0.5' $ +; ,'tau = 1.0' $ +; ,'tau = 2.0' $ +; ,'tau = 5.0' $ +; ] +;tttitles = tttitle +;shorttt = tttitle Index: /trunk/mesoscale/PLOT/SPEC/LES/test.pro =================================================================== --- /trunk/mesoscale/PLOT/SPEC/LES/test.pro (revision 113) +++ /trunk/mesoscale/PLOT/SPEC/LES/test.pro (revision 113) @@ -0,0 +1,123 @@ +pro test + + +getcdf, $ +; file='./wrfout_d01_9999-01-01_05:08:21', $ + file='./wrfout_d01_9999-01-01_06:10:00', $ +; file='./wrfout_d01_9999-01-01_07:11:40', $ + charvar='W', $ + invar=w +w1 = max ( w, DIMENSION=3 ) +wm1 = min ( TEMPORARY(w), DIMENSION=3 ) + +;getcdf, $ +; file='./wrfout_d01_9999-01-01_06:10:00', $ +; charvar='W', $ +; invar=w +;w2 = max ( w, DIMENSION=3 ) +;wm2 = min ( TEMPORARY(w), DIMENSION=3 ) +; +;getcdf, $ +; file='./wrfout_d01_9999-01-01_07:11:40', $ +; charvar='W', $ +; invar=w +;w3 = max ( w, DIMENSION=3 ) +;wm3 = min ( TEMPORARY(w), DIMENSION=3 ) +; +;getcdf, $ +; file='./wrfout_d01_9999-01-01_04:06:40', $ +; charvar='W', $ +; invar=w +;w4 = max ( w, DIMENSION=3 ) +;wm4 = min ( TEMPORARY(w), DIMENSION=3 ) + +getcdf, $ + file='./wrfout_d01_9999-01-01_06:10:00', $ + charvar='PSFC', $ + invar=psfc + +gros=50 ;20;5;10 +;gros=5 ;; pour le cas hill +var1 = psfc[ where( w1 gt 10.) ] +var2 = psfc +var3 = psfc[ where( wm1 lt -5.) ] +;var4 = psfc[ where( w1 gt 13.) ] ;; bien mais deja montre avec courbes max +r1 = ( float(histogram(floor(var1*gros))) ) ;/ n_elements(var1) * 100. +r2 = ( float(histogram(floor(var2*gros))) ) ;/ n_elements(var2) * 100. +r3 = ( float(histogram(floor(var3*gros))) ) ;/ n_elements(var3) * 100. +r1 = 100. * r1 / max(r1) +r2 = 100. * r2 / max(r2) +r3 = 100. * r3 / max(r3) +t1 = min(var1) + findgen(n_elements(r1))/gros +t2 = min(var2) + findgen(n_elements(r2))/gros +t3 = min(var3) + findgen(n_elements(r3))/gros +PS_START, file='stats_w.ps' + !P.Charsize = 1.2 + !p.charthick = 2.0 + !p.thick = 2.0 + !x.thick = 2.0 + !y.thick = 2.0 +!p.psym = 10 +plot, t2, r2, xrange=[min(t2)+(max(t2)-min(t2))/2., max(t2)], xtickinterval=1, xtitle='Surface pressure (Pa)', ytitle='Relative quantity (ratio to max value in %)';, yrange=[0,4] +oplot, t3, r3, psym=2 +oplot, t1, r1, psym=5 +PS_END, /PNG + +getcdf, $ + file='./wrfout_d01_9999-01-01_00:00:00', $ + charvar='PSFC', $ + invar=psfc + +xx=findgen(n_elements(psfc(*,0,0)))*50. +yy=findgen(n_elements(psfc(0,*,0)))*50. + +PS_START, file='psfc.ps' + !P.Charsize = 1.2 + !p.charthick = 2.0 + !p.thick = 2.0 + !x.thick = 2.0 + !y.thick = 2.0 +levu=[685.,687.,689.,691.,693.,695.,697.,699.,701.,703.,705.,707.,709.,711.,713.,715.,717.,719.] + +psfc(0,0,1)=684. +psfc(0,1,1)=718. + +loadct, 4 +contour, reform(psfc(*,*,1)), xx, yy, /cell_fill, nlevels=60, xtitle='x dimension (m)', ytitle='y dimension (m)', /isotropic, max_value=718., min_value=684. +contour, reform(psfc(*,*,1)), xx, yy, lev=levu, c_labels=findgen(n_elements(levu))*0+1, /overplot +PS_END, /PNG + +stop + + +;m = max(psfc) +;psfc = -11000. * alog ( psfc / m ) + +PS_START, file='stats_w.ps' + !P.Charsize = 1.2 + !p.charthick = 2.0 + !p.thick = 2.0 + !x.thick = 2.0 + !y.thick = 2.0 +plot, psfc, w1, psym=3, title='Maximum of vertical wind amplitude along the vertical (between LT 11 and 15)', xtitle='Surface pressure (Pa)', ytitle='Maximum vertical wind (m s!U-1!N)', yrange=[5.,20.], xrange=[685.,720.] +oplot, psfc, w2, psym=3 +oplot, psfc, w3, psym=3 +oplot, psfc, w4, psym=3 +PS_END, /PNG + +PS_START, file='stats_w2.ps' + !P.Charsize = 1.2 + !p.charthick = 2.0 + !p.thick = 2.0 + !x.thick = 2.0 + !y.thick = 2.0 +plot, psfc, wm1, psym=3, title='Maximum of vertical wind amplitude along the vertical (between LT 11 and 15)', xtitle='Surface pressure (Pa)', ytitle='Maximum vertical wind (m s!U-1!N)', yrange=[-12.,-4.], xrange=[685.,720.] +oplot, psfc, wm2, psym=3 +oplot, psfc, wm3, psym=3 +oplot, psfc, wm4, psym=3 + +PS_END, /PNG + + + +end