Index: trunk/mesoscale/PLOT/SPEC/LES/compturb.pro
===================================================================
--- trunk/mesoscale/PLOT/SPEC/LES/compturb.pro	(revision 85)
+++ trunk/mesoscale/PLOT/SPEC/LES/compturb.pro	(revision 113)
@@ -10,5 +10,5 @@
 ;smoothampl=4.
 smoothampl=0.
-;smoothampl=5.
+smoothampl=5.
 ;;;;;;;
 ;
@@ -41,219 +41,219 @@
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 
-;; pour article, decommenter ici et commenter l'appel au deuxieme fichier
+;;; pour article, decommenter ici et commenter l'appel au deuxieme fichier
+;;;com
+;;!p.multi = [0,2,1]
+;
+;        path = paths(0) & print, path 
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;;;
+;;; COEFFICIENT POUR LOIS DE SIMILITUDE
+;;;
+;;@coeffhf
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zefield       =  a_wt_bot 
+;zey           =  a_h
+;set_name      =  saveplot+'comp_A_HFLUX.ps'
+;set_title     =  "Similarity: vert. eddy heat flux on Mars"
+;set_titlex    =  "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N"
+;set_titley    =  'Dimensionless height z/z!Di!N'
+;set_subtitle  =  ''
+;set_xrange    =  [-0.5,1.]
+;set_yrange    =  [0.,1.5]
+;set_tickx     =  0.25
+;set_ticky     =  0.1
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zesim = 1 - 1.2 * a_h
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
+;;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )]
+;;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )]
+;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;!P.Charsize = 1.2
+;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin
+;for ll = 1, n_elements(localtimes)-1 do begin
+;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
+;endfor
+;;wherex = set_xrange(0)
+;;wherey = set_yrange(1)
+;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zecolor = 0. & loadct, 33
+;for i=1, n_elements(paths)-1 do begin
+;;        zecolor = 50. + zecolor
+;zecolor = 0. & loadct, 0
+;        path = paths(i) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        zefield = a_wt_bot
+;        zey = a_h
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        for ll = 0, n_elements(localtimes)-1 do begin
+;          user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;          if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
+;        endfor
+;        ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;endfor
+;loadct, 0
+;;oplot, [0.,0.], set_yrange, linestyle=0, color=0
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
+;	;oplot, set_xrange, [0.,0.], linestyle=0, color=0
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_End, /PNG
 ;;com
