Changeset 2387 in lmdz_wrf for trunk/tools/module_ForDiagnostics.f90

Timestamp:

Mar 11, 2019, 11:40:11 AM (6 years ago)

Author:

lfita

Message:

Adding:

• `tws': Wet Bulb temperature after Stull, 2011

File:

• trunk/tools/module_ForDiagnostics.f90 (modified) (3 diffs)

trunk/tools/module_ForDiagnostics.f90

@@ -35 +35 @@
 ! compute_psl_ecmwf: Compute sea level pressure using ECMWF method following Mats Hamrud and Philippe Courtier [Pa]
 ! compute_range_faces: Subroutine to compute faces [uphill, valleys, downhill] of a mountain range along a given face
+! compute_tws_RK[1/2/3/4]: Subroutine to compute Wet Bulb temperature of 1/2/3/4D series of values
 ! compute_vertint1D: Subroutine to vertically integrate a 1D variable in any vertical coordinates
 ! compute_zint4D: Subroutine to vertically integrate a 4D variable in any vertical coordinates
@@ -845 +846 @@
 ! Subroutine to compute fog: qcloud + qice /= 0.
 ! And visibility following Kunkel, B. A., (1984): Parameterization of droplet terminal velocity and 
-!   extinction coefficient in fog models. J. Climate Appl. Meteor., 23, 34–41.
+!   extinction coefficient in fog models. J. Climate Appl. Meteor., 23, 34-41.
 
     IMPLICIT NONE
@@ -1357 +1358 @@
   END SUBROUTINE compute_Koeppen_Geiger_climates
 
+  SUBROUTINE compute_tws_RK1(d1, tas, hurs, tws)
+! Subroutine to compute Wet Bulb temperature of 1D series of values using equation after:
+!    Stull, R. (2011), J. Appl. Meteor. Climatol. 50(11):2267-2269. doi: 10.1175/JAMC-D-11-0143.1
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1
+    REAL(r_k), DIMENSION(d1), INTENT(in)                 :: tas, hurs
+    REAL(r_k), DIMENSION(d1), INTENT(out)                :: tws
+ 
+! Local
+    INTEGER                                              :: it
+
+!!!!!!! Variables
+! tas: 2-m air temperature [K]
+! hurs: 2-m relative humidity [1]
+
+    fname = 'compute_tws_RK1'
+
+    DO it=1, d1
+      tws(it) = var_tws_S11(tas(it), hurs(it))
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_tws_RK1
+
+  SUBROUTINE compute_tws_RK2(d1, d2, tas, hurs, tws)
+! Subroutine to compute Wet Bulb temperature of 2D series of values using equation after:
+!    Stull, R. (2011), J. Appl. Meteor. Climatol. 50(11):2267-2269. doi: 10.1175/JAMC-D-11-0143.1
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2
+    REAL(r_k), DIMENSION(d1,d2), INTENT(in)              :: tas, hurs
+    REAL(r_k), DIMENSION(d1,d2), INTENT(out)             :: tws
+ 
+! Local
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! tas: 2-m air temperature [K]
+! hurs: 2-m relative humidity [1]
+
+    fname = 'compute_tws_RK2'
+
+    DO i=1, d1
+      DO j=1, d2
+        tws(i,j) = var_tws_S11(tas(i,j), hurs(i,j))
+      END DO
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_tws_RK2
+
+  SUBROUTINE compute_tws_RK3(d1, d2, d3, tas, hurs, tws)
+! Subroutine to compute Wet Bulb temperature of 3D series of values using equation after:
+!    Stull, R. (2011), J. Appl. Meteor. Climatol. 50(11):2267-2269. doi: 10.1175/JAMC-D-11-0143.1
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2, d3
+    REAL(r_k), DIMENSION(d1,d2,d3), INTENT(in)           :: tas, hurs
+    REAL(r_k), DIMENSION(d1,d2,d3), INTENT(out)          :: tws
+ 
+! Local
+    INTEGER                                              :: i, j, k
+
+!!!!!!! Variables
+! tas: 2-m air temperature [K]
+! hurs: 2-m relative humidity [1]
+
+    fname = 'compute_tws_RK3'
+
+    DO i=1, d1
+      DO j=1, d2
+        DO k=1, d3
+          tws(i,j,k) = var_tws_S11(tas(i,j,k), hurs(i,j,k))
+        END DO
+      END DO
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_tws_RK3
+
+  SUBROUTINE compute_tws_RK4(d1, d2, d3, d4, tas, hurs, tws)
+! Subroutine to compute Wet Bulb temperature of 4D series of values using equation after:
+!    Stull, R. (2011), J. Appl. Meteor. Climatol. 50(11):2267-2269. doi: 10.1175/JAMC-D-11-0143.1
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2, d3, d4
+    REAL(r_k), DIMENSION(d1,d2,d3,d4), INTENT(in)        :: tas, hurs
+    REAL(r_k), DIMENSION(d1,d2,d3,d4), INTENT(out)       :: tws
+ 
+! Local
+    INTEGER                                              :: i,j,k,l
+
+!!!!!!! Variables
+! tas: 2-m air temperature [K]
+! hurs: 2-m relative humidity [1]
+
+    fname = 'compute_tws_RK4'
+
+    DO i=1, d1
+      DO j=1, d2
+       DO k=1, d3
+         DO l=1, d4
+           tws(i,j,k,l) = var_tws_S11(tas(i,j,k,l), hurs(i,j,k,l))
+         END DO
+        END DO
+      END DO
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_tws_RK4
+
 END MODULE module_ForDiagnostics