Changeset 613 for trunk/MESOSCALE/LMD_MM_MARS/SIMU

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/MESOSCALE/LMD_MM_MARS/SIMU

Files:

• DEF/DUST_STORM_paper (added)
• DEF/DUST_STORM_paper/callphys.def (added)
• DEF/DUST_STORM_paper/localstorm.def (added)
• DEF/DUST_STORM_paper/namelist.input (added)
• DEF/DUST_STORM_paper/namelist.wps (added)
• DEF/DUST_STORM_paper/run.def (added)
• DEF/DUST_STORM_paper/stress.def (added)
• RUN/vert_level/levspe.pro (modified) (2 diffs)

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