Changeset 3687 for LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90

Timestamp:

May 27, 2020, 4:52:16 PM (5 years ago)

Author:

lguez

Message:

Modify sensible heat due to rain sent to the ocean

Modify the sensible heat flux due to rain which is sent to the
ocean. Replace the computation of sens_prec_liq in procedure
calcul_fluxs by a call to sens_heat_rain. Set sens_prec_sol in
procedure calcul_fluxs to 0 because, for now, sens_heat_rain is
supposed to account for both rain and snow.

For the call to sens_heat_rain in procedure calcul_fluxs, we need
an additional dummy argument rhoa of calcul_fluxs. Add dummy
argument rhoa to ocean_cpl_noice, ocean_forced_noice,
ocean_forced_ice and ocean_cpl_ice because we need to pass it down
to calcul_fluxs.

Change the dimension of sens_prec_liq and sens_prec_sol in
procedures calcul_fluxs, ocean_cpl_noice, ocean_cpl_ice,
ocean_forced_noice, ocean_forced_ice, cpl_send_ocean_fields and
cpl_send_seaice_fields to knon.

In procedures ocean_forced_noice and ocean_cpl_noice, promote
sens_prec_liq from local variable to dummy argument because we need
it in surf_ocean. Remove useless initialization of sens_prec_liq
and sens_prec_sol in ocean_cpl_noice, ocean_cpl_ice,
ocean_forced_ice and ocean_forced_noice: they are intent out in
calcul_fluxs.

Remove variable rf of module phys_output_var_mod, we use
sens_prec_liq instead. Remove local variable yrf of procedure
pbl_surface. rf and yrf appeared in pbl_surface only to be output.
Remove variable o_rf of module phys_output_ctrlout_mod. Remove
dummy argument rf of procedure surf_ocean.

Do not call sens_heat_rain in surf_ocean since we now call it
from calcul_fluxs. Move the computation of rhoa in surf_ocean
before the calls to ocean_cpl_noice and ocean_forced_noice.

Add the computation of rhoa in surf_seaice, to pass it down to
ocean_cpl_ice and ocean_forced_ice.

If activate_ocean_skin == 1 then the results are changed because the
call to sens_heat_rain in calcul_fluxs now uses the SST from the
current time step of physics. On this point, the present revision
reverses revision [3463].

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90 (modified) (8 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90

@@ -20 +20 @@
        z0m, z0h, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & 
        tsurf_new, dflux_s, dflux_l, lmt_bils, &
-       flux_u1, flux_v1, t_int, s_int, ds_ns, dt_ns, dter, dser, tkt, tks, rf, &
-       taur)
+       flux_u1, flux_v1, t_int, s_int, ds_ns, dt_ns, dter, dser, tkt, tks, taur)
 
     use albedo, only: alboc, alboc_cd
@@ -32 +31 @@
     USE indice_sol_mod, ONLY : nbsrf, is_oce
     USE limit_read_mod
-    use sens_heat_rain_m, only: sens_heat_rain
     use config_ocean_skin_m, only: activate_ocean_skin
     !
@@ -122 +120 @@
     ! épaisseur (m) de la couche de diffusion de masse (microlayer)
 
-    REAL, intent(out):: rf(:) ! (knon)
-    ! sensible heat flux at the surface due to rainfall, in  W m-2
-    ! positive upward
-
     REAL, intent(out):: taur(:) ! (knon)
     ! momentum flux due to rain, in Pa
@@ -139 +133 @@
     CHARACTER(len=20),PARAMETER :: modname="surf_ocean"
     real rhoa(knon) ! density of moist air  (kg / m3)
-    real xlv(knon) ! chaleur latente d'évaporation (J / kg)
-    real precip_tot(knon) ! rain + snow
+    REAL sens_prec_liq(knon)
 
     real sss(klon)
@@ -175 +168 @@
     ENDIF
 
-
+    rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon)))
     !******************************************************************************
     ! Switch according to type of ocean (couple, slab or forced)
@@ -191 +184 @@
             radsol, snow, agesno, &
             qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-            tsurf_new, dflux_s, dflux_l, sss, s_int)
+            tsurf_new, dflux_s, dflux_l, sens_prec_liq, sss, s_int, rhoa)
 
     CASE('slab')
@@ -214 +207 @@
             radsol, snow, agesno, &
             qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-            tsurf_new, dflux_s, dflux_l)
+            tsurf_new, dflux_s, dflux_l, sens_prec_liq, rhoa)
     END SELECT
 
@@ -319 +312 @@
 
     if (activate_ocean_skin >= 1) then
-       rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon)))
-       xlv = rlvtt
-       precip_tot = precip_rain(:knon) + precip_snow(:knon)
-       rf =  sens_heat_rain(precip_tot, temp_air(:knon), spechum(:knon), rhoa, &
-            xlv, tsurf_in(:knon), ps(:knon))
        if (type_ocean /= 'couple') sss(:knon) = 35.
        call bulk_flux(tkt, tks, taur, dter, dser, t_int, s_int, ds_ns, dt_ns, &
             u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = sss(:knon), &
-            rain = precip_tot, hf = - fluxsens(:knon), hlb = - fluxlat(:knon), &
+            rain = precip_rain(:knon) + precip_snow(:knon), &
+            hf = - fluxsens(:knon), hlb = - fluxlat(:knon), &
             rnl = - lwnet(:knon), &
             tau = sqrt(flux_u1(:knon)**2 + flux_v1(:knon)**2), rhoa = rhoa, &
-            xlv = xlv, rf = rf, dtime = dtime, rns = swnet(:knon))
+            xlv = [(rlvtt, i = 1, knon)], rf = - sens_prec_liq, dtime = dtime, &
+            rns = swnet(:knon))
        if (activate_ocean_skin == 2) tsurf_new(:knon) = t_int
     end if