Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin

Timestamp:

Jul 29, 2024, 11:01:04 PM (5 months ago)

Author:

abarral

Message:

Put YOMCST.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/sens_heat_rain_m.F90 (modified) (2 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/sens_heat_rain_m.F90

@@ -25 +25 @@
 
 #ifdef IN_LMDZ
-    include "YOMCST.h"
-    ! for eps_w
+    USE lmdz_yomcst, ONLY: eps_w
 #endif
 
-    REAL, INTENT(IN):: rain ! rain mass flux, in kg m-2 s-1
-    REAL, INTENT(IN):: t ! air temperature at 10 m, in K
-    REAL, INTENT(IN):: q ! specific humidity at 10 m
-    REAL, INTENT(IN):: rhoa ! density of moist air  (kg / m3)
-    REAL, INTENT(IN):: xlv ! latent heat of evaporation (J / kg)
-    REAL, INTENT(IN):: t_int ! interface temperature, in K
-    REAL, INTENT(IN):: p ! surface pressure, in Pa
+    REAL, INTENT(IN) :: rain ! rain mass flux, in kg m-2 s-1
+    REAL, INTENT(IN) :: t ! air temperature at 10 m, in K
+    REAL, INTENT(IN) :: q ! specific humidity at 10 m
+    REAL, INTENT(IN) :: rhoa ! density of moist air  (kg / m3)
+    REAL, INTENT(IN) :: xlv ! latent heat of evaporation (J / kg)
+    REAL, INTENT(IN) :: t_int ! interface temperature, in K
+    REAL, INTENT(IN) :: p ! surface pressure, in Pa
 
     ! Local:
-    
+
     REAL es ! saturation pressure of wator vapor, in Pa
     REAL alfac ! wet bulb factor
@@ -58 +57 @@
     t_celsius = t - 273.15
     dtmp = (1. + 3.309e-3 * t_celsius - 1.44e-6 * t_celsius**2) * 0.02411 &
-         / (rhoa * cpa)
+            / (rhoa * cpa)
 
     ! Gosnell 1995 k0991, equation (11):
-    alfac =  1. / (1. + (wetc * xlv * dwat) / (cpa * dtmp))
+    alfac = 1. / (1. + (wetc * xlv * dwat) / (cpa * dtmp))
 
     ! Gosnell 1995 k0991, equation (12):
-    sens_heat_rain =  rain * alfac * cpw * (t_int - t + (q_int - q) * xlv / cpa)
+    sens_heat_rain = rain * alfac * cpw * (t_int - t + (q_int - q) * xlv / cpa)
 
   END FUNCTION sens_heat_rain
-  
+
 END MODULE sens_heat_rain_m