Changeset 1181 for trunk

Timestamp:

Feb 18, 2014, 11:19:52 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE. moving manual and notes into main repository.

Location:

trunk

Files:

• MESOSCALE/LMD_MM_MARS/SRC/DEV (added)
• MESOSCALE/LMD_MM_MARS/SRC/DEV/BENCH (moved) (moved from trunk/MESOSCALE_DEV/BENCH)
• MESOSCALE/LMD_MM_MARS/SRC/DEV/NOTES.txt (moved) (moved from trunk/MESOSCALE_DEV/NOTES.txt) (1 diff)
• MESOSCALE/LMD_MM_MARS/SRC/SCRIPTS/MPI (moved) (moved from trunk/MESOSCALE_DEV/MPI)
• MESOSCALE/LMD_MM_MARS/SRC/SCRIPTS/NETCDF (moved) (moved from trunk/MESOSCALE_DEV/NETCDF)
• MESOSCALE/MANUAL (moved) (moved from trunk/MESOSCALE_DEV/MANUAL)
• MESOSCALE/MANUAL/LMD_MM_MARS_User_Manual_previous.pdf (deleted)
• MESOSCALE/MANUAL/WP11_1_Mesoscale_Model_Design.pdf (deleted)
• MESOSCALE_DEV/PLOT/SPEC/LES/getturb.pro (modified) (9 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/NOTES.txt

@@ +1 @@
+
+
+bad synchronization physics / dynamics with LES
+--> si on a fait un update, il faut reinstaller les SOURCES du LES.
+
 
 si le systeme est sensible il vaut mieux r_aspect=1 et nsplit = 4 que

trunk/MESOSCALE_DEV/PLOT/SPEC/LES/getturb.pro

@@ -20 +20 @@
 ;p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191.
 ;print, 'ATTENTION ATTENTION R/cp !!!!', r_cp 
+
+; INTERCOMP INTERCOMP
+r_cp = 192./770.
 
 ;
@@ -32 +35 @@
 endif else begin
    ;PREF_SET, 'IDL_PATH', '/home/spiga/Save/SOURCES/IDL/fsc_psconfig:<IDL_DEFAULT>', /COMMIT
-   PREF_SET, 'IDL_PATH', '/home/spiga/SVN/trunk/mesoscale/PLOT/MINIMAL:<IDL_DEFAULT>', /COMMIT
+   ;PREF_SET, 'IDL_PATH', '/home/spiga/SVN/trunk/mesoscale/PLOT/MINIMAL:<IDL_DEFAULT>', /COMMIT
+   PREF_SET, 'IDL_PATH', '/home/spiga/MODELES/MESOSCALE_DEV/PLOT/MINIMAL:<IDL_DEFAULT>', /COMMIT
 endelse
 
@@ -172 +176 @@
 
 ;@include_ustar.pro
-
 ;@include_hfx.pro
-
 ;@include_w.pro
-
 ;@include_tprime.pro
 
@@ -374 +375 @@
 ;;if (saveps eq 'true') then PS_End, /PNG
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  localtime
-;zey           =  pbl_height
-;set_name      =  'HEIGHT.ps'
-;set_title     =  'Boundary layer height (km)'
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Boundary layer height (km)'
-;set_subtitle  =  '' ;'Criterion is : static stability > 1.5 K.m!U-1!N'
-;set_xrange    =  [11.,17.] ;[8.,17.]
-;set_yrange    =  [0.,8.]
-;set_tickx     =  1.
-;set_ticky     =  1.
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
-;;oplot, zefield, zey, psym=5
-;if (saveps eq 'true') then PS_End, /PNG
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
+
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+zefield       =  localtime
+zey           =  pbl_height
+set_name      =  'HEIGHT.ps'
+set_title     =  'Boundary layer height (km)'
+set_titlex    =  'Local Time (h)'
+set_titley    =  'Boundary layer height (km)'
+set_subtitle  =  '' ;'Criterion is : static stability > 1.5 K.m!U-1!N'
+set_xrange    =  [7.,19.] ;[11.,17.] ;[8.,17.]
+set_yrange    =  [0.,9.]
+set_tickx     =  1.
+set_ticky     =  1.
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+if (saveps eq 'true') then PS_Start, FILENAME=set_name
+!P.Charsize = 1.2
+plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin
+;oplot, zefield, zey, psym=5
+if (saveps eq 'true') then PS_End, /PNG
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;zefield       =  localtime
@@ -437 +438 @@
 ;;if (saveps eq 'true') then PS_End, /PNG
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  transpose(tke) 
-;zex           =  localtime 
-;zey           =  h
-;set_name      =  'TKE.ps'
-;set_title     =  "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Altitude above surface (km)'
-;set_subtitle  =  '' ;'Mean over the simulation domain'
-;set_xrange    =  [11.,17.] ;[8.,17.]
-;set_yrange    =  [0.,8.]  ;; [0.,200.]
-;set_tickx     =  1.
-;set_ticky     =  1. ;; 50.
-;minval        =  0.
-;maxval        =  20.
-;nlev          =  maxval-minval
-;pal           =  22
-;rrr           =  'no'
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;;; 0. levels
-;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
-;;;; 1. background
-;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 
-;  ;; 2. color field
-;  loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 
-;              contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
-;;; 3. contour field
-;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)
-;;; 4. choose output
-;if (saveps eq 'true') then PS_End, /PNG
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
+
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+zefield       =  transpose(tke) 
+zex           =  localtime 
+zey           =  h
+set_name      =  'TKE.ps'
+set_title     =  "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<u'!U2!N>+<v'!U2!N>+<w'!U2!N>]
+set_titlex    =  'Local Time (h)'
+set_titley    =  'Altitude above surface (km)'
+set_subtitle  =  '' ;'Mean over the simulation domain'
+set_xrange    =  [7.,19.] ;[8.,17.]
+set_yrange    =  [0.,9.]  ;; [0.,200.]
+set_tickx     =  1.
+set_ticky     =  1. ;; 50.
+minval        =  0.
+maxval        =  20.
+nlev          =  maxval-minval
+pal           =  22
+rrr           =  'no'
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+if (saveps eq 'true') then PS_Start, FILENAME=set_name
+!P.Charsize = 1.2
+;; 0. levels
+lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
+;;; 1. background
+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 
+  ;; 2. color field
+  loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) 
+              contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
+;; 3. contour field
+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)
+;; 4. choose output
+if (saveps eq 'true') then PS_End, /PNG
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;zefield       =  2.*transpose(ztke)
@@ -478 +479 @@
 ;;maxval        =  12.
 ;;pal           =  0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  transpose(ztke)
