Changeset 113
- Timestamp:
- Apr 18, 2011, 8:19:21 AM (14 years ago)
- Location:
- trunk/mesoscale/PLOT/SPEC/LES
- Files:
-
- 1 added
- 3 edited
-
compturb.pro (modified) (12 diffs)
-
getturb.pro (modified) (12 diffs)
-
report_tasi.inc (modified) (1 diff)
-
test.pro (added)
Legend:
- Unmodified
- Added
- Removed
-
trunk/mesoscale/PLOT/SPEC/LES/compturb.pro
r85 r113 10 10 ;smoothampl=4. 11 11 smoothampl=0. 12 ;smoothampl=5.12 smoothampl=5. 13 13 ;;;;;;; 14 14 ; … … 41 41 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 42 42 43 ;; pour article, decommenter ici et commenter l'appel au deuxieme fichier 43 ;;; pour article, decommenter ici et commenter l'appel au deuxieme fichier 44 ;;;com 45 ;;!p.multi = [0,2,1] 46 ; 47 ; path = paths(0) & print, path 48 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 49 ; file = path+'/compturb.dat' & restore, filename=file 50 ; file = path+'/getturb.dat' & restore, filename=file 51 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 52 ; 53 ;;; 54 ;;; COEFFICIENT POUR LOIS DE SIMILITUDE 55 ;;; 56 ;;@coeffhf 57 ; 58 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 59 ;zefield = a_wt_bot 60 ;zey = a_h 61 ;set_name = saveplot+'comp_A_HFLUX.ps' 62 ;set_title = "Similarity: vert. eddy heat flux on Mars" 63 ;set_titlex = "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N" 64 ;set_titley = 'Dimensionless height z/z!Di!N' 65 ;set_subtitle = '' 66 ;set_xrange = [-0.5,1.] 67 ;set_yrange = [0.,1.5] 68 ;set_tickx = 0.25 69 ;set_ticky = 0.1 70 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 71 ;zesim = 1 - 1.2 * a_h 72 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 73 ;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 74 ;;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )] 75 ;;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )] 76 ;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] 77 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 78 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 79 ;!P.Charsize = 1.2 80 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 81 ;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 82 ;for ll = 1, n_elements(localtimes)-1 do begin 83 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 84 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin 85 ;endfor 86 ;;wherex = set_xrange(0) 87 ;;wherey = set_yrange(1) 88 ;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path 89 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 90 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 91 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 92 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 93 ;zecolor = 0. & loadct, 33 94 ;for i=1, n_elements(paths)-1 do begin 95 ;; zecolor = 50. + zecolor 96 ;zecolor = 0. & loadct, 0 97 ; path = paths(i) & print, path 98 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 99 ; file = path+'/compturb.dat' & restore, filename=file 100 ; file = path+'/getturb.dat' & restore, filename=file 101 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 102 ; zefield = a_wt_bot 103 ; zey = a_h 104 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 105 ; for ll = 0, n_elements(localtimes)-1 do begin 106 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 107 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor 108 ; endfor 109 ; ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor 110 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 111 ;endfor 112 ;loadct, 0 113 ;;oplot, [0.,0.], set_yrange, linestyle=0, color=0 114 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 115 ; ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP 116 ; ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 117 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 118 ;if (saveps eq 'true') then PS_End, /PNG 44 119 ;;com 45 ;!p.multi = [0,2,1] 46 47 path = paths(0) & print, path 48 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 49 file = path+'/compturb.dat' & restore, filename=file 50 file = path+'/getturb.dat' & restore, filename=file 51 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 52 53 ;; 54 ;; COEFFICIENT POUR LOIS DE SIMILITUDE 55 ;; 56 ;@coeffhf 57 58 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 59 zefield = a_wt_bot 60 zey = a_h 61 set_name = saveplot+'comp_A_HFLUX.ps' 62 set_title = "Similarity: vert. eddy heat flux on Mars" 63 set_titlex = "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N" 64 set_titley = 'Dimensionless height z/z!Di!N' 65 set_subtitle = '' 66 set_xrange = [-0.5,1.] 67 set_yrange = [0.,1.5] 68 set_tickx = 0.25 69 set_ticky = 0.1 70 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 71 zesim = 1 - 1.2 * a_h 72 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 73 localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 74 ;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )] 75 ;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )] 76 ;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] 77 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 78 if (saveps eq 'true') then PS_Start, FILENAME=set_name 79 !P.Charsize = 1.2 80 user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 81 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 82 for ll = 1, n_elements(localtimes)-1 do begin 83 user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 84 if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin 85 endfor 86 ;wherex = set_xrange(0) 87 ;wherey = set_yrange(1) 88 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path 89 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 90 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 91 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 92 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 93 zecolor = 0. & loadct, 33 94 for i=1, n_elements(paths)-1 do begin 95 ; zecolor = 50. + zecolor 96 zecolor = 0. & loadct, 0 97 path = paths(i) & print, path 98 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 99 file = path+'/compturb.dat' & restore, filename=file 100 file = path+'/getturb.dat' & restore, filename=file 101 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 102 zefield = a_wt_bot 103 zey = a_h 104 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 105 for ll = 0, n_elements(localtimes)-1 do begin 106 user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 107 if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor 108 endfor 109 ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor 110 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 111 endfor 112 loadct, 0 113 ;oplot, [0.,0.], set_yrange, linestyle=0, color=0 114 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 115 ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP 116 ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 117 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 118 if (saveps eq 'true') then PS_End, /PNG 119 ;com 120 121 122 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 123 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 124 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 125 126 path = paths(0) & print, path 127 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 128 file = path+'/compturb.dat' & restore, filename=file 129 file = path+'/getturb.dat' & restore, filename=file 130 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 131 132 ;; 133 ;; COEFFICIENT POUR LOIS DE SIMILITUDE 134 ;; 135 ;@coeffw 136 137 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 138 zefield = a_vel_bot 139 zey = a_h 140 set_name = saveplot+'comp_A_VEL.ps' 141 set_title = "Similarity: vert. vel. variance on Mars" 142 set_titlex = "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N" 143 set_titley = 'Dimensionless height z/z!Di!N' 144 set_subtitle = '' 145 set_xrange = [0.,1.2] 146 set_yrange = [0.,1.5] 147 set_tickx = 0.25 ;0.1 148 set_ticky = 0.1 149 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 150 zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 151 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 152 localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 153 ;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] 154 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 155 ;com 156 if (saveps eq 'true') then PS_Start, FILENAME=set_name 157 !P.Charsize = 1.2 158 user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 159 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 160 for ll = 1, n_elements(localtimes)-1 do begin 161 user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 162 if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin 163 endfor 164 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 165 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 166 zecolor = 0. & loadct, 33 167 for i=1, n_elements(paths)-1 do begin 168 zecolor = 0. & loadct, 0 169 ; zecolor = 50. + zecolor 170 path = paths(i) & print, path 171 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 172 file = path+'/compturb.dat' & restore, filename=file 173 file = path+'/getturb.dat' & restore, filename=file 174 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 175 zefield = a_vel_bot 176 zey = a_h 177 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 178 for ll = 0, n_elements(localtimes)-1 do begin 179 user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 180 if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor 181 endfor 182 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 183 endfor 184 loadct, 0 185 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 186 ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP 187 ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 188 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 189 if (saveps eq 'true') then PS_End, /PNG 190 120 ; 121 ; 122 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 123 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 124 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 125 ; 126 ; path = paths(0) & print, path 127 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 128 ; file = path+'/compturb.dat' & restore, filename=file 129 ; file = path+'/getturb.dat' & restore, filename=file 130 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 131 ; 132 ;;; 133 ;;; COEFFICIENT POUR LOIS DE SIMILITUDE 134 ;;; 135 ;;@coeffw 136 ; 137 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 138 ;zefield = a_vel_bot 139 ;zey = a_h 140 ;set_name = saveplot+'comp_A_VEL.ps' 141 ;set_title = "Similarity: vert. vel. variance on Mars" 142 ;set_titlex = "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N" 143 ;set_titley = 'Dimensionless height z/z!Di!N' 144 ;set_subtitle = '' 145 ;set_xrange = [0.,1.2] 146 ;set_yrange = [0.,1.5] 147 ;set_tickx = 0.25 ;0.1 148 ;set_ticky = 0.1 149 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 150 ;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 151 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 152 ;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 153 ;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )] 154 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 191 155 ;;com 192 ;stop 193 194 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 195 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 196 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 197 198 path = paths(0) & print, path 199 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 200 file = path+'/compturb.dat' & restore, filename=file 201 file = path+'/getturb.dat' & restore, filename=file 202 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 203 204 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 205 zefield = localtime 206 zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE) 207 set_name = saveplot+'comp_W_STAR.ps' 208 set_title = '';"Free convection velocity scale W!D*!N (m.s!U-1!N)" 209 set_titlex = 'Local Time (h)' 210 set_titley = 'Free convection velocity scale (m s!U-1!N)' 211 set_subtitle = '' 212 set_xrange = [08.,17.] 213 set_yrange = [00.,06.] 214 set_tickx = 1.0 215 set_ticky = 0.5 216 wherex = 8.5;11.0 217 wherey = 5.5;1.75 218 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 219 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 220 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 221 if (saveps eq 'true') then PS_Start, FILENAME=set_name 222 !P.Charsize = 1.2 223 altlin=0 224 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 225 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path 226 plots, wherex + 2.*dddx/1.1, wherey 227 plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 228 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 229 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 156 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 157 ;!P.Charsize = 1.2 158 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 159 ;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 160 ;for ll = 1, n_elements(localtimes)-1 do begin 161 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 162 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin 163 ;endfor 164 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 165 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 230 166 ;zecolor = 0. & loadct, 33 231 for i=1, n_elements(paths)-1 do begin 232 ;zecolor = float(255)/n_elements(paths) + zecolor 233 zecolor = 0 234 path = paths(i) & print, path 235 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 236 file = path+'/compturb.dat' & restore, filename=file 237 file = path+'/getturb.dat' & restore, filename=file 238 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 239 zefield = localtime 240 zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot 241 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 242 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 243 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 244 oplot, zefield, zey, color=zecolor, linestyle=altlin+i 245 xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 246 plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 247 plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 248 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 249 endfor 250 loadct, 0 251 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 252 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 253 if (saveps eq 'true') then PS_End, /PNG 254 255 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 256 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 257 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 167 ;for i=1, n_elements(paths)-1 do begin 168 ;zecolor = 0. & loadct, 0 169 ;; zecolor = 50. + zecolor 170 ; path = paths(i) & print, path 171 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 172 ; file = path+'/compturb.dat' & restore, filename=file 173 ; file = path+'/getturb.dat' & restore, filename=file 174 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 175 ; zefield = a_vel_bot 176 ; zey = a_h 177 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 178 ; for ll = 0, n_elements(localtimes)-1 do begin 179 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 180 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor 181 ; endfor 182 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 183 ;endfor 184 ;loadct, 0 185 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 186 ; ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP 187 ; ;oplot, set_xrange, [0.,0.], linestyle=0, color=0 188 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 189 ;if (saveps eq 'true') then PS_End, /PNG 190 ; 191 ;;;com 192 ;;stop 193 ; 194 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 195 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 196 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 197 ; 198 ; path = paths(0) & print, path 199 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 200 ; file = path+'/compturb.dat' & restore, filename=file 201 ; file = path+'/getturb.dat' & restore, filename=file 202 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 203 ; 204 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 205 ;zefield = localtime 206 ;zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE) 207 ;set_name = saveplot+'comp_W_STAR.ps' 208 ;set_title = '';"Free convection velocity scale W!D*!N (m.s!U-1!N)" 209 ;set_titlex = 'Local Time (h)' 210 ;set_titley = 'Free convection velocity scale (m s!U-1!N)' 211 ;set_subtitle = '' 212 ;set_xrange = [08.,17.] 213 ;set_yrange = [00.,06.] 214 ;set_tickx = 1.0 215 ;set_ticky = 0.5 216 ;wherex = 8.5;11.0 217 ;wherey = 5.5;1.75 218 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 219 ;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 220 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 221 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 222 ;!P.Charsize = 1.2 223 ;altlin=0 224 ;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 225 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path 226 ; plots, wherex + 2.*dddx/1.1, wherey 227 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 228 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 229 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 230 ;;zecolor = 0. & loadct, 33 231 ;for i=1, n_elements(paths)-1 do begin 232 ; ;zecolor = float(255)/n_elements(paths) + zecolor 233 ; zecolor = 0 234 ; path = paths(i) & print, path 235 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 236 ; file = path+'/compturb.dat' & restore, filename=file 237 ; file = path+'/getturb.dat' & restore, filename=file 238 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 239 ; zefield = localtime 240 ; zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot 241 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 242 ; 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 243 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 244 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 245 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 246 ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 247 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 248 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 249 ;endfor 250 ;loadct, 0 251 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 252 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 253 ;if (saveps eq 'true') then PS_End, /PNG 254 ; 255 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 256 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 257 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 258 258 259 259 path = paths(0) & print, path … … 271 271 set_titley = 'Boundary layer height (km)' 272 272 set_subtitle = '' 273 ;set_xrange = [09.,18.]274 set_yrange = [00.,0 7.]273 set_xrange = [08.,18.] 274 set_yrange = [00.,08.] ;[00.,07.] 275 275 set_tickx = 1. 276 276 set_ticky = 1. 277 277 wherex = 09.;13.25 ;12.75 278 wherey = 06. ;2.5 ;2.0278 wherey = 06.5;2.5 ;2.0 279 279 ; set_xrange = [13.5,18.] & set_yrange = [06.4,08.0] & set_tickx = 0.5 & set_ticky = 0.2 280 280 ; wherex=15.30 … … 321 321 if (saveps eq 'true') then PS_End, /PNG 322 322 323 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 324 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 325 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 326 ;; 327 ;; path = paths(0) & print, path 328 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 329 ;; file = path+'/compturb.dat' & restore, filename=file 330 ;; file = path+'/getturb.dat' & restore, filename=file 331 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 332 ;; 333 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 334 ;;zefield = localtime 335 ;;zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 336 ;;set_name = saveplot+'comp_dHEIGHT.ps' 337 ;;set_title = "Boundary layer height growth dz!Di!N/dt (km/hour)" 338 ;;set_titlex = 'Local Time (h)' 339 ;;set_titley = 'Boundary layer height growth (km/hour)' 340 ;;set_subtitle = '' 341 ;;set_xrange = [9.,18.] 342 ;;set_yrange = [-0.4,1.6] 343 ;;set_tickx = 1. 344 ;;set_ticky = 0.2 345 ;;wherex = 10.30 346 ;;wherey = 0.1 347 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 348 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 349 ;;!P.Charsize = 1.2 350 ;;altlin = 0 351 ;;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 352 ;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path 353 ;; plots, wherex + 2.*dddx/1.1, wherey 354 ;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 355 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 356 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 357 ;;;zecolor = 0. & loadct, 33 358 ;;for i=1, n_elements(paths)-1 do begin 359 ;; ;zecolor = float(255)/n_elements(paths) + zecolor 360 ;; zecolor = 0 361 ;; path = paths(i) & print, path 362 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 363 ;; file = path+'/compturb.dat' & restore, filename=file 364 ;; file = path+'/getturb.dat' & restore, filename=file 365 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 366 ;; zefield = localtime 367 ;; zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 368 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 369 ;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 370 ;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor 371 ;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 372 ;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 373 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 374 ;;endfor 375 ;;loadct, 0 376 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 377 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 378 ;;if (saveps eq 'true') then PS_End, /PNG 379 ; 323 380 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 324 381 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; … … 333 390 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 334 391 ;zefield = localtime 335 ;zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 336 ;set_name = saveplot+'comp_dHEIGHT.ps' 337 ;set_title = "Boundary layer height growth dz!Di!N/dt (km/hour)" 392 ;zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) 393 ;zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 394 ;set_name = saveplot+'comp_WMAX.ps' 395 ;set_title = '';"Maximum and minimum vertical velocity (m.s!U-1!N)" 338 396 ;set_titlex = 'Local Time (h)' 339 ;set_titley = ' Boundary layer height growth (km/hour)'397 ;set_titley = 'Maximum and minimum vertical velocity (m s!U-1!N)' 340 398 ;set_subtitle = '' 341 ; set_xrange = [9.,18.]342 ;set_yrange = [- 0.4,1.6]399 ;;set_xrange = [09.0,18.0] 400 ;set_yrange = [-12.0,18.0];[-15.0,20.0] 343 401 ;set_tickx = 1. 344 ;set_ticky = 0.2345 ;wherex = 10.30346 ;wherey = 0.1402 ;set_ticky = 2. 403 ;wherex = 8.5;11.00 404 ;wherey = 16.5;5.5 347 405 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 348 406 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name … … 350 408 ;altlin = 0 351 409 ;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 352 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path 410 ;oplot, zefield, zey2, color=0, linestyle=altlin 411 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) 353 412 ; plots, wherex + 2.*dddx/1.1, wherey 354 413 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 355 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 356 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 414 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 415 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 357 416 ;;zecolor = 0. & loadct, 33 358 417 ;for i=1, n_elements(paths)-1 do begin 359 418 ; ;zecolor = float(255)/n_elements(paths) + zecolor 360 ; zecolor = 0419 ; zecolor = 0. 361 420 ; path = paths(i) & print, path 362 421 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; … … 365 424 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 366 425 ; zefield = localtime 367 ; zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 426 ; zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) 427 ; zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 368 428 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 369 429 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 370 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor 430 ; oplot, zefield, zey2, color=zecolor, linestyle=altlin+i 431 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 371 432 ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 372 433 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue … … 377 438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 378 439 ;if (saveps eq 'true') then PS_End, /PNG 379 380 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 381 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 382 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 383 384 path = paths(0) & print, path 385 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 386 file = path+'/compturb.dat' & restore, filename=file 387 file = path+'/getturb.dat' & restore, filename=file 388 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 389 390 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 391 zefield = localtime 392 zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) 393 zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 394 set_name = saveplot+'comp_WMAX.ps' 395 set_title = '';"Maximum and minimum vertical velocity (m.s!U-1!N)" 396 set_titlex = 'Local Time (h)' 397 set_titley = 'Maximum and minimum vertical velocity (m s!U-1!N)' 398 set_subtitle = '' 399 ;set_xrange = [09.0,18.0] 400 set_yrange = [-12.0,18.0];[-15.0,20.0] 401 set_tickx = 1. 402 set_ticky = 2. 403 wherex = 8.5;11.00 404 wherey = 16.5;5.5 405 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 406 if (saveps eq 'true') then PS_Start, FILENAME=set_name 407 !P.Charsize = 1.2 408 altlin = 0 409 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 410 oplot, zefield, zey2, color=0, linestyle=altlin 411 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) 412 plots, wherex + 2.*dddx/1.1, wherey 413 plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 414 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 415 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 416 ;zecolor = 0. & loadct, 33 417 for i=1, n_elements(paths)-1 do begin 418 ;zecolor = float(255)/n_elements(paths) + zecolor 419 zecolor = 0. 420 path = paths(i) & print, path 421 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 422 file = path+'/compturb.dat' & restore, filename=file 423 file = path+'/getturb.dat' & restore, filename=file 424 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 425 zefield = localtime 426 zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE) 427 zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 428 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 429 oplot, zefield, zey, color=zecolor, linestyle=altlin+i 430 oplot, zefield, zey2, color=zecolor, linestyle=altlin+i 431 xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 432 plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 433 plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 434 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 435 endfor 436 loadct, 0 437 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 439 if (saveps eq 'true') then PS_End, /PNG 440 441 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 442 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 443 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 444 445 path = paths(0) & print, path 446 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 447 file = path+'/compturb.dat' & restore, filename=file 448 file = path+'/getturb.dat' & restore, filename=file 449 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 450 451 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 452 zefield = localtime 453 zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 454 set_name = saveplot+'comp_TKE.ps' 455 set_title = '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)" 456 set_titlex = 'Local Time (h)' 457 set_titley = 'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)' 458 set_subtitle = '' 459 ;set_xrange = [09.0,18.0] 460 set_yrange = [00.0,18.0];[00.0,26.0] 461 set_tickx = 1. 462 set_ticky = 2. 463 wherex = 8.5;11.70 464 wherey = 16.;5.4 465 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 466 if (saveps eq 'true') then PS_Start, FILENAME=set_name 467 !P.Charsize = 1.2 468 altlin = 0 469 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 470 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0) 471 plots, wherex + 2.*dddx/1.1, wherey 472 plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 473 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 475 ;zecolor = 0. & loadct, 33 476 for i=1, n_elements(paths)-1 do begin 477 ;zecolor = float(255)/n_elements(paths) + zecolor 478 zecolor = 0. 479 path = paths(i) & print, path 480 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 481 file = path+'/compturb.dat' & restore, filename=file 482 file = path+'/getturb.dat' & restore, filename=file 483 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 484 zefield = localtime 485 zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 486 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 487 oplot, zefield, zey, color=zecolor, linestyle=altlin+i 488 xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitles(i), color=zecolor 489 plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 490 plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 491 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 492 493 endfor 494 loadct, 0 495 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 496 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 497 if (saveps eq 'true') then PS_End, /PNG 498 499 500 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 501 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 502 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 503 504 path = paths(0) & print, path 505 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 506 file = path+'/compturb.dat' & restore, filename=file 507 file = path+'/getturb.dat' & restore, filename=file 508 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 509 510 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 511 zefield = localtime 512 zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 513 set_name = saveplot+'comp_HFBOT.ps' 514 set_title = '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)" 515 set_titlex = 'Local Time (h)' 516 set_titley = 'Maximum vertical eddy heat flux (K m s!U-1!N)' 517 set_subtitle = '' 518 ;set_xrange = [09.0,18.0] 519 set_yrange = [00.0,02.0];[00.0,02.8] 520 set_tickx = 1. 521 set_ticky = 0.2 522 wherex = 8.5;10.30 523 wherey = 1.8;0.7 524 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 525 if (saveps eq 'true') then PS_Start, FILENAME=set_name 526 !P.Charsize = 1.2 527 altlin = 0 528 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 529 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) 530 plots, wherex + 2.*dddx/1.1, wherey 531 plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 532 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 533 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 534 ;zecolor = 0. & loadct, 33 535 for i=1, n_elements(paths)-1 do begin 536 ;zecolor = float(255)/n_elements(paths) + zecolor 537 zecolor = 0. 538 path = paths(i) & print, path 539 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 540 file = path+'/compturb.dat' & restore, filename=file 541 file = path+'/getturb.dat' & restore, filename=file 542 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 543 zefield = localtime 544 zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 545 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 546 oplot, zefield, zey, color=zecolor, linestyle=altlin+i 547 xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 548 plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 549 plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 550 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 551 552 endfor 553 loadct, 0 554 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 555 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 556 if (saveps eq 'true') then PS_End, /PNG 557 558 440 ; 559 441 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 560 442 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; … … 569 451 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 570 452 ;zefield = localtime 571 ;zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 572 ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 573 ;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 574 ;set_name = saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps' 575 ;set_title = "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)" 453 ;zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 454 ;set_name = saveplot+'comp_TKE.ps' 455 ;set_title = '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)" 576 456 ;set_titlex = 'Local Time (h)' 577 ;set_titley = ' Potential temperature time derivative (K.hour!U-1!N)'457 ;set_titley = 'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)' 578 458 ;set_subtitle = '' 579 ; set_xrange = [ 09., 18.]580 ;set_yrange = [ -01., 06.]459 ;;set_xrange = [09.0,18.0] 460 ;set_yrange = [00.0,18.0];[00.0,26.0] 581 461 ;set_tickx = 1. 582 ;set_ticky = 1.583 ;wherex = 13.25584 ;wherey = 5.5462 ;set_ticky = 2. 463 ;wherex = 8.5;11.70 464 ;wherey = 16.;5.4 585 465 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 586 466 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 467 ;!P.Charsize = 1.2 468 ;altlin = 0 469 ;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 470 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0) 471 ; plots, wherex + 2.*dddx/1.1, wherey 472 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 473 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 475 ;;zecolor = 0. & loadct, 33 476 ;for i=1, n_elements(paths)-1 do begin 477 ; ;zecolor = float(255)/n_elements(paths) + zecolor 478 ; zecolor = 0. 479 ; path = paths(i) & print, path 480 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 481 ; file = path+'/compturb.dat' & restore, filename=file 482 ; file = path+'/getturb.dat' & restore, filename=file 483 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 484 ; zefield = localtime 485 ; zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 486 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 487 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 488 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitles(i), color=zecolor 489 ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 490 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 491 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 492 ; 493 ;endfor 494 ;loadct, 0 495 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 496 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 497 ;if (saveps eq 'true') then PS_End, /PNG 498 ; 499 ; 500 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 501 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 502 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 503 ; 504 ; path = paths(0) & print, path 505 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 506 ; file = path+'/compturb.dat' & restore, filename=file 507 ; file = path+'/getturb.dat' & restore, filename=file 508 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 509 ; 510 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 511 ;zefield = localtime 512 ;zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 513 ;set_name = saveplot+'comp_HFBOT.ps' 514 ;set_title = '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)" 515 ;set_titlex = 'Local Time (h)' 516 ;set_titley = 'Maximum vertical eddy heat flux (K m s!U-1!N)' 517 ;set_subtitle = '' 518 ;;set_xrange = [09.0,18.0] 519 ;set_yrange = [00.0,02.0];[00.0,02.8] 520 ;set_tickx = 1. 521 ;set_ticky = 0.2 522 ;wherex = 8.5;10.30 523 ;wherey = 1.8;0.7 524 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 525 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 587 526 ;!P.Charsize = 1.2 588 527 ;altlin = 0 … … 591 530 ; plots, wherex + 2.*dddx/1.1, wherey 592 531 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 532 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 533 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 534 ;;zecolor = 0. & loadct, 33 535 ;for i=1, n_elements(paths)-1 do begin 536 ; ;zecolor = float(255)/n_elements(paths) + zecolor 537 ; zecolor = 0. 538 ; path = paths(i) & print, path 539 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 540 ; file = path+'/compturb.dat' & restore, filename=file 541 ; file = path+'/getturb.dat' & restore, filename=file 542 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 543 ; zefield = localtime 544 ; zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 545 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 546 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 547 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 548 ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 549 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 550 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 551 ; 552 ;endfor 553 ;loadct, 0 554 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 555 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 556 ;if (saveps eq 'true') then PS_End, /PNG 557 ; 558 ; 559 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 560 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 561 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 562 ;; 563 ;; path = paths(0) & print, path 564 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 565 ;; file = path+'/compturb.dat' & restore, filename=file 566 ;; file = path+'/getturb.dat' & restore, filename=file 567 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 568 ;; 569 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 570 ;;zefield = localtime 571 ;;zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 572 ;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 573 ;;;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 574 ;;set_name = saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps' 575 ;;set_title = "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)" 576 ;;set_titlex = 'Local Time (h)' 577 ;;set_titley = 'Potential temperature time derivative (K.hour!U-1!N)' 578 ;;set_subtitle = '' 579 ;;set_xrange = [ 09., 18.] 580 ;;set_yrange = [-01., 06.] 581 ;;set_tickx = 1. 582 ;;set_ticky = 1. 583 ;;wherex = 13.25 584 ;;wherey = 5.5 585 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 586 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 587 ;;!P.Charsize = 1.2 588 ;;altlin = 0 589 ;;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 590 ;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) 591 ;; plots, wherex + 2.*dddx/1.1, wherey 592 ;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 594 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 595 ;;;zecolor = 0. & loadct, 33 596 ;;for i=1, n_elements(paths)-1 do begin 597 ;; ;zecolor = float(255)/n_elements(paths) + zecolor 598 ;; zecolor = 0. 599 ;; path = paths(i) & print, path 600 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 601 ;; file = path+'/compturb.dat' & restore, filename=file 602 ;; file = path+'/getturb.dat' & restore, filename=file 603 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 604 ;; zefield = localtime 605 ;; zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 606 ;; ;zey = DERIV(localtime, t_bot) 607 ;; zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 608 ;; ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 609 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 610 ;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 611 ;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor 612 ;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 613 ;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 614 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 615 ;;endfor 616 ;;loadct, 0 617 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 618 ;;;oplot, set_xrange, [0,0], linestyle=2 619 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 620 ;;if (saveps eq 'true') then PS_End, /PNG 621 ; 622 ; 623 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 624 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 625 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 626 ;; 627 ;; path = paths(0) & print, path 628 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 629 ;; file = path+'/compturb.dat' & restore, filename=file 630 ;; file = path+'/getturb.dat' & restore, filename=file 631 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 632 ;; 633 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 634 ;;zefield = localtime 635 ;;zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top 636 ;;set_name = saveplot+'comp_TTOP.ps' 637 ;;set_title = "Temperature at the top of the boundary layer (K)" 638 ;;set_titlex = 'Local Time (h)' 639 ;;set_titley = 'Temperature at the top of the boundary layer (K)' 640 ;;set_subtitle = '' 641 ;;set_xrange = [09.,18.] 642 ;;set_yrange = [195.,225.] 643 ;;set_tickx = 1. 644 ;;set_ticky = 5. 645 ;;wherex = 13.00 646 ;;wherey = 223.5 647 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 648 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 649 ;;!P.Charsize = 1.2 650 ;;altlin = 0 651 ;;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 652 ;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path 653 ;; plots, wherex + 2.*dddx/1.1, wherey 654 ;; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 655 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 656 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 657 ;;;zecolor = 0. & loadct, 33 658 ;;for i=1, n_elements(paths)-1 do begin 659 ;; ;zecolor = float(255)/n_elements(paths) + zecolor 660 ;; zecolor = 0. 661 ;; path = paths(i) & print, path 662 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 663 ;; file = path+'/compturb.dat' & restore, filename=file 664 ;; file = path+'/getturb.dat' & restore, filename=file 665 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 666 ;; zefield = localtime 667 ;; zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top 668 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 669 ;; oplot, zefield, zey, color=zecolor, linestyle=altlin+i 670 ;; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor 671 ;; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy 672 ;; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 673 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 674 ;;endfor 675 ;;loadct, 0 676 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 677 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 678 ;;if (saveps eq 'true') then PS_End, /PNG 679 ; 680 ; 681 ;;;;;;;;;;;;;;;; 682 ;;;;;;;;;;;;;;;; 683 ;lev = 4 684 ;;;;;;;;;;;;;;;; 685 ;;;;;;;;;;;;;;;; 686 ; 687 ; path = paths(0) & print, path 688 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 689 ; file = path+'/compturb.dat' & restore, filename=file 690 ; file = path+'/addturb.dat' & restore, filename=file 691 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 692 ; 693 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 694 ;zefield = localtime 695 ;zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) 696 ;set_name = saveplot+'comp_VELMAX.ps' 697 ;set_title = '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)" 698 ;set_titlex = 'Local Time (h)' 699 ;set_titley = "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)" 700 ;set_subtitle = '' 701 ;;set_xrange = [09.0,18.0] 702 ;set_yrange = [00.0,14.0] 703 ;set_tickx = 1. 704 ;set_ticky = 2. 705 ;wherex = 8.5;10.30 706 ;wherey = 12.;0.7 707 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 708 ;print, h(lev) 709 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 710 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 711 ;!P.Charsize = 1.2 712 ;altlin = 0 713 ;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 714 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0) 715 ; plots, wherex + 2.*dddx/1.1, wherey 716 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue 717 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 594 718 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 595 719 ;;zecolor = 0. & loadct, 33 … … 600 724 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 601 725 ; file = path+'/compturb.dat' & restore, filename=file 602 ; file = path+'/ getturb.dat' & restore, filename=file726 ; file = path+'/addturb.dat' & restore, filename=file 603 727 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 604 728 ; zefield = localtime 605 ; zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 606 ; ;zey = DERIV(localtime, t_bot) 607 ; zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 608 ; ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 729 ; zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE) 609 730 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 610 731 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i … … 613 734 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue 614 735 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 736 ; 615 737 ;endfor 616 738 ;loadct, 0 617 739 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 618 ;;oplot, set_xrange, [0,0], linestyle=2619 740 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 620 741 ;if (saveps eq 'true') then PS_End, /PNG 621 622 623 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;624 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;625 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;626 ;627 ; path = paths(0) & print, path628 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;629 ; file = path+'/compturb.dat' & restore, filename=file630 ; file = path+'/getturb.dat' & restore, filename=file631 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;632 ;633 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;634 ;zefield = localtime635 ;zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top636 ;set_name = saveplot+'comp_TTOP.ps'637 ;set_title = "Temperature at the top of the boundary layer (K)"638 ;set_titlex = 'Local Time (h)'639 ;set_titley = 'Temperature at the top of the boundary layer (K)'640 ;set_subtitle = ''641 ;set_xrange = [09.,18.]642 ;set_yrange = [195.,225.]643 ;set_tickx = 1.644 ;set_ticky = 5.645 ;wherex = 13.00646 ;wherey = 223.5647 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;648 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name649 ;!P.Charsize = 1.2650 ;altlin = 0651 ;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=altlin652 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path653 ; plots, wherex + 2.*dddx/1.1, wherey654 ; plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue655 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;656 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;657 ;;zecolor = 0. & loadct, 33658 ;for i=1, n_elements(paths)-1 do begin659 ; ;zecolor = float(255)/n_elements(paths) + zecolor660 ; zecolor = 0.661 ; path = paths(i) & print, path662 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;663 ; file = path+'/compturb.dat' & restore, filename=file664 ; file = path+'/getturb.dat' & restore, filename=file665 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;666 ; zefield = localtime667 ; zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top668 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;669 ; oplot, zefield, zey, color=zecolor, linestyle=altlin+i670 ; xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, path, color=zecolor671 ; plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy672 ; plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue673 ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;674 ;endfor675 ;loadct, 0676 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;677 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;678 ;if (saveps eq 'true') then PS_End, /PNG679 680 681 ;;;;;;;;;;;;;;;682 ;;;;;;;;;;;;;;;683 lev = 4684 ;;;;;;;;;;;;;;;685 ;;;;;;;;;;;;;;;686 687 path = paths(0) & print, path688 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;689 file = path+'/compturb.dat' & restore, filename=file690 file = path+'/addturb.dat' & restore, filename=file691 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;692 693 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;694 zefield = localtime695 zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)696 set_name = saveplot+'comp_VELMAX.ps'697 set_title = '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)"698 set_titlex = 'Local Time (h)'699 set_titley = "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)"700 set_subtitle = ''701 ;set_xrange = [09.0,18.0]702 set_yrange = [00.0,14.0]703 set_tickx = 1.704 set_ticky = 2.705 wherex = 8.5;10.30706 wherey = 12.;0.7707 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;708 print, h(lev)709 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;710 if (saveps eq 'true') then PS_Start, FILENAME=set_name711 !P.Charsize = 1.2712 altlin = 0713 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=altlin714 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)715 plots, wherex + 2.*dddx/1.1, wherey716 plots, wherex + 2.*dddx/20., wherey, linestyle=altlin, /continue717 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;718 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;719 ;zecolor = 0. & loadct, 33720 for i=1, n_elements(paths)-1 do begin721 ;zecolor = float(255)/n_elements(paths) + zecolor722 zecolor = 0.723 path = paths(i) & print, path724 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;725 file = path+'/compturb.dat' & restore, filename=file726 file = path+'/addturb.dat' & restore, filename=file727 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;728 zefield = localtime729 zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)730 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;731 oplot, zefield, zey, color=zecolor, linestyle=altlin+i732 xyouts, wherex + 2.*dddx, wherey - float(i)*dddy, tttitle(i), color=zecolor733 plots, wherex + 2.*dddx/1.1, wherey - float(i)*dddy734 plots, wherex + 2.*dddx/20., wherey - float(i)*dddy, linestyle=altlin+i, /continue735 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;736 737 endfor738 loadct, 0739 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;740 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;741 if (saveps eq 'true') then PS_End, /PNG742 742 743 743 -
trunk/mesoscale/PLOT/SPEC/LES/getturb.pro
r96 r113 17 17 ;; constantes 18 18 ;; 19 ;p0=610. & t0=220. & r_cp=1/4.4 & grav=3.72 & R=192.20 p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191.21 print, 'ATTENTION ATTENTION R/cp !!!!', r_cp19 p0=610. & t0=220. & r_cp=1.0/4.4 & grav=3.72 & R=192. 20 ;p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191. 21 ;print, 'ATTENTION ATTENTION R/cp !!!!', r_cp 22 22 23 23 ; … … 171 171 OPENR, 23, 'input_more' & READF, 23, hgtu, tsurfu & CLOSE, 23 172 172 173 @include_ustar.pro173 ;@include_ustar.pro 174 174 175 175 ;@include_hfx.pro … … 344 344 ;*******************; 345 345 346 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 347 ;zefield = a_vel_bot 348 ;zey = a_h 349 ;set_name = 'A_VEL.ps' 350 ;set_title = "Scaled vertical velocity variance (m.s!U-1!N)" 351 ;set_titlex = 'Scaled vertical velocity variance (m.s!U-1!N)' 352 ;set_titley = 'Scaled height (km)' 353 ;set_subtitle = '' 354 ;set_xrange = [0.,1.2] 355 ;set_yrange = [0.,1.2] 356 ;set_tickx = 0.1 357 ;set_ticky = 0.1 358 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 359 ;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 360 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 361 ;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 362 ;;localtimes = localtime 363 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 364 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 365 ;!P.Charsize = 1.2 366 ;altlin=0 367 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 368 ;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 369 ;for ll = 1, n_elements(localtimes)-1 do begin 370 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 371 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 372 ;endfor 373 ;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0 374 ;if (saveps eq 'true') then PS_End, /PNG 375 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 376 377 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 378 zefield = localtime 379 zey = pbl_height 380 set_name = 'HEIGHT.ps' 381 set_title = 'Boundary layer height (km)' 382 set_titlex = 'Local Time (h)' 383 set_titley = 'Boundary layer height (km)' 384 set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' 385 set_xrange = [8.,17.] 386 set_yrange = [0.,7.] 387 set_tickx = 1. 388 set_ticky = 1. 389 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 390 if (saveps eq 'true') then PS_Start, FILENAME=set_name 391 !P.Charsize = 1.2 392 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 393 ;oplot, zefield, zey, psym=5 394 if (saveps eq 'true') then PS_End, /PNG 395 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 396 346 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 347 ;;zefield = a_vel_bot 348 ;;zey = a_h 349 ;;set_name = 'A_VEL.ps' 350 ;;set_title = "Scaled vertical velocity variance (m.s!U-1!N)" 351 ;;set_titlex = 'Scaled vertical velocity variance (m.s!U-1!N)' 352 ;;set_titley = 'Scaled height (km)' 353 ;;set_subtitle = '' 354 ;;set_xrange = [0.,1.2] 355 ;;set_yrange = [0.,1.2] 356 ;;set_tickx = 0.1 357 ;;set_ticky = 0.1 358 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 359 ;;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2 360 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 361 ;;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )] 362 ;;;localtimes = localtime 363 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 364 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 365 ;;!P.Charsize = 1.2 366 ;;altlin=0 367 ;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 368 ;;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 369 ;;for ll = 1, n_elements(localtimes)-1 do begin 370 ;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 371 ;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 372 ;;endfor 373 ;;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0 374 ;;if (saveps eq 'true') then PS_End, /PNG 375 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 376 ; 397 377 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 398 378 ;zefield = localtime 399 ;zey = w_star_bot400 ;set_name = ' W_STAR.ps'401 ;set_title = "Free convection velocity scale (m.s!U-1!N)"379 ;zey = pbl_height 380 ;set_name = 'HEIGHT.ps' 381 ;set_title = 'Boundary layer height (km)' 402 382 ;set_titlex = 'Local Time (h)' 403 ;set_titley = ' Free convection velocity scale (m.s!U-1!N)'404 ;set_subtitle = '' 405 ;set_xrange = [ 8.,18.]406 ;set_yrange = [0., 6.]383 ;set_titley = 'Boundary layer height (km)' 384 ;set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' 385 ;set_xrange = [11.,17.] ;[8.,17.] 386 ;set_yrange = [0.,8.] 407 387 ;set_tickx = 1. 408 ;set_ticky = 0.5388 ;set_ticky = 1. 409 389 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 410 390 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name … … 414 394 ;if (saveps eq 'true') then PS_End, /PNG 415 395 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 416 417 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 418 ;zefield = localtime 419 ;zey = hfbot 420 ;zey2 = hfsurf 421 ;zey3 = hfbotmax 422 ;set_name = 'HFBOT.ps' 423 ;set_title = "Heat flux at the top of the radiative layer (K.m.s!U-1!N)" 424 ;set_titlex = 'Local Time (h)' 425 ;set_titley = 'Heat flux at the top of the radiative layer (K.m.s!U-1!N)' 426 ;set_subtitle = '' 427 ;set_xrange = [8.,18.] 428 ;set_yrange = [0.,2.5] 429 ;set_tickx = 1. 430 ;set_ticky = 0.5 431 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 432 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 433 ;!P.Charsize = 1.2 434 ;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 435 ;oplot, zefield, zey2, linestyle=1 436 ;oplot, zefield, zey3, linestyle=2 437 ;if (saveps eq 'true') then PS_End, /PNG 438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 439 440 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 441 zefield = transpose(tke) 442 zex = localtime 443 zey = h 444 set_name = 'TKE.ps' 445 set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>] 446 set_titlex = 'Local Time (h)' 447 set_titley = 'Altitude above surface (km)' 448 set_subtitle = '' ;'Mean over the simulation domain' 449 set_xrange = [8.,17.] 450 set_yrange = [0.,7.] ;; [0.,200.] 451 set_tickx = 1. 452 set_ticky = 1. ;; 50. 453 minval = 0. 454 maxval = 15. 455 nlev = maxval-minval 456 pal = 22 457 rrr = 'no' 458 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 459 if (saveps eq 'true') then PS_Start, FILENAME=set_name 460 !P.Charsize = 1.2 461 ;; 0. levels 462 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 463 ;;; 1. background 464 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 465 ;; 2. color field 466 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 467 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 468 ;; 3. contour field 469 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) 470 ;; 4. choose output 471 if (saveps eq 'true') then PS_End, /PNG 472 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 473 474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 475 ;zefield = 2.*transpose(ztke) 476 ;set_title = "Vertical wind variance (m!U2!N.s!U-2!N)" 477 ;minval = -6. 478 ;maxval = 12. 479 ;pal = 0 480 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 481 zefield = transpose(ztke) 482 zex = localtime 483 zey = h 484 set_name = 'zTKE.ps' 485 set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>] 486 set_titlex = 'Local Time (h)' 487 set_titley = 'Altitude above surface (km)' 488 set_subtitle = '' ;'Mean over the simulation domain' 489 set_xrange = [8.,17.] 490 set_yrange = [0.,7.] 491 set_tickx = 1. 492 set_ticky = 1. 493 minval = 0. 494 maxval = 8. 495 nlev = maxval-minval 496 pal = 22 497 rrr = 'no' 498 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 499 if (saveps eq 'true') then PS_Start, FILENAME=set_name 500 !P.Charsize = 1.2 501 ;; 0. levels 502 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 503 ;;; 1. background 504 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 505 ;; 2. color field 506 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 507 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 508 ;; 3. contour field 509 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) 510 ;; 4. choose output 511 if (saveps eq 'true') then PS_End, /PNG 512 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 513 514 ;restore, filename='addturb.dat' 515 ;modvar = SMOOTH(TEMPORARY(modvar), [0,smoothampl], /EDGE_TRUNCATE) 516 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 517 ;zefield = transpose(modvar) 518 ;zex = localtime 396 ; 397 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 398 ;;zefield = localtime 399 ;;zey = w_star_bot 400 ;;set_name = 'W_STAR.ps' 401 ;;set_title = "Free convection velocity scale (m.s!U-1!N)" 402 ;;set_titlex = 'Local Time (h)' 403 ;;set_titley = 'Free convection velocity scale (m.s!U-1!N)' 404 ;;set_subtitle = '' 405 ;;set_xrange = [8.,18.] 406 ;;set_yrange = [0.,6.] 407 ;;set_tickx = 1. 408 ;;set_ticky = 0.5 409 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 410 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 411 ;;!P.Charsize = 1.2 412 ;;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 413 ;;;oplot, zefield, zey, psym=5 414 ;;if (saveps eq 'true') then PS_End, /PNG 415 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 416 ; 417 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 418 ;;zefield = localtime 419 ;;zey = hfbot 420 ;;zey2 = hfsurf 421 ;;zey3 = hfbotmax 422 ;;set_name = 'HFBOT.ps' 423 ;;set_title = "Heat flux at the top of the radiative layer (K.m.s!U-1!N)" 424 ;;set_titlex = 'Local Time (h)' 425 ;;set_titley = 'Heat flux at the top of the radiative layer (K.m.s!U-1!N)' 426 ;;set_subtitle = '' 427 ;;set_xrange = [8.,18.] 428 ;;set_yrange = [0.,2.5] 429 ;;set_tickx = 1. 430 ;;set_ticky = 0.5 431 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 432 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 433 ;;!P.Charsize = 1.2 434 ;;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 435 ;;oplot, zefield, zey2, linestyle=1 436 ;;oplot, zefield, zey3, linestyle=2 437 ;;if (saveps eq 'true') then PS_End, /PNG 438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 439 ; 440 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 441 ;zefield = transpose(tke) 442 ;zex = localtime 519 443 ;zey = h 520 ;set_name = ' modvar.ps'521 ;set_title = " Horizontal wind speed variance (m!U2!N.s!U-2!N)"444 ;set_name = 'TKE.ps' 445 ;set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>] 522 446 ;set_titlex = 'Local Time (h)' 523 447 ;set_titley = 'Altitude above surface (km)' 524 448 ;set_subtitle = '' ;'Mean over the simulation domain' 525 ;set_xrange = [7.,17.] 449 ;set_xrange = [11.,17.] ;[8.,17.] 450 ;set_yrange = [0.,8.] ;; [0.,200.] 451 ;set_tickx = 1. 452 ;set_ticky = 1. ;; 50. 453 ;minval = 0. 454 ;maxval = 20. 455 ;nlev = maxval-minval 456 ;pal = 22 457 ;rrr = 'no' 458 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 459 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 460 ;!P.Charsize = 1.2 461 ;;; 0. levels 462 ;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 463 ;;;; 1. background 464 ;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 465 ; ;; 2. color field 466 ; loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 467 ; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 468 ;;; 3. contour field 469 ;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) 470 ;;; 4. choose output 471 ;if (saveps eq 'true') then PS_End, /PNG 472 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 473 ; 474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 475 ;;zefield = 2.*transpose(ztke) 476 ;;set_title = "Vertical wind variance (m!U2!N.s!U-2!N)" 477 ;;minval = -6. 478 ;;maxval = 12. 479 ;;pal = 0 480 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 481 ;zefield = transpose(ztke) 482 ;zex = localtime 483 ;zey = h 484 ;set_name = 'zTKE.ps' 485 ;set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>] 486 ;set_titlex = 'Local Time (h)' 487 ;set_titley = 'Altitude above surface (km)' 488 ;set_subtitle = '' ;'Mean over the simulation domain' 489 ;set_xrange = [8.,17.] 526 490 ;set_yrange = [0.,7.] 527 491 ;set_tickx = 1. 528 492 ;set_ticky = 1. 529 493 ;minval = 0. 530 ;maxval = 10.494 ;maxval = 8. 531 495 ;nlev = maxval-minval 532 496 ;pal = 22 … … 547 511 ;if (saveps eq 'true') then PS_End, /PNG 548 512 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 549 550 551 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 552 zefield = transpose(wt) 553 zex = localtime 554 zey = h 555 set_name = 'HF.ps' 556 set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'> 557 set_titlex = 'Local Time (h)' 558 set_titley = 'Altitude above surface (km)' 559 set_subtitle = '' ;'Mean over the simulation domain' 560 set_xrange = [8.,17.] 561 set_yrange = [0.,7.] ;; [0.,200.] 562 set_tickx = 1. 563 set_ticky = 1. ;; 50. 564 minval = -2. 565 maxval = 2. 566 nlev = floor(maxval-minval)*10 567 pal = 33 568 rrr = 'no' 569 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 570 ;minval = -0.8 571 ;maxval = 1.2 572 ;pal = 0 573 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 574 if (saveps eq 'true') then PS_Start, FILENAME=set_name 575 !P.Charsize = 1.2 576 ;; 0. levels 577 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 578 ;;; 1. background 579 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 580 ;; 2. color field 581 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 582 ;;;-------------------------------------------------------------------------------------------------------------------------------- 583 ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone 584 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 585 nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) 586 TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b 587 ;;;-------------------------------------------------------------------------------------------------------------------------------- 588 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 589 ;; 3. contour field 590 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) 591 ;; 4. choose output 592 if (saveps eq 'true') then PS_End, /PNG 593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 594 595 ;;restore, filename='tpot_profB' 596 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 597 ;zefield = t 598 ;zey = ht 599 ;set_name = 'T.ps' 600 ;set_title = "Potential Temperature (K)" 601 ;set_titlex = 'Potential Temperature (K)' 513 ; 514 ;;restore, filename='addturb.dat' 515 ;;modvar = SMOOTH(TEMPORARY(modvar), [0,smoothampl], /EDGE_TRUNCATE) 516 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 517 ;;zefield = transpose(modvar) 518 ;;zex = localtime 519 ;;zey = h 520 ;;set_name = 'modvar.ps' 521 ;;set_title = "Horizontal wind speed variance (m!U2!N.s!U-2!N)" 522 ;;set_titlex = 'Local Time (h)' 523 ;;set_titley = 'Altitude above surface (km)' 524 ;;set_subtitle = '' ;'Mean over the simulation domain' 525 ;;set_xrange = [7.,17.] 526 ;;set_yrange = [0.,7.] 527 ;;set_tickx = 1. 528 ;;set_ticky = 1. 529 ;;minval = 0. 530 ;;maxval = 10. 531 ;;nlev = maxval-minval 532 ;;pal = 22 533 ;;rrr = 'no' 534 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 535 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 536 ;;!P.Charsize = 1.2 537 ;;;; 0. levels 538 ;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 539 ;;;;; 1. background 540 ;;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 541 ;;;; 2. color field 542 ;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 543 ;; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 544 ;;;; 3. contour field 545 ;;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) 546 ;;;; 4. choose output 547 ;;if (saveps eq 'true') then PS_End, /PNG 548 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 549 ; 550 ; 551 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 552 ;zefield = transpose(wt) 553 ;zex = localtime 554 ;zey = h 555 ;set_name = 'HF.ps' 556 ;set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'> 557 ;set_titlex = 'Local Time (h)' 602 558 ;set_titley = 'Altitude above surface (km)' 603 559 ;set_subtitle = '' ;'Mean over the simulation domain' 604 ;set_xrange = [190.,230.] ;[min(t),max(t)] 605 ;set_yrange = [0.,7.] 606 ;set_tickx = 5. 607 ;set_ticky = 1. 608 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 609 ;;localtimes = [9,10,11,12,13,14,15,16,17,18] 610 ;localtimes = [7,9,11,13,15,17] 611 ;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] 612 ;;localtimes = [15.0] 560 ;set_xrange = [8.,17.] 561 ;set_yrange = [0.,8.] ;; [0.,200.] 562 ;set_tickx = 1. 563 ;set_ticky = 1. ;; 50. 564 ;minval = -1.5 ;-2. 565 ;maxval = 2.5 ;2. 566 ;nlev = floor(maxval-minval)*10 567 ;pal = 33 568 ;rrr = 'no' 569 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 570 ;;minval = -0.8 571 ;;maxval = 1.2 572 ;;pal = 0 613 573 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 614 574 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 615 575 ;!P.Charsize = 1.2 576 ;;; 0. levels 577 ;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 578 ;;;; 1. background 579 ;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 580 ;;; 2. color field 581 ;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 582 ; ;;;-------------------------------------------------------------------------------------------------------------------------------- 583 ; ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone 584 ; 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 585 ; nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) 586 ; TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b 587 ; ;;;-------------------------------------------------------------------------------------------------------------------------------- 588 ; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 589 ;;; 3. contour field 590 ;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) 591 ;;; 4. choose output 592 ;if (saveps eq 'true') then PS_End, /PNG 593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 594 ; 595 ;;;restore, filename='tpot_profB' 596 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 597 ;;zefield = t 598 ;;zey = ht 599 ;;set_name = 'T.ps' 600 ;;set_title = "Potential Temperature (K)" 601 ;;set_titlex = 'Potential Temperature (K)' 602 ;;set_titley = 'Altitude above surface (km)' 603 ;;set_subtitle = '' ;'Mean over the simulation domain' 604 ;;set_xrange = [190.,230.] ;[min(t),max(t)] 605 ;;set_yrange = [0.,7.] 606 ;;set_tickx = 5. 607 ;;set_ticky = 1. 608 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 609 ;;;localtimes = [9,10,11,12,13,14,15,16,17,18] 610 ;;localtimes = [7,9,11,13,15,17] 611 ;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] 612 ;;;localtimes = [15.0] 613 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 614 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 615 ;;!P.Charsize = 1.2 616 ;;;altlin=0 & loadct, 0 617 ;;altlin=1 & loadct, 0 618 ;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 619 ;;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 620 ;;;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 621 ;;for ll = 1, n_elements(localtimes)-1 do begin 622 ;; CASE altlin OF 623 ;; 0: altlin=1 624 ;; 1: altlin=0 625 ;; ENDCASE 626 ;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 627 ;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 628 ;;; if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin 629 ;;endfor 630 ;;;oplot, ro_tpot, ro_column, psym=7 631 ;;;xyouts, 192.0, 0.25, '09:00' 632 ;;;xyouts, 205.5, 0.25, '11:00' 633 ;;;xyouts, 214.5, 0.25, '13:00' 634 ;;;xyouts, 219.5, 0.25, '17:00' 635 ;;;xyouts, 220.5, 2.30, '17:00 [RO]' 636 ;;if (saveps eq 'true') then PS_End, /PNG 637 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 638 ; 639 ; 640 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 641 ;;zefield = stst 642 ;;zey = ht 643 ;;set_name = 'STST.ps' 644 ;;set_title = "Static stability (K.m!U-1!N)" 645 ;;set_titlex = 'Static stability (K.m!U-1!N)' 646 ;;set_titley = 'Altitude above surface (km)' 647 ;;set_subtitle = 'Mean over the simulation domain' 648 ;;set_xrange = [-5.,5.] 649 ;;set_yrange = [0.,7.] 650 ;;set_tickx = 1. 651 ;;set_ticky = 1. 652 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 653 ;;localtimes = [9,11,13,15,17] 654 ;;localtimes = [12,13,14,15,16] 655 ;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] 656 ;;;localtimes = [15.0] 657 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 658 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 659 ;;!P.Charsize = 1.2 616 660 ;;altlin=0 & loadct, 0 617 ;altlin=1 & loadct, 0 618 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 619 ;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 620 ;;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 621 ;for ll = 1, n_elements(localtimes)-1 do begin 622 ; CASE altlin OF 623 ; 0: altlin=1 624 ; 1: altlin=0 625 ; ENDCASE 626 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 627 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 628 ;; if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin 629 ;endfor 630 ;;oplot, ro_tpot, ro_column, psym=7 631 ;;xyouts, 192.0, 0.25, '09:00' 632 ;;xyouts, 205.5, 0.25, '11:00' 633 ;;xyouts, 214.5, 0.25, '13:00' 634 ;;xyouts, 219.5, 0.25, '17:00' 635 ;;xyouts, 220.5, 2.30, '17:00 [RO]' 636 ;if (saveps eq 'true') then PS_End, /PNG 637 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 638 639 640 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 641 ;zefield = stst 642 ;zey = ht 643 ;set_name = 'STST.ps' 644 ;set_title = "Static stability (K.m!U-1!N)" 645 ;set_titlex = 'Static stability (K.m!U-1!N)' 646 ;set_titley = 'Altitude above surface (km)' 647 ;set_subtitle = 'Mean over the simulation domain' 648 ;set_xrange = [-5.,5.] 649 ;set_yrange = [0.,7.] 650 ;set_tickx = 1. 651 ;set_ticky = 1. 652 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 653 ;localtimes = [9,11,13,15,17] 654 ;localtimes = [12,13,14,15,16] 655 ;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0] 656 ;;localtimes = [15.0] 657 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 658 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 659 ;!P.Charsize = 1.2 660 ;altlin=0 & loadct, 0 661 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 662 ;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 663 ;oplot, zefield(*,nntt), zey(*,nntt), psym=5 664 ;for ll = 1, n_elements(localtimes)-1 do begin 665 ; CASE altlin OF 666 ; 0: altlin=1 667 ; 1: altlin=0 668 ; ENDCASE 669 ; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 670 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 671 ; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5 672 ;endfor 673 ;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 674 ;if (saveps eq 'true') then PS_End, /PNG 675 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 676 677 678 679 if (n_elements(xtke) eq 0) then stop 680 681 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 682 ;;zefield = 2.*transpose(xtke) 683 ;;minval = -4 684 ;;maxval = 8. ;12. 685 ;;pal = 0 686 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 661 ;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 662 ;;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 663 ;;oplot, zefield(*,nntt), zey(*,nntt), psym=5 664 ;;for ll = 1, n_elements(localtimes)-1 do begin 665 ;; CASE altlin OF 666 ;; 0: altlin=1 667 ;; 1: altlin=0 668 ;; ENDCASE 669 ;; user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0) 670 ;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin 671 ;; if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5 672 ;;endfor 673 ;;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 674 ;;if (saveps eq 'true') then PS_End, /PNG 675 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 676 ; 677 ; 678 ; 679 ;if (n_elements(xtke) eq 0) then stop 680 ; 681 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 682 ;;;zefield = 2.*transpose(xtke) 683 ;;;minval = -4 684 ;;;maxval = 8. ;12. 685 ;;;pal = 0 686 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 687 ;;zex = localtime 688 ;;zey = h 689 ;;set_name = 'xTKE.ps' 690 ;;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>] 691 ;;set_titlex = 'Local Time (h)' 692 ;;set_titley = 'Altitude above surface (km)' 693 ;;set_subtitle = '';'Mean over the simulation domain' 694 ;;set_xrange = [8.,17.] 695 ;;set_yrange = [0.,7.] 696 ;;set_tickx = 1. 697 ;;set_ticky = 1. 698 ;;minval = 0. 699 ;;maxval = 8. 700 ;;nlev = maxval-minval 701 ;;pal = 22 702 ;;rrr = 'no' 703 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 704 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 705 ;;!P.Charsize = 1.2 706 ;;;; 0. levels 707 ;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 708 ;;;;; 1. background 709 ;;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 710 ;;;; 2. color field 711 ;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 712 ;; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 713 ;;;; 3. contour field 714 ;;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) 715 ;;;; 4. choose output 716 ;;if (saveps eq 'true') then PS_End, /PNG 717 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 718 ; 719 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 720 ;zefield = transpose(ytke) + transpose(xtke) 687 721 ;zex = localtime 688 722 ;zey = h 689 ;set_name = ' xTKE.ps'690 ;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[< u'!U2!N>]723 ;set_name = 'hTKE.ps' 724 ;set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<v'!U2!N>] 691 725 ;set_titlex = 'Local Time (h)' 692 726 ;set_titley = 'Altitude above surface (km)' 693 ;set_subtitle = '' ;'Mean over the simulation domain'727 ;set_subtitle = '' ;'Mean over the simulation domain' 694 728 ;set_xrange = [8.,17.] 695 ;set_yrange = [0., 7.]729 ;set_yrange = [0.,1.] ;[0.,7.] ;; [0.,200.] 696 730 ;set_tickx = 1. 697 ;set_ticky = 1.731 ;set_ticky = 0.1 ;1. ;50. 698 732 ;minval = 0. 699 ;maxval = 8.733 ;maxval = 12. ;10. 700 734 ;nlev = maxval-minval 701 735 ;pal = 22 … … 716 750 ;if (saveps eq 'true') then PS_End, /PNG 717 751 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 718 719 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 720 zefield = transpose(ytke) + transpose(xtke) 721 zex = localtime 722 zey = h 723 set_name = 'hTKE.ps' 724 set_title = "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<v'!U2!N>] 725 set_titlex = 'Local Time (h)' 726 set_titley = 'Altitude above surface (km)' 727 set_subtitle = '' ;'Mean over the simulation domain' 728 set_xrange = [8.,17.] 729 set_yrange = [0.,1.] ;[0.,7.] ;; [0.,200.] 730 set_tickx = 1. 731 set_ticky = 0.1 ;1. ;50. 732 minval = 0. 733 maxval = 12. ;10. 734 nlev = maxval-minval 735 pal = 22 736 rrr = 'no' 737 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 738 if (saveps eq 'true') then PS_Start, FILENAME=set_name 739 !P.Charsize = 1.2 740 ;; 0. levels 741 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 742 ;;; 1. background 743 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 744 ;; 2. color field 745 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 746 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 747 ;; 3. contour field 748 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) 749 ;; 4. choose output 750 if (saveps eq 'true') then PS_End, /PNG 751 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 752 753 if (n_elements(wmax) eq 0) then stop 752 ; 753 ;if (n_elements(wmax) eq 0) then stop 754 754 755 755 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; … … 762 762 set_titley = 'Altitude above surface (km)' 763 763 set_subtitle = '' 764 set_xrange = [ 8.,18.]765 set_yrange = [0., 7.]764 set_xrange = [11.,17.] ;[8.,18.] 765 set_yrange = [0.,8.] ;[0.,8.] 766 766 set_tickx = 1. 767 767 set_ticky = 1. 768 768 minval = 0. 769 maxval = 1 6. ;15.;14.;13.769 maxval = 18. ;16. ;15.;14.;13. 770 770 nlev = maxval-minval 771 771 pal = 22 … … 776 776 ;; 0. levels 777 777 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 778 levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)] 778 779 ;;; 1. background 779 780 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 783 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 783 784 ;; 3. contour field 784 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 = (levLT 0.0)785 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) 785 786 ;; 4. choose output 786 787 if (saveps eq 'true') then PS_End, /PNG … … 796 797 set_titley = 'Altitude above surface (km)' 797 798 set_subtitle = '' 798 set_xrange = [ 8.,18.]799 set_yrange = [0., 7.]799 set_xrange = [11.,17.] ;[8.,18.] 800 set_yrange = [0.,8.] 800 801 set_tickx = 1. 801 802 set_ticky = 1. 802 803 minval = 0. 803 maxval = 9.;12.804 maxval = 10.;12. 804 805 nlev = maxval-minval 805 806 pal = 22 … … 810 811 ;; 0. levels 811 812 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 813 levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)] 812 814 ;;; 1. background 813 815 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 818 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 817 819 ;; 3. contour field 818 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 = (levLT 0.0)820 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) 819 821 ;; 4. choose output 820 822 if (saveps eq 'true') then PS_End, /PNG 821 823 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 822 823 824 restore, filename='addturb.dat' 825 velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) 826 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 827 zefield = transpose(velmax) 828 zex = localtime 829 zey = h 830 set_name = 'velmax.ps' 831 set_title = "Maximum horizontal wind speed (m.s!U-1!N)" 832 set_titlex = 'Local Time (h)' 833 set_titley = 'Altitude above surface (km)' 834 set_subtitle = '' ;'Mean over the simulation domain' 835 set_xrange = [8.,17.] 836 set_yrange = [0.,4.] 837 set_tickx = 1. 838 set_ticky = 1. 839 minval = 0. 840 maxval = 30. ;12. 841 nlev = maxval-minval 842 pal = 22 843 rrr = 'no' 844 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 845 if (saveps eq 'true') then PS_Start, FILENAME=set_name 846 !P.Charsize = 1.2 847 ;; 0. levels 848 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 849 ;;; 1. background 850 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 851 ;;; 2. color field 852 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 853 contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 854 ;; 3. contour field 855 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) 856 ;; 4. choose output 857 if (saveps eq 'true') then PS_End, /PNG 858 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 859 860 861 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 862 ;zefield = localtime 863 ;zey = depressions 864 ;set_name = 'DEPR.ps' 865 ;set_title = "Percentage of depressions in the domain" 824 ; 825 ; 826 ;restore, filename='addturb.dat' 827 ;velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) 828 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 829 ;zefield = transpose(velmax) 830 ;zex = localtime 831 ;zey = h 832 ;set_name = 'velmax.ps' 833 ;set_title = "Maximum horizontal wind speed (m.s!U-1!N)" 866 834 ;set_titlex = 'Local Time (h)' 867 ;set_titley = ' Percentage'868 ;set_subtitle = '' 869 ;set_xrange = [ 8.,18.]870 ;set_yrange = [0., 3.]835 ;set_titley = 'Altitude above surface (km)' 836 ;set_subtitle = '' ;'Mean over the simulation domain' 837 ;set_xrange = [11.,17.] ;[8.,17.] 838 ;set_yrange = [0.,4.] 871 839 ;set_tickx = 1. 872 840 ;set_ticky = 1. 841 ;minval = 0. 842 ;maxval = 30. ;12. 843 ;nlev = maxval-minval 844 ;pal = 22 845 ;rrr = 'no' 873 846 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 874 847 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name 875 848 ;!P.Charsize = 1.2 876 ;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 877 ;oplot, zefield, zey, psym=5 878 ; oplot, localtime, abs(psmin), linestyle=1 849 ;;; 0. levels 850 ;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] 851 ;;;; 1. background 852 ;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 853 ;;;; 2. color field 854 ;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 855 ; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill 856 ;;; 3. contour field 857 ;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) 858 ;;; 4. choose output 879 859 ;if (saveps eq 'true') then PS_End, /PNG 880 860 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 861 ; 862 ; 863 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 864 ;;zefield = localtime 865 ;;zey = depressions 866 ;;set_name = 'DEPR.ps' 867 ;;set_title = "Percentage of depressions in the domain" 868 ;;set_titlex = 'Local Time (h)' 869 ;;set_titley = 'Percentage' 870 ;;set_subtitle = '' 871 ;;set_xrange = [8.,18.] 872 ;;set_yrange = [0.,3.] 873 ;;set_tickx = 1. 874 ;;set_ticky = 1. 875 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 876 ;;if (saveps eq 'true') then PS_Start, FILENAME=set_name 877 ;;!P.Charsize = 1.2 878 ;;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 879 ;;oplot, zefield, zey, psym=5 880 ;; oplot, localtime, abs(psmin), linestyle=1 881 ;;if (saveps eq 'true') then PS_End, /PNG 882 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 881 883 882 884 stop -
trunk/mesoscale/PLOT/SPEC/LES/report_tasi.inc
r85 r113 1 1 2 saveplot = './' 2 cases = [' MERIDIANI_tau050', 'MERIDIANI_tau100', 'MERIDIANI_tau200', 'MERIDIANI_tau500']3 cases = ['baseline', 'hill', 'little_hill', 'realistic_extreme'] 3 4 paths = cases 4 tttitle = ['tau = 0.5' $ 5 ,'tau = 1.0' $ 6 ,'tau = 2.0' $ 7 ,'tau = 5.0' $ 8 ] 5 tttitle = cases 9 6 tttitles = tttitle 10 7 shorttt = tttitle 8 9 10 ;saveplot = './' 11 ;cases = ['MERIDIANI_tau050', 'MERIDIANI_tau100', 'MERIDIANI_tau200', 'MERIDIANI_tau500'] 12 ;paths = cases 13 ;tttitle = ['tau = 0.5' $ 14 ; ,'tau = 1.0' $ 15 ; ,'tau = 2.0' $ 16 ; ,'tau = 5.0' $ 17 ; ] 18 ;tttitles = tttitle 19 ;shorttt = tttitle
Note: See TracChangeset
for help on using the changeset viewer.
