Changeset 1238 in lmdz_wrf for trunk

Timestamp:

Oct 27, 2016, 10:37:16 AM (9 years ago)

Author:

lfita

Message:

Commiting with some tests due to 'segmentation fault'!

Location:

trunk/tools

Files:

• model_graphics.py (modified) (5 diffs)
• nc_var_tools.py (modified) (5 diffs)

trunk/tools/model_graphics.py

@@ -914 +914 @@
     otrackf.close()
 
+    quit(-1)
     # Joining variable files
     if not os.path.isfile(fileon):
@@ -1158 +1159 @@
             # pinterp: pressure interpolation (to the given $plevels)
             vals=plevels + ',1,1'
+            print fname + '; Lluis vals:', vals, 'plevels:', plevels
             if dofile:
                 try:
@@ -3807 +3809 @@
         op = vnop.split('_')[1]
         values = allcompvar[vnop]
-        print fname + ': Lluis first values:', values
         VnOpS = varnoper(vn, op, opersurnames)
 
@@ -4393 +4394 @@
         if addfiles:
             if dbg:
-                print '  ' + main + ": adding variables to copmute removing " +      \
+                print '  ' + main + ": adding variables to compute removing " +      \
                   "file '" + varcompf + "' ..."
             sub.call('rm ' + varcompf +' >& /dev/null', shell=True)
@@ -4410 +4411 @@
                     sout = sub.call(ins, shell=True)
                 files1h = gen.files_folder(iwdir,fh)
+                print sub.call('ls ' + iwdir + '/' + fh + '*', shell=True)
+                print main + '; Lluis file1h:', files1h, 'iwdir:', iwdir, 'fh:', fh
                 if len(files1h) < 1:
                     print errmsg

trunk/tools/nc_var_tools.py

@@ -17078 +17078 @@
         ovar2 = onc.variables['PB']
 
-        pres = (ovar1[:] + ovar2[:]).astype('float64')
+        pres0 = ovar1[:] + ovar2[:]
+        pres = pres0.astype('float64')
 
         dimx = pres.shape[3]
@@ -17171 +17172 @@
     # temperature
     if modname == 'WRF':
-        Rd = 287.04
-        Cp = 7.*Rd/2.
-        RCP = Rd/Cp
-        p0 = 100000.
-        ovar1 = onc.variables['T']
-        temp = ((ovar1[:]+300.)*(pres[:]/p0)**RCP).astype('float64')
+        Rd = np.float64(287.04)
+        Cp = np.float64(7.*Rd/2.)
+        RCP = np.float64(Rd/Cp)
+        p0 = np.float64(100000.)
+        ovar10 = onc.variables['T']
+        var10 = ovar10[:]
+        ovar1 = (var10).astype('float64')
+        temp0 = (ovar1[:]+300.)*(pres[:]/p0)**RCP
+        temp = temp0.astype('float64')
     elif modname == 'LMDZ' or modname == 'cfLMDZ':
         ovar1 = onc.variables['temp']
@@ -17226 +17230 @@
     # Creation of pressure variable dimension
     newvar = onewnc.createVariable('pres', 'f8', ('pres'))
-    newvar[:] = interplevs
+    newvar[:] = interplevs[:]
     basicvardef(newvar, 'pressure', 'Pressure', 'Pa')
     attr = set_attribute(newvar, 'positive', 'down')
@@ -17295 +17299 @@
                 newvarattr[attrn] = attrv
 
+        print fname + '; Lluis varin shape:', varin.shape
         if varin is not None:
             varint = varin.transpose()
@@ -17303 +17308 @@
             qvt = qv.transpose()
 
+            print fname + '; Lluis going to the subroutine....'
+            print fname + '; Lluis shapes: varint:', varint.shape,'prest:',prest.shape,'psfct:',psfct.shape,\
+              'hgtt:', hgtt.shape,'tempt:',tempt.shape,'qvt:',qvt.shape
             varinterp = fin.module_forinterpolate.interp( data_in=varint,            \
               pres_field=prest, interp_levels=interplevs, psfc=psfct,                \
               ter=hgtt, tk=tempt, qv=qvt, linlog=linloginterp,                       \
-              extrapolate=extrap, geopt=isgeop, missing=gen.fillValueF)
+              extrapolate=extrap, geopt=isgeop, missing=gen.fillValueF, ix=dimx,     \
+              iy=dimy, iz=dimz, it=dimt, num_metgrid_levels=len(interplevs))
+            print fname + '; Lluis after the subroutine'
 
             newvar = onewnc.createVariable(CFvn, 'f4', tuple(newMODdims),            \