-;zex           =  localtime
-;zey           =  h
-;set_name      =  'zTKE.ps'
-;set_title     =  "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>]
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Altitude above surface (km)'
-;set_subtitle  =  '' ;'Mean over the simulation domain'
-;set_xrange    =  [8.,17.]
-;set_yrange    =  [0.,7.]
-;set_tickx     =  1.
-;set_ticky     =  1.
-;minval        =  0.
-;maxval        =  8.
-;nlev          =  maxval-minval
-;pal           =  22
-;rrr           =  'no'
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;;; 0. levels
-;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
-;;;; 1. background
-;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
-;;; 2. color field
-;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
-;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
-;;; 3. contour field
-;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)
-;;; 4. choose output
-;if (saveps eq 'true') then PS_End, /PNG
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+zefield       =  transpose(ztke)
+zex           =  localtime
+zey           =  h
+set_name      =  'zTKE.ps'
+set_title     =  "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[<w'!U2!N>]
+set_titlex    =  'Local Time (h)'
+set_titley    =  'Altitude above surface (km)'
+set_subtitle  =  '' ;'Mean over the simulation domain'
+set_xrange    =  [7.,19.]
+set_yrange    =  [0.,9.]
+set_tickx     =  1.
+set_ticky     =  1.
+minval        =  0.
+maxval        =  8.
+nlev          =  maxval-minval
+pal           =  22
+rrr           =  'no'
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+if (saveps eq 'true') then PS_Start, FILENAME=set_name
+!P.Charsize = 1.2
+;; 0. levels
+lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
+;;; 1. background
+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
+;; 2. color field
+loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
+            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
+;; 3. contour field
+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)
+;; 4. choose output
+if (saveps eq 'true') then PS_End, /PNG
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
 ;;restore, filename='addturb.dat'
 ;;modvar = SMOOTH(TEMPORARY(modvar),  [0,smoothampl], /EDGE_TRUNCATE)
