Changeset 1472 in lmdz_wrf

Timestamp:

Mar 28, 2017, 11:19:04 PM (8 years ago)

Author:

lfita

Message:

Improving `netcdf_fold_concatenation_HMT':

• Including list of files
• Including PyNCplot global attributes

File:

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

trunk/tools/nc_var_tools.py

@@ -7229 +7229 @@
     fname='compute_opervartimes'
 
-    ofile = 'opervartimes_' + varname + '.nc'
-
     if values == 'h':
         print fname + '_____________________________________________________________'
-        print opervartimes.__doc__
+        print compute_opervartimes.__doc__
         quit()
+
+    ofile = 'opervartimes_' + varname + '.nc'
 
     if not os.path.isfile(ncfile):
@@ -10604 +10604 @@
     """
     import subprocess as sub
-    fname='netcdf_fold_concatenate_HMT'
+    fname='netcdf_fold_concatenation_HMT'
 
     ofile = 'netcdf_fold_concatenated_HMT.nc'
@@ -10623 +10623 @@
     if Nfiles == 0:
         print errormsg
-        print '  ' + fname + ": there are no files as '" + fold + '/' + ncfile + "*' !!"
+        print '  ' +fname+ ": there are no files as '" + fold + '/' + '*'.join(HMT) +\
+          "*' !!"
         quit(-1)
     print confiles
@@ -10631 +10632 @@
     for filen in confiles:
         print 'charging: ',filen
-        ncobjs.append(NetCDFFile(fold + '/' + filen, 'r'))
+        if filen.find(fold) != -1:
+            ncobjs.append(NetCDFFile(filen, 'r'))
+        else:
+            ncobjs.append(NetCDFFile(fold + '/' + filen, 'r'))
 
 # Looking for the new length of the concatenation dimension
@@ -10656 +10660 @@
     else:
         newnc.createDimension(condim, totcondim)
+    newnc.createDimension('file', Nfiles)
+    newnc.createDimension('Lstring', 256)
+
+    # dim-variables
+    newvar = newnc.createVariable('file','c',('file','Lstring'))
+    newvals = writing_str_nc(newvar, confiles, 256)
+    if confiles[0].find(fold) != -1:
+        str="concatenated files following dimension '" + condim + "'"
+    else:
+        str="concatenated files from folder '" + fold + "' following dimension '" +  \
+          condim + "'"
+    vardef = basicvardef(newvar, 'files', str, '-')
+
 # Fake variable for the dimension
     print 'condim:',condim
@@ -10735 +10752 @@
 
     newnc.sync()
+    add_global_PyNCplot(newnc, main, fname, '1.1')
+
     newnc.close()