-;!p.multi = [0,2,1]
-
-        path = paths(0) & print, path 
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;
-;; COEFFICIENT POUR LOIS DE SIMILITUDE
-;;
-;@coeffhf
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  a_wt_bot 
-zey           =  a_h
-set_name      =  saveplot+'comp_A_HFLUX.ps'
-set_title     =  "Similarity: vert. eddy heat flux on Mars"
-set_titlex    =  "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N"
-set_titley    =  'Dimensionless height z/z!Di!N'
-set_subtitle  =  ''
-set_xrange    =  [-0.5,1.]
-set_yrange    =  [0.,1.5]
-set_tickx     =  0.25
-set_ticky     =  0.1
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zesim = 1 - 1.2 * a_h
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
-;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )]
-;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )]
-;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name
-!P.Charsize = 1.2
-user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin
-for ll = 1, n_elements(localtimes)-1 do begin
-  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
-endfor
-;wherex = set_xrange(0)
-;wherey = set_yrange(1)
-;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-;        zecolor = 50. + zecolor
-zecolor = 0. & loadct, 0
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = a_wt_bot
-        zey = a_h
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        for ll = 0, n_elements(localtimes)-1 do begin
-          user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-          if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
-        endfor
-        ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-endfor
-loadct, 0
-;oplot, [0.,0.], set_yrange, linestyle=0, color=0
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
-	;oplot, set_xrange, [0.,0.], linestyle=0, color=0
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
-;com
-
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-	path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;
-;; COEFFICIENT POUR LOIS DE SIMILITUDE
-;;
-;@coeffw
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  a_vel_bot 
-zey           =  a_h  
-set_name      =  saveplot+'comp_A_VEL.ps'
-set_title     =  "Similarity: vert. vel. variance on Mars"
-set_titlex    =  "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N"
-set_titley    =  'Dimensionless height z/z!Di!N'
-set_subtitle  =  ''
-set_xrange    =  [0.,1.2]
-set_yrange    =  [0.,1.5]
-set_tickx     =  0.25 ;0.1
-set_ticky     =  0.1
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
-;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;com
-if (saveps eq 'true') then PS_Start, FILENAME=set_name
-!P.Charsize = 1.2
-user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin
-for ll = 1, n_elements(localtimes)-1 do begin
-  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
-endfor
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-zecolor = 0. & loadct, 0
-;	zecolor = 50. + zecolor
-	path = paths(i) & print, path
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	file = path+'/compturb.dat' & restore, filename=file
-	file = path+'/getturb.dat'  & restore, filename=file
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	zefield = a_vel_bot 
-	zey = a_h
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	for ll = 0, n_elements(localtimes)-1 do begin
-	  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
-	  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
-	endfor
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-endfor
-loadct, 0
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
-	;oplot, set_xrange, [0.,0.], linestyle=0, color=0
-	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
-
+;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;	path = paths(0) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;;;
+;;; COEFFICIENT POUR LOIS DE SIMILITUDE
+;;;
+;;@coeffw
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zefield       =  a_vel_bot 
+;zey           =  a_h  
+;set_name      =  saveplot+'comp_A_VEL.ps'
+;set_title     =  "Similarity: vert. vel. variance on Mars"
+;set_titlex    =  "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N"
+;set_titley    =  'Dimensionless height z/z!Di!N'
+;set_subtitle  =  ''
+;set_xrange    =  [0.,1.2]
+;set_yrange    =  [0.,1.5]
+;set_tickx     =  0.25 ;0.1
+;set_ticky     =  0.1
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
+;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;com
-;stop
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-        path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  localtime
-zey           =  SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE)
-set_name      =  saveplot+'comp_W_STAR.ps'
-set_title     =  '';"Free convection velocity scale W!D*!N (m.s!U-1!N)"
-set_titlex    =  'Local Time (h)'
-set_titley    =  'Free convection velocity scale (m s!U-1!N)'
-set_subtitle  =  ''
-set_xrange    =  [08.,17.]
-set_yrange    =  [00.,06.]
-set_tickx     =  1.0
-set_ticky     =  0.5
-wherex	      =  8.5;11.0
-wherey	      =  5.5;1.75
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name
-!P.Charsize = 1.2
-altlin=0
-plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path
-        plots,  wherex + 2.*dddx/1.1,  wherey
-        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;!P.Charsize = 1.2
+;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin
+;for ll = 1, n_elements(localtimes)-1 do begin
+;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
+;endfor
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-        ;zecolor = float(255)/n_elements(paths) + zecolor
-	zecolor = 0
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = localtime
-        zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-	fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
-        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-endfor
-loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;for i=1, n_elements(paths)-1 do begin
+;zecolor = 0. & loadct, 0
+;;	zecolor = 50. + zecolor
+;	path = paths(i) & print, path
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	file = path+'/compturb.dat' & restore, filename=file
+;	file = path+'/getturb.dat'  & restore, filename=file
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	zefield = a_vel_bot 
+;	zey = a_h
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	for ll = 0, n_elements(localtimes)-1 do begin
+;	  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
+;	  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
+;	endfor
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;endfor
+;loadct, 0
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
+;	;oplot, set_xrange, [0.,0.], linestyle=0, color=0
+;	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_End, /PNG
+;
+;;;com
+;;stop
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;        path = paths(0) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zefield       =  localtime
+;zey           =  SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE)
+;set_name      =  saveplot+'comp_W_STAR.ps'
+;set_title     =  '';"Free convection velocity scale W!D*!N (m.s!U-1!N)"
+;set_titlex    =  'Local Time (h)'
+;set_titley    =  'Free convection velocity scale (m s!U-1!N)'
+;set_subtitle  =  ''
+;set_xrange    =  [08.,17.]
+;set_yrange    =  [00.,06.]
+;set_tickx     =  1.0
+;set_ticky     =  0.5
+;wherex	      =  8.5;11.0
+;wherey	      =  5.5;1.75
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;!P.Charsize = 1.2
+;altlin=0
+;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path
+;        plots,  wherex + 2.*dddx/1.1,  wherey
+;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zecolor = 0. & loadct, 33
+;for i=1, n_elements(paths)-1 do begin
+;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;	zecolor = 0
+;        path = paths(i) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        zefield = localtime
+;        zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;	fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
+;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;endfor
+;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_End, /PNG
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 
         path = paths(0) & print, path
@@ -271,10 +271,10 @@
 set_titley    =  'Boundary layer height (km)'
 set_subtitle  =  ''
-;set_xrange    =  [09.,18.]
-set_yrange    =  [00.,07.]
+set_xrange    =  [08.,18.]
+set_yrange    =  [00.,08.] ;[00.,07.]
 set_tickx     =  1. 
 set_ticky     =  1.
 wherex        =  09.;13.25 ;12.75
-wherey        =  06.;2.5 ;2.0
+wherey        =  06.5;2.5 ;2.0
 ;	set_xrange    =  [13.5,18.] & set_yrange    =  [06.4,08.0] & set_tickx     =  0.5 & set_ticky     =  0.2 
 ;       wherex=15.30
@@ -321,4 +321,61 @@
 if (saveps eq 'true') then PS_End, /PNG
 
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;        path = paths(0) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zefield       =  localtime
+;;zey           =  smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)  
+;;set_name      =  saveplot+'comp_dHEIGHT.ps'
+;;set_title     =  "Boundary layer height growth dz!Di!N/dt (km/hour)"
+;;set_titlex    =  'Local Time (h)'
+;;set_titley    =  'Boundary layer height growth (km/hour)'
+;;set_subtitle  =  ''
+;;set_xrange    =  [9.,18.]
+;;set_yrange    =  [-0.4,1.6] 
+;;set_tickx     =  1.
+;;set_ticky     =  0.2
+;;wherex        =  10.30
+;;wherey        =  0.1
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;;!P.Charsize = 1.2
+;;altlin = 0
+;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
+;;        plots,  wherex + 2.*dddx/1.1,  wherey
+;;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;zecolor = 0. & loadct, 33
+;;for i=1, n_elements(paths)-1 do begin
+;;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;;	zecolor = 0
+;;        path = paths(i) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        zefield = localtime
+;;        zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
+;;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;endfor
+;;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_End, /PNG
+;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
@@ -333,16 +390,17 @@
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;zefield       =  localtime
-;zey           =  smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)  
-;set_name      =  saveplot+'comp_dHEIGHT.ps'
-;set_title     =  "Boundary layer height growth dz!Di!N/dt (km/hour)"
+;zey           =  SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
+;zey2          =  SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
+;set_name      =  saveplot+'comp_WMAX.ps'
+;set_title     =  '';"Maximum and minimum vertical velocity (m.s!U-1!N)"
 ;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Boundary layer height growth (km/hour)'
+;set_titley    =  'Maximum and minimum vertical velocity (m s!U-1!N)'
 ;set_subtitle  =  ''
-;set_xrange    =  [9.,18.]
-;set_yrange    =  [-0.4,1.6] 
+;;set_xrange    =  [09.0,18.0]
+;set_yrange    =  [-12.0,18.0];[-15.0,20.0]
 ;set_tickx     =  1.
-;set_ticky     =  0.2
-;wherex        =  10.30
-;wherey        =  0.1
+;set_ticky     =  2.
+;wherex        =  8.5;11.00
+;wherey        =  16.5;5.5
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name
@@ -350,13 +408,14 @@
 ;altlin = 0
 ;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
+;oplot, zefield, zey2, color=0, linestyle=altlin
+;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
 ;        plots,  wherex + 2.*dddx/1.1,  wherey
 ;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;zecolor = 0. & loadct, 33
 ;for i=1, n_elements(paths)-1 do begin
 ;        ;zecolor = float(255)/n_elements(paths) + zecolor
-;	zecolor = 0
+;        zecolor = 0.
 ;        path = paths(i) & print, path
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
@@ -365,8 +424,10 @@
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;        zefield = localtime
-;        zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 
+;        zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
+;	zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
+;        oplot, zefield, zey2, color=zecolor, linestyle=altlin+i
+;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
 ;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
 ;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
