Changeset 1048 in lmdz_wrf

Timestamp:

Aug 26, 2016, 9:21:12 PM (9 years ago)

Author:

lfita

Message:

Up to first differences with LMDZ

File:

• trunk/tools/model_graphics.py (modified) (11 diffs)

trunk/tools/model_graphics.py

@@ -556 +556 @@
                 try:
                     with gen.Capturing() as output:
-                        ncvar.DataSetSection_multivars(values, cf, vs)
+                        ncvar.DataSetSection_multivars(values, idir + '/' + cf, vs)
                 except:
                     print errmsg
-                    print 'ncvar.DataSetSection_multivars('+values+', '+cf+', '+vs+')'
+                    print 'ncvar.DataSetSection_multivars('+values+', '+idir+'/'+cf+ \
+                      ', '+vs+')'
                     for sout in output: print sout
                     quit(-1)
@@ -568 +569 @@
 
                 # Keeping track of the operations
-                pyins = pyH + "/nc_var.py -f "+cf+" -o DataSetSection_multivars -v "+\
-                  vs + " -S '" + values + "'"
+                pyins = pyH + "/nc_var.py -f " + idir + '/' + cf +                   \
+                  " -o DataSetSection_multivars -v " + vs + " -S '" + values + "'"
                 otrackf.write('\n')
                 otrackf.write('# ' + CFvarn + " " + modvar + '\n')
-                otrackf.write(pyins + '\n')
+                otrackf.write('python ' + pyins + '\n')
 
                 # CF renaming of variable
@@ -584 +585 @@
                 diagc = '@'.join(diagvar[1:])
             
-                values = '-f ' + cf + " -d '" + dims + "' -v '" + diagn + '|' +      \
-                  diagc + "'"
+                values = '-f ' + idir + '/' + cf + " -d '" + dims + "' -v '" +       \
+                  diagn + '|' + diagc + "'"
                 try:
                     with gen.Capturing() as output:
@@ -624 +625 @@
                         for sout in output: print sout
 
-            # adding CF lon,lat,time in WRF files
-            if minf.name == 'WRF':
+            # CFification of files
+            if minf.name == 'WRF' or minf.name == 'WRF_LMDZ':
                 values = vdnx + ':' + vdny + ':'+ Tref + ':' + Tunits
-                ncvar.WRF_toCF(values, ifilen) 
+                try:
+                    with gen.Capturing() as output:
+                        ncvar.WRF_toCF(values, ifilen) 
+                except:
+                    print errmsg
+                    print 'WRF_toCF('+ values +', ' + ifilen + ')'
+                    for sout in output: print sout
+                    quit(-1)
+                if db:
+                    for sout in output: print sout
+            elif minf.name == 'LMDZ':
+                ncvar.LMDZ_toCF(ifilen) 
+
+                try:
+                    with gen.Capturing() as output:
+                        ncvar.LMDZ_toCF(ifilen) 
+                except:
+                    print errmsg
+                    print 'LMDZ_toCF('+ ifilen + ')'
+                    for sout in output: print sout
+                    quit(-1)
+                if db:
+                    for sout in output: print sout
+            else:
+                print errmsg
+                print '  ' + fname + ": no CFification for model '" +minf.name+ "' !!"
+                print "    available ones: 'WRF', 'WRF_LMDZ', 'LMDZ'"
+                quit(-1)
 
         ifile = ifile + 1
@@ -2940 +2968 @@
         # End of avoiding to repeat all the experiment search
 
-        print "  For experiment '"+exp+"' is required to compute:", Nvar, "variables"
+        print "    For experiment '" + exp + "' is required to compute:", Nvar,      \
+          "variables"
 
         if dbg:
@@ -2955 +2984 @@
 # Figures of variables direct from files
 ##
-        print "  " + main + ": Plotting direct figures ..."
+        print "    Plotting direct figures ..."
         dirfigf = owdir + '/directplotsdraw.inf'
         if figscratch:
@@ -2976 +3005 @@
         # End of avoiding to repeat all the plots search
 
-        print "  For experiment '"+exp+"' is required to plot:", Nplt, "plots"
+        print "    For experiment '" + exp + "' is required to plot:", Nplt, "plots"
 
         if dbg:
@@ -2992 +3021 @@
 # Experiments loop
 ##
+    print '  ******* ****** ***** **** *** ** * ** *** **** ***** ****** *******'
+    print "  ** '" + mod + "': Inter experiments differences  "
+    print '  ******* ****** ***** **** *** ** * ** *** **** ***** ****** *******'
+    print '    experiments:', exps
     # There are two kind of differences:
     #   DIFFop: differences between operations of each given variable. 
@@ -3002 +3035 @@
 
     Nexps = len(exps)
-    for exp1 in exps[0:Nexps-1]:
-        for exp2 in exps[1:Nexps]:
+    for iexp1 in range(0,Nexps-1):
+        exp1 = exps[iexp1]
+        for iexp2 in range(iexp1+1,Nexps):
+            exp2 = exps[iexp2]
             Sexps = exp2 + '-' + exp1
             print '  ' + Sexps + '...'
@@ -3043 +3078 @@
         # End of avoiding to repeat all the experiment search
 
-        print "  For experiments '"+exp2+"'-'"+exp1+"' is required to compute:",     \
+        print "    For experiments '"+exp2+"'-'"+exp1+"' is required to compute:",   \
           Nvar, "differences"
 
@@ -3094 +3129 @@
           dbg)
 
-    quit()
+#    quit()
 # end of mods loop