Changeset 2561 in lmdz_wrf

Timestamp:

May 28, 2019, 2:56:00 PM (6 years ago)

Author:

lfita

Message:

removing spurious print in module_scientific.f90' and in nc_var_tools.py'

Location:

trunk/tools

Files:

• module_scientific.f90 (modified) (1 diff)
• nc_var_tools.py (modified) (11 diffs)

trunk/tools/module_scientific.f90

@@ -2629 +2629 @@
   ! Filling with the points of all the space with .TRUE.
   Npts = COUNT(boolmat)
-  PRINT *, dx, dy, 'Lluis Npts:', Npts, 'Nppt:', Nppt
 
   IF (ALLOCATED(points)) DEALLOCATE(points)

trunk/tools/nc_var_tools.py

@@ -29391 +29391 @@
                 getvarybndst = ydimvarbnds2D.transpose()
 
-                print 'ref:', dxref, dyref
-                print 'Lluis looking for:', reflon, 'lat.', reflat
-                for j in range(dyref):
-                    print 'Lluis ref lat.', yslice2D[:,j,0]
-                print 'Lluis get lon.', getlon[0,:]
-                print 'Lluis get lat.', getlat[:,0]
-                for i in range(dxget):
-                    print i,'Lluis get lat.', xdimvarbnds2D[:,0,i]
-                for j in range(dyget):
-                    print j,'Lluis get lat.', ydimvarbnds2D[:,j,0]
-
                 Ngridsint, gridsint, areas2Dt, areast, percenst =                    \
                   fsci.module_scientific.grid_spacepercen(xcavals=reflont,           \
@@ -29644 +29633 @@
                 newdim = onewnc.createDimension('coord',2)
 
+            print '    Lluis dims ', dn+'Ngrid', 'Srgrid', Srgrid
+
             newvar = onewnc.createVariable(dn+'Ngrid','i', tuple(Srgrid))
             if Ngridsin.shape[1] == 1:
@@ -29773 +29764 @@
     dxyB = osliceinB.shape[1]
 
+    print 'Lluis dims osliceNB', osliceNB.dimensions, 'osliceNA', osliceNA.dimensions
+
+
     Srgrid = list(osliceNB.dimensions) + list(osliceNA.dimensions)
     Sigrid = ['coord', dn+'gridin'] + Srgrid
@@ -29804 +29798 @@
     inpointsB = inpB.transpose()
 
+    print 'Lluis here 00'
+
     # Remembering that it is python (C-like...)
     inpointsA = inpointsA-1
@@ -29819 +29815 @@
     if not gen.searchInlist(onewnc.dimensions, dn+'gridin'):
         newdim = onewnc.createDimension(dn+'gridin', maxNpointsAB)
+
+    print 'Lluis file creation'
+
+    print 'Lluis shapes: newvar', newvar.shape, 'NpointsAB', NpointsAB.shape, 'Srgrid', Srgrid
 
     newvar = onewnc.createVariable(dn+'Ngrid','i', tuple(Srgrid))
@@ -29826 +29826 @@
     newvar.setncattr('coordinates',' '.join(Srgrid[::-1]))
     
+    print 'Lluis file creation 1'
+
     innewvar = onewnc.createVariable(dn+'gridin', 'i', tuple(Sigrid),                \
       fill_value=gen.fillValueI)
@@ -29832 +29834 @@
     innewvar.setncattr('coordinates',' '.join(Srgrid[::-1]))
 
+    print 'Lluis file creation 2'
+
     aanewvar = onewnc.createVariable(dn+'area','f',tuple(Srgrid),                    \
       fill_value=gen.fillValueF)
@@ -29838 +29842 @@
     aanewvar.setncattr('coordinates',' '.join(Srgrid[::-1]))
 
+    print 'Lluis file creation 3'
+
     anewvar = onewnc.createVariable(dn+'gridarea','f',tuple(Spgrid),                 \
       fill_value=gen.fillValueF)
@@ -29844 +29850 @@
     anewvar.setncattr('coordinates',' '.join(Spgrid[::-1]))
 
+    print 'Lluis file creation 4'
+
     pnewvar = onewnc.createVariable(dn+'gridpercen','f',tuple(Spgrid),               \
       fill_value=gen.fillValueF)
@@ -29850 +29858 @@
     pnewvar.setncattr('coordinates',' '.join(Spgrid[::-1]))
 
-    jjA = 10
-    iiA = 2
-    jjB = 0
-    iiB = 4
+    print 'Lluis prior to... dyA, dxA, dyB, dxB', dyA, dxA, dyB, dxB
+
     for jA in range(dyA):
         for iA in range(dxA):
@@ -29890 +29896 @@
     #onewnc.close()
     #quit()
+
+    print 'Lluis here !!'
 
     if Nnewslcs >= 3: