Changeset 1025 in lmdz_wrf for trunk

Timestamp:

Aug 16, 2016, 8:53:49 PM (9 years ago)

Author:

lfita

Message:

Getting almost there for the first figure!

Location:

trunk/tools

Files:

• model_graphics.py (modified) (16 diffs)
• model_graphics_template.dat (modified) (5 diffs)

trunk/tools/model_graphics.py

@@ -7 +7 @@
 import nc_var_tools as ncvar
 import generic_tools as gen
+import drawing_tools as drw
 import time as tim
 #  To avoid errors like: '/bin/sh: 1: Syntax error: Bad fd number'
@@ -289 +290 @@
             #   [operation]: operation to compute
             #   [fheader]: header of file with the required variables
-            #   [vmodel]: 
-            #   [vdiag]:
+            #   [vmodel]: Variable from model output
+            #   [vdiag]: Varible as diagnostic from different variables of the model output
             #   [globalP]: Wether vertical interpolation is:
             #     'global': for all file
@@ -320 +321 @@
 
     # Outwritting the varcompute to avoid next time (if it is not filescratch!)
-    objf = open(owdir + '/varcompute.inf', 'w')
+    objf = open(odir + '/varcompute.inf', 'w')
     objf.write('files: ' + strfiles + '\n')
     objf.write('testfiles: ' + strtestfiles + '\n')
     objf.write('varcompute: ' + Svarcompute + '\n')
     objf.write('itotv: ' + str(ivop) + '\n')
-    objf.close()
+    objf.close() 
 
     return allvarcomp, ivop
@@ -365 +366 @@
     return Files, TestFiles, allvarcomp, ivop
 
