Changeset 1777 in lmdz_wrf for trunk/tools/module_ForDiagnosticsVars.f90

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/module_ForDiagnosticsVars.f90 (modified) (6 diffs)

trunk/tools/module_ForDiagnosticsVars.f90

@@ -1027 +1027 @@
   END SUBROUTINE var_zmla_generic
 
-  SUBROUTINE var_zwind(d1, u, v, z, u10, v10, sa, ca, topo, newz, unewz, vnewz)
+  SUBROUTINE var_zwind(d1, u, v, z, u10, v10, sa, ca, newz, unewz, vnewz)
 ! Subroutine to extrapolate the wind at a given height following the 'power law' methodology
 !    wss[newz] = wss[z1]*(newz/z1)**alpha
@@ -1039 +1039 @@
     INTEGER, INTENT(in)                                  :: d1
     REAL(r_k), DIMENSION(d1), INTENT(in)                 :: u,v,z
-    REAL(r_k), INTENT(in)                                :: u10, v10, topo, sa, ca, newz
+    REAL(r_k), INTENT(in)                                :: u10, v10, sa, ca, newz
     REAL(r_k), INTENT(out)                               :: unewz, vnewz
 
@@ -1049 +1049 @@
 !!!!!!! Variables
 ! u,v: vertical wind components [ms-1]
-! z: height above surface [m]
+! z: height above surface on half-mass levels [m]
 ! u10,v10: 10-m wind components [ms-1]
-! topo: topographical height [m]
 ! sa, ca: local sine and cosine of map rotation [1.]
 ! newz: desired height above grpund of extrapolation
@@ -1058 +1057 @@
     fname = 'var_zwind'
 
-    PRINT *,' ilev zaground newz z[ilev+1] z[ilev+2] _______'
+    !PRINT *,' ilev zaground newz z[ilev+1] z[ilev+2] _______'
     IF (z(1) < newz ) THEN
       DO inear = 1,d1-2
         zaground = z(inear+2)
-        PRINT *, inear, z(inear), newz, z(inear+1), z(inear+2)
+        !PRINT *, inear, z(inear), newz, z(inear+1), z(inear+2)
         IF ( zaground >= newz) EXIT
       END DO
     ELSE
-      PRINT *, 1, z(1), newz, z(2), z(3), ' z(1) > newz'
+      !PRINT *, 1, z(1), newz, z(2), z(3), ' z(1) > newz'
       inear = d1 - 2
     END IF
@@ -1083 +1082 @@
        v2(1) = u(inear+1)
        v2(2) = v(inear+1)
-       zz(1) = z(inear) - topo
-       zz(2) = z(inear+1) - topo
+       zz(1) = z(inear)
+       zz(2) = z(inear+1)
     END IF
 
     ! Computing for each component
     alpha = (LOG(ABS(v2))-LOG(ABS(v1)))/(LOG(zz(2))-LOG(zz(1)))
-    PRINT *,' Computing with v1:', v1, ' ms-1 v2:', v2, ' ms-1'
-    PRINT *,' z1:', zz(1), 'm z2:', zz(2), ' m'
-    PRINT *,' alhpa u:', alpha(1), ' alpha 2:', alpha(2)
+    !PRINT *,' Computing with v1:', v1, ' ms-1 v2:', v2, ' ms-1'
+    !PRINT *,' z1:', zz(1), 'm z2:', zz(2), ' m'
+    !PRINT *,' alhpa u:', alpha(1), ' alpha 2:', alpha(2)
     
     uvnewz = v1*(newz/zz(1))**alpha
@@ -1098 +1097 @@
     vnewz = uvnewz(1)*sa + uvnewz(2)*ca
 
-    PRINT *,'  result vz:', uvnewz
+    !PRINT *,'  result vz:', uvnewz
 
     !STOP