Changeset 1193 in lmdz_wrf

Timestamp:

Oct 17, 2016, 6:08:12 PM (9 years ago)

Author:

lfita

Message:

Working version with differences (not understood yet, why it was not working...)

File:

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

trunk/tools/model_graphics.py

@@ -1058 +1058 @@
     Fopers = ''
     if db: print "    computing statistics of variable:'", CFvarn, "' operation '" + \
-      Opers + "'...'"
+      Opers + "' using file '" + ifilen + "...'"
     for op in opers:
 
@@ -3808 +3808 @@
         op = vnop.split('_')[1]
         values = allcompvar[vnop]
-        VnOpS = varnoper(vn, op, opsurnames)
+        VnOpS = varnoper(vn, op, opersurnames)
 
         # Only are reprojected that files which do contain 'lon' and 'lat', 
@@ -3839 +3839 @@
             sout = sub.call('rm ' + odir+'/'+ofileop + ' >& /dev/null', shell=True)
 
+        # List of files to re-project (first without operation)
         ifiles = [odir+'/'+ifileorig, odir+'/'+ifileop]
+        # List of output files after re-projection (first without operation)
         ofiles = [odir+'/'+ofileorig, odir+'/'+ofileop]
+
+        print fname + '; Lluis vn:', vn,'op:',op,'VnOpS:',VnOpS
+        print fname + '; Lluis ifiles:', ifiles
+        print fname + '; Lluis ofiles:', ofiles
+
 
         for iif in range(2):
@@ -3853 +3860 @@
                 print '  ' + fname + ": operation '" + op + "' can not be " +        \
                   "reprojected due to the lack of 'lon' or 'lat' skipping it !!"
-                print "    NOTE: difference can only be computed as 'variable' " +   \
-                  "difference"
-                print "      As operating along the differences from the original" + \
-                  " values from each model"
+                print "    NOTE: difference can only be computed as:"
+                print "      - 1st: compute variable's differences between models"
+                print "      - 2nd: Perform the operations"
 
                 # Building up resources to compute re-projected statistics
@@ -3892 +3898 @@
                         quit(-1)
                 else:
-                    values = 'lon,lat,' + config['RefProjfile'] + ':'.join(CFvardims)
+                    # Name of the variable to reproject
+                    if iif == 0:
+                        reprjvn = vn
+                    else: 
+                        reprjvn = VnOpS
+              
+                    values = 'lon,lat,' + config['RefProjfile'] + ',lon,lat,' +      \
+                      reprjop + ',' + ':'.join(CFvardims)
+                    ins= 'ncvar.reproject('+values + ', ' + ifile + ', '+ reprjvn + ')'
                     try:
                         with gen.Capturing() as output:
-                            sout=ncvar.reproject(values, ifile, VnOpS)
-                            ins = 'ncvar.reproject(' + values + ', ' + ifile + ', '+ \
-                              VnOpS + ')'
+                            sout=ncvar.reproject(values, ifile, reprjvn)
                     except:
                         print errmsg
@@ -3903 +3915 @@
                         for s1out in output: print s1out
                         quit(-1)
+
+                    sub.call('mv reproject.nc ' + ofile, shell=True)
 
                 if debug: