Ignore:
Timestamp:
Apr 18, 2011, 8:19:21 AM (14 years ago)
Author:
aslmd
Message:

LMD_LES_MARS: qq modifications pour les routines graphiques

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/PLOT/SPEC/LES/compturb.pro

    r85 r113  
    1010;smoothampl=4.
    1111smoothampl=0.
    12 ;smoothampl=5.
     12smoothampl=5.
    1313;;;;;;;
    1414;
    … …  
    4141;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    4242
    43 ;; pour article, decommenter ici et commenter l'appel au deuxieme fichier
     43;;; pour article, decommenter ici et commenter l'appel au deuxieme fichier
     44;;;com
     45;;!p.multi = [0,2,1]
     46;
     47;        path = paths(0) & print, path
     48;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     49;        file = path+'/compturb.dat' & restore, filename=file
     50;        file = path+'/getturb.dat'  & restore, filename=file
     51;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     52;
     53;;;
     54;;; COEFFICIENT POUR LOIS DE SIMILITUDE
     55;;;
     56;;@coeffhf
     57;
     58;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     59;zefield       =  a_wt_bot
     60;zey           =  a_h
     61;set_name      =  saveplot+'comp_A_HFLUX.ps'
     62;set_title     =  "Similarity: vert. eddy heat flux on Mars"
     63;set_titlex    =  "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N"
     64;set_titley    =  'Dimensionless height z/z!Di!N'
     65;set_subtitle  =  ''
     66;set_xrange    =  [-0.5,1.]
     67;set_yrange    =  [0.,1.5]
     68;set_tickx     =  0.25
     69;set_ticky     =  0.1
     70;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     71;zesim = 1 - 1.2 * a_h
     72;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     73;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
     74;;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )]
     75;;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )]
     76;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
     77;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     78;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     79;!P.Charsize = 1.2
     80;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     81;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin
     82;for ll = 1, n_elements(localtimes)-1 do begin
     83;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     84;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
     85;endfor
     86;;wherex = set_xrange(0)
     87;;wherey = set_yrange(1)
     88;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path
     89;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     90;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     91;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     92;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     93;zecolor = 0. & loadct, 33
     94;for i=1, n_elements(paths)-1 do begin
     95;;        zecolor = 50. + zecolor
     96;zecolor = 0. & loadct, 0
     97;        path = paths(i) & print, path
     98;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     99;        file = path+'/compturb.dat' & restore, filename=file
     100;        file = path+'/getturb.dat'  & restore, filename=file
     101;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     102;        zefield = a_wt_bot
     103;        zey = a_h
     104;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     105;        for ll = 0, n_elements(localtimes)-1 do begin
     106;          user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     107;          if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
     108;        endfor
     109;        ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor
     110;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     111;endfor
     112;loadct, 0
     113;;oplot, [0.,0.], set_yrange, linestyle=0, color=0
     114;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     115;       ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
     116;       ;oplot, set_xrange, [0.,0.], linestyle=0, color=0
     117;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     118;if (saveps eq 'true') then PS_End, /PNG
    44119;;com
    45 ;!p.multi = [0,2,1]
    46 
    47         path = paths(0) & print, path
    48         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    49         file = path+'/compturb.dat' & restore, filename=file
    50         file = path+'/getturb.dat'  & restore, filename=file
    51         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    52 
    53 ;;
    54 ;; COEFFICIENT POUR LOIS DE SIMILITUDE
    55 ;;
    56 ;@coeffhf
    57 
    58 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    59 zefield       =  a_wt_bot
    60 zey           =  a_h
    61 set_name      =  saveplot+'comp_A_HFLUX.ps'
    62 set_title     =  "Similarity: vert. eddy heat flux on Mars"
    63 set_titlex    =  "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N"
    64 set_titley    =  'Dimensionless height z/z!Di!N'
    65 set_subtitle  =  ''
    66 set_xrange    =  [-0.5,1.]
    67 set_yrange    =  [0.,1.5]
    68 set_tickx     =  0.25
    69 set_ticky     =  0.1
    70 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    71 zesim = 1 - 1.2 * a_h
    72 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    73 localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
    74 ;;localtimes = localtime[where( (localtime ge 13) and (localtime le 16) )]
    75 ;;localtimes = localtime[where( (localtime ge 12) and (localtime le 15) )]
    76 ;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
    77 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    78 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    79 !P.Charsize = 1.2
    80 user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    81 plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.12, 0.22, 0.42, 1.0] ;color=0, linestyle=altlin
    82 for ll = 1, n_elements(localtimes)-1 do begin
    83   user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    84   if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
    85 endfor
    86 ;wherex = set_xrange(0)
    87 ;wherey = set_yrange(1)
    88 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + dddx, wherey, path
    89 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    90 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    91 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    92 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    93 zecolor = 0. & loadct, 33
    94 for i=1, n_elements(paths)-1 do begin
    95 ;        zecolor = 50. + zecolor
    96 zecolor = 0. & loadct, 0
    97         path = paths(i) & print, path
    98         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    99         file = path+'/compturb.dat' & restore, filename=file
    100         file = path+'/getturb.dat'  & restore, filename=file
    101         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    102         zefield = a_wt_bot
    103         zey = a_h
    104         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    105         for ll = 0, n_elements(localtimes)-1 do begin
    106           user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    107           if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
    108         endfor
    109         ;xyouts, wherex + dddx, wherey - float(i)*dddy, path, color=zecolor
    110         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    111 endfor
    112 loadct, 0
    113 ;oplot, [0.,0.], set_yrange, linestyle=0, color=0
    114         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    115         ;ye, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
    116         ;oplot, set_xrange, [0.,0.], linestyle=0, color=0
    117         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    118 if (saveps eq 'true') then PS_End, /PNG
    119 ;com
    120 
    121 
    122 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    123 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    124 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    125 
    126         path = paths(0) & print, path
    127         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    128         file = path+'/compturb.dat' & restore, filename=file
    129         file = path+'/getturb.dat'  & restore, filename=file
    130         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    131 
    132 ;;
    133 ;; COEFFICIENT POUR LOIS DE SIMILITUDE
    134 ;;
    135 ;@coeffw
    136 
    137 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    138 zefield       =  a_vel_bot
    139 zey           =  a_h 
    140 set_name      =  saveplot+'comp_A_VEL.ps'
    141 set_title     =  "Similarity: vert. vel. variance on Mars"
    142 set_titlex    =  "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N"
    143 set_titley    =  'Dimensionless height z/z!Di!N'
    144 set_subtitle  =  ''
    145 set_xrange    =  [0.,1.2]
    146 set_yrange    =  [0.,1.5]
    147 set_tickx     =  0.25 ;0.1
    148 set_ticky     =  0.1
    149 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    150 zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
    151 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    152 localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
    153 ;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
    154 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    155 ;com
    156 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    157 !P.Charsize = 1.2
    158 user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    159 plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin
    160 for ll = 1, n_elements(localtimes)-1 do begin
    161   user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    162   if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
    163 endfor
    164 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    165 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    166 zecolor = 0. & loadct, 33
    167 for i=1, n_elements(paths)-1 do begin
    168 zecolor = 0. & loadct, 0
    169 ;       zecolor = 50. + zecolor
    170         path = paths(i) & print, path
    171         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    172         file = path+'/compturb.dat' & restore, filename=file
    173         file = path+'/getturb.dat'  & restore, filename=file
    174         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    175         zefield = a_vel_bot
    176         zey = a_h
    177         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    178         for ll = 0, n_elements(localtimes)-1 do begin
    179           user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    180           if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
    181         endfor
    182         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    183 endfor
    184 loadct, 0
    185         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    186         ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
    187         ;oplot, set_xrange, [0.,0.], linestyle=0, color=0
    188         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    189 if (saveps eq 'true') then PS_End, /PNG
    190 
     120;
     121;
     122;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     123;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     124;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     125;
     126;       path = paths(0) & print, path
     127;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     128;        file = path+'/compturb.dat' & restore, filename=file
     129;        file = path+'/getturb.dat'  & restore, filename=file
     130;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     131;
     132;;;
     133;;; COEFFICIENT POUR LOIS DE SIMILITUDE
     134;;;
     135;;@coeffw
     136;
     137;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     138;zefield       =  a_vel_bot
     139;zey           =  a_h 
     140;set_name      =  saveplot+'comp_A_VEL.ps'
     141;set_title     =  "Similarity: vert. vel. variance on Mars"
     142;set_titlex    =  "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N"
     143;set_titley    =  'Dimensionless height z/z!Di!N'
     144;set_subtitle  =  ''
     145;set_xrange    =  [0.,1.2]
     146;set_yrange    =  [0.,1.5]
     147;set_tickx     =  0.25 ;0.1
     148;set_ticky     =  0.1
     149;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     150;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
     151;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     152;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
     153;;localtimes = localtime[where( (localtime ge 13) and (localtime le 15) )]
     154;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    191155;;com
    192 ;stop
    193 
    194 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    195 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    196 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    197 
    198         path = paths(0) & print, path
    199         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    200         file = path+'/compturb.dat' & restore, filename=file
    201         file = path+'/getturb.dat'  & restore, filename=file
    202         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    203 
    204 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    205 zefield       =  localtime
    206 zey           =  SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE)
    207 set_name      =  saveplot+'comp_W_STAR.ps'
    208 set_title     =  '';"Free convection velocity scale W!D*!N (m.s!U-1!N)"
    209 set_titlex    =  'Local Time (h)'
    210 set_titley    =  'Free convection velocity scale (m s!U-1!N)'
    211 set_subtitle  =  ''
    212 set_xrange    =  [08.,17.]
    213 set_yrange    =  [00.,06.]
    214 set_tickx     =  1.0
    215 set_ticky     =  0.5
    216 wherex        =  8.5;11.0
    217 wherey        =  5.5;1.75
    218 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    219 fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
    220 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    221 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    222 !P.Charsize = 1.2
    223 altlin=0
    224 plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    225 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path
    226         plots,  wherex + 2.*dddx/1.1,  wherey
    227         plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    228 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    229 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     156;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     157;!P.Charsize = 1.2
     158;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     159;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, psym=3;, /ISOTROPIC;, POS=[0.62, 0.22, 0.92, 1.0] ;color=0, linestyle=altlin
     160;for ll = 1, n_elements(localtimes)-1 do begin
     161;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     162;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3 ;linestyle=altlin
     163;endfor
     164;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     165;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    230166;zecolor = 0. & loadct, 33
    231 for i=1, n_elements(paths)-1 do begin
    232         ;zecolor = float(255)/n_elements(paths) + zecolor
    233         zecolor = 0
    234         path = paths(i) & print, path
    235         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    236         file = path+'/compturb.dat' & restore, filename=file
    237         file = path+'/getturb.dat'  & restore, filename=file
    238         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    239         zefield = localtime
    240         zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot
    241         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    242         fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
    243         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    244         oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    245         xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
    246         plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
    247         plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    248         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    249 endfor
    250 loadct, 0
    251 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    252 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    253 if (saveps eq 'true') then PS_End, /PNG
    254 
    255 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    256 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    257 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     167;for i=1, n_elements(paths)-1 do begin
     168;zecolor = 0. & loadct, 0
     169;;      zecolor = 50. + zecolor
     170;       path = paths(i) & print, path
     171;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     172;       file = path+'/compturb.dat' & restore, filename=file
     173;       file = path+'/getturb.dat'  & restore, filename=file
     174;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     175;       zefield = a_vel_bot
     176;       zey = a_h
     177;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     178;       for ll = 0, n_elements(localtimes)-1 do begin
     179;         user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     180;         if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=3, color=zecolor
     181;       endfor
     182;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     183;endfor
     184;loadct, 0
     185;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     186;       ;ye2, ryeah, tyeah & oplot, ryeah, tyeah, linestyle=0, color=255, /NOCLIP
     187;       ;oplot, set_xrange, [0.,0.], linestyle=0, color=0
     188;       ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     189;if (saveps eq 'true') then PS_End, /PNG
     190;
     191;;;com
     192;;stop
     193;
     194;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     195;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     196;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     197;
     198;        path = paths(0) & print, path
     199;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     200;        file = path+'/compturb.dat' & restore, filename=file
     201;        file = path+'/getturb.dat'  & restore, filename=file
     202;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     203;
     204;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     205;zefield       =  localtime
     206;zey           =  SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE)
     207;set_name      =  saveplot+'comp_W_STAR.ps'
     208;set_title     =  '';"Free convection velocity scale W!D*!N (m.s!U-1!N)"
     209;set_titlex    =  'Local Time (h)'
     210;set_titley    =  'Free convection velocity scale (m s!U-1!N)'
     211;set_subtitle  =  ''
     212;set_xrange    =  [08.,17.]
     213;set_yrange    =  [00.,06.]
     214;set_tickx     =  1.0
     215;set_ticky     =  0.5
     216;wherex       =  8.5;11.0
     217;wherey       =  5.5;1.75
     218;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     219;fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
     220;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     221;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     222;!P.Charsize = 1.2
     223;altlin=0
     224;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     225;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0);path
     226;        plots,  wherex + 2.*dddx/1.1,  wherey
     227;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
     228;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     229;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     230;;zecolor = 0. & loadct, 33
     231;for i=1, n_elements(paths)-1 do begin
     232;        ;zecolor = float(255)/n_elements(paths) + zecolor
     233;       zecolor = 0
     234;        path = paths(i) & print, path
     235;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     236;        file = path+'/compturb.dat' & restore, filename=file
     237;        file = path+'/getturb.dat'  & restore, filename=file
     238;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     239;        zefield = localtime
     240;        zey = SMOOTH(w_star_bot, [smoothampl], /EDGE_TRUNCATE);w_star_bot
     241;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     242;       fixfix = where((w_star_bot le 0.) and (localtime ge 15.)) & if (fixfix(0) ne -1) then zey(fixfix(0):n_elements(zey)-1) = !Values.F_NAN
     243;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     244;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
     245;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
     246;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
     247;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
     248;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     249;endfor
     250;loadct, 0
     251;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     252;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     253;if (saveps eq 'true') then PS_End, /PNG
     254;
     255;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     256;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     257;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    258258
    259259        path = paths(0) & print, path
    … …  
    271271set_titley    =  'Boundary layer height (km)'
    272272set_subtitle  =  ''
    273 ;set_xrange    =  [09.,18.]
    274 set_yrange    =  [00.,07.]
     273set_xrange    =  [08.,18.]
     274set_yrange    =  [00.,08.] ;[00.,07.]
    275275set_tickx     =  1.
    276276set_ticky     =  1.
    277277wherex        =  09.;13.25 ;12.75
    278 wherey        =  06.;2.5 ;2.0
     278wherey        =  06.5;2.5 ;2.0
    279279;       set_xrange    =  [13.5,18.] & set_yrange    =  [06.4,08.0] & set_tickx     =  0.5 & set_ticky     =  0.2
    280280;       wherex=15.30
    … …  
    321321if (saveps eq 'true') then PS_End, /PNG
    322322
     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;
    323380;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    324381;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    333390;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    334391;zefield       =  localtime
    335 ;zey           =  smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE) 
    336 ;set_name      =  saveplot+'comp_dHEIGHT.ps'
    337 ;set_title     =  "Boundary layer height growth dz!Di!N/dt (km/hour)"
     392;zey           =  SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
     393;zey2          =  SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
     394;set_name      =  saveplot+'comp_WMAX.ps'
     395;set_title     =  '';"Maximum and minimum vertical velocity (m.s!U-1!N)"
    338396;set_titlex    =  'Local Time (h)'
    339 ;set_titley    =  'Boundary layer height growth (km/hour)'
     397;set_titley    =  'Maximum and minimum vertical velocity (m s!U-1!N)'
    340398;set_subtitle  =  ''
    341 ;set_xrange    =  [9.,18.]
    342 ;set_yrange    =  [-0.4,1.6]
     399;;set_xrange    =  [09.0,18.0]
     400;set_yrange    =  [-12.0,18.0];[-15.0,20.0]
    343401;set_tickx     =  1.
    344 ;set_ticky     =  0.2
    345 ;wherex        =  10.30
    346 ;wherey        =  0.1
     402;set_ticky     =  2.
     403;wherex        =  8.5;11.00
     404;wherey        =  16.5;5.5
    347405;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    348406;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    … …  
    350408;altlin = 0
    351409;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    352 ;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
     410;oplot, zefield, zey2, color=0, linestyle=altlin
     411;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
    353412;        plots,  wherex + 2.*dddx/1.1,  wherey
    354413;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    355 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    356 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     414;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     415;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    357416;;zecolor = 0. & loadct, 33
    358417;for i=1, n_elements(paths)-1 do begin
    359418;        ;zecolor = float(255)/n_elements(paths) + zecolor
    360 ;       zecolor = 0
     419;        zecolor = 0.
    361420;        path = paths(i) & print, path
    362421;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    365424;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    366425;        zefield = localtime
    367 ;        zey = smooth(DERIV(localtime, smooth(pbl_height,smoothampl,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     426;        zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
     427;       zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
    368428;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    369429;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    370 ;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
     430;        oplot, zefield, zey2, color=zecolor, linestyle=altlin+i
     431;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
    371432;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
    372433;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    … …  
    377438;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    378439;if (saveps eq 'true') then PS_End, /PNG
    379 
    380 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    381 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    382 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    383 
    384         path = paths(0) & print, path
    385         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    386         file = path+'/compturb.dat' & restore, filename=file
    387         file = path+'/getturb.dat'  & restore, filename=file
    388         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    389 
    390 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    391 zefield       =  localtime
    392 zey           =  SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
    393 zey2          =  SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
    394 set_name      =  saveplot+'comp_WMAX.ps'
    395 set_title     =  '';"Maximum and minimum vertical velocity (m.s!U-1!N)"
    396 set_titlex    =  'Local Time (h)'
    397 set_titley    =  'Maximum and minimum vertical velocity (m s!U-1!N)'
    398 set_subtitle  =  ''
    399 ;set_xrange    =  [09.0,18.0]
    400 set_yrange    =  [-12.0,18.0];[-15.0,20.0]
    401 set_tickx     =  1.
    402 set_ticky     =  2.
    403 wherex        =  8.5;11.00
    404 wherey        =  16.5;5.5
    405 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    406 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    407 !P.Charsize = 1.2
    408 altlin = 0
    409 plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    410 oplot, zefield, zey2, color=0, linestyle=altlin
    411 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
    412         plots,  wherex + 2.*dddx/1.1,  wherey
    413         plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    414 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    415 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    416 ;zecolor = 0. & loadct, 33
    417 for i=1, n_elements(paths)-1 do begin
    418         ;zecolor = float(255)/n_elements(paths) + zecolor
    419         zecolor = 0.
    420         path = paths(i) & print, path
    421         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    422         file = path+'/compturb.dat' & restore, filename=file
    423         file = path+'/getturb.dat'  & restore, filename=file
    424         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    425         zefield = localtime
    426         zey = SMOOTH(wmaxmax, [smoothampl], /EDGE_TRUNCATE)
    427         zey2 = SMOOTH(wminmin, [smoothampl], /EDGE_TRUNCATE)
    428         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    429         oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    430         oplot, zefield, zey2, color=zecolor, linestyle=altlin+i
    431         xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
    432         plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
    433         plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    434         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    435 endfor
    436 loadct, 0
    437 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    439 if (saveps eq 'true') then PS_End, /PNG
    440 
    441 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    442 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    443 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    444 
    445         path = paths(0) & print, path
    446         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    447         file = path+'/compturb.dat' & restore, filename=file
    448         file = path+'/getturb.dat'  & restore, filename=file
    449         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    450 
    451 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    452 zefield       =  localtime
    453 zey           =  SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
    454 set_name      =  saveplot+'comp_TKE.ps'
    455 set_title     =  '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)"
    456 set_titlex    =  'Local Time (h)'
    457 set_titley    =  'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)'
    458 set_subtitle  =  ''
    459 ;set_xrange    =  [09.0,18.0]
    460 set_yrange    =  [00.0,18.0];[00.0,26.0]
    461 set_tickx     =  1.
    462 set_ticky     =  2.
    463 wherex        =  8.5;11.70
    464 wherey        =  16.;5.4
    465 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    466 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    467 !P.Charsize = 1.2
    468 altlin = 0
    469 plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    470 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/25. & xyouts, wherex + 2.*dddx, wherey, tttitles(0)
    471         plots,  wherex + 2.*dddx/1.1,  wherey
    472         plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    473 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    475 ;zecolor = 0. & loadct, 33
    476 for i=1, n_elements(paths)-1 do begin
    477         ;zecolor = float(255)/n_elements(paths) + zecolor
    478         zecolor = 0.
    479         path = paths(i) & print, path
    480         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    481         file = path+'/compturb.dat' & restore, filename=file
    482         file = path+'/getturb.dat'  & restore, filename=file
    483         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    484         zefield = localtime
    485         zey = SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
    486         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    487         oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    488         xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitles(i), color=zecolor
    489         plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
    490         plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    491         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    492 
    493 endfor
    494 loadct, 0
    495 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    496 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    497 if (saveps eq 'true') then PS_End, /PNG
    498 
    499 
    500 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    501 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    502 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    503 
    504         path = paths(0) & print, path
    505         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    506         file = path+'/compturb.dat' & restore, filename=file
    507         file = path+'/getturb.dat'  & restore, filename=file
    508         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    509 
    510 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    511 zefield       =  localtime
    512 zey           =  SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
    513 set_name      =  saveplot+'comp_HFBOT.ps'
    514 set_title     =  '';"Maximum heat flux <w'T'>!Dmax!N (K.m.s!U-1!N)"
    515 set_titlex    =  'Local Time (h)'
    516 set_titley    =  'Maximum vertical eddy heat flux (K m s!U-1!N)'
    517 set_subtitle  =  ''
    518 ;set_xrange    =  [09.0,18.0]
    519 set_yrange    =  [00.0,02.0];[00.0,02.8]
    520 set_tickx     =  1.
    521 set_ticky     =  0.2
    522 wherex        =  8.5;10.30
    523 wherey        =  1.8;0.7
    524 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    525 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    526 !P.Charsize = 1.2
    527 altlin = 0
    528 plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    529 dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
    530         plots,  wherex + 2.*dddx/1.1,  wherey
    531         plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    532 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    533 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    534 ;zecolor = 0. & loadct, 33
    535 for i=1, n_elements(paths)-1 do begin
    536         ;zecolor = float(255)/n_elements(paths) + zecolor
    537         zecolor = 0.
    538         path = paths(i) & print, path
    539         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    540         file = path+'/compturb.dat' & restore, filename=file
    541         file = path+'/getturb.dat'  & restore, filename=file
    542         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    543         zefield = localtime
    544         zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
    545         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    546         oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    547         xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
    548         plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
    549         plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    550         ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    551 
    552 endfor
    553 loadct, 0
    554 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    555 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    556 if (saveps eq 'true') then PS_End, /PNG
    557 
    558 
     440;
    559441;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    560442;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    569451;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    570452;zefield       =  localtime
    571 ;zey           =  t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
    572 ;zey          =  SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
    573 ;;zey           =  SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
    574 ;set_name      =  saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps'
    575 ;set_title     =  "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)"
     453;zey           =  SMOOTH(tkemax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
     454;set_name      =  saveplot+'comp_TKE.ps'
     455;set_title     =  '';"Maximum turbulent kinetic energy (m!U2!N.s!U-2!N)" ; 0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]!Dmax!N (m!U2!N.s!U-2!N)"
    576456;set_titlex    =  'Local Time (h)'
    577 ;set_titley    =  'Potential temperature time derivative (K.hour!U-1!N)'
     457;set_titley    =  'Maximum turbulent kinetic energy (m!U2!N s!U-2!N)'
    578458;set_subtitle  =  ''
    579 ;set_xrange    =  [ 09., 18.]
    580 ;set_yrange    =  [-01., 06.]
     459;;set_xrange    =  [09.0,18.0]
     460;set_yrange    =  [00.0,18.0];[00.0,26.0]
    581461;set_tickx     =  1.
    582 ;set_ticky     =  1.
    583 ;wherex        =  13.25
    584 ;wherey        =  5.5
     462;set_ticky     =  2.
     463;wherex        =  8.5;11.70
     464;wherey        =  16.;5.4
    585465;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    586466;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
    587526;!P.Charsize = 1.2
    588527;altlin = 0
    … …  
    591530;        plots,  wherex + 2.*dddx/1.1,  wherey
    592531;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
    593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     532;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     533;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     534;;zecolor = 0. & loadct, 33
     535;for i=1, n_elements(paths)-1 do begin
     536;        ;zecolor = float(255)/n_elements(paths) + zecolor
     537;       zecolor = 0.
     538;        path = paths(i) & print, path
     539;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     540;        file = path+'/compturb.dat' & restore, filename=file
     541;        file = path+'/getturb.dat'  & restore, filename=file
     542;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     543;        zefield = localtime
     544;        zey = SMOOTH(hfbotmax, [smoothampl], /EDGE_TRUNCATE);hfbotmax
     545;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     546;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
     547;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
     548;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
     549;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
     550;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     551;
     552;endfor
     553;loadct, 0
     554;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     555;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     556;if (saveps eq 'true') then PS_End, /PNG
     557;
     558;
     559;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     560;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     561;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     562;;
     563;;        path = paths(0) & print, path
     564;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     565;;        file = path+'/compturb.dat' & restore, filename=file
     566;;        file = path+'/getturb.dat'  & restore, filename=file
     567;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     568;;
     569;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     570;;zefield       =  localtime
     571;;zey           =  t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
     572;;zey         =  SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     573;;;zey           =  SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     574;;set_name      =  saveplot+'comp_TdtBOT.ps';saveplot+'comp_TBOT.ps'
     575;;set_title     =  "Potential temperature time derivative dT!Db!N/dt (K.hour!U-1!N)"
     576;;set_titlex    =  'Local Time (h)'
     577;;set_titley    =  'Potential temperature time derivative (K.hour!U-1!N)'
     578;;set_subtitle  =  ''
     579;;set_xrange    =  [ 09., 18.]
     580;;set_yrange    =  [-01., 06.]
     581;;set_tickx     =  1.
     582;;set_ticky     =  1.
     583;;wherex        =  13.25
     584;;wherey        =  5.5
     585;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     586;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     587;;!P.Charsize = 1.2
     588;;altlin = 0
     589;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     590;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
     591;;        plots,  wherex + 2.*dddx/1.1,  wherey
     592;;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
     593;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     594;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     595;;;zecolor = 0. & loadct, 33
     596;;for i=1, n_elements(paths)-1 do begin
     597;;        ;zecolor = float(255)/n_elements(paths) + zecolor
     598;;        zecolor = 0.
     599;;        path = paths(i) & print, path
     600;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     601;;        file = path+'/compturb.dat' & restore, filename=file
     602;;        file = path+'/getturb.dat'  & restore, filename=file
     603;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     604;;        zefield = localtime
     605;;        zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
     606;;              ;zey = DERIV(localtime, t_bot)
     607;;              zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     608;;              ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     609;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     610;;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
     611;;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, tttitle(i), color=zecolor
     612;;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
     613;;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
     614;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     615;;endfor
     616;;loadct, 0
     617;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     618;;;oplot, set_xrange, [0,0], linestyle=2
     619;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     620;;if (saveps eq 'true') then PS_End, /PNG
     621;
     622;
     623;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     624;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     625;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     626;;
     627;;        path = paths(0) & print, path
     628;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     629;;        file = path+'/compturb.dat' & restore, filename=file
     630;;        file = path+'/getturb.dat'  & restore, filename=file
     631;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     632;;
     633;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     634;;zefield       =  localtime
     635;;zey           =  SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top
     636;;set_name      =  saveplot+'comp_TTOP.ps'
     637;;set_title     =  "Temperature at the top of the boundary layer (K)"
     638;;set_titlex    =  'Local Time (h)'
     639;;set_titley    =  'Temperature at the top of the boundary layer (K)'
     640;;set_subtitle  =  ''
     641;;set_xrange    =  [09.,18.]
     642;;set_yrange    =  [195.,225.]
     643;;set_tickx     =  1.
     644;;set_ticky     =  5.
     645;;wherex        =  13.00
     646;;wherey        =  223.5
     647;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     648;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     649;;!P.Charsize = 1.2
     650;;altlin = 0
     651;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     652;;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, path
     653;;        plots,  wherex + 2.*dddx/1.1,  wherey
     654;;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
     655;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     656;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     657;;;zecolor = 0. & loadct, 33
     658;;for i=1, n_elements(paths)-1 do begin
     659;;        ;zecolor = float(255)/n_elements(paths) + zecolor
     660;;      zecolor = 0.
     661;;        path = paths(i) & print, path
     662;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     663;;        file = path+'/compturb.dat' & restore, filename=file
     664;;        file = path+'/getturb.dat'  & restore, filename=file
     665;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     666;;        zefield = localtime
     667;;        zey = SMOOTH(tt_top, [smoothampl], /EDGE_TRUNCATE);tt_top
     668;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     669;;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
     670;;        xyouts, wherex + 2.*dddx,      wherey - float(i)*dddy, path, color=zecolor
     671;;        plots,  wherex + 2.*dddx/1.1,  wherey - float(i)*dddy
     672;;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
     673;;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     674;;endfor
     675;;loadct, 0
     676;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     677;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     678;;if (saveps eq 'true') then PS_End, /PNG
     679;
     680;
     681;;;;;;;;;;;;;;;;
     682;;;;;;;;;;;;;;;;
     683;lev = 4
     684;;;;;;;;;;;;;;;;
     685;;;;;;;;;;;;;;;;
     686;
     687;        path = paths(0) & print, path
     688;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     689;        file = path+'/compturb.dat' & restore, filename=file
     690;        file = path+'/addturb.dat'  & restore, filename=file
     691;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     692;
     693;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     694;zefield       =  localtime
     695;zey          =  SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
     696;set_name      =  saveplot+'comp_VELMAX.ps'
     697;set_title     =  '';"Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m.s!U-1!N)"
     698;set_titlex    =  'Local Time (h)'
     699;set_titley    =  "Maximum horizontal velocity "+string(h(lev)*1000.,'(I0)')+"m above ground (m s!U-1!N)"
     700;set_subtitle  =  ''
     701;;set_xrange    =  [09.0,18.0]
     702;set_yrange    =  [00.0,14.0]
     703;set_tickx     =  1.
     704;set_ticky     =  2.
     705;wherex        =  8.5;10.30
     706;wherey        =  12.;0.7
     707;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     708;print, h(lev)
     709;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     710;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     711;!P.Charsize = 1.2
     712;altlin = 0
     713;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     714;dddx = (set_xrange(1) - set_xrange(0))/20. & dddy = (set_yrange(1) - set_yrange(0))/20. & xyouts, wherex + 2.*dddx, wherey, tttitle(0)
     715;        plots,  wherex + 2.*dddx/1.1,  wherey
     716;        plots,  wherex + 2.*dddx/20.,  wherey, linestyle=altlin, /continue
     717;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    594718;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    595719;;zecolor = 0. & loadct, 33
    … …  
    600724;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    601725;        file = path+'/compturb.dat' & restore, filename=file
    602 ;        file = path+'/getturb.dat'  & restore, filename=file
     726;        file = path+'/addturb.dat'  & restore, filename=file
    603727;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    604728;        zefield = localtime
    605 ;        zey = t_bot;SMOOTH(t_bot, [smoothampl], /EDGE_TRUNCATE);tt_top
    606 ;               ;zey = DERIV(localtime, t_bot)
    607 ;               zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot,smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
    608 ;               ;zey = SMOOTH(DERIV(localtime, SMOOTH(t_bot/((p0 / p_bot)^r_cp),smoothampl/2.,/EDGE_TRUNCATE)),smoothampl,/EDGE_TRUNCATE)
     729;       zey = SMOOTH(reform(velmax(lev,*)), [smoothampl], /EDGE_TRUNCATE)
    609730;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    610731;        oplot, zefield, zey, color=zecolor, linestyle=altlin+i
    … …  
    613734;        plots,  wherex + 2.*dddx/20.,  wherey - float(i)*dddy, linestyle=altlin+i, /continue
    614735;        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     736;
    615737;endfor
    616738;loadct, 0
    617739;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    618 ;;oplot, set_xrange, [0,0], linestyle=2
    619740;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    620741;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
    742742
    743743
Note: See TracChangeset for help on using the changeset viewer.