Changeset 1043 in lmdz_wrf

Timestamp:

Aug 22, 2016, 2:14:28 PM (9 years ago)

Author:

lfita

Message:

Adding `tvar'

Location:

trunk/tools

Files:

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

trunk/tools/model_graphics.py

@@ -682 +682 @@
         tstd: temporal standard deviation values
         turb: Taylor's turbulence decomposition value
+        tvar: temporal variance values (Taylor's turbulence)
         xmean: x-axis mean values
         ymean: y-axis mean values
@@ -956 +957 @@
 #            vals='turbulence|'${CFvarn}
 #            dims='time@time,'${dnz}'@'${vdnz}',lat@lat,lon@lon,pres@pres'
-#            pyout=`python ${pyHOME}/diagnostics.py -d ${dims} -v ${vals} -f ${cfiles}`
+#            pyins = pyH + "/diagnostics.py -d " + dims + " -v '" +vals + \
+#                "' -f " prevfiles
+#            if dofile:
+#                try:
+#                    with gen.Capturing() as output:
+#                        sout = sub.call('python ' + pyins, shell=True)
+#                except:
+#                    print errmsg
+#                    print pyins
+#                    for sout in output: print sout
+#                    quit(-1)
+
+#                if db:
+#                    for sout in output: print sout
+
+#                sout = sub.call('mv diagnostics.nc '+fileon, shell=True)
 
 #            # Keeping the operations
@@ -966 +982 @@
 
 #            varkeep=':'${CFvarn}'turb'
+
+        elif op == 'tvar':
+            # temporal variance values (Taylor's turbulence)
+            vals='time|-1,time,var,' + ':'.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_var.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 == 'xmean':
@@ -2046 +2093 @@
         tstd: temporal standard deviation values
         turb: Taylor's turbulence decomposition value
+        tvar: temporal variance values (Taylor's turbulence)
         xmean: x-axis mean values
         ymean: y-axis mean values
@@ -2251 +2299 @@
 #            varkeep=':'${CFvarn}'turb'
 
+        elif op == 'tvar':
+            # temporal variance values (Taylor's turbulence)
+            vals='time|-1,time,var,' + ':'.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_var.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 == 'xmean':
             # x-axis mean values

trunk/tools/model_graphics_template.dat

@@ -93 +93 @@
 #   tstd: temporal standard deviation values
 #   turb: Taylor's turbulence decomposition value
+#   tvar: temporal variance values (Taylor's turbulence)
 #   xmean: x-axis mean values
 #   ymean: y-axis mean values
@@ -106 +107 @@
 VAR_tmean = uas:vas:hfss:hfls
 VAR_tstd = tas:wss
+VAR_tvar = tas:wss
 VAR_xmean = uas:vas:hfss:hfls
 #VAR_last+pinterp+xmean = ua:va:hur:hus:turb
@@ -128 +130 @@
 DIFFOP_tmean = hfss:hfls
 DIFFOP_tstd = tas:wss
+DIFFOP_tvar = tas:wss
 DIFFOP_pinterp+tmean+xmean = ua:va
 
@@ -134 +137 @@
 DIFFVAR_tmean = hfss:hfls
 DIFFVAR_tstd = tas:wss
+DIFFVAR_tvar = tas:wss
 DIFFVAR_pinterp+tmean+xmean = ua:va
 
@@ -142 +146 @@
 opsur_mean = Lmean:lmean:tmean:xmean:ymean
 opsur_std = tstd
+opsur_var = tvar
 
 # 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:tstd|tstd: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:tvar|tvar:xmean|xmean:ymean|ymean:zsum|zsum
 
 # Pressure levels
@@ -170 +175 @@
 # DIRPLT_: Figures which can be made directly without another model/experiment output
 
-DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:tas|tstd#wss|tstd
+DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:tas|tstd#wss|tstd:wss|tvar#tas|tvar
 DIRPLT_shadconthovmsfc = hfls|xmean#hfss|xmean
 DIRPLT_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean
@@ -178 +183 @@
 #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:tas|tstd#wss|tstd
+PLOTDIFFOP_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean:tas|tstd#wss|tstd:wss|tvar#tas|tvar
 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:tas|tstd#wss|tstd
+PLOTDIFFVAR_shadcont2Dsfc = hfls|last#hfss|last:hfss|tmean#hfls|tmean:tas|tstd#wss|tstd:wss|tvar#tas|tvar
 PLOTDIFFVAR_shadconthovmsfc = uas|xmean#vas|xmean
 PLOTDIFFVAR_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean
@@ -197 +202 @@
 #        [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: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
+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:tas|tstd|shadcont2Dsfc|0.|100.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|100.|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:tas|tstd|shadcont2Dsfc|0.|100.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|100.|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:tas|tstd|shadcont2Dsfc|0.|100.|Reds|%g|black:wss|tstd|shadcont2Dsfc|0.|100.|Blues|%g|black
 
 # Figures output kind