Changeset 114 for trunk/mesoscale/PLOT/SPEC/SERIES

Timestamp:

Apr 22, 2011, 8:14:18 PM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: update graphic tools for GW, water cycle + generic graphic stuff (ps_start, map_latlon) + update notes

Location:

trunk/mesoscale/PLOT/SPEC/SERIES

Files:

• mawd.pro (modified) (4 diffs)
• mawd_inc.pro (modified) (3 diffs)

trunk/mesoscale/PLOT/SPEC/SERIES/mawd.pro

@@ -2 +2 @@
 
 ;;;;;;;
-;f_user='/donnees/aslmd/SVN/trunk/mesoscale/TMPDIR/GCMINI/mawd.nc'
-;f_user='/tmp7/aslmd/mawd.nc'
-f_user='/d5/aslmd/mawd.nc'
+experiment='a035'
+;experiment='a040'
+;experiment='a045'
+;;f_user='/donnees/aslmd/SVN/trunk/mesoscale/TMPDIR/GCMINI/mawd.nc'
+;;f_user='/tmp7/aslmd/mawd.nc'
+f_user='/d5/aslmd/WATERCYCLE/mawd_'+experiment+'.nc'
 v_user = 'mtot'
 x_user = 'Time'
@@ -26 +29 @@
 SPAWN, '\rm param_plot.idl ; cp mawd_inc.pro param_plot.idl'
 
-PS_START, file='mawd.ps'
-  !P.Charsize = 1.2
-  !p.charthick = 2.0
-  !p.thick = 2.0
-  !x.thick = 2.0
-  !y.thick = 2.0
+PS_START, file='mawd_'+experiment+'.ps'
+;  !P.Charsize = 1.2
+;  !p.charthick = 2.0
+;  !p.thick = 2.0
+;  !x.thick = 2.0
+;  !y.thick = 2.0
 
 cdfid = ncdf_open(f_user)
@@ -41 +44 @@
 ncdf_varget, cdfid, varid, yy
 
+varid=ncdf_varid(cdfid,x_user)
+ncdf_varget, cdfid, varid, xx
+
 if (tes eq 'no') then begin
   ;;; en precip-microns
   champ = champ * 1.e6 / 917.
-  ;;; prendre Time depuis un autre fichier
-  cdfid = ncdf_open('/d5/aslmd/gcm_LT14_a035.nc')
-  varid=ncdf_varid(cdfid,x_user)
-  ncdf_varget, cdfid, varid, xx
   ;;;; entre 0 et 360
   ;xx = xx MOD 360
+  xx = xx - 360.*3.
 endif else begin
-  varid=ncdf_varid(cdfid,x_user)
-  ncdf_varget, cdfid, varid, xx
   xx = xx - 360.*2.
 endelse
@@ -76 +77 @@
 endif
 
+what_I_plot = smooth(what_I_plot,[10,1])
+
 map_latlon, $
         what_I_plot, $                          ; 2D field

trunk/mesoscale/PLOT/SPEC/SERIES/mawd_inc.pro

@@ -8 +8 @@
         format          =       '(I0)' 
         colors          =       32
-        pal             =       22;4;33              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
+        pal             =       33;22;4;33              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
 
 ;;;;;*************************************** FILL LIMITS
         minfield_init   =       0.
-        maxfield_init   =       65.
-        ndiv            =       13
+        maxfield_init   =       100. ;120. ;65.
+        ndiv            =       10 ;12 ;13
 
         ;;;;;*************************************** LIMIT TRICKS
@@ -20 +20 @@
                 lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max
                 lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
-                ;lim_blank = 0.1 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
+                lim_blank = 1.0 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
 
 ;;;;;*************************************** CONTOUR
@@ -38 +38 @@
 ;;;;; MULTIYEAR
         isotropic       =       'false'
-        windowx         =       [1500.,4000.]
-        intervalx       =       500.0
+        windowx         =       [0.,2880.] ;[0.,3600.] ;[1500.,4000.]
+        intervalx       =       360.0
         windowy         =       [-90.,90.]
         intervaly       =       30.0
-        lev             =       [100.]
+        lev             =       [600.] ;[100.]