Changeset 2620 in lmdz_wrf for trunk/tools

Timestamp:

Jun 21, 2019, 2:50:32 PM (6 years ago)

Author:

lfita

Message:

Fixing computation of zmla using U. Wyoming sounding data

File:

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

trunk/tools/diagnostics.py

@@ -2144 +2144 @@
         var1 = ncobj.variables[depvars[1]][:]
         dimz = var0.shape[1]
-        var2 = ncobj.variables[depvars[2]][:]
+        var2 = ncobj.variables[depvars[2]][:]/9.8
         var3 = ncobj.getncattr('Station_elevation')
-
-        # Wrong order of dimensions !
-        var0 = var0.transpose()
-        var1 = var1.transpose()
-        var2 = var2.transpose()
 
         diagout, diagoutd, diagoutvd = diag.Forcompute_zmla_gen(var0,var1,var2,var3, \
           dnames,dvnames)
+
+        # No need to remove topography height
+        diagout = diagout + var3
 
         # Removing the nonChecking variable-dimensions from the initial list