Changeset 2254 for LMDZ5/trunk/libf/phylmd/ocean_forced_mod.F90

Timestamp:

Apr 3, 2015, 12:16:36 PM (9 years ago)

Author:

fhourdin

Message:

Modification du calcul des flux air/mer
1) Introduction d'un facteur f_qsat_oce=0.98 devant qsat dans le calcul
de l'évaporation sur océan pour tenir compte de la moindre évaporation
de l'eau salée.
2) Introduction d'une différentiation entre z0 pour le sensible, z0h,
et le latent, z0q, imposé constant z0q=f_z0qh_oce*z0h

Modification of air/sea fluxes computation
1) Introduction of a correcting factor f_qsat_oce=0.98 on qsat
to account for the weaker evaporation of salty water.
2) Introduction of z0q=f_z0qh_oce*z0h

File:

• LMDZ5/trunk/libf/phylmd/ocean_forced_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/ocean_forced_mod.F90

@@ -13 +13 @@
   SUBROUTINE ocean_forced_noice( &
        itime, dtime, jour, knon, knindex, &
-       p1lay, cdragh, cdragm, precip_rain, precip_snow, &
+       p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, &
        temp_air, spechum, &
        AcoefH, AcoefQ, BcoefH, BcoefQ, &
@@ -33 +33 @@
     USE indice_sol_mod
     INCLUDE "YOMCST.h"
+    INCLUDE "clesphys.h"
+
 
 ! Input arguments
@@ -40 +42 @@
     REAL, INTENT(IN)                         :: dtime
     REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
-    REAL, DIMENSION(klon), INTENT(IN)        :: cdragh, cdragm
+    REAL, DIMENSION(klon), INTENT(IN)        :: cdragh, cdragq, cdragm
     REAL, DIMENSION(klon), INTENT(IN)        :: precip_rain, precip_snow
     REAL, DIMENSION(klon), INTENT(IN)        :: temp_air, spechum
@@ -109 +111 @@
 ! Calcul de tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l and qsurf
     CALL calcul_fluxs(knon, is_oce, dtime, &
-         tsurf_lim, p1lay, cal, beta, cdragh, ps, &
+         tsurf_lim, p1lay, cal, beta, cdragh, cdragq, ps, &
          precip_rain, precip_snow, snow, qsurf,  &
          radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
-         AcoefH, AcoefQ, BcoefH, BcoefQ, &
+         f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, &
          tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
 
@@ -231 +233 @@
     v1_lay(:) = v1(:) - v0(:)
     CALL calcul_fluxs(knon, is_sic, dtime, &
-         tsurf_tmp, p1lay, cal, beta, cdragh, ps, &
+         tsurf_tmp, p1lay, cal, beta, cdragh, cdragh, ps, &
          precip_rain, precip_snow, snow, qsurf,  &
          radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
-         AcoefH, AcoefQ, BcoefH, BcoefQ, &
+         f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, &
          tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)