Changeset 2417 in lmdz_wrf

Timestamp:

Mar 25, 2019, 4:05:16 PM (6 years ago)

Author:

lfita

Message:

Adding fix for Nmaxgrisin == 1 in `usefile_compute_slices_stats_areaweighted'

File:

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

trunk/tools/nc_var_tools.py

@@ -30682 +30682 @@
 
     # removing monotones. Fortran only accepts rank <= 7 arrays !!
+    Nmaxgridsinone = False
+    if sin.shape[1] == 1: Nmaxgridsinone = True
+
     sN, Ndims = gen.remove_monotones(sN, osN.dimensions)
     sin, indims = gen.remove_monotones(sin, osin.dimensions)
@@ -30689 +30692 @@
     dimsslice = list(Ndims)
     shapeslice = list(sN.shape)
+
+    # Due to the removal of monotones we are losing the cases when Nmaxgrisin == 1
+    if Nmaxgridsinone:
+        print warnmsg
+        print '  ' + fname + ": maximum number of grid points per slice == 1 !!"
+        print '    adjusting accordingly certain matrices ...' 
+        shapein = list(sin.shape)
+        Lshape = len(shapein)
+        reshapein = [shapein[0]] + [1] + shapein[1:Lshape+1]
+
+        sin0 = sin + 0
+
+        newsin = np.zeros(tuple(reshapein), dtype=int)
+        for i in range(shapein[0]):
+            newsin[i,0,...] = sin[i,...]
+
+        sin = None
+        sin = newsin+ 0
 
     sNt = sN.transpose()