Changeset 1041 in lmdz_wrf for trunk

Timestamp:

Aug 22, 2016, 12:00:16 PM (9 years ago)

Author:

lfita

Message:

Addding `tstd': temporal standard deviation

Location:

trunk/tools

Files:

• model_graphics.py (modified) (9 diffs)
• model_graphics_template.dat (modified) (9 diffs)

trunk/tools/model_graphics.py

@@ -29 +29 @@
     mandatorykeys = ['pyHOME', 'cdoHOME', 'scratch', 'filescratch', 'figscratch',    \
       'diffscratch', 'figdiffscratch', 'addfiles', 'addfigures', 'adddiffs',         \
-      'debug', 'models', 'ifold', 'ofold', 'warnmsg', 'errmsg', 'titleoperations',   \
-      'kindfig', 'CFreftime', 'CFunitstime', 'mapval', 'timekind', 'timefmt',        \
-      'timelabel' ]
+      'adddifffigures', 'debug', 'models', 'ifold', 'ofold', 'warnmsg', 'errmsg',    \
+      'titleoperations', 'kindfig', 'CFreftime', 'CFunitstime', 'mapval',            \
+      'timekind', 'timefmt', 'timelabel' ]
 
     # List with the optional keys
@@ -148 +148 @@
         adddiffs = False
 
+    if config['adddifffigures'] == 'true':
+        adddifffigs = True
+    else:
+        adddifffigs = False
+
     if config['debug'] == 'true':
         debug = True
@@ -153 +158 @@
         debug = False
 
-    return scr, filescr, figscr, difscr, figdifscr, addfils, addfigs, adddiffs, debug
+    return scr, filescr, figscr, difscr, figdifscr, addfils, addfigs, adddiffs, adddifffigs, debug
 
 def exp_headers(mod,config):
@@ -675 +680 @@
         pinterp: pressure interpolation (to the given $plevels)
         tmean: temporal mean values
+        tstd: temporal standard deviation values
         turb: Taylor's turbulence decomposition value
         xmean: x-axis mean values
@@ -915 +921 @@
             varkeep.append('timestats')
 
+        elif op == 'tstd':
+            # temporal standard deviation values
+            vals='time|-1,time,std,' + ':'.join(varnCFs) + ':' + vdnz
+            dims = gen.dictvar_listS(varnCFs,CFdimvardict,',','@') + dnz+'@'+vdnz
+            if dofile:
+                try:
+                    with gen.Capturing() as output:
+                        ncvar.file_oper_alongdims(vals,prevfile,CFvarnp)
+                except:
+                    print errmsg
+                    print 'file_oper_alongdims('+vals+',', prevfile, ','+CFvarnp+')'
+                    for sout in output: print sout
+                    quit(-1)
+
+                if db:
+                    for sout in output: print sout
+
+                sout = sub.call('mv file_oper_alongdims_std.nc '+fileon, shell=True)
+
+                # Keeping the operations
+                pyins=pyH + "/nc_var.py -o file_oper_alongdims -S '" + vals +        \
+                  "' -f " + prevfile + " -v " + CFvarnp
+                otrackf.write("\n")
+                otrackf.write("# " + CFvarn + " " + Fopers + "\n")
+                otrackf.write(pyins + "\n")
+
+            # removing dimension variable-dimension 'time'
+            varnCFs.remove('time')
+
+            varkeep.append('timestats')
+
 #        elif op == 'turb':
 #            # turbulence values
@@ -2007 +2044 @@
         pinterp: pressure interpolation (to the given $plevels)
         tmean: temporal mean values
+        tstd: temporal standard deviation values
         turb: Taylor's turbulence decomposition value
         xmean: x-axis mean values
@@ -2167 +2205 @@
             varkeep.append('timestats')
 
+        elif op == 'tstd':
+            # temporal standard deviation values
+            vals='time|-1,time,std,' + ':'.join(varnCFs) + ':' + vdnz
+            dims = gen.dictvar_listS(varnCFs,CFdimvardict,',','@') + dnz+'@'+vdnz
+            if dofile:
+                try:
+                    with gen.Capturing() as output:
+                        ncvar.file_oper_alongdims(vals,prevfile,CFvarnp)
+                except:
+                    print errmsg
+                    print 'file_oper_alongdims('+vals+',', prevfile, ','+CFvarnp+')'
+                    for sout in output: print sout
+                    quit(-1)
+
+                if db:
+                    for sout in output: print sout
+
+                sout = sub.call('mv file_oper_alongdims_std.nc '+fileon, shell=True)
+
+                # Keeping the operations
+                pyins=pyH + "/nc_var.py -o file_oper_alongdims -S '" + vals +        \
+                  "' -f " + prevfile + " -v " + CFvarnp
+                otrackf.write("\n")
+                otrackf.write("# " + varn + " " + Fopers + "\n")
+                otrackf.write(pyins + "\n")
+
+            # removing dimension variable-dimension 'time'
+            varnCFs.remove('time')
+
+            varkeep.append('timestats')
+
 #        elif op == 'turb':
 #            # turbulence values
@@ -2680 +2749 @@
 # scratches
 scratch, filescratch, figscratch, diffscratch, figdiffscratch, addfiles, addfigures, \
-  adddiffs, dbg = scratches(cnf)
+  adddiffs, adddifffigures, dbg = scratches(cnf)
 
 # Getting models
@@ -2921 +2990 @@
             objf.close()
 
-        if addfigures:
+        if adddifffigures:
             sout = sub.call('rm ' + dirfigf + ' >& /dev/null', shell=True)
 

trunk/tools/model_graphics_template.dat

@@ -36 +36 @@
 
 # Have new files been added and do they need to be processed?
-addfiles = false
+addfiles = true
 
 # Have new figures been added and do they need to be drawn?
-addfigures = false
+addfigures = true
 
 # Have new differences been added and do they need to be processed and drawn?
-adddiffs = false
+adddiffs = true
+
+# Have new difference figures been added and do they need to be drawn?
+adddifffigures = true
 
 # Debug
@@ -88 +91 @@
 #   pinterp: pressure interpolation (to the given $plevels, see below)
 #   tmean: temporal mean values
+#   tstd: temporal standard deviation values
 #   turb: Taylor's turbulence decomposition value
 #   xmean: x-axis mean values
@@ -101 +105 @@
 VAR_last = uas:vas:hfss:hfls
 VAR_tmean = uas:vas:hfss:hfls
+VAR_tstd = tas:wss
 VAR_xmean = uas:vas:hfss:hfls
 #VAR_last+pinterp+xmean = ua:va:hur:hus:turb
@@ -122 +127 @@
 DIFFOP_xmean = uas:vas
 DIFFOP_tmean = hfss:hfls
+DIFFOP_tstd = tas:wss
 DIFFOP_pinterp+tmean+xmean = ua:va
 
@@ -127 +133 @@
 DIFFVAR_xmean = uas:vas
 DIFFVAR_tmean = hfss:hfls
+DIFFVAR_tstd = tas:wss
 DIFFVAR_pinterp+tmean+xmean = ua:va
 
@@ -134 +141 @@
 #      once 'tmean', or 'xmean' are computed, variable [varn] will become [varn]mean
 opsur_mean = Lmean:lmean:tmean:xmean:ymean
+opsur_std = tstd
 
 # Text as title for each operation
 #   as ':' separated list of [op]|[explanation '!' as spaces]
-titleoperations = acc|tacc:diff|diff:direct|!:last|last:Lmean|latmean:Lsec|latsec:lmean|lonmean:lsec|lonsec:pinterp|pinterp:tmean|tmean:turb|turbulence:xmean|xmean:ymean|ymean:zsum|zsum
+titleoperations = acc|tacc:diff|diff:direct|!:last|last:Lmean|latmean:Lsec|latsec:lmean|lonmean:lsec|lonsec:pinterp|pinterp:tmean|tmean:tstd|tstd:turb|turbulence:xmean|xmean:ymean|ymean:zsum|zsum
 
 # Pressure levels
@@ -162 +170 @@
 # DIRPLT_: Figures which can be made directly without another model/experiment output
 
-DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last
+DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:tas|tstd#wss|tstd
 DIRPLT_shadconthovmsfc = hfls|xmean#hfss|xmean
 DIRPLT_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean
@@ -170 +178 @@
 #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
 
-PLOTDIFFOP_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean
+PLOTDIFFOP_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean:tas|tstd#wss|tstd
 PLOTDIFFOP_shadconthovmsfc = uas|xmean#vas|xmean
 PLOTDIFFOP_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean
 
-PLOTDIFFVAR_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean
+PLOTDIFFVAR_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean:tas|tstd#wss|tstd
 PLOTDIFFVAR_shadconthovmsfc = uas|xmean#vas|xmean
 PLOTDIFFVAR_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean
@@ -189 +197 @@
 #        [cntformat]: format of the contour labels
 #        [colorcnt]: color for the countor lines
-specificvarplot = cll|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:clh|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:husturbmean|pinterp@turb@xmean|shadcont2Dzsec|0.|1.e-6|Blues|%g|black
-specificdiffopplot = hfls|tmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:hfss|tmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:hfls|last|shadcont2Dsfc|-800.|0.|seismic|%g|black:hfss|last|shadcont2Dsfc|-400.|10.|seismic|%g|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black
-specificdiffvarplot = hfls|tmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:hfss|tmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:hfls|last|shadcont2Dsfc|-800.|0.|seismic|%g|black:hfss|last|shadcont2Dsfc|-400.|10.|seismic|%g|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black
+specificvarplot = cll|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:clh|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:husturbmean|pinterp@turb@xmean|shadcont2Dzsec|0.|1.e-6|Blues|%g|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|10.|Blues|%g|black
+specificdiffopplot = hfls|tmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:hfss|tmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:hfls|last|shadcont2Dsfc|-800.|0.|seismic|%g|black:hfss|last|shadcont2Dsfc|-400.|10.|seismic|%g|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|10.|Blues|%g|black
+specificdiffvarplot = hfls|tmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:hfss|tmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:hfls|last|shadcont2Dsfc|-800.|0.|seismic|%g|black:hfss|last|shadcont2Dsfc|-400.|10.|seismic|%g|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|10.|Blues|%g|black
 
 # Figures output kind