Changeset 2538 for LMDZ5/trunk/libf/phylmd/calcul_fluxs_mod.F90

Timestamp:

Jun 3, 2016, 4:12:16 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Computation of heat fluxes associated with solid and liquid precipitations
over ocean and seaice. Quantities are sent to the coupler
LF

File:

• LMDZ5/trunk/libf/phylmd/calcul_fluxs_mod.F90 (modified) (6 diffs, 1 prop)

LMDZ5/trunk/libf/phylmd/calcul_fluxs_mod.F90

@@ +1 @@
+!
+! $Id$
 !
 MODULE calcul_fluxs_mod
@@ -9 +11 @@
        radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, gustiness, &
        fqsat, petAcoef, peqAcoef, petBcoef, peqBcoef, &
-       tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
+       tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
+       sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol)
+ 
     
     USE dimphy, ONLY : klon
@@ -49 +53 @@
 !   dflux_s      derivee du flux de chaleur sensible / Ts
 !   dflux_l      derivee du flux de chaleur latente  / Ts
-!
+!   sens_prec_liq flux sensible lié aux echanges de precipitations liquides
+!   sens_prec_sol                                    precipitations solides
+!   lat_prec_liq  flux latent lié aux echanges de precipitations liquides
+!   lat_prec_sol                                  precipitations solides
 
     INCLUDE "YOETHF.h"
@@ -77 +84 @@
     REAL, DIMENSION(klon), INTENT(OUT)   :: tsurf_new, evap, fluxsens, fluxlat
     REAL, DIMENSION(klon), INTENT(OUT)   :: dflux_s, dflux_l
+    REAL, DIMENSION(klon), OPTIONAL      :: sens_prec_liq, sens_prec_sol
+    REAL, DIMENSION(klon), OPTIONAL      :: lat_prec_liq, lat_prec_sol
 
 ! Variables locales
@@ -130 +139 @@
     dflux_s = 0.
     dflux_l = 0.
+    if (PRESENT(sens_prec_liq)) sens_prec_liq = 0.
+    if (PRESENT(sens_prec_sol)) sens_prec_sol = 0.
+    if (PRESENT(lat_prec_liq)) lat_prec_liq = 0.
+    if (PRESENT(lat_prec_sol)) lat_prec_sol = 0.
 !
 ! zx_qs = qsat en kg/kg
@@ -243 +256 @@
 !!$     &  run_off(i) = run_off(i) + max(qsol(i) - max_eau_sol, 0.0)
 !!$    qsol(i) = min(qsol(i), max_eau_sol) 
+!
+! calcul de l'enthalpie des precipitations liquides et solides
+!
+!       if (PRESENT(enth_prec_liq))                   &
+!       enth_prec_liq(i) = rcw * (t1lay(i) - tsurf(i)) * &
+!                          precip_rain(i)
+!       if (PRESENT(enth_prec_sol))                  &
+!       enth_prec_sol(i) = rcs * (t1lay(i) - tsurf(i)) * &
+!                          precip_snow(i)
+! On calcule par rapport a T=0
+       if (PRESENT(sens_prec_liq))                   &
+         sens_prec_liq(i) = rcw * (t1lay(i) - RTT) * &
+                          precip_rain(i)
+       if (PRESENT(sens_prec_sol))                   &
+         sens_prec_sol(i) = rcs * (t1lay(i) - RTT) * &
+                          precip_snow(i)
+       if (PRESENT(lat_prec_liq))                    &
+         lat_prec_liq(i) =  precip_rain(i) * (RLVTT - RLVTT) 
+       if (PRESENT(lat_prec_sol))                    &
+         lat_prec_sol(i) = precip_snow(i) * (RLSTT - RLVTT)
     ENDDO
+
+    
+!       if (PRESENT(sens_prec_liq))                  &
+!          WRITE(*,*)' calculs_fluxs sens_prec_liq (min, max)', &
+!          MINVAL(sens_prec_liq(1:knon)), MAXVAL(sens_prec_liq(1:knon))
+!       if (PRESENT(sens_prec_sol))                  &
+!          WRITE(*,*)' calculs_fluxs sens_prec_sol (min, max)', &
+!          MINVAL(sens_prec_sol(1:knon)), MAXVAL(sens_prec_sol(1:knon))
+ 
 !
 !****************************************************************************************