Changeset 613

Timestamp:

Apr 4, 2012, 4:36:50 PM (13 years ago)

Author:

aslmd

Message:

UTIL PYTHON: an example script for pp.py (whole ExoMars? report). minor corrections. MESOSCALE : dust storm def files

Location:

trunk

Files:

• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/callphys.def (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/localstorm.def (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/namelist.input (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/namelist.wps (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/run.def (added)
• MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_paper/stress.def (added)
• MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro (modified) (2 diffs)
• UTIL/PYTHON/myplot.py (modified) (1 diff)
• UTIL/PYTHON/planetoplot.py (modified) (1 diff)
• UTIL/PYTHON/rapport_exomars.sh (added)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro

@@ -12 +12 @@
 ;
 !p.multi=[0,2,2]
+;!p.multi=[0,1,4]
+
+!p.charthick = 2.0
+!p.thick = 5.0 
+!x.thick = 2.0 
+!y.thick = 2.0
+
+
+
 set_plot, 'ps'
-device, file='plot.ps', /landscape
-tinv=5      ;; intervalle plot
+device, file='plot.ps' ;, /landscape
+
+tinv=10     ;; intervalle plot
 hache=10.   ;; scale height
 psurf=610.  ;; pression surface -- ne change pas les eta levels
@@ -66 +76 @@
 ; PLOT
 ;
-plot, x, etas, title='ETA '+param, xtickinterval=tinv 
-!psym=7 & oplot, x, etas
-!psym=0
-plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv
-!psym=7 & oplot, x, pseudo 
-!psym=0
-plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv
-!psym=7 & oplot, x, pi 
-!psym=0
+;plot, x, etas, title='ETA '+param, xtickinterval=tinv 
+plot, x, etas, title='ETA LEVELS', xtickinterval=tinv, psym=3
+;!psym=3 & oplot, x, etas
+;!psym=0
+plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv, psym=3
+;!psym=3 & oplot, x, pseudo 
+;!psym=0
+plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv, psym=3
+;!psym=3 & oplot, x, pi 
+;!psym=0
 ;
 diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2)
 ;
-plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv
-!psym=7 & oplot, x, diff
-!psym=0
+;plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv
+;plot, x, ecart, title='d ALTITUDE (km)', xtickinterval=tinv, psym=1
+plot, x, diff, title='d ALTITUDE (km)', xtickinterval=tinv, psym=3
+;!psym=3 & oplot, x, diff
+;!psym=0
 
 

trunk/UTIL/PYTHON/myplot.py

@@ -318 +318 @@
         subv = 3
         subh = 1
-        fig.subplots_adjust(wspace = 0.5)
+        fig.subplots_adjust(hspace = 0.5)
         rcParams['font.size'] = int( rcParams['font.size'] * 1. / 2. )
     elif   numplot == 4:

trunk/UTIL/PYTHON/planetoplot.py

@@ -244 +244 @@
               elif axtime == "sol":   time = soltab
           else:
-              if "Times" in nc.variables:   time = count + np.arange(0,len(nc.variables["Times"]),1)
-              elif "Time" in nc.variables:  time = count + np.arange(0,len(nc.variables["Time"]),1)
-              else:                         time = count + np.arange(0,1,1)
+              if "Times" in nc.dimensions:   time = count + np.arange(0,len(nc.dimensions["Times"]),1)
+              elif "Time" in nc.dimensions:  time = count + np.arange(0,len(nc.dimensions["Time"]),1)
+              else:                          time = count + np.arange(0,1,1)
               if ope in ["cat"]:
                  if nnn > 0:  count = time[-1] + 1  ## so that a cat is possible with simple subscripts