Changeset 287

Timestamp:

Sep 8, 2011, 7:14:13 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: minor updates to parameter files, plotting python module, launching and post-processing

Location:

trunk/MESOSCALE/LMD_MM_MARS

Files:

• SIMU/DEF/NEWPHYS_TRACERS/namelist.input (modified) (1 diff)
• SIMU/RUN/launch.gnome (modified) (2 diffs)
• SRC/POSTPROC/api.F90 (modified) (1 diff)
• SRC/POSTPROC/namelist.api.minim (added)
• SRC/PYTHON/myplot.py (modified) (2 diffs)
• SRC/PYTHON/winds.py (modified) (12 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input

@@ -44 +44 @@
  tiso = 100.
  pd_scalar = T    !!! tres IMPORTANT, sinon valeurs negatives 
+ diff_6th_opt = 0
  !!!
  !!! OPTIONAL

trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome

@@ -34 +34 @@
 cp      $LOADL_STEP_INITDIR/*.exe       .
 cp      $LOADL_STEP_INITDIR/namelist*   .
-cp      $LOADL_STEP_INITDIR/wrfi*       .       ### executer sur la frontale
-                                                ### au prealable !!!
+#cp     $LOADL_STEP_INITDIR/wrfi*       .       ### executer sur la frontale
+#                                               ### au prealable !!!
+ln -sf  $LOADL_STEP_INITDIR/wrfi*       .
 #######################################################################
 if [[ "${real}" == "0" ]]
@@ -46 +47 @@
   ##mv rsl.* ideal/
 else
-  cp    $LOADL_STEP_INITDIR/wrfb* .
+#  cp $LOADL_STEP_INITDIR/wrfb* .
+ln -sf $LOADL_STEP_INITDIR/wrfb* .
 fi
 #######################################################################

trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/api.F90

@@ -69 +69 @@
       output_name     = ' '
       interp_levels   = -99999.
-      process         = 'all'
+      process         = 'list' !!'all'
       
 

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py

@@ -487 +487 @@
              "anomaly":      "%.1f",\
              "W":            "%.1f",\
+             "WMAX_TH":      "%.1f",\
+             "QSURFICE":     "%.0f",\
                     }
     if whichvar not in fmtvar:
@@ -507 +509 @@
              "VMR_ICE":      "Blues",\
              "W":            "jet",\
+             "WMAX_TH":      "spectral",\
              "anomaly":      "RdBu_r",\
+             "QSURFICE":     "hot_r",\
                      }
 #W --> spectral ou jet

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/winds.py

@@ -31 +31 @@
            var2=None,\
            ndiv=10,\
-           first=1):
+           first=1,\
+           mult=1.,\
+           zetitle="fill"):
 
     ####################################################################################################################
@@ -48 +50 @@
     import numpy as np
     from numpy.core.defchararray import find
-
-    #if save == 'eps': rcParams['backend'] = 'PS'
 
     ######################
@@ -87 +87 @@
     ### Time loop for plotting device
     found_lct = False
+    nplot = 1
     itime = first
-    nplot = 1
     error = False
     if itstep is None and numplot > 0: itstep = int(24./numplot)
@@ -126 +126 @@
                if var2 == 'HGT':  zelevels = np.arange(-10000.,30000.,2000.)
                contour( x, y, what_I_contour, zelevels, colors='k', linewidths = 0.33 ) #colors='w' )# , alpha=0.5)
+           else:
+               errormess("There is an error in reducing field !")
 
        #### Shaded plot
        if var:
            what_I_plot, error = reducefield( getfield(nc,var), d4=itime, d3=nvert )
+           what_I_plot = what_I_plot*mult
            if not error: 
                fvar = var
@@ -136 +139 @@
                    what_I_plot = 100. * ((what_I_plot / smooth(what_I_plot,12)) - 1.)
                    fvar = 'anomaly'
+               #if mult != 1:     
+               #    fvar = str(mult) + "*" + var
                ###
                if typefile in ['mesoapi','meso']:    what_I_plot = dumpbdy(what_I_plot,6)
@@ -143 +148 @@
                if not tile:
                    if not hole: what_I_plot = bounds(what_I_plot,zevmin,zevmax)
-                   zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)
+                   #zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)
+                   zelevels = np.linspace(zevmin,zevmax)
                    contourf( x, y, what_I_plot, zelevels, cmap = palette )
                else:
@@ -154 +160 @@
                                            extend='neither',spacing='proportional')
                                            # both min max neither 
-           
+           else:
+               errormess("There is an error in reducing field !")
+ 
        ### Vector plot
        if winds:
@@ -179 +187 @@
             if addchar:  plottitle = plottitle + addchar + "_LT"+str(ltst)
             else:        plottitle = plottitle + "_LT"+str(ltst)
+       if mult != 1:           plottitle = str(mult) + "*" + plottitle
+       if zetitle != "fill":   plottitle = zetitle
        ptitle( plottitle )
        itime += itstep
@@ -186 +196 @@
     ### Save the figure in a file in the data folder or an user-defined folder
     if typefile in ['meso','mesoapi']:   prefix = getprefix(nc)   
-    elif typefile in ['gcm','gcmex']:            prefix = 'LMD_GCM_'
+    elif typefile in ['gcm','gcmex']:    prefix = 'LMD_GCM_'
     else:                                prefix = ''
     ###
@@ -208 +218 @@
             print "save mode not supported. using gui instead."
             show()
-    else:             print "Local time not found"
+    else:   print "Local time not found"
 
     ###############
@@ -263 +273 @@
     parser.add_option('-S', '--save',   action='store',dest='save',      type="string",  default="gui", help='save mode (png,eps,svg,pdf or gui)(def=gui)')
     parser.add_option('-a', '--anomaly',action='store_true',dest='anomaly',              default=False, help='compute and plot relative anomaly in %')
-    parser.add_option('-w', '--with',   action='store',dest='var2',     type="string",   default=None,  help='variable contoured')
-    parser.add_option('--div',          action='store',dest='ndiv',     type="int",      default=10,    help='number of divisions in colorbar (def: 10)')
-    parser.add_option('-F', '--first',  action='store',dest='first',    type="int",      default=1,     help='first subscript to plot (def: 1)')
+    parser.add_option('-w', '--with',   action='store',dest='var2',      type="string",  default=None,  help='variable contoured')
+    parser.add_option('--div',          action='store',dest='ndiv',      type="int",     default=10,    help='number of divisions in colorbar (def: 10)')
+    parser.add_option('-F', '--first',  action='store',dest='first',     type="int",     default=1,     help='first subscript to plot (def: 1)')
+    parser.add_option('--mult',         action='store',dest='mult',      type="float",   default=1.,    help='a multiplicative factor to plotted field')
+    parser.add_option('--title',        action='store',dest='zetitle',   type="string",  default="fill",help='customize the whole title')
     #parser.add_option('-V', action='store', dest='comb',        type="float",   default=None,  help='a defined combination of variables')
     (opt,args) = parser.parse_args()
@@ -340 +352 @@
                 tile=opt.tile,zoom=opt.zoom,display=opt.display,\
                 itstep=opt.itstep,hole=opt.hole,save=opt.save,\
-                anomaly=opt.anomaly,var2=opt.var2,ndiv=opt.ndiv,first=opt.first)
+                anomaly=opt.anomaly,var2=opt.var2,ndiv=opt.ndiv,first=opt.first,\
+                mult=opt.mult,zetitle=opt.zetitle)
             print 'Done: '+name
             system("rm -f to_be_erased")