Changeset 2305 in lmdz_wrf

Timestamp:

Jan 30, 2019, 9:41:55 PM (6 years ago)

Author:

lfita

Message:

Trying to fix. HGT values are transpose!!! (Due to Fortran retrieving indices?)

File:

• trunk/tools/nc_var_tools.py (modified) (1 diff)

trunk/tools/nc_var_tools.py

@@ -28524 +28524 @@
 
     sliceNAt = sliceNA.transpose()
+    sliceNBt = sliceNB.transpose()
+
     sliceinAt = sliceinA.transpose()
-    sliceNBt = sliceNB.transpose()
     sliceinBt = sliceinB.transpose()
+    ## We are in Fortran thus (dx,dy) and 1st:1 !!!!
+    #sliceinAt0 = sliceinA.transpose()
+    #sliceinBt0 = sliceinB.transpose()
+    #sliceinAt = np.ones((sliceinAt0.shape), dtype=np.float)*gen.fillValueF
+    #sliceinBt = np.ones((sliceinBt0.shape), dtype=np.float)*gen.fillValueF
+    #sliceinAt[...,0] = sliceinAt0[...,1]+1
+    #sliceinAt[...,1] = sliceinAt0[...,0]+1
+    #sliceinBt[...,0] = sliceinBt0[...,1]+1
+    #sliceinBt[...,1] = sliceinBt0[...,0]+1
 
     NpointsABt, pointsABt, inpA, inpB = fsci.module_scientific.coincident_gridsin2d( \