Changeset 1026 in lmdz_wrf for trunk

Timestamp:

Aug 17, 2016, 9:22:22 AM (9 years ago)

Author:

lfita

Message:

Working version with figures shadcont2Dsfc' and shadconthovmsfc'

Location:

trunk/tools

Files:

• model_graphics.py (modified) (10 diffs)
• model_graphics_template.dat (modified) (4 diffs)

trunk/tools/model_graphics.py

@@ -1140 +1140 @@
 
 def draw_plot(kplot, CFvplot, fplot, vplot, dplot, pplot, finame, tfig, kfig, mapval,\
-  tunits, od, pyH, fscr, db):
+  tvals, od, pyH, fscr, db):
     """ Function to draw a plot
       kplot= kind of plot
@@ -1152 +1152 @@
       kfig= kind of the figure
       mapval= value of the map
-      tunits= units of time
+      tvals= list with time values
+        tunits: units of time ('!' for spaces)
+        timekind: kind of time ticks in the plot
+        timefmt: format of time ticks in the plot
+        timelabel: label of time-axis in the plot
       od= output directory
       pyH= python HOME
@@ -1169 +1173 @@
     # Tracking file with all the figures
     trkobjf = open('all_figures.inf', 'a')
+
+    # Time-characteristics
+    tunits = tvals[0]
+    timekind = tvals[1]
+    timefmt = tvals[2]
+    timelabel = tvals[3]
 
     if not os.path.isfile(finame):
@@ -1201 +1211 @@
             cntstdn = CFvplot[1]
             figfs = ','.join(fplot)
+
             shad = pplot[0]
             srange = str(shad[0]) + ',' + str(shad[1])
@@ -1206 +1217 @@
             cnt = pplot[1]
             crange = str(cnt[0]) + ',' + str(cnt[1])
-            cfmt = cnt[2]
-            cline = cnt[3]
+            cfmt = cnt[3]
+            cline = cnt[4]
 
             graphvals = ','.join(CFvplot)
@@ -1220 +1231 @@
             try:
                 with gen.Capturing() as output:
-                    sout = sub.call(plotins, shell=True)
+                    sout0 = sub.call(plotins, shell=True)
             except:
                 print 'drawing.draw_2D_shad_cont(' + graphvals + ',' + figfs + ',' + \
                   fvarS + ')'
+                print sout0
                 for sout in output: print sout
                 quit(-1)
@@ -1233 +1245 @@
             trkobjf.write("#" + tfig.replace('!',' ') + '\n')
             trkobjf.write(plotins + '\n')
+
+        elif kplot == 'shadconthovmsfc':
+            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[3]
+            cline = cnt[4]
+            # It is assumed that if the space variable is 'lon': is desired a 
+            #   (lon, time) plot it it is 'lat': then (time, lat) plot
+            if gen.searchInlist(dplot,'lon'):
+                spacedim ='lon'
+                figmid = 'longitudinal|evolution|of'
+                figtit = '|'.join(tfig.split('!')[0:2]) + '|' + figmid + '|' +       \
+                  '|'.join(tfig.split('!')[2:])
+                reverse= 'None'
+                dims= spacedim+'|-1,time|-1;'+spacedim+'|-1,time|-1;'+spacedim+';time'
+            else:
+                spacedim='lat'
+                figmid = 'meridional|evolution|of'
+                figtit = '|'.join(tfig.split('!')[0:2]) + '|' + figmid + '|' +       \
+                  '|'.join(tfig.split('!')[2:])
+                reverse='None'
+                dims= spacedim+'|-1,time|-1;'+spacedim+'|-1,time|-1;time;'+spacedim
+
+            graphvals = ','.join(CFvplot)
+            graphvals = graphvals + ';'+ dims +';'+ cbar +  \
+              ';fixsigc,' + cline + ';' + cfmt + ';' + srange + ';' + crange +',9;'+ \
+              figtit + ';' + kfig + ';' + reverse + ';' + 'time|' + tunits + '|'+ \
+              timekind + '|' + timefmt + '|' + timelabel
+
+            fvarS = ','.join(vplot)
+
+            drwins = 'draw_2D_shad_cont_time'
+            plotins = 'python ' + pyH + '/drawing.py -f ' + figfs +' -o ' + drwins + \
+              " -S '" + graphvals + "' -v " + fvarS
+
+            try:
+                with gen.Capturing() as output:
+                    sout0 = sub.call(plotins, shell=True)
+            except:
+                print 'draw_2D_shad_cont_time(' + graphvals + ',' + figfs + ',' +    \
+                  fvarS + ')'
+                print sout0
+                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
@@ -1280 +1352 @@
         opersurnames[opn] = vls
 
+    # time values
     # Units time for the plots
     rd = config['CFreftime']
-    timeunits = config['CFunitstime'] + '!since!' + rd[0:4] + '-' + rd[4:6] + '-' +  \
+    tunits = 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
+    # time ticks kind
+    tkind = config['timekind']
+    # time ticks format
+    tfmt = config['timefmt']
+    # time axis label
+    tlab = config['timelabel']
+    timevals = [tunits, tkind, tfmt, tlab]
+
+    # Dictionary of plot specificities
     #   [minval]: minimum value
     #   [maxval]: minimum value
@@ -1405 +1485 @@
                 # Title of the figure:
                 if ivp == 0:
-                    titfigure = mod + '!' + exp + '!' + vn 
+                    titfigure = mod + '!' + exp + "!'" + vn + "'"
                 else:
-                    titfigure = titfigure + '!&!' + vn
-                for op1 in op.split('+'):
+                    titfigure = titfigure + "!&!'" + vn + "'"
+                opvals = op.split('+')
+                for op1 in opvals:
                     if not opexplained.has_key(op1):
                         print errormsg
@@ -1414 +1495 @@
                         print '    provided:', opexplained.keys()
                         quit(-1)
-                    titfigure = titfigure + '!' + opexplained[op1]
+                    if op1 == opvals[0]:
+                        titfigure = titfigure + '$_{[' + opexplained[op1]
+                        if len(opvals) == 1: titfigure = titfigure + ']}$'
+                    elif op1 == opvals[len(opvals)-1]:
+                        titfigure = titfigure + '\!' + opexplained[op1] + ']}$'
+                    else:
+                        titfigure = titfigure + '\!' + opexplained[op1]
 
                 ivp = ivp + 1
             # End of variable-operation
+            figname = figname + '.' + kindfigure
 
             draw_plot(kplot, CFvarsplot, filesplot, varsplot, dimsplot, pictplot,    \
-              figname, titfigure, kindfigure, mapvalue, timeunits, odir, pyHOME,     \
+              figname, titfigure, kindfigure, mapvalue, timevals, odir, pyHOME,      \
               figscr, debug)
 

trunk/tools/model_graphics_template.dat

@@ -33 +33 @@
 
 # Have new figures been added and do they need to be drawn?
-addfigures = true
+addfigures = false
 
 # Debug
-debug = true
+debug = false
 
 # Folder with the files
@@ -112 +112 @@
 # 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
+titleoperations = acc|tacc:diff|diff:direct|!:last|last:Lmean|latmean:Lsec|latsec:lmean|lonmean:lsec|lonsec:pinterp|pressure!interpolation:tmean|temporal!mean:turb|turbulence:xmean|xmean:ymean|ymean:zsum|zsum
 
 # Pressure levels
@@ -138 +138 @@
 
 DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:wss|last#tas|last
-#DIRPLT_shadconthovmsfc = wss|xmean#tas|xmean
+DIRPLT_shadconthovmsfc = wss|xmean#tas|xmean
 #DIRPLT_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean:hus|pinterp+tmean+xmean#ta|pinterp+tmean+xmean
 
@@ -177 +177 @@
 #  [timefmt]; desired format for the time ticks (combination of LaTeX and 'linux date formats')
 #  [timelabel]; label at the graph ('!' for spaces)
-
 timekind = exct,6,h
 timefmt = $%d^{%H}$