Changeset 1777 in lmdz_wrf for trunk/tools/diag_tools.py

Timestamp:

Feb 16, 2018, 3:17:02 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `WRFz': as height above surface unstaggering WRFgeop
• Fixing zwind computation

File:

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

trunk/tools/diag_tools.py

@@ -911 +911 @@
 
 
-def Forcompute_zwind(ua, va, zsl, uas, vas, hgt, sina, cosa, zval, dimns, dimvns):
+def Forcompute_zwind(ua, va, z, uas, vas, sina, cosa, zval, dimns, dimvns):
     """ Function to compute the wind at a given height following the power law method
     Forcompute_zwind(ua, va, zsl, uas, vas, hgt, sina, cosa, zval, dimns, dimvns)
       [ua]= x-component of unstaggered 3D wind (assuming [[t],z,y,x]) [ms-1]
       [va]= y-component of unstaggered 3D wind (assuming [[t],z,y,x]) [ms-1]
-      [zsl]= height above sea level [m]
+      [z]= height above surface [m]
       [uas]= x-component of unstaggered 10 m wind (assuming [[t],z,y,x]) [ms-1]
       [vas]= y-component of unstaggered 10 m wind (assuming [[t],z,y,x]) [ms-1]
-      [hgt]= topographical height (assuming [m]
       [sina]= local sine of map rotation [1.]
       [cosa]= local cosine of map rotation [1.]
@@ -944 +943 @@
 
         pvar1, pvar2= fdin.module_fordiagnostics.compute_zwind4d(ua=ua.transpose(),  \
-          va=va[:].transpose(), zsl=zsl[:].transpose(), uas=uas.transpose(),         \
-          vas=vas.transpose(), hgt=hgt.transpose(), sina=sina.transpose(),           \
-          cosa=cosa.transpose(), zextrap=zval, d1=dx, d2=dy, d3=dz, d4=dt)
+          va=va[:].transpose(), z=z[:].transpose(), uas=uas.transpose(),             \
+          vas=vas.transpose(), sina=sina.transpose(), cosa=cosa.transpose(),         \
+          zextrap=zval, d1=dx, d2=dy, d3=dz, d4=dt)
         var1 = pvar1.transpose()
         var2 = pvar2.transpose()