@@ -377,184 +438,5 @@
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;if (saveps eq 'true') then PS_End, /PNG
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-        path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  localtime
-zey           =  SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
-zey2          =  SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
-set_name      =  saveplot+'comp_WMAX.ps'
-set_title     =  '';"Maximum and minimum vertical velocity (m.s!U-1!N)"
-set_titlex    =  'Local Time (h)'
-set_titley    =  'Maximum and minimum vertical velocity (m s!U-1!N)'
-set_subtitle  =  ''
-;set_xrange    =  [09.0,18.0]
-set_yrange    =  [-12.0,18.0];[-15.0,20.0]
-set_tickx     =  1.
-set_ticky     =  2.
-wherex        =  8.5;11.00
-wherey        =  16.5;5.5
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name
-!P.Charsize = 1.2
-altlin = 0
-plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-oplot, zefield, zey2, color=0, linestyle=altlin
-dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
-        plots,  wherex + 2.*dddx/1.1,  wherey
-        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-        ;zecolor = float(255)/n_elements(paths) + zecolor
-        zecolor = 0.
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = localtime
-        zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
-	zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE) 
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-        oplot, zefield, zey2, color=zecolor, linestyle=altlin+i
-        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
-        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-endfor
-loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-        path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  localtime
-zey           =  SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
-set_name      =  saveplot+'comp_TKE.ps'
-set_title     =  '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)"
-set_titlex    =  'Local Time (h)'
-set_titley    =  'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)'
-set_subtitle  =  ''
-;set_xrange    =  [09.0,18.0]
-set_yrange    =  [00.0,18.0];[00.0,26.0]
-set_tickx     =  1.
-set_ticky     =  2.
-wherex        =  8.5;11.70
-wherey        =  16.;5.4
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name
-!P.Charsize = 1.2
-altlin = 0
-plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0)
-        plots,  wherex + 2.*dddx/1.1,  wherey
-        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-        ;zecolor = float(255)/n_elements(paths) + zecolor
-        zecolor = 0.
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = localtime
-        zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitles(i), color=zecolor
-        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-endfor
-loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
-
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-        path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  localtime
-zey           =  SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
-set_name      =  saveplot+'comp_HFBOT.ps'
-set_title     =  '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)"
-set_titlex    =  'Local Time (h)'
-set_titley    =  'Maximum vertical eddy heat flux (K m s!U-1!N)'
-set_subtitle  =  ''
-;set_xrange    =  [09.0,18.0]
-set_yrange    =  [00.0,02.0];[00.0,02.8]
-set_tickx     =  1.
-set_ticky     =  0.2
-wherex        =  8.5;10.30
-wherey        =  1.8;0.7
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name 
-!P.Charsize = 1.2
-altlin = 0
-plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
-        plots,  wherex + 2.*dddx/1.1,  wherey
-        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-        ;zecolor = float(255)/n_elements(paths) + zecolor
-	zecolor = 0.
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/getturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = localtime
-        zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
-        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-endfor
-loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG 
-
-
+;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
@@ -569,20 +451,77 @@
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;zefield       =  localtime
-;zey           =  t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
-;zey	      =  SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
-;;zey           =  SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
-;set_name      =  saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps'
-;set_title     =  "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)"
+;zey           =  SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
+;set_name      =  saveplot+'comp_TKE.ps'
+;set_title     =  '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)"
 ;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Potential temperature time derivative (K.hour!U-1!N)'
+;set_titley    =  'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)'
 ;set_subtitle  =  ''
-;set_xrange    =  [ 09., 18.] 
-;set_yrange    =  [-01., 06.] 
+;;set_xrange    =  [09.0,18.0]
+;set_yrange    =  [00.0,18.0];[00.0,26.0]
 ;set_tickx     =  1.
