source: trunk/mesoscale/PLOT/SPEC/LES/compturb.pro @ 93

Last change on this file since 93 was 85, checked in by aslmd, 14 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
File size: 46.0 KB
RevLine 
[85]1pro compturb
2saveps='true' & p0=610. & t0=220. & r_cp=1/4.4 & grav=3.72 & R=192.
3
4;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
5;; grep frames_per_outfile */*/namelist.input
6history_interval_s = 100.
7smoothampl=3700/history_interval_s
8;smoothampl=10.
9;smoothampl=20.
10;smoothampl=4.
11smoothampl=0.
12;smoothampl=5.
13;;;;;;;
14;
15; ATTENTION IL FAUT CHANGER LES VARIABLES A DEUX ENDROITS
16;
17;;;;;;;
18;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
19
20
21;;;
22;;; SET HERE WHAT YOU WANT TO COMPARE
23;;;
24;@article_qjrms.inc
25@report_tasi.inc
26;;;
27;;;
28;;;
29
30
31;;;
32;;; TWEAKS
33;;;
34;SPAWN, 'cp compturb.pro '+saveplot
35;SPAWN, 'gqview '+saveplot+' &'
36;goto, quick
37
38
39;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
40;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
41;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
42
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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
59zefield       =  a_wt_bot
60zey           =  a_h
61set_name      =  saveplot+'comp_A_HFLUX.ps'
62set_title     =  "Similarity: vert. eddy heat flux on Mars"
63set_titlex    =  "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N"
64set_titley    =  'Dimensionless height z/z!Di!N'
65set_subtitle  =  ''
66set_xrange    =  [-0.5,1.]
67set_yrange    =  [0.,1.5]
68set_tickx     =  0.25
69set_ticky     =  0.1
70;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
71zesim = 1 - 1.2 * a_h
72;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
73localtimes = 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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
78if (saveps eq 'true') then PS_Start, FILENAME=set_name
79!P.Charsize = 1.2
80user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
81plot, 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
82for 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
85endfor
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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
93zecolor = 0. & loadct, 33
94for i=1, n_elements(paths)-1 do begin
95;        zecolor = 50. + zecolor
96zecolor = 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
111endfor
112loadct, 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
118if (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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
138zefield       =  a_vel_bot
139zey           =  a_h 
140set_name      =  saveplot+'comp_A_VEL.ps'
141set_title     =  "Similarity: vert. vel. variance on Mars"
142set_titlex    =  "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N"
143set_titley    =  'Dimensionless height z/z!Di!N'
144set_subtitle  =  ''
145set_xrange    =  [0.,1.2]
146set_yrange    =  [0.,1.5]
147set_tickx     =  0.25 ;0.1
148set_ticky     =  0.1
149;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
150zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
151;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
152localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
153;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
154;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
155;com
156if (saveps eq 'true') then PS_Start, FILENAME=set_name
157!P.Charsize = 1.2
158user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
159plot, 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
160for 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
163endfor
164;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
165;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
166zecolor = 0. & loadct, 33
167for i=1, n_elements(paths)-1 do begin
168zecolor = 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
183endfor
184loadct, 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
189if (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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
205zefield       =  localtime
206zey           =  SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE)
207set_name      =  saveplot+'comp_W_STAR.ps'
208set_title     =  '';"Free convection velocity scale W!D*!N (m.s!U-1!N)"
209set_titlex    =  'Local Time (h)'
210set_titley    =  'Free convection velocity scale (m s!U-1!N)'
211set_subtitle  =  ''
212set_xrange    =  [08.,17.]
213set_yrange    =  [00.,06.]
214set_tickx     =  1.0
215set_ticky     =  0.5
216wherex        =  8.5;11.0
217wherey        =  5.5;1.75
218;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
219fixfix = 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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
221if (saveps eq 'true') then PS_Start, FILENAME=set_name
222!P.Charsize = 1.2
223altlin=0
224plot, 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
225dddx = (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
231for 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
249endfor
250loadct, 0
251;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
252;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
253if (saveps eq 'true') then PS_End, /PNG
254
255;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
256;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
257;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
258
259        path = paths(0) & print, path
260        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
261        file = path+'/compturb.dat' & restore, filename=file
262        file = path+'/getturb.dat'  & restore, filename=file
263        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
264
265;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
266zefield       =  localtime
267zey           =  SMOOTH(pbl_height, [smoothampl], /EDGE_TRUNCATE)
268set_name      =  saveplot+'comp_HEIGHT.ps'
269set_title     =  '';"Boundary layer height z!Di!N (km)"; [definition : static stability > 1.5 K.m!U-1!N] "
270set_titlex    =  'Local Time (h)'
271set_titley    =  'Boundary layer height (km)'
272set_subtitle  =  ''
273;set_xrange    =  [09.,18.]
274set_yrange    =  [00.,07.]
275set_tickx     =  1.
276set_ticky     =  1.
277wherex        =  09.;13.25 ;12.75
278wherey        =  06.;2.5 ;2.0
279;       set_xrange    =  [13.5,18.] & set_yrange    =  [06.4,08.0] & set_tickx     =  0.5 & set_ticky     =  0.2
280;       wherex=15.30
281;       wherey=06.90
282;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
283if (saveps eq 'true') then PS_Start, FILENAME=set_name
284!P.Charsize = 1.2
285altlin=0
286plot, 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
287dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
288        plots,  wherex + 2.*dddx/1.1,  wherey
289        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
290;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
291;w = where((localtime ge 17.) and (localtime le 17.5)) & if (w(0) ne -1) then print, 'BL height 17:00-17:30 : ', mean(pbl_height[w])
292w = where((localtime ge 14.) and (localtime le 14.5)) & if (w(0) ne -1) then print, 'BL height 14:00-14:30 : ', mean(pbl_height[w])
293;print, max(SMOOTH(tke, [0,37], /EDGE_TRUNCATE))
294;print, max(SMOOTH(wmax, [0,37], /EDGE_TRUNCATE))
295;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
296;zecolor = 0. & loadct, 33
297for i=1, n_elements(paths)-1 do begin
298        ;zecolor = float(255)/n_elements(paths) + zecolor
299        zecolor=0
300        path = paths(i) & print, path
301        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
302        file = path+'/compturb.dat' & restore, filename=file
303        file = path+'/getturb.dat'  & restore, filename=file
304        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
305        zefield = localtime
306        zey = SMOOTH(pbl_height, [smoothampl], /EDGE_TRUNCATE);pbl_height
307        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
308        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
309        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
310        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
311        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
312        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
313        ;w = where((localtime ge 17.) and (localtime le 17.5)) & if (w(0) ne -1) then print, 'BL height 17:00-17:30 : ', mean(pbl_height[w])
314        w = where((localtime ge 14.) and (localtime le 14.5)) & if (w(0) ne -1) then print, 'BL height 14:00-14:30 : ', mean(pbl_height[w])
315        ;print, max(SMOOTH(tke, [0,37], /EDGE_TRUNCATE))
316        ;print, max(SMOOTH(wmax, [0,37], /EDGE_TRUNCATE))
317endfor
318loadct, 0
319;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
320;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
321if (saveps eq 'true') then PS_End, /PNG
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
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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
391zefield       =  localtime
392zey           =  SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
393zey2          =  SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
394set_name      =  saveplot+'comp_WMAX.ps'
395set_title     =  '';"Maximum and minimum vertical velocity (m.s!U-1!N)"
396set_titlex    =  'Local Time (h)'
397set_titley    =  'Maximum and minimum vertical velocity (m s!U-1!N)'
398set_subtitle  =  ''
399;set_xrange    =  [09.0,18.0]
400set_yrange    =  [-12.0,18.0];[-15.0,20.0]
401set_tickx     =  1.
402set_ticky     =  2.
403wherex        =  8.5;11.00
404wherey        =  16.5;5.5
405;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
406if (saveps eq 'true') then PS_Start, FILENAME=set_name
407!P.Charsize = 1.2
408altlin = 0
409plot, 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
410oplot, zefield, zey2, color=0, linestyle=altlin
411dddx = (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
417for 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        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
435endfor
436loadct, 0
437;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
438;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
439if (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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
452zefield       =  localtime
453zey           =  SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
454set_name      =  saveplot+'comp_TKE.ps'
455set_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)"
456set_titlex    =  'Local Time (h)'
457set_titley    =  'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)'
458set_subtitle  =  ''
459;set_xrange    =  [09.0,18.0]
460set_yrange    =  [00.0,18.0];[00.0,26.0]
461set_tickx     =  1.
462set_ticky     =  2.
463wherex        =  8.5;11.70
464wherey        =  16.;5.4
465;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
466if (saveps eq 'true') then PS_Start, FILENAME=set_name
467!P.Charsize = 1.2
468altlin = 0
469plot, 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
470dddx = (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
476for 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
493endfor
494loadct, 0
495;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
496;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
497if (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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
511zefield       =  localtime
512zey           =  SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
513set_name      =  saveplot+'comp_HFBOT.ps'
514set_title     =  '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)"
515set_titlex    =  'Local Time (h)'
516set_titley    =  'Maximum vertical eddy heat flux (K m s!U-1!N)'
517set_subtitle  =  ''
518;set_xrange    =  [09.0,18.0]
519set_yrange    =  [00.0,02.0];[00.0,02.8]
520set_tickx     =  1.
521set_ticky     =  0.2
522wherex        =  8.5;10.30
523wherey        =  1.8;0.7
524;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
525if (saveps eq 'true') then PS_Start, FILENAME=set_name
526!P.Charsize = 1.2
527altlin = 0
528plot, 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
529dddx = (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
535for 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
552endfor
553loadct, 0
554;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
555;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
556if (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;;;;;;;;;;;;;;;
683lev = 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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
694zefield       =  localtime
695zey           =  SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
696set_name      =  saveplot+'comp_VELMAX.ps'
697set_title     =  '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)"
698set_titlex    =  'Local Time (h)'
699set_titley    =  "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)"
700set_subtitle  =  ''
701;set_xrange    =  [09.0,18.0]
702set_yrange    =  [00.0,14.0]
703set_tickx     =  1.
704set_ticky     =  2.
705wherex        =  8.5;10.30
706wherey        =  12.;0.7
707;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
708print, h(lev)
709;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
710if (saveps eq 'true') then PS_Start, FILENAME=set_name
711!P.Charsize = 1.2
712altlin = 0
713plot, 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
714dddx = (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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
718;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
719;zecolor = 0. & loadct, 33
720for i=1, n_elements(paths)-1 do begin
721        ;zecolor = float(255)/n_elements(paths) + zecolor
722        zecolor = 0.
723        path = paths(i) & print, path
724        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
725        file = path+'/compturb.dat' & restore, filename=file
726        file = path+'/addturb.dat'  & restore, filename=file
727        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
728        zefield = localtime
729        zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
730        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
731        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
732        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
733        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
734        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
735        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
736
737endfor
738loadct, 0
739;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
740;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
741if (saveps eq 'true') then PS_End, /PNG
742
743
744
745
746
747
748;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
749;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
750;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
751
752
753goto, no_therm
754quick:
755;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
756cp = R / r_cp
757;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
758;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
759x  = fltarr(n_elements(paths)) & y  = fltarr(n_elements(paths))
760x2 = fltarr(n_elements(paths)) & y2 = fltarr(n_elements(paths))
761for i=0, n_elements(paths)-1 do begin
762        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
763        path = paths(i) & print, path
764        file = path+'/compturb.dat' & restore, filename=file
765        file = path+'/getturb.dat'  & restore, filename=file
766        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
767        yeye = 17.
768        ;yeye = 15.
769        ;yeye = 13.
770        ;yeye = 09.
771        w = where(localtime eq yeye) ;& print, w
772
773
774                ;yeye = 08.
775                ;w2 = where(localtime eq yeye) ;& print, w
776        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
777        x(i) = (p0 / p_bot[w])^r_cp
778        y(i) = t_bot[w];-t_bot[w2] & print, y(i)
779        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
780        y2(i) = z_top[w]
781        ;y2(i) = pbl_height[w]
782        x2(i) = t_bot[w]
783        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
784
785;w2 = where( abs(h-0.01) eq min(abs(h-0.01)) )
786;w2(0)=0
787;w2 = where( abs(h-0.1) eq min(abs(h-0.1)) )
788;w2 = where( abs(h-1.5) eq min(abs(h-1.5)) )
789;w2 = where( abs(h-1.) eq min(abs(h-1.)) )
790;print, w2(0), w(0)
791;x(i) = (p0 / pt[w2(0),w(0)])^r_cp
792;y(i) = t[w2(0),w(0)]
793;print, h[w2]
794;stop
795
796;ye=smooth(DERIV(localtime, smooth(t_bot,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
797;ye=DERIV(localtime, t_bot)
798;w3 = where(localtime eq 13.) ;;11. flat!
799;y(i) = ye[w3]
800
801endfor
802print, x2
803print, y2
804fit  = LINFIT(x ,y )
805print, fit
806fit2 = LINFIT(x2,y2)
807;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
808;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
809;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
810file='./data_hinson' & header='' & nlines_header=0 & ncol = 7
811;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
812nlines = FILE_LINES(file)-nlines_header & data=FLTARR(ncol,nlines)
813OPENR, lun, file, /GET_LUN & READF, lun, data & CLOSE, lun & FREE_LUN, lun
814;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
815alt = reform(data(5,*)) - reform(data(3,*))             ;; altitude
816xdata = (exp(alt/10.))^r_cp                             ;; p/p0 with assumed scale height
817ydata = reform(data(6,*))                               ;; potential temperature
818;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
819ydata2 = reform(data(5,*));reform(data(3,*));reform(data(5,*))                                  ;; altitude of the top of the boundary layer zt
820xdata2 = reform(data(6,*))                              ;; potential temperature
821;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
822data = 0.
823fitdata  = LINFIT(xdata ,ydata )
824fitdata2 = LINFIT(xdata2,ydata2)
825;fitdata2 = POLY_FIT(xdata2,ydata2,2)
826;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
827;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
828;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
829
830;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
831set_name      =  saveplot+'comp_EXN.ps'
832set_title     =  ""
833set_titlex    =  'Base Exner function P!U-1!N=(p!D0!N/p!Db!N)!UR/cp!N'
834set_titley    =  'Base potential temperature T!Db!N (K)'
835set_subtitle  =  ''
836set_xrange    =  [0.90,1.25]
837set_yrange    =  [210.,265.]
838;set_xrange    =  [0.90,1.10]
839;set_yrange    =  [210.,250.]
840        ;set_xrange    =  [0.80,1.20]
841        ;set_yrange    =  [2.,5.]
842        ;set_xrange    =  [0.80,1.20]
843        ;set_yrange    =  [180.,230.]
844set_tickx     =  0.05
845set_ticky     =  5.
846wherex        =  1.00
847wherey        =  220.   
848;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
849if (saveps eq 'true') then PS_Start, FILENAME=set_name
850!P.Charsize = 1.2
851altlin=0
852dddy = (set_yrange(1) - set_yrange(0))/20.
853dddx = (set_xrange(1) - set_xrange(0))/20.
854
855!p.multi=[0,2,1]
856
857set_title='LES results. T!Db!N = '+string(fit(1), '(F5.1)')+' P!U-1!N + '+string(fit(0), '(F5.1)')
858                plot, $
859set_xrange, fit[1]*set_xrange + fit[0], linestyle=0, $
860                xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, $
861                ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, $
862                 title=set_title, subtitle=set_subtitle, color=0
863for i=0, n_elements(paths)-1 do xyouts, x(i), y(i), shorttt(i)
864;               xyouts, $
865;wherex, wherey, 'squares : data. T!Db!N = '+string(fitdata(1), '(F6.2)')+' P!U-1!N + '+string(fitdata(0), '(F6.2)')
866;                xyouts, $
867;wherex, wherey-dddy, 'triangles : LES. T!Db!N = '+string(fit(1), '(F6.2)')+' P!U-1!N + '+string(fit(0), '(F6.2)')
868
869set_title='RO data. T!Db!N = '+string(fitdata(1), '(F5.1)')+' P!U-1!N + '+string(fitdata(0), '(F5.1)')
870                plot, $
871set_xrange, fitdata[1]*set_xrange + fitdata[0], $
872                xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, $
873                ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, $
874                 title=set_title, subtitle=set_subtitle, color=0, linestyle=1
875                oplot, $
876xdata, ydata, psym=6
877;                oplot, $
878;x, y, psym=5
879
880;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
881;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
882if (saveps eq 'true') then PS_End, /PNG
883
884
885;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
886fit     = fit2
887x       = x2
888y       = y2
889fitdata = fitdata2
890xdata   = xdata2
891ydata   = ydata2 
892;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
893set_name      =  saveplot+'comp_ZTOP.ps'
894set_title     =  "Absolute altitude of the BL top z!Dt!N (km)"
895set_titlex    =  'Base potential temperature T!Db!N (K)'
896set_titley    =  'Absolute altitude of the BL top (km)'
897set_subtitle  =  ''
898set_xrange    =  [150.,280.] ;[200.,265.]
899set_yrange    =  [-8.,25.]   ;[-4.,20.]
900set_tickx     =  10.
901set_ticky     =  2.
902wherex        =  205.
903wherey        =  17.0
904;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
905if (saveps eq 'true') then PS_Start, FILENAME=set_name
906!P.Charsize = 1.2
907altlin=0
908dddy = (set_yrange(1) - set_yrange(0))/20.
909dddx = (set_xrange(1) - set_xrange(0))/20.
910;yeye = set_xrange[0] + (set_xrange[1]-set_xrange[0])*findgen(100.)/99.
911                plot, $
912set_xrange, fitdata[1]*set_xrange + fitdata[0], $
913;yeye, fitdata[2]*yeye^2 + fitdata[1]*yeye + fitdata[0], $
914                xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, $
915                ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, $
916                 title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
917                oplot, $
918xdata, ydata, psym=6
919                oplot, $
920x, y, psym=5
921                oplot, $
922set_xrange, fit[1]*set_xrange + fit[0], linestyle=1
923                xyouts, $
924wherex, wherey, 'squares : data.  z!Dt!N  =  ('+strtrim(string(fitdata(1)),1)+')  T!Db!N + ('+strtrim(string(fitdata(0)),1)+')'
925                xyouts, $
926wherex, wherey-dddy, 'triangles : model. z!Dt!N  =  ('+strtrim(string(fit(1)),1)+')  T!Db!N + ('+strtrim(string(fit(0)),1)+')'
927;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
928;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
929if (saveps eq 'true') then PS_End, /PNG
930
931
932;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
933;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
934;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
935
936
937no_therm:
938
939end
Note: See TracBrowser for help on using the repository browser.