+def pinterpS(gP):
+    """ For that variables which require vertical interpolation 'p' suffix to the file header is added
+      gP= kind of pinterp: 'global', 'local', 'none'
+    """
+    fname = 'pinterpS'
+
+    if gP != 'none':
+        gPS = 'p'
+    else:
+        gPS = ''
+
+    return gPS
+
 def compute_variable(minf, idir, usefiles, odir, cvar, gP, scr, pyH, Tref, Tunits, db):
     """ Function to compute a variable
@@ -401 +415 @@
     # For that variables which require vertical interpolation 'p' suffix to the 
     #   file header is added
-    if gP != 'none':
-        SgP = 'p'
-    else:
-        SgP = ''
+    SgP = pinterpS(gP)
 
     # Getting in working dir
@@ -435 +446 @@
             if modvar is not None and diagvar is None:
                 # model variable
-                values = modvar + ',0,-1,-1'
-                vs = modvar + ',' + vdnx + ',' + vdny + ',' + vdnz + ',' + vdnt
+                values = modvar[0] + ',0,-1,-1'
+                vs = modvar[0] + ',' + vdnx + ',' + vdny + ',' + vdnz + ',' + vdnt
                 with gen.Capturing() as output:
                     ncvar.DataSetSection_multivars(values, cf, vs)
@@ -448 +459 @@
                   vs + " -S '" + values + "'"
                 otrackf.write('\n')
-                otrackf.write('# ' + CFvarn + " " + modvar + '\n')
+                otrackf.write('# ' + CFvarn + " " + modvar[0] + '\n')
                 otrackf.write(pyins + '\n')
 
                 # CF renaming of variable
-                ncvar.chvarname(CFvarn,ifilen,modvar)
+                ncvar.chvarname(CFvarn,ifilen,modvar[0])
             else:
                 # diagnostic variable
@@ -587 +598 @@
     Tunits = config['CFunitstime']
     #  opsur = operation surnames
+
     opsur = gen.get_specdictionary_HMT(config, H='opsur_',M='',T='')
     Opsurs = {}
@@ -593 +605 @@
         vals = opsur[opk]
         Opsurs[opn] = vals.split(':')
-    if db:
-        print '  ' + fname + ' operation surnames _______'
-        gen.printing_dictionary(Opsurs)
+#    if db:
+#        print '  ' + fname + ' operation surnames _______'
+#        gen.printing_dictionary(Opsurs)
 
     # Getting in working dir
@@ -976 +988 @@
     return
 
+def get_plots_var(vc,pltk,db):
+    """ Function to provide the plots to make for each variable from pltk ('DIRPLT_', or 'DIFFPLT_') dictionary
+      vc= pltk dictionary, dictionary with all the parameters started with pltk
+      pltk= kind of plots. Values started by 'DIRPLT_', or 'DIFFPLT_' in model configuration
+    """
+    fname = 'get_plots_var'
+
+    LH = len(pltk)
+    # list of plots
+    doplots = {}
+    # plots by variable
+    plotsvar = {}
+    # individual plots by variable (not repeated)
+    indivplotsvar = {}
+
+    # Variables with a plot which requires vertical interpolation
+    varplotp = []
+
+    ipl = 0
+    for plot in vc.keys():
+        Lplot = len(plot)
+        pln = plot[LH:Lplot+1]
+        varops = vc[plot].split(':')
+        for vn in varops:
+            VOn = vn.split('#')
+            # Variables and operations in plot
+            vplt = []
+            oplt = []
+            voplt = []
+            for vop in VOn:
+                vv = vop.split('|')[0]
+                oo = vop.split('|')[1]
+                vplt.append(vv)
+                oplt.append(oo)
+                if not gen.searchInlist(voplt,vop): voplt.append(vop)
+            vpltS = '#'.join(vplt)
+            opltS = '#'.join(oplt)
+            
+            if not doplots.has_key(pln):
+                doplots[pln] = [vn]
+            else:
+                plots = doplots[pln]
+                plots.append(vn)
+                doplots[pln] = plots
+
+            if not plotsvar.has_key(vn):
+                plotsvar[vn] = [pln]
+            else:
+                plots = plotsvar[vn]
+                plots.append(pln)
+                plotsvar[vn] = plots
+
+            for VOv in voplt:
+                if not indivplotsvar.has_key(VOv):
+                    indivplotsvar[VOv] = [pln]
+                else:
+                    plots = indivplotsvar[VOv]
+                    indivplots.append(pln)
+                    indivplotsvar[VOv] = plots
+
+        if vc[plot].find('pinterp') != -1:
+            if not gen.searchInlist(varplotp,vn): varplotp.append(pln)
+
+        ipl = ipl + 1
+
+    if db:
+        print '  plots to make _______'
+        gen.printing_dictionary(doplots)
+        print '  plots by variable|operation _______'
+        gen.printing_dictionary(plotsvar)
+        print '  individual plots by variable|operation _______'
+        gen.printing_dictionary(indivplotsvar)
+        print '  plots with the vertical interpolation _______'
+        print '    #', varplotp
+
+    return doplots, plotsvar, indivplotsvar, varplotp
+
+def plots_listconstruct(config, odir, debug):
+    """ Function to create the list of plots to draw
+      config= Configuration of the experiment
+      odir= output experiment folder
+    """
+    fname='plots_listconstruct'
+  
+    dirplot = gen.get_specdictionary_HMT(config, H='DIRPLT_')
+
+    if debug:
+        print '  direct plots to draw ________'
+        gen.printing_dictionary(dirplot)
+
+    # Getting plots by variable
+    plots, varplot, indivarplot, plotp = get_plots_var(dirplot,'DIRPLT_',debug)    
+
+    Svarplot = gen.dictKeysVals_stringList(plots,cV=':')
+
+    Nplots = 0
+    for pl in Svarplot.split(','):
+        Nplots = Nplots + len(pl.split(':'))
+
+    # Outwritting the varcompute to avoid next time (if it is not filescratch!)
+    objf = open(odir + '/directplotsdraw.inf', 'w')
+    objf.write('plots: ' + Svarplot + '\n')
+    objf.write('itotp: ' + str(Nplots) + '\n')
+    objf.close() 
+
+    return plots, Nplots
+
+def read_plot_file(figfile):
+    """ Function to read the file with the information about the  plots to draw
+      figfile= file with the information
+    """
+    fname = 'read_plot_file'
+
+    if not os.path.isfile(figfile):
+        print errormsg
+        print '  ' + fname + ": figures file '" + figfile + "' does not exist !!"
+        quit(-1)
+
+    objf = open(figfile, 'r')
+    for line in objf:
+        if line[0:1] != '#' and len(line) > 1:
+            values = line.split(' ')
+            if values[0] == 'plots:':
+                plots = gen.stringList_dictKeysVals(values[1],cV=':')
+            elif values[0] == 'itotp:':
+                Nplots = int(values[1])
+
+    objf.close()
+
+    return plots, Nplots
+
+def varnoper(vn, oper, opsurdict):
+    """ Function to provide name of the variable after operation as result of addition of the 'surnames'
+      vn= variable name
+      oper= '+' separated list of operations
+      opsurdict= dictionary with the surname of operations
+    >>> OpSurDict = {'mean': ['tmean', 'xmean', 'ymean']}
+    >>> varnoper('tas', 'pinterp+tmean+xmean', OpSurDict)
+    tasmeanmean
+    """
+    fname = 'varnoper'
+
+    newvn = vn
+
+    for op in oper.split('+'):
+        surname = gen.dictionary_key_list(opsurdict,op)
+        if surname is not None:
+            newvn = newvn + surname
+
+    return newvn
+
+def draw_plot(kplot, CFvplot, fplot, vplot, dplot, pplot, finame, tfig, kfig, mapval,\
+  tunits, od, pyH, fscr, db):
+    """ Function to draw a plot
+      kplot= kind of plot
+      CFvplot= list of the CFnames of the variables
+      fplot= list with the files in the plot
+      vplot= list with the name of the variables in each file
+      dplot= list of the dimensions in each file
+      pplot= list of pictoric characteristics for each variable
+      finame= name of the figure
+      tfig= title of the figure
+      kfig= kind of the figure
+      mapval= value of the map
+      tunits= units of time
+      od= output directory
+      pyH= python HOME
+      fscr= whether figure should be done from scratch
+    """
+    fname = 'draw_plot'
+
+    if fscr:
+        sout = sub.call('rm ' + finame + ' >& /dev/null', shell=True)
+
+    os.chdir(od)
+
+    # CF variable-dimensions
+    CFvardims = {'lon':'lon', 'lat': 'lat', 'pres': 'pres', 'time': 'time'}
+
+    # Tracking file with all the figures
+    trkobjf = open('all_figures.inf', 'a')
+
+    if not os.path.isfile(finame):
+        if db:
+            print "    drawing '" + finame + "' ..."
+            print '      CFvars:', CFvplot
+            print '      files:', fplot
+            print '      in file vars:', vplot
+            print '      dims:', dplot
+            print '      pictoric values:', pplot
+            print '      figure name:', finame
+            print '      title:', tfig.replace('!',' ')
+            print '      kind:', kfig
+            print '      map:', mapval
+            print '      time units:', tunits
+
+        if kplot == 'diffmap2Dsfc':
+            print "  " + fname + ": kind of plot '" + kplot + "' not ready !!"
+            quit(-1)
+        elif kplot == 'diffmap2Dz':
+            print "  " + fname + ": kind of plot '" + kplot + "' not ready !!"
+            quit(-1)
+        elif kplot == 'map2Dsfc':
+            print "  " + fname + ": kind of plot '" + kplot + "' not ready !!"
+            quit(-1)
+        elif kplot == 'map3D':
+            print "  " + fname + ": kind of plot '" + kplot + "' not ready !!"
+            quit(-1)
+        elif kplot == 'shadcont2Dsfc':
+            figtit = tfig.replace('!','|')
+            shdstdn = CFvplot[0]
+            cntstdn = CFvplot[1]
+            figfs = ','.join(fplot)
+            shad = pplot[0]
+            srange = str(shad[0]) + ',' + str(shad[1])
+            cbar = shad[2]
+            cnt = pplot[1]
+            crange = str(cnt[0]) + ',' + str(cnt[1])
+            cfmt = cnt[2]
+            cline = cnt[3]
+
+            graphvals = ','.join(CFvplot)
+            graphvals = graphvals + ':lon|-1,lat|-1:lon|-1,lat|-1:lon:lat:' + cbar + \
+              ':fixc,' + cline + ':' + cfmt + ':' + srange + ':' + crange + ',9:' +  \
+              figtit + ':' + kfig + ':None:' + mapval
+
+            fvarS = ','.join(vplot)
+
+            plotins = "python "+pyH+"/drawing.py -f "+figfs+" -o draw_2D_shad_cont "+\
+              "-S '" + graphvals + "' -v " + fvarS
+            try:
+                with gen.Capturing() as output:
+                    sout = sub.call(plotins, shell=True)
+            except:
+                print 'drawing.draw_2D_shad_cont(' + graphvals + ',' + figfs + ',' + \
+                  fvarS + ')'
+                for sout in output: print sout
+                quit(-1)
+
+            sout = sub.call('mv 2Dfields_shadow-contour.'+kfig+' '+finame, shell=True)
+
+            # keeping all figures
+            trkobjf.write('\n')
+            trkobjf.write("#" + tfig.replace('!',' ') + '\n')
+            trkobjf.write(plotins + '\n')
+        else:
+            print errmsg
+            print "  " + fname + ": plot kind '" + kplot + "' not ready !!"
+            quit(-1)
+
+    trkobjf.close()
+
+    return
+
+def draw_plots(config, plots, mod, exp, odir, allvarcomp, figscr, debug):
+    """ Function to draw all plots
+      config= Configuration of the experiment
+      plots= dictionary with the plots
+      mod= model of the plots
+      exp= experiment of the plots
+      odir= output experiment folder
+      allvarcomp= dictionary with all the variables to compute and their information
+      figscr= whether figures should be done from the scratch or not
+
+    * Plot as 
+      {[kplot]} = [varn1]|[op1]#[varn2]|[op2]#[...[varnN]|[opN]], ...
+        [kplot] ___ 
+          diffmap2Dsfc: 2D map of surface differences values of 1 variable
+          diffmap2Dz: 2D map of 3D differences values of 1 variable
+          map2Dsfc: 2D map of surface values of 1 variable
+          map3D: 2D map of 3D values of 1 variable
+          shadconthovmsfc: Hovmoeller diagrams of 2 variable at the surface in shadow and the other in contourn
+          shadcont2Dsfc: 2D map of shadow (1st variable) and countour (2nd variable) [stvar1]#[stvar2]
+          shadcont2Dzsec: 2D map of vertical section of 2 variables one in shadow and the other in contourn
+        [varn] 
+          variable
+        [op]
+          '+' separated list of operations
+        in figures with more than 1 variable, use '#' to separate the [varn]|[op]
+    """ 
+    fname = 'draw_plots'
+
+    os.chdir(odir)
+
+    # Dictionary with the operations with surnames for the operated variable
+    opersurnames = {}
+    opsur = gen.get_specdictionary_HMT(config, H='opsur_',M='',T='')
+    for opsr in opsur.keys():
+        opn = opsr.split('_')[1]
+        vls = opsur[opsr].split(':')
+        opersurnames[opn] = vls
+
+    # Units time for the plots
+    rd = config['CFreftime']
+    timeunits = config['CFunitstime'] + '!since!' + rd[0:4] + '-' + rd[4:6] + '-' +  \
+      rd[6:8] + '!' + rd[8:10] + ':' + rd[10:12] + ':' + rd[12:14]
+
+    # Dictinoary of plot specificities
+    #   [minval]: minimum value
+    #   [maxval]: minimum value
+    #   [colorbar]: name of the colorbar (from matplotlib) to use
+    #   [cntformat]: format of the contour labels
+    #   [colorcnt]: color for the countor lines
+    plotspecifics = {}
+    plotspecs = config['specificvarplot'].split(':')
+    for pltspc in plotspecs:
+        pltvls = pltspc.split('|')
+        vn = pltvls[0]
+        op = pltvls[1]
+        fn = pltvls[2]
+        plotspecifics[fn + '_' + vn + '_' + op] = pltvls[3:]
+    if debug:
+        print 'Specific values for plots _______'
+        gen.printing_dictionary(plotspecifics)
+
+    # Kind of figures
+    kindfigure = config['kindfig']
+
+    # Map value
+    mapvalue = config['mapval']
+
+    # pythone scripts HOME
+    pyHOME = config['pyHOME']
+
+    # Title-text of operations
+    opexplained = {}
+    optits = config['titleoperations'].split(':')
+    for optit in optits:
+        opn = optit.split('|')[0]
+        opt = optit.split('|')[1]
+        opexplained[opn] = opt
+    if debug:
+        print 'Titles for operations  _______'
+        gen.printing_dictionary(opexplained)
+
+    for kplot in plots.keys():
+        varsplt = plots[kplot]
+        for varplt in varsplt:
+            if debug:
+                print "  printing '" + kplot + "' var ':" + varplt + "'..."
+            varops = varplt.split('#')
+
+            # CF variables in plot
+            CFvarsplot = []
+            # Files in plot
+            filesplot = []
+            # Variables in plot within the files
+            varsplot = []
+            # Dims in figure
+            dimsplot = []
+            # pictoric values in figure
+            pictplot = []
+            # Name of the figure
+            figname = ''
+            # Title of the figure
+            titfigure = ''
+
+            ivp = 0
+            for varop in varops:
+                vn = varop.split('|')[0]
+                op = varop.split('|')[1]
+
+                # CF variables in plot
+                CFvarsplot.append(vn)
+ 
+                vnopS = vn + '_' + op
+                if not allvarcomp.has_key(vnopS):
+                    print errormsg
+                    print '  ' + fname + ": no entry for variable-operation '" +     \
+                      vnopS + "' !!"
+                vopvals = allvarcomp[vnopS]
+                headf = vopvals[0]
+                globalP = vopvals[3]
+                gP = pinterpS(globalP)
+
+                filen = vn + '_' + headf + gP + '_' + op + '.nc'
+                filesplot.append(filen)
+                # Do we have processed the given variable?
+                if not os.path.isfile(filen):
+                    print warnmsg
+                    print "  " + fname + ": there is no file for variable '" + varop \
+                      + "' skiping it !!"
+                    break
+
+                # Name of the variable inside the file
+                vnsur = varnoper(vn, op, opersurnames)
+                varsplot.append(vnsur)
+
+                # Dimensions in file
+                try:
+                    with gen.Capturing() as output:
+                        dims = ncvar.idims(filen)
+                except:
+                    print 'ncvar.idims('+filen+')'
+                    for sout in output: print sout
+                    quit(-1)
+
+                dimsplot.append(dims)
+
+                # pictoric values for the figure
+                Sfivaop = kplot + '_' + vn + '_' + op
+                if plotspecifics.has_key(Sfivaop):
+                    pictvals = plotspecifics[Sfivaop]
+                else:
+                    Vvals = gen.variables_values(vn)
+                    pictvals = [Vvals[2], Vvals[3], Vvals[6], '%g', 'black']
+
+                pictplot.append(pictvals)
+
+                # Header of the name of the figure
+                if ivp == 0:
+                    figname = kplot +'_'+ vn + '_' + headf + '_' + op.replace('+','-')
+                else:
+                    figname = figname + '-'+vn+'_'+headf+'_'+op.replace('+','-')
+
+                # Title of the figure:
+                if ivp == 0:
+                    titfigure = mod + '!' + exp + '!' + vn 
+                else:
+                    titfigure = titfigure + '!&!' + vn
+                for op1 in op.split('+'):
+                    if not opexplained.has_key(op1):
+                        print errormsg
+                        print '  '+fname+": no explanation for operation '"+op1+"' !!"
+                        print '    provided:', opexplained.keys()
+                        quit(-1)
+                    titfigure = titfigure + '!' + opexplained[op1]
+
+                ivp = ivp + 1
+            # End of variable-operation
+
+            draw_plot(kplot, CFvarsplot, filesplot, varsplot, dimsplot, pictplot,    \
+              figname, titfigure, kindfigure, mapvalue, timeunits, odir, pyHOME,     \
+              figscr, debug)
+
+        # End of variables-operations
+
+    # End of kind of plots
+
+    return
+
 # Files with information about the configuration of the script
 inffiles = ['varcompute.inf', 'all_computevars.inf', 'all_statsvars.inf']
@@ -996 +1449 @@
 for mod in mods:
     print mod
+    if scratch:
+        sout = sub.call('rm -rf '+cnf['ofold']+'/'+mod+' > /dev/null', shell=True)
+
     # Get experiments and headers of model
     exps, fheaders = exp_headers(mod,cnf)
@@ -1049 +1505 @@
             objf.close()
 
+        varcompf = owdir + '/varcompute.inf'
         if addfiles:
-            sub.call('rm ' + owdir +'/varcompute.inf >& /dev/null', shell=True)
-
-        varcompf = owdir + '/varcompute.inf'
+            sub.call('rm ' + varcompf +' >& /dev/null', shell=True)
+
         if not os.path.isfile(varcompf):
             # Does input folder has header files?
@@ -1065 +1521 @@
                     ins = 'rm '+ owdir+'/*_' + fh + '*.nc >& /dev/null'
                     sout = sub.call(ins, shell=True)
-                    files1h = gen.files_folder(iwdir,fh)
-                    if len(files1h) < 1:
-                        print errmsg
-                        print '  ' + main + ": folder '" + iwdir + "' does not " +   \
-                          "contain files '" + fh + "*' !!"
-                        quit(-1)
-                    files[fh] = files1h
-                    testfiles[fh] = files1h[0]
+                files1h = gen.files_folder(iwdir,fh)
+                if len(files1h) < 1:
+                    print errmsg
+                    print '  ' + main + ": folder '" + iwdir + "' does not " +   \
+                      "contain files '" + fh + "*' !!"
+                    quit(-1)
+                files[fh] = files1h
+                testfiles[fh] = files1h[0]
 
             if dbg:
@@ -1086 +1542 @@
             print warnmsg
             print '  ' + main + ": getting variables to compute already from file !!"
-            files, testfiles, allcompvar, Nvar = read_varcomp_file(owdir +           \
-              '/varcompute.inf')
+            files, testfiles, allcompvar, Nvar = read_varcomp_file(varcompf)
 
         # End of avoiding to repeat all the experiment search
@@ -1094 +1549 @@
 
         if dbg:
-            print 'Variables to compute _______', type(allcompvar)
+            print 'Variables to compute _______'
             gen.printing_dictionary(allcompvar)
 
@@ -1103 +1558 @@
         compute_vars(cnf, Modinf, iwdir, owdir, files, allcompvar, filescratch, dbg)
 
+##
+# Figures of variables direct from files
+##
+        print "  " + main + ": Plotting direct figures ..."
+        dirfigf = owdir + '/directplotsdraw.inf'
+        if figscratch:
+            sout = sub.call('rm ' + dirfigf + ' >& /dev/null', shell=True)
+
+            objf = open(owdir+'/all_figures.inf','w')
+            objf.write("## Drawing of all figures \n")
+            objf.close()
+
+        if addfigures:
+            sout = sub.call('rm ' + dirfigf + ' >& /dev/null', shell=True)
+
+        if not os.path.isfile(dirfigf):
+            listplots, Nplt = plots_listconstruct(cnf, owdir, dbg)
+        else:
+            print warnmsg
+            print '  ' + main + ": getting plots to draw already from file !!"
+            listplots, Nplt = read_plot_file(dirfigf)
+
+        # End of avoiding to repeat all the plots search
+
+        print "  For experiment '"+exp+"' is required to plot:", Nplt, "plots"
+
+        if dbg:
+            print 'Plots to draw _______'
+            gen.printing_dictionary(listplots)
+
+        draw_plots(cnf, listplots, mod, exp, owdir, allcompvar, figscratch, dbg)
+
         quit()
-
     # end of experiments loop
+
 # end of mods loop
 

trunk/tools/model_graphics_template.dat

@@ -30 +30 @@
 
 # Have new files been added and do they need to be processed?
-addfiles = true
+addfiles = false
 
 # Have new figures been added and do they need to be drawn?
-addfigures = false
+addfigures = true
 
 # Debug
@@ -86 +86 @@
 #   in foudre values from: ${HOME}/UNSW-CCRC-WRF/tools/postprocess/GMS-UC/WRF4G/util/postprocess/wrfncxnj/wrfncxnj.table
 #   must exist on the table $pyHOME/variables_values.dat
-#   varcombo = [varn];[vark1]@...@[varkN]: variable as consecutive kinds of operation
+#
 # Variables to compute (as possible multiple consecutive combination of operations must start by 'VAR_' separated by '+'
 #   VAR_[calc1]+[calc2] = tas:wss
@@ -108 +108 @@
 #    opsur_mean = tmean:xmean
 #      once 'tmean', or 'xmean' are computed, variable [varn] will become [varn]mean
-#   NOTE: surnames are attaced only once 
 opsur_mean = Lmean:lmean:tmean:xmean:ymean
+
+# Text as title for each operation
+#   as ':' separated list of [op]|[explanation '!' as spaces]
+titleoperations = acc|temporal!accumulated:diff|difference:direct|!:last|last!temporal!value:Lmean|latitudinal!mean:Lsec|latitudinal!section:lmean|longitudinal!mean:lsec|longitudinal!section:pinterp|pressure!interpolation:tmean|temporal!mean:turb|Taylor's!turbulence:xmean|x-axis!mean:ymean|y-axis!mean:zsum|vertical!aggregated
 
 # Pressure levels
@@ -130 +133 @@
 #   shadcont2Dzsec: 2D map of vertical section of 2 variables one in shadow and the other in contourn
 #
-#   ':' separated list of statitsics variable values are given as: [var]|[kind(including combo values)]
+#   ':' separated list of statitsics variable values are given as: [var]|['+' separated list of operations]
 #     in figures with more than 1 variable, use '#' to separate them
-#drawplots = shadcont2Dsfc:shadconthovmsfc:sahdcont2Dzsec
-drawplots = shadcont2Dzsec
-drawdiffplots = shadcont2Dsfc:shadconthovmsfc:sahdcont2Dzsec
+# DIRPLT_: Figures which can be made directly without another model/experiment output
 
-# Figures which can be made directly without another model/experiment output
-directplots = shadconthovmsfc:shadcont2Dsfc:shadcont2Dzsec
+DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:wss|last#tas|last
+#DIRPLT_shadconthovmsfc = wss|xmean#tas|xmean
+#DIRPLT_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean:hus|pinterp+tmean+xmean#ta|pinterp+tmean+xmean
 
-pltshadcont2Dsfc = huss|tmean#tas|tmean:vas|tmean#uas|tmean:ps|tmean#wss|tmean:pr|tmean#rsds|tmean:prw|tmean#clt|tmean:hfls|tmean#hfss|tmean:evspsbl|tmean#bils|tmean:zmaxth|tmean#stherm|tmean:wakeh|tmean#wakes|tmean:zmaxth|tmean#stherm|tmean:prls|tmean#prc|tmean:clt|tmean#cll|tmean:clm|tmean#clh|tmean:hurs|tmean#zmla|tmean:huss|last#tas|last:vas|last#uas|last:ps|last#wss|tmean:pr|last#rsds|last:prw|last#clt|last:hfls|last#hfss|last:evspsbl|last#bils|last:zmaxth|last#stherm|last:wakeh|last#wakes|last:zmaxth|last#stherm|last:prls|last#prc|last:clt|last#cll|last:clm|last#clh|last:hurs|last#zmla|last
-pltshadconthovmsfc = hurs|xmean#tas|xmean:wss|xmean#ps|xmean
-pltshadcont2Dzsec = ua|pinterp@last@xmean#va|pinterp@last@xmean:hur|pinterp@last@xmean#ta|pinterp@last@xmean:ws|pinterp@turb@xmean#hus|pinterp@turb@xmean
+#pltshadcont2Dsfc = huss|tmean#tas|tmean:vas|tmean#uas|tmean:ps|tmean#wss|tmean:pr|tmean#rsds|tmean:prw|tmean#clt|tmean:hfls|tmean#hfss|tmean:evspsbl|tmean#bils|tmean:zmaxth|tmean#stherm|tmean:wakeh|tmean#wakes|tmean:zmaxth|tmean#stherm|tmean:prls|tmean#prc|tmean:clt|tmean#cll|tmean:clm|tmean#clh|tmean:hurs|tmean#zmla|tmean:huss|last#tas|last:vas|last#uas|last:ps|last#wss|tmean:pr|last#rsds|last:prw|last#clt|last:hfls|last#hfss|last:evspsbl|last#bils|last:zmaxth|last#stherm|last:wakeh|last#wakes|last:zmaxth|last#stherm|last:prls|last#prc|last:clt|last#cll|last:clm|last#clh|last:hurs|last#zmla|last
+#pltshadconthovmsfc = hurs|xmean#tas|xmean:wss|xmean#ps|xmean
+#pltshadcont2Dzsec = ua|pinterp@last@xmean#va|pinterp@last@xmean:hur|pinterp@last@xmean#ta|pinterp@last@xmean:ws|pinterp@turb@xmean#hus|pinterp@turb@xmean
 
 pltdiffshadcont2Dsfc = tas|tmean#wss|tmean:wss|last#tas|last
@@ -149 +151 @@
 # Specific variable-statistisc graphic parameters (if no value is given, they will be retrieved from 'variables_values.dat')
 #  ':' separated list of specific parameters for plotting as:
-#      [varname]|[statistics]|[minval]|[maxval]|[colorbar]|[cntformat]|[colorcnt]
+#      [varname]|[operation]|[minval]|[maxval]|[colorbar]|[cntformat]|[colorcnt]
 #        [varname]: name of the variable
-#        [statistics]: statistics (also combo)
+#        [operation]: '+' separated list of operations
 #        [figure]: figure to be used
 #        [minval]: minimum value