source: trunk/UTIL/IDLplot/SPEC/LES/compturb.pro @ 1306

Last change on this file since 1306 was 113, checked in by aslmd, 15 years ago

LMD_LES_MARS: qq modifications pour les routines graphiques
File size: 46.7 KB
Line 
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.
12smoothampl=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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
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;
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
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  =  ''
273set_xrange    =  [08.,18.]
274set_yrange    =  [00.,08.] ;[00.,07.]
275set_tickx     =  1.
276set_ticky     =  1.
277wherex        =  09.;13.25 ;12.75
278wherey        =  06.5;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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
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;
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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
718;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
719;;zecolor = 0. & loadct, 33
720;for 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;
737;endfor
738;loadct, 0
739;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
740;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
741;if (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.