-;set_ticky     =  1. 
-;wherex        =  13.25 
-;wherey        =  5.5
+;set_ticky     =  2.
+;wherex        =  8.5;11.70
+;wherey        =  16.;5.4
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;!P.Charsize = 1.2
+;altlin = 0
+;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0)
+;        plots,  wherex + 2.*dddx/1.1,  wherey
+;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zecolor = 0. & loadct, 33
+;for i=1, n_elements(paths)-1 do begin
+;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;        zecolor = 0.
+;        path = paths(i) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        zefield = localtime
+;        zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitles(i), color=zecolor
+;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;endfor
+;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_End, /PNG
+;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;        path = paths(0) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zefield       =  localtime
+;zey           =  SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
+;set_name      =  saveplot+'comp_HFBOT.ps'
+;set_title     =  '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)"
+;set_titlex    =  'Local Time (h)'
+;set_titley    =  'Maximum vertical eddy heat flux (K m s!U-1!N)'
+;set_subtitle  =  ''
+;;set_xrange    =  [09.0,18.0]
+;set_yrange    =  [00.0,02.0];[00.0,02.8]
+;set_tickx     =  1.
+;set_ticky     =  0.2
+;wherex        =  8.5;10.30
+;wherey        =  1.8;0.7
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_Start, FILENAME=set_name 
 ;!P.Charsize = 1.2
 ;altlin = 0
@@ -591,5 +530,190 @@
 ;        plots,  wherex + 2.*dddx/1.1,  wherey
 ;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zecolor = 0. & loadct, 33
