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/getturb.pro

    r96 r113  
    1717;; constantes
    1818;;
    19 ;p0=610. & t0=220. & r_cp=1/4.4 & grav=3.72 & R=192.
    20 p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191.
    21 print, 'ATTENTION ATTENTION R/cp !!!!', r_cp
     19p0=610. & t0=220. & r_cp=1.0/4.4 & grav=3.72 & R=192.
     20;p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191.
     21;print, 'ATTENTION ATTENTION R/cp !!!!', r_cp
    2222
    2323;
    … …  
    171171OPENR, 23, 'input_more' & READF, 23, hgtu, tsurfu & CLOSE, 23
    172172
    173 @include_ustar.pro
     173;@include_ustar.pro
    174174
    175175;@include_hfx.pro
    … …  
    344344;*******************;
    345345
    346 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    347 ;zefield       =  a_vel_bot
    348 ;zey           =  a_h 
    349 ;set_name      =  'A_VEL.ps'
    350 ;set_title     =  "Scaled vertical velocity variance (m.s!U-1!N)"
    351 ;set_titlex    =  'Scaled vertical velocity variance (m.s!U-1!N)'
    352 ;set_titley    =  'Scaled height (km)'
    353 ;set_subtitle  =  ''
    354 ;set_xrange    =  [0.,1.2]
    355 ;set_yrange    =  [0.,1.2]
    356 ;set_tickx     =  0.1
    357 ;set_ticky     =  0.1
    358 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    359 ;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
    360 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    361 ;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
    362 ;;localtimes = localtime
    363 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    364 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    365 ;!P.Charsize = 1.2
    366 ;altlin=0
    367 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    368 ;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin, /ISOTROPIC
    369 ;for ll = 1, n_elements(localtimes)-1 do begin
    370 ;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    371 ;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
    372 ;endfor
    373 ;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0
    374 ;if (saveps eq 'true') then PS_End, /PNG
    375 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    376 
    377 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    378 zefield       =  localtime
    379 zey           =  pbl_height
    380 set_name      =  'HEIGHT.ps'
    381 set_title     =  'Boundary layer height (km)'
    382 set_titlex    =  'Local Time (h)'
    383 set_titley    =  'Boundary layer height (km)'
    384 set_subtitle  =  '' ;'Criterion is : static stability > 1.5 K.m!U-1!N'
    385 set_xrange    =  [8.,17.]
    386 set_yrange    =  [0.,7.]
    387 set_tickx     =  1.
    388 set_ticky     =  1.
    389 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    390 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    391 !P.Charsize = 1.2
    392 plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    393 ;oplot, zefield, zey, psym=5
    394 if (saveps eq 'true') then PS_End, /PNG
    395 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    396 
     346;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     347;;zefield       =  a_vel_bot
     348;;zey           =  a_h 
     349;;set_name      =  'A_VEL.ps'
     350;;set_title     =  "Scaled vertical velocity variance (m.s!U-1!N)"
     351;;set_titlex    =  'Scaled vertical velocity variance (m.s!U-1!N)'
     352;;set_titley    =  'Scaled height (km)'
     353;;set_subtitle  =  ''
     354;;set_xrange    =  [0.,1.2]
     355;;set_yrange    =  [0.,1.2]
     356;;set_tickx     =  0.1
     357;;set_ticky     =  0.1
     358;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     359;;zesim = 1.8 * a_h^(2./3.) * ( 1. - 0.8 * a_h )^2
     360;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     361;;localtimes = localtime[where( (localtime ge 11) and (localtime le 16) )]
     362;;;localtimes = localtime
     363;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     364;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     365;;!P.Charsize = 1.2
     366;;altlin=0
     367;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     368;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin, /ISOTROPIC
     369;;for ll = 1, n_elements(localtimes)-1 do begin
     370;;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     371;;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
     372;;endfor
     373;;loadct, 33 & oplot, zesim, zey, psym=5, color=255 & loadct, 0
     374;;if (saveps eq 'true') then PS_End, /PNG
     375;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     376;
    397377;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    398378;zefield       =  localtime
    399 ;zey           =  w_star_bot
    400 ;set_name      =  'W_STAR.ps'
    401 ;set_title     =  "Free convection velocity scale (m.s!U-1!N)"
     379;zey           =  pbl_height
     380;set_name      =  'HEIGHT.ps'
     381;set_title     =  'Boundary layer height (km)'
    402382;set_titlex    =  'Local Time (h)'
    403 ;set_titley    =  'Free convection velocity scale (m.s!U-1!N)'
    404 ;set_subtitle  =  ''
    405 ;set_xrange    =  [8.,18.]
    406 ;set_yrange    =  [0.,6.]
     383;set_titley    =  'Boundary layer height (km)'
     384;set_subtitle  =  '' ;'Criterion is : static stability > 1.5 K.m!U-1!N'
     385;set_xrange    =  [11.,17.] ;[8.,17.]
     386;set_yrange    =  [0.,8.]
    407387;set_tickx     =  1.
    408 ;set_ticky     =  0.5
     388;set_ticky     =  1.
    409389;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    410390;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    … …  
    414394;if (saveps eq 'true') then PS_End, /PNG
    415395;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    416 
    417 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    418 ;zefield       =  localtime
    419 ;zey           =  hfbot
    420 ;zey2         =  hfsurf
    421 ;zey3          =  hfbotmax
    422 ;set_name      =  'HFBOT.ps'
    423 ;set_title     =  "Heat flux at the top of the radiative layer (K.m.s!U-1!N)"
    424 ;set_titlex    =  'Local Time (h)'
    425 ;set_titley    =  'Heat flux at the top of the radiative layer (K.m.s!U-1!N)'
    426 ;set_subtitle  =  ''
    427 ;set_xrange    =  [8.,18.]
    428 ;set_yrange    =  [0.,2.5]
    429 ;set_tickx     =  1.
    430 ;set_ticky     =  0.5
    431 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    432 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    433 ;!P.Charsize = 1.2
    434 ;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    435 ;oplot, zefield, zey2, linestyle=1
    436 ;oplot, zefield, zey3, linestyle=2
    437 ;if (saveps eq 'true') then PS_End, /PNG
    438 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    439 
    440 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    441 zefield       =  transpose(tke)
    442 zex           =  localtime
    443 zey           =  h
    444 set_name      =  'TKE.ps'
    445 set_title     =  "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]
    446 set_titlex    =  'Local Time (h)'
    447 set_titley    =  'Altitude above surface (km)'
    448 set_subtitle  =  '' ;'Mean over the simulation domain'
    449 set_xrange    =  [8.,17.]
    450 set_yrange    =  [0.,7.]  ;; [0.,200.]
    451 set_tickx     =  1.
    452 set_ticky     =  1. ;; 50.
    453 minval        =  0.
    454 maxval        =  15.
    455 nlev          =  maxval-minval
    456 pal           =  22
    457 rrr           =  'no'
    458 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    459 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    460 !P.Charsize = 1.2
    461 ;; 0. levels
    462 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
    463 ;;; 1. background
    464 loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    465   ;; 2. color field
    466   loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
    467               contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    468 ;; 3. contour field
    469 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
    470 ;; 4. choose output
    471 if (saveps eq 'true') then PS_End, /PNG
    472 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    473 
    474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    475 ;zefield       =  2.*transpose(ztke)
    476 ;set_title     =  "Vertical wind variance (m!U2!N.s!U-2!N)"
    477 ;minval        =  -6.
    478 ;maxval        =  12.
    479 ;pal           =  0
    480 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    481 zefield       =  transpose(ztke)
    482 zex           =  localtime
    483 zey           =  h
    484 set_name      =  'zTKE.ps'
    485 set_title     =  "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>]
    486 set_titlex    =  'Local Time (h)'
    487 set_titley    =  'Altitude above surface (km)'
    488 set_subtitle  =  '' ;'Mean over the simulation domain'
    489 set_xrange    =  [8.,17.]
    490 set_yrange    =  [0.,7.]
    491 set_tickx     =  1.
    492 set_ticky     =  1.
    493 minval        =  0.
    494 maxval        =  8.
    495 nlev          =  maxval-minval
    496 pal           =  22
    497 rrr           =  'no'
    498 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    499 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    500 !P.Charsize = 1.2
    501 ;; 0. levels
    502 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
    503 ;;; 1. background
    504 loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    505 ;; 2. color field
    506 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
    507             contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    508 ;; 3. contour field
    509 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
    510 ;; 4. choose output
    511 if (saveps eq 'true') then PS_End, /PNG
    512 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    513 
    514 ;restore, filename='addturb.dat'
    515 ;modvar = SMOOTH(TEMPORARY(modvar),  [0,smoothampl], /EDGE_TRUNCATE)
    516 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    517 ;zefield       =  transpose(modvar)
    518 ;zex           =  localtime
     396;
     397;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     398;;zefield       =  localtime
     399;;zey           =  w_star_bot
     400;;set_name      =  'W_STAR.ps'
     401;;set_title     =  "Free convection velocity scale (m.s!U-1!N)"
     402;;set_titlex    =  'Local Time (h)'
     403;;set_titley    =  'Free convection velocity scale (m.s!U-1!N)'
     404;;set_subtitle  =  ''
     405;;set_xrange    =  [8.,18.]
     406;;set_yrange    =  [0.,6.]
     407;;set_tickx     =  1.
     408;;set_ticky     =  0.5
     409;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     410;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     411;;!P.Charsize = 1.2
     412;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     413;;;oplot, zefield, zey, psym=5
     414;;if (saveps eq 'true') then PS_End, /PNG
     415;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     416;
     417;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     418;;zefield       =  localtime
     419;;zey           =  hfbot
     420;;zey2        =  hfsurf
     421;;zey3          =  hfbotmax
     422;;set_name      =  'HFBOT.ps'
     423;;set_title     =  "Heat flux at the top of the radiative layer (K.m.s!U-1!N)"
     424;;set_titlex    =  'Local Time (h)'
     425;;set_titley    =  'Heat flux at the top of the radiative layer (K.m.s!U-1!N)'
     426;;set_subtitle  =  ''
     427;;set_xrange    =  [8.,18.]
     428;;set_yrange    =  [0.,2.5]
     429;;set_tickx     =  1.
     430;;set_ticky     =  0.5
     431;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     432;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     433;;!P.Charsize = 1.2
     434;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     435;;oplot, zefield, zey2, linestyle=1
     436;;oplot, zefield, zey3, linestyle=2
     437;;if (saveps eq 'true') then PS_End, /PNG
     438;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     439;
     440;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     441;zefield       =  transpose(tke)
     442;zex           =  localtime
    519443;zey           =  h
    520 ;set_name      =  'modvar.ps'
    521 ;set_title     =  "Horizontal wind speed variance (m!U2!N.s!U-2!N)"
     444;set_name      =  'TKE.ps'
     445;set_title     =  "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]
    522446;set_titlex    =  'Local Time (h)'
    523447;set_titley    =  'Altitude above surface (km)'
    524448;set_subtitle  =  '' ;'Mean over the simulation domain'
    525 ;set_xrange    =  [7.,17.]
     449;set_xrange    =  [11.,17.] ;[8.,17.]
     450;set_yrange    =  [0.,8.]  ;; [0.,200.]
     451;set_tickx     =  1.
     452;set_ticky     =  1. ;; 50.
     453;minval        =  0.
     454;maxval        =  20.
     455;nlev          =  maxval-minval
     456;pal           =  22
     457;rrr           =  'no'
     458;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     459;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     460;!P.Charsize = 1.2
     461;;; 0. levels
     462;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     463;;;; 1. background
     464;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
     465;  ;; 2. color field
     466;  loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
     467;              contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
     468;;; 3. contour field
     469;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     470;;; 4. choose output
     471;if (saveps eq 'true') then PS_End, /PNG
     472;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     473;
     474;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     475;;zefield       =  2.*transpose(ztke)
     476;;set_title     =  "Vertical wind variance (m!U2!N.s!U-2!N)"
     477;;minval        =  -6.
     478;;maxval        =  12.
     479;;pal           =  0
     480;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     481;zefield       =  transpose(ztke)
     482;zex           =  localtime
     483;zey           =  h
     484;set_name      =  'zTKE.ps'
     485;set_title     =  "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>]
     486;set_titlex    =  'Local Time (h)'
     487;set_titley    =  'Altitude above surface (km)'
     488;set_subtitle  =  '' ;'Mean over the simulation domain'
     489;set_xrange    =  [8.,17.]
    526490;set_yrange    =  [0.,7.]
    527491;set_tickx     =  1.
    528492;set_ticky     =  1.
    529493;minval        =  0.
    530 ;maxval        =  10.
     494;maxval        =  8.
    531495;nlev          =  maxval-minval
    532496;pal           =  22
    … …  
    547511;if (saveps eq 'true') then PS_End, /PNG
    548512;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    549 
    550 
    551 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    552 zefield       =  transpose(wt)
    553 zex           =  localtime
    554 zey           =  h
    555 set_name      =  'HF.ps'
    556 set_title     =  "Turbulent Heat Flux (K.m.s!U-1!N)"  ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'>
    557 set_titlex    =  'Local Time (h)'
    558 set_titley    =  'Altitude above surface (km)'
    559 set_subtitle  =  '' ;'Mean over the simulation domain'
    560 set_xrange    =  [8.,17.]
    561 set_yrange    =  [0.,7.]  ;; [0.,200.]
    562 set_tickx     =  1.
    563 set_ticky     =  1. ;; 50.
    564 minval        =  -2.
    565 maxval        =  2.
    566 nlev          =  floor(maxval-minval)*10
    567 pal           =  33 
    568 rrr           =  'no'
    569 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    570 ;minval        =  -0.8
    571 ;maxval        =  1.2
    572 ;pal           =  0
    573 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    574 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    575 !P.Charsize = 1.2
    576 ;; 0. levels
    577 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
    578 ;;; 1. background
    579 loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    580 ;; 2. color field
    581 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
    582             ;;;--------------------------------------------------------------------------------------------------------------------------------
    583             ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone
    584             ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5
    585             nu = nd + (nu-nd)/2  ;; otherwise the interval is too large (because we removed 0)
    586             TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b
    587             ;;;--------------------------------------------------------------------------------------------------------------------------------
    588             contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    589 ;; 3. contour field
    590 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
    591 ;; 4. choose output
    592 if (saveps eq 'true') then PS_End, /PNG
    593 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    594 
    595 ;;restore, filename='tpot_profB'
    596 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    597 ;zefield       =  t
    598 ;zey           =  ht
    599 ;set_name      =  'T.ps'
    600 ;set_title     =  "Potential Temperature (K)"
    601 ;set_titlex    =  'Potential Temperature (K)'
     513;
     514;;restore, filename='addturb.dat'
     515;;modvar = SMOOTH(TEMPORARY(modvar),  [0,smoothampl], /EDGE_TRUNCATE)
     516;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     517;;zefield       =  transpose(modvar)
     518;;zex           =  localtime
     519;;zey           =  h
     520;;set_name      =  'modvar.ps'
     521;;set_title     =  "Horizontal wind speed variance (m!U2!N.s!U-2!N)"
     522;;set_titlex    =  'Local Time (h)'
     523;;set_titley    =  'Altitude above surface (km)'
     524;;set_subtitle  =  '' ;'Mean over the simulation domain'
     525;;set_xrange    =  [7.,17.]
     526;;set_yrange    =  [0.,7.]
     527;;set_tickx     =  1.
     528;;set_ticky     =  1.
     529;;minval        =  0.
     530;;maxval        =  10.
     531;;nlev          =  maxval-minval
     532;;pal           =  22
     533;;rrr           =  'no'
     534;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     535;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     536;;!P.Charsize = 1.2
     537;;;; 0. levels
     538;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     539;;;;; 1. background
     540;;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
     541;;;; 2. color field
     542;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
     543;;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
     544;;;; 3. contour field
     545;;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     546;;;; 4. choose output
     547;;if (saveps eq 'true') then PS_End, /PNG
     548;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     549;
     550;
     551;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     552;zefield       =  transpose(wt)
     553;zex           =  localtime
     554;zey           =  h
     555;set_name      =  'HF.ps'
     556;set_title     =  "Turbulent Heat Flux (K.m.s!U-1!N)"  ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'>
     557;set_titlex    =  'Local Time (h)'
    602558;set_titley    =  'Altitude above surface (km)'
    603559;set_subtitle  =  '' ;'Mean over the simulation domain'
    604 ;set_xrange    =  [190.,230.] ;[min(t),max(t)]
    605 ;set_yrange    =  [0.,7.]
    606 ;set_tickx     =  5.
    607 ;set_ticky     =  1.
    608 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    609 ;;localtimes = [9,10,11,12,13,14,15,16,17,18]
    610 ;localtimes = [7,9,11,13,15,17]
    611 ;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0]
    612 ;;localtimes = [15.0]
     560;set_xrange    =  [8.,17.]
     561;set_yrange    =  [0.,8.]  ;; [0.,200.]
     562;set_tickx     =  1.
     563;set_ticky     =  1. ;; 50.
     564;minval        =  -1.5 ;-2.
     565;maxval        =  2.5  ;2.
     566;nlev          =  floor(maxval-minval)*10
     567;pal           =  33 
     568;rrr           =  'no'
     569;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     570;;minval        =  -0.8
     571;;maxval        =  1.2
     572;;pal           =  0
    613573;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    614574;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    615575;!P.Charsize = 1.2
     576;;; 0. levels
     577;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     578;;;; 1. background
     579;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
     580;;; 2. color field
     581;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
     582;            ;;;--------------------------------------------------------------------------------------------------------------------------------
     583;            ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone
     584;            ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5
     585;            nu = nd + (nu-nd)/2  ;; otherwise the interval is too large (because we removed 0)
     586;            TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b
     587;            ;;;--------------------------------------------------------------------------------------------------------------------------------
     588;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
     589;;; 3. contour field
     590;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     591;;; 4. choose output
     592;if (saveps eq 'true') then PS_End, /PNG
     593;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     594;
     595;;;restore, filename='tpot_profB'
     596;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     597;;zefield       =  t
     598;;zey           =  ht
     599;;set_name      =  'T.ps'
     600;;set_title     =  "Potential Temperature (K)"
     601;;set_titlex    =  'Potential Temperature (K)'
     602;;set_titley    =  'Altitude above surface (km)'
     603;;set_subtitle  =  '' ;'Mean over the simulation domain'
     604;;set_xrange    =  [190.,230.] ;[min(t),max(t)]
     605;;set_yrange    =  [0.,7.]
     606;;set_tickx     =  5.
     607;;set_ticky     =  1.
     608;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     609;;;localtimes = [9,10,11,12,13,14,15,16,17,18]
     610;;localtimes = [7,9,11,13,15,17]
     611;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0]
     612;;;localtimes = [15.0]
     613;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     614;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     615;;!P.Charsize = 1.2
     616;;;altlin=0 & loadct, 0
     617;;altlin=1 & loadct, 0
     618;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     619;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     620;;;plot, zefield(*,nntt), les_column, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     621;;for ll = 1, n_elements(localtimes)-1 do begin
     622;;  CASE altlin OF
     623;;  0: altlin=1
     624;;  1: altlin=0
     625;;  ENDCASE 
     626;;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     627;;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
     628;;;  if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin
     629;;endfor
     630;;;oplot, ro_tpot, ro_column, psym=7
     631;;;xyouts, 192.0, 0.25, '09:00'
     632;;;xyouts, 205.5, 0.25, '11:00'
     633;;;xyouts, 214.5, 0.25, '13:00'
     634;;;xyouts, 219.5, 0.25, '17:00'
     635;;;xyouts, 220.5, 2.30, '17:00 [RO]'
     636;;if (saveps eq 'true') then PS_End, /PNG
     637;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     638;
     639;
     640;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     641;;zefield       =  stst
     642;;zey           =  ht
     643;;set_name      =  'STST.ps'
     644;;set_title     =  "Static stability (K.m!U-1!N)"
     645;;set_titlex    =  'Static stability (K.m!U-1!N)'
     646;;set_titley    =  'Altitude above surface (km)'
     647;;set_subtitle  =  'Mean over the simulation domain'
     648;;set_xrange    =  [-5.,5.]
     649;;set_yrange    =  [0.,7.]
     650;;set_tickx     =  1.
     651;;set_ticky     =  1.
     652;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     653;;localtimes = [9,11,13,15,17]
     654;;localtimes = [12,13,14,15,16]
     655;;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0]
     656;;;localtimes = [15.0]
     657;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     658;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     659;;!P.Charsize = 1.2
    616660;;altlin=0 & loadct, 0
    617 ;altlin=1 & loadct, 0
    618 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    619 ;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    620 ;;plot, zefield(*,nntt), les_column, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    621 ;for ll = 1, n_elements(localtimes)-1 do begin
    622 ;  CASE altlin OF
    623 ;  0: altlin=1
    624 ;  1: altlin=0
    625 ;  ENDCASE 
    626 ;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    627 ;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
    628 ;;  if (nntt ne -1) then oplot, zefield(*,nntt), les_column, linestyle=altlin
    629 ;endfor
    630 ;;oplot, ro_tpot, ro_column, psym=7
    631 ;;xyouts, 192.0, 0.25, '09:00'
    632 ;;xyouts, 205.5, 0.25, '11:00'
    633 ;;xyouts, 214.5, 0.25, '13:00'
    634 ;;xyouts, 219.5, 0.25, '17:00'
    635 ;;xyouts, 220.5, 2.30, '17:00 [RO]'
    636 ;if (saveps eq 'true') then PS_End, /PNG
    637 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    638 
    639 
    640 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    641 ;zefield       =  stst
    642 ;zey           =  ht
    643 ;set_name      =  'STST.ps'
    644 ;set_title     =  "Static stability (K.m!U-1!N)"
    645 ;set_titlex    =  'Static stability (K.m!U-1!N)'
    646 ;set_titley    =  'Altitude above surface (km)'
    647 ;set_subtitle  =  'Mean over the simulation domain'
    648 ;set_xrange    =  [-5.,5.]
    649 ;set_yrange    =  [0.,7.]
    650 ;set_tickx     =  1.
    651 ;set_ticky     =  1.
    652 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    653 ;localtimes = [9,11,13,15,17]
    654 ;localtimes = [12,13,14,15,16]
    655 ;;localtimes = [17.0,17.1,17.2,17.3,17.4,17.5,17.6,17.7,17.8,17.9,18.0]
    656 ;;localtimes = [15.0]
    657 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    658 ;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    659 ;!P.Charsize = 1.2
    660 ;altlin=0 & loadct, 0
    661 ;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    662 ;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    663 ;oplot, zefield(*,nntt), zey(*,nntt), psym=5
    664 ;for ll = 1, n_elements(localtimes)-1 do begin
    665 ;  CASE altlin OF
    666 ;  0: altlin=1
    667 ;  1: altlin=0
    668 ;  ENDCASE
    669 ;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
    670 ;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
    671 ;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5
    672 ;endfor
    673 ;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 
    674 ;if (saveps eq 'true') then PS_End, /PNG
    675 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    676 
    677 
    678 
    679 if (n_elements(xtke) eq 0) then stop
    680 
    681 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    682 ;;zefield       =  2.*transpose(xtke)
    683 ;;minval        =  -4
    684 ;;maxval        =  8. ;12.
    685 ;;pal           =  0
    686 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     661;;user_lt=localtimes[0] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     662;;plot, zefield(*,nntt), zey(*,nntt), xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     663;;oplot, zefield(*,nntt), zey(*,nntt), psym=5
     664;;for ll = 1, n_elements(localtimes)-1 do begin
     665;;  CASE altlin OF
     666;;  0: altlin=1
     667;;  1: altlin=0
     668;;  ENDCASE
     669;;  user_lt=localtimes[ll] & yeah=where(abs(localtime-user_lt) eq (min(abs(localtime-user_lt)))) & nntt = yeah(0)
     670;;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), linestyle=altlin
     671;;  if (nntt ne -1) then oplot, zefield(*,nntt), zey(*,nntt), psym=5
     672;;endfor
     673;;oplot, 0.*zefield(*,nntt) + 1.5, zey(*,nntt), linestyle=2 
     674;;if (saveps eq 'true') then PS_End, /PNG
     675;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     676;
     677;
     678;
     679;if (n_elements(xtke) eq 0) then stop
     680;
     681;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     682;;;zefield       =  2.*transpose(xtke)
     683;;;minval        =  -4
     684;;;maxval        =  8. ;12.
     685;;;pal           =  0
     686;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     687;;zex           =  localtime
     688;;zey           =  h
     689;;set_name      =  'xTKE.ps'
     690;;set_title     =  "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>]
     691;;set_titlex    =  'Local Time (h)'
     692;;set_titley    =  'Altitude above surface (km)'
     693;;set_subtitle  =  '';'Mean over the simulation domain'
     694;;set_xrange    =  [8.,17.]
     695;;set_yrange    =  [0.,7.]
     696;;set_tickx     =  1.
     697;;set_ticky     =  1.
     698;;minval        =  0.
     699;;maxval        =  8.
     700;;nlev          =  maxval-minval
     701;;pal           =  22
     702;;rrr           =  'no'
     703;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     704;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     705;;!P.Charsize = 1.2
     706;;;; 0. levels
     707;;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     708;;;;; 1. background
     709;;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
     710;;;; 2. color field
     711;;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
     712;;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
     713;;;; 3. contour field
     714;;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     715;;;; 4. choose output
     716;;if (saveps eq 'true') then PS_End, /PNG
     717;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     718;
     719;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     720;zefield       =  transpose(ytke) + transpose(xtke)
    687721;zex           =  localtime
    688722;zey           =  h
    689 ;set_name      =  'xTKE.ps'
    690 ;set_title     =  "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>]
     723;set_name      =  'hTKE.ps'
     724;set_title     =  "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<v'!U2!N>]
    691725;set_titlex    =  'Local Time (h)'
    692726;set_titley    =  'Altitude above surface (km)'
    693 ;set_subtitle  =  '';'Mean over the simulation domain'
     727;set_subtitle  =  '' ;'Mean over the simulation domain'
    694728;set_xrange    =  [8.,17.]
    695 ;set_yrange    =  [0.,7.]
     729;set_yrange    =  [0.,1.] ;[0.,7.] ;; [0.,200.]
    696730;set_tickx     =  1.
    697 ;set_ticky     =  1.
     731;set_ticky     =  0.1 ;1. ;50.
    698732;minval        =  0.
    699 ;maxval        =  8.
     733;maxval        =  12. ;10.
    700734;nlev          =  maxval-minval
    701735;pal           =  22
    … …  
    716750;if (saveps eq 'true') then PS_End, /PNG
    717751;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    718 
    719 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    720 zefield       =  transpose(ytke) + transpose(xtke)
    721 zex           =  localtime
    722 zey           =  h
    723 set_name      =  'hTKE.ps'
    724 set_title     =  "Horizontal Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<v'!U2!N>]
    725 set_titlex    =  'Local Time (h)'
    726 set_titley    =  'Altitude above surface (km)'
    727 set_subtitle  =  '' ;'Mean over the simulation domain'
    728 set_xrange    =  [8.,17.]
    729 set_yrange    =  [0.,1.] ;[0.,7.] ;; [0.,200.]
    730 set_tickx     =  1.
    731 set_ticky     =  0.1 ;1. ;50.
    732 minval        =  0.
    733 maxval        =  12. ;10.
    734 nlev          =  maxval-minval
    735 pal           =  22
    736 rrr           =  'no'
    737 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    738 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    739 !P.Charsize = 1.2
    740 ;; 0. levels
    741 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
    742 ;;; 1. background
    743 loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    744 ;; 2. color field
    745 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
    746             contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    747 ;; 3. contour field
    748 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
    749 ;; 4. choose output
    750 if (saveps eq 'true') then PS_End, /PNG
    751 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    752 
    753 if (n_elements(wmax) eq 0) then stop
     752;
     753;if (n_elements(wmax) eq 0) then stop
    754754
    755755;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    762762set_titley    =  'Altitude above surface (km)'
    763763set_subtitle  =  ''
    764 set_xrange    =  [8.,18.]
    765 set_yrange    =  [0.,7.]
     764set_xrange    =  [11.,17.] ;[8.,18.]
     765set_yrange    =  [0.,8.] ;[0.,8.]
    766766set_tickx     =  1.
    767767set_ticky     =  1.
    768768minval        =  0.
    769 maxval        =  16. ;15.;14.;13.
     769maxval        =  18. ;16. ;15.;14.;13.
    770770nlev          =  maxval-minval
    771771pal           =  22
    … …  
    776776;; 0. levels
    777777lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     778levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)]
    778779;;; 1. background
    779780loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    … …  
    782783              contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    783784;; 3. contour field
    784 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     785loadct, 0 & contour, zefield, zex, zey, levels=levhalf, c_labels=findgen(n_elements(levhalf))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (levhalf LT 0.0)
    785786;; 4. choose output
    786787if (saveps eq 'true') then PS_End, /PNG
    … …  
    796797set_titley    =  'Altitude above surface (km)'
    797798set_subtitle  =  ''
    798 set_xrange    =  [8.,18.]
    799 set_yrange    =  [0.,7.]
     799set_xrange    =  [11.,17.] ;[8.,18.]
     800set_yrange    =  [0.,8.]
    800801set_tickx     =  1.
    801802set_ticky     =  1.
    802803minval        =  0.
    803 maxval        =  9.;12.
     804maxval        =  10.;12.
    804805nlev          =  maxval-minval
    805806pal           =  22
    … …  
    810811;; 0. levels
    811812lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     813levhalf = minval + (maxval-minval)*findgen((nlev/2)+1)/float(nlev/2) & if (minval ne 0.) then levhalf = levhalf[where(levhalf ne 0.)]
    812814;;; 1. background
    813815loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    … …  
    816818            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    817819;; 3. contour field
    818 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     820loadct, 0 & contour, zefield, zex, zey, levels=levhalf, c_labels=findgen(n_elements(levhalf))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (levhalf LT 0.0)
    819821;; 4. choose output
    820822if (saveps eq 'true') then PS_End, /PNG
    821823;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    822 
    823 
    824 restore, filename='addturb.dat'
    825 velmax = SMOOTH(TEMPORARY(velmax),  [0,smoothampl], /EDGE_TRUNCATE)
    826 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    827 zefield       =  transpose(velmax)
    828 zex           =  localtime
    829 zey           =  h
    830 set_name      =  'velmax.ps'
    831 set_title     =  "Maximum horizontal wind speed (m.s!U-1!N)"
    832 set_titlex    =  'Local Time (h)'
    833 set_titley    =  'Altitude above surface (km)'
    834 set_subtitle  =  '' ;'Mean over the simulation domain'
    835 set_xrange    =  [8.,17.]
    836 set_yrange    =  [0.,4.]
    837 set_tickx     =  1.
    838 set_ticky     =  1.
    839 minval        =  0.
    840 maxval        =  30. ;12.
    841 nlev          =  maxval-minval
    842 pal           =  22
    843 rrr           =  'no'
    844 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    845 if (saveps eq 'true') then PS_Start, FILENAME=set_name
    846 !P.Charsize = 1.2
    847 ;; 0. levels
    848 lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
    849 ;;; 1. background
    850 loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
    851 ;;; 2. color field
    852 loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
    853             contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
    854 ;; 3. contour field
    855 loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
    856 ;; 4. choose output
    857 if (saveps eq 'true') then PS_End, /PNG
    858 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    859 
    860 
    861 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    862 ;zefield       =  localtime
    863 ;zey           =  depressions
    864 ;set_name      =  'DEPR.ps'
    865 ;set_title     =  "Percentage of depressions in the domain"
     824;
     825;
     826;restore, filename='addturb.dat'
     827;velmax = SMOOTH(TEMPORARY(velmax),  [0,smoothampl], /EDGE_TRUNCATE)
     828;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     829;zefield       =  transpose(velmax)
     830;zex           =  localtime
     831;zey           =  h
     832;set_name      =  'velmax.ps'
     833;set_title     =  "Maximum horizontal wind speed (m.s!U-1!N)"
    866834;set_titlex    =  'Local Time (h)'
    867 ;set_titley    =  'Percentage'
    868 ;set_subtitle  =  ''
    869 ;set_xrange    =  [8.,18.]
    870 ;set_yrange    =  [0.,3.]
     835;set_titley    =  'Altitude above surface (km)'
     836;set_subtitle  =  '' ;'Mean over the simulation domain'
     837;set_xrange    =  [11.,17.] ;[8.,17.]
     838;set_yrange    =  [0.,4.]
    871839;set_tickx     =  1.
    872840;set_ticky     =  1.
     841;minval        =  0.
     842;maxval        =  30. ;12.
     843;nlev          =  maxval-minval
     844;pal           =  22
     845;rrr           =  'no'
    873846;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    874847;if (saveps eq 'true') then PS_Start, FILENAME=set_name
    875848;!P.Charsize = 1.2
    876 ;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
    877 ;oplot, zefield, zey, psym=5
    878 ;       oplot, localtime, abs(psmin), linestyle=1
     849;;; 0. levels
     850;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
     851;;;; 1. background
     852;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0
     853;;;; 2. color field
     854;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
     855;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
     856;;; 3. contour field
     857;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0)
     858;;; 4. choose output
    879859;if (saveps eq 'true') then PS_End, /PNG
    880860;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     861;
     862;
     863;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     864;;zefield       =  localtime
     865;;zey           =  depressions
     866;;set_name      =  'DEPR.ps'
     867;;set_title     =  "Percentage of depressions in the domain"
     868;;set_titlex    =  'Local Time (h)'
     869;;set_titley    =  'Percentage'
     870;;set_subtitle  =  ''
     871;;set_xrange    =  [8.,18.]
     872;;set_yrange    =  [0.,3.]
     873;;set_tickx     =  1.
     874;;set_ticky     =  1.
     875;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     876;;if (saveps eq 'true') then PS_Start, FILENAME=set_name
     877;;!P.Charsize = 1.2
     878;;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
     879;;oplot, zefield, zey, psym=5
     880;;      oplot, localtime, abs(psmin), linestyle=1
     881;;if (saveps eq 'true') then PS_End, /PNG
     882;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    881883
    882884stop
Note: See TracChangeset for help on using the changeset viewer.