@@ -547 +548 @@
 ;;if (saveps eq 'true') then PS_End, /PNG
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  transpose(wt)
-;zex           =  localtime
-;zey           =  h
-;set_name      =  'HF.ps'
-;set_title     =  "Turbulent Heat Flux (K.m.s!U-1!N)"  ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'>
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Altitude above surface (km)'
-;set_subtitle  =  '' ;'Mean over the simulation domain'
-;set_xrange    =  [8.,17.]
-;set_yrange    =  [0.,8.]  ;; [0.,200.]
-;set_tickx     =  1.
-;set_ticky     =  1. ;; 50.
-;minval        =  -1.5 ;-2.
-;maxval        =  2.5  ;2.
-;nlev          =  floor(maxval-minval)*10
-;pal           =  33  
-;rrr           =  'no'
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;minval        =  -0.8
-;;maxval        =  1.2
-;;pal           =  0
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;;; 0. levels
-;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
-;;;; 1. background
-;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
-;;; 2. color field
-;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
-;            ;;;--------------------------------------------------------------------------------------------------------------------------------
-;            ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone
-;            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
-;            nu = nd + (nu-nd)/2  ;; otherwise the interval is too large (because we removed 0)
-;            TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b
-;            ;;;--------------------------------------------------------------------------------------------------------------------------------
-;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
-;;; 3. contour field
-;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)
-;;; 4. choose output
-;if (saveps eq 'true') then PS_End, /PNG
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
+
+
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+zefield       =  transpose(wt)
+zex           =  localtime
+zey           =  h
+set_name      =  'HF.ps'
+set_title     =  "Turbulent Heat Flux (K.m.s!U-1!N)"  ;"Vertical Eddy Heat Flux <w'T'>" ;<w'!7h!3'>
+set_titlex    =  'Local Time (h)'
+set_titley    =  'Altitude above surface (km)'
+set_subtitle  =  '' ;'Mean over the simulation domain'
+set_xrange    =  [7.,19.]
+set_yrange    =  [0.,9.]  ;; [0.,200.]
+set_tickx     =  1.
+set_ticky     =  1. ;; 50.
+minval        =  -1.5 ;-2.
+maxval        =  2.5  ;2.
+nlev          =  floor(maxval-minval)*10
+pal           =  33  
+rrr           =  'no'
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;minval        =  -0.8
+;maxval        =  1.2
+;pal           =  0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+if (saveps eq 'true') then PS_Start, FILENAME=set_name
+!P.Charsize = 1.2
+;; 0. levels
+lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
+;;; 1. background
+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
+;; 2. color field
+loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
+            ;;;--------------------------------------------------------------------------------------------------------------------------------
+            ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone
+            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
+            nu = nd + (nu-nd)/2  ;; otherwise the interval is too large (because we removed 0)
+            TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b
+            ;;;--------------------------------------------------------------------------------------------------------------------------------
+            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
+;; 3. contour field
+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)
+;; 4. choose output
+if (saveps eq 'true') then PS_End, /PNG
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
 ;;;restore, filename='tpot_profB'
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
@@ -750 +751 @@
 ;if (saveps eq 'true') then PS_End, /PNG
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;if (n_elements(wmax) eq 0) then stop
+
+if (n_elements(wmax) eq 0) then stop
 
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
@@ -822 +823 @@
 if (saveps eq 'true') then PS_End, /PNG
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;
-;restore, filename='addturb.dat'
-;velmax = SMOOTH(TEMPORARY(velmax),  [0,smoothampl], /EDGE_TRUNCATE)
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;zefield       =  transpose(velmax)
-;zex           =  localtime
-;zey           =  h
-;set_name      =  'velmax.ps'
-;set_title     =  "Maximum horizontal wind speed (m.s!U-1!N)"
-;set_titlex    =  'Local Time (h)'
-;set_titley    =  'Altitude above surface (km)'
-;set_subtitle  =  '' ;'Mean over the simulation domain'
-;set_xrange    =  [11.,17.] ;[8.,17.]
-;set_yrange    =  [0.,4.]
-;set_tickx     =  1.
-;set_ticky     =  1.
-;minval        =  0.
-;maxval        =  30. ;12.
-;nlev          =  maxval-minval
-;pal           =  22
-;rrr           =  'no'
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;if (saveps eq 'true') then PS_Start, FILENAME=set_name
-;!P.Charsize = 1.2
-;;; 0. levels 
-;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
-;;;; 1. background
-;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
-;;;; 2. color field
-;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
-;            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
-;;; 3. contour field
-;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)
-;;; 4. choose output
-;if (saveps eq 'true') then PS_End, /PNG
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;
-;
+
+
+restore, filename='addturb.dat'
+velmax = SMOOTH(TEMPORARY(velmax),  [0,smoothampl], /EDGE_TRUNCATE)
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+zefield       =  transpose(velmax)
+zex           =  localtime
+zey           =  h
+set_name      =  'velmax.ps'
+set_title     =  "Maximum horizontal wind speed (m.s!U-1!N)"
+set_titlex    =  'Local Time (h)'
+set_titley    =  'Altitude above surface (km)'
+set_subtitle  =  '' ;'Mean over the simulation domain'
+set_xrange    =  [07.,19.] ;[8.,17.]
+set_yrange    =  [0.,4.]
+set_tickx     =  1.
+set_ticky     =  1.
+minval        =  0.
+maxval        =  30. ;12.
+nlev          =  maxval-minval
+pal           =  22
+rrr           =  'no'
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+if (saveps eq 'true') then PS_Start, FILENAME=set_name
+!P.Charsize = 1.2
+;; 0. levels 
+lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)]
+;;; 1. background
+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
+;;; 2. color field
+loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b)
+            contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill
+;; 3. contour field
+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)
+;; 4. choose output
+if (saveps eq 'true') then PS_End, /PNG
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
+
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;zefield       =  localtime