+;for i=1, n_elements(paths)-1 do begin
+;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;	zecolor = 0.
+;        path = paths(i) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/getturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        zefield = localtime
+;        zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax 
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
+;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;endfor
+;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_End, /PNG 
+;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;        path = paths(0) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zefield       =  localtime
+;;zey           =  t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
+;;zey	      =  SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
+;;;zey           =  SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
+;;set_name      =  saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps'
+;;set_title     =  "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)"
+;;set_titlex    =  'Local Time (h)'
+;;set_titley    =  'Potential temperature time derivative (K.hour!U-1!N)'
+;;set_subtitle  =  ''
+;;set_xrange    =  [ 09., 18.] 
+;;set_yrange    =  [-01., 06.] 
+;;set_tickx     =  1.
+;;set_ticky     =  1. 
+;;wherex        =  13.25 
+;;wherey        =  5.5
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;;!P.Charsize = 1.2
+;;altlin = 0
+;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
+;;        plots,  wherex + 2.*dddx/1.1,  wherey
+;;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;zecolor = 0. & loadct, 33
+;;for i=1, n_elements(paths)-1 do begin
+;;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;;        zecolor = 0.
+;;        path = paths(i) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        zefield = localtime
+;;        zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 
+;;		;zey = DERIV(localtime, t_bot)
+;;		zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
+;;		;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
+;;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;endfor
+;;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;oplot, set_xrange, [0,0], linestyle=2
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_End, /PNG
+;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;        path = paths(0) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;zefield       =  localtime
+;;zey           =  SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top
+;;set_name      =  saveplot+'comp_TTOP.ps'
+;;set_title     =  "Temperature at the top of the boundary layer (K)"
+;;set_titlex    =  'Local Time (h)'
+;;set_titley    =  'Temperature at the top of the boundary layer (K)'
+;;set_subtitle  =  ''
+;;set_xrange    =  [09.,18.]
+;;set_yrange    =  [195.,225.]
+;;set_tickx     =  1.
+;;set_ticky     =  5.
+;;wherex        =  13.00
+;;wherey        =  223.5
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
+;;!P.Charsize = 1.2
+;;altlin = 0
+;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
+;;        plots,  wherex + 2.*dddx/1.1,  wherey
+;;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;zecolor = 0. & loadct, 33
+;;for i=1, n_elements(paths)-1 do begin
+;;        ;zecolor = float(255)/n_elements(paths) + zecolor
+;;	zecolor = 0.
+;;        path = paths(i) & print, path
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        file = path+'/compturb.dat' & restore, filename=file
+;;        file = path+'/getturb.dat'  & restore, filename=file
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        zefield = localtime
+;;        zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top 
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
+;;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
+;;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
+;;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
+;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;endfor
+;;loadct, 0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;;if (saveps eq 'true') then PS_End, /PNG
+;
+;
+;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;
+;lev = 4
+;;;;;;;;;;;;;;;;
+;;;;;;;;;;;;;;;;
+;
+;        path = paths(0) & print, path
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;        file = path+'/compturb.dat' & restore, filename=file
+;        file = path+'/addturb.dat'  & restore, filename=file
+;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;zefield       =  localtime
+;zey 	      =  SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
+;set_name      =  saveplot+'comp_VELMAX.ps'
+;set_title     =  '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)"
+;set_titlex    =  'Local Time (h)'
+;set_titley    =  "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)"
+;set_subtitle  =  ''
+;;set_xrange    =  [09.0,18.0]
+;set_yrange    =  [00.0,14.0]
+;set_tickx     =  1.
+;set_ticky     =  2.
+;wherex        =  8.5;10.30
+;wherey        =  12.;0.7
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;print, h(lev)
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;if (saveps eq 'true') then PS_Start, FILENAME=set_name 
+;!P.Charsize = 1.2
+;altlin = 0
+;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
+;        plots,  wherex + 2.*dddx/1.1,  wherey
+;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;zecolor = 0. & loadct, 33
@@ -600,11 +724,8 @@
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;        file = path+'/compturb.dat' & restore, filename=file
-;        file = path+'/getturb.dat'  & restore, filename=file
+;        file = path+'/addturb.dat'  & restore, filename=file
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;        zefield = localtime
-;        zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top 
-;		;zey = DERIV(localtime, t_bot)
-;		zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
-;		;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
+;	zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
@@ -613,131 +734,10 @@
 ;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
 ;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;
 ;endfor
 ;loadct, 0
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;oplot, set_xrange, [0,0], linestyle=2
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;if (saveps eq 'true') then PS_End, /PNG
-
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;        path = paths(0) & print, path
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;        file = path+'/compturb.dat' & restore, filename=file
-;        file = path+'/getturb.dat'  & restore, filename=file
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  localtime
-;zey           =  SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top
-;set_name      =  saveplot+'comp_TTOP.ps'
-;set_title     =  "Temperature at the top of the boundary layer (K)"
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Temperature at the top of the boundary layer (K)'
-;set_subtitle  =  ''
-;set_xrange    =  [09.,18.]
-;set_yrange    =  [195.,225.]
-;set_tickx     =  1.
-;set_ticky     =  5.
-;wherex        =  13.00
-;wherey        =  223.5
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;altlin = 0
-;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
-;        plots,  wherex + 2.*dddx/1.1,  wherey
-;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;zecolor = 0. & loadct, 33
-;for i=1, n_elements(paths)-1 do begin
-;        ;zecolor = float(255)/n_elements(paths) + zecolor
-;	zecolor = 0.
-;        path = paths(i) & print, path
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;        file = path+'/compturb.dat' & restore, filename=file
-;        file = path+'/getturb.dat'  & restore, filename=file
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;        zefield = localtime
-;        zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top 
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
-;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;endfor
-;loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_End, /PNG
-
-
-;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;
-lev = 4
-;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;
-
-        path = paths(0) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/addturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-zefield       =  localtime
-zey 	      =  SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
-set_name      =  saveplot+'comp_VELMAX.ps'
-set_title     =  '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)"
-set_titlex    =  'Local Time (h)'
-set_titley    =  "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)"
-set_subtitle  =  ''
-;set_xrange    =  [09.0,18.0]
-set_yrange    =  [00.0,14.0]
-set_tickx     =  1.
-set_ticky     =  2.
-wherex        =  8.5;10.30
-wherey        =  12.;0.7
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-print, h(lev)
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_Start, FILENAME=set_name 
-!P.Charsize = 1.2
-altlin = 0
-plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
-        plots,  wherex + 2.*dddx/1.1,  wherey
-        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zecolor = 0. & loadct, 33
-for i=1, n_elements(paths)-1 do begin
-        ;zecolor = float(255)/n_elements(paths) + zecolor
-        zecolor = 0.
-        path = paths(i) & print, path
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        file = path+'/compturb.dat' & restore, filename=file
-        file = path+'/addturb.dat'  & restore, filename=file
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        zefield = localtime
-	zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
-        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
-        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
-        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
-        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
-endfor
-loadct, 0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-if (saveps eq 'true') then PS_End, /PNG
 
 
