Ignore:
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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro

    r11 r613  
    1212;
    1313!p.multi=[0,2,2]
     14;!p.multi=[0,1,4]
     15
     16!p.charthick = 2.0
     17!p.thick = 5.0
     18!x.thick = 2.0
     19!y.thick = 2.0
     20
     21
     22
    1423set_plot, 'ps'
    15 device, file='plot.ps', /landscape
    16 tinv=5      ;; intervalle plot
     24device, file='plot.ps' ;, /landscape
     25
     26tinv=10     ;; intervalle plot
    1727hache=10.   ;; scale height
    1828psurf=610.  ;; pression surface -- ne change pas les eta levels
     
    6676; PLOT
    6777;
    68 plot, x, etas, title='ETA '+param, xtickinterval=tinv
    69 !psym=7 & oplot, x, etas
    70 !psym=0
    71 plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv
    72 !psym=7 & oplot, x, pseudo
    73 !psym=0
    74 plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv
    75 !psym=7 & oplot, x, pi
    76 !psym=0
     78;plot, x, etas, title='ETA '+param, xtickinterval=tinv
     79plot, x, etas, title='ETA LEVELS', xtickinterval=tinv, psym=3
     80;!psym=3 & oplot, x, etas
     81;!psym=0
     82plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv, psym=3
     83;!psym=3 & oplot, x, pseudo
     84;!psym=0
     85plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv, psym=3
     86;!psym=3 & oplot, x, pi
     87;!psym=0
    7788;
    7889diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2)
    7990;
    80 plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv
    81 !psym=7 & oplot, x, diff
    82 !psym=0
     91;plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv
     92;plot, x, ecart, title='d ALTITUDE (km)', xtickinterval=tinv, psym=1
     93plot, x, diff, title='d ALTITUDE (km)', xtickinterval=tinv, psym=3
     94;!psym=3 & oplot, x, diff
     95;!psym=0
    8396
    8497
Note: See TracChangeset for help on using the changeset viewer.