Changeset 3628 for LMDZ6/branches/Ocean_skin/libf/phylmd/ocean_cpl_mod.F90

Timestamp:

Feb 4, 2020, 10:36:32 PM (5 years ago)

Author:

lguez

Message:

If the ocean skin parameterization is working actively
(activate_ocean_skin == 2) and we are coupled to the ocean then send
ocean-air interface salinity to the ocean. New dummy argument s_int
of procedures ocean_cpl_noice and cpl_send_ocean_fields. We can
only send interface salinity from the previous time-step since
communication with the ocean is before the call to bulk_flux. So make
s_int a state variable: move s_int from phys_output_var_mod to
phys_state_var_mod. Still, we only read s_int from startphy,
define it before the call to surf_ocean and write it to restartphy
if activate_ocean_skin == 2 and type_ocean == 'couple'. In
procedure pbl_surface, for clarity, move the definition of output
variables t_int, dter, dser, tkt, tks, rf, taur to missing_val to
after the call to surf_ocean, with the definition of s_int,
ds_ns, dt_ns to missing_val. This does not change anything for
t_int, dter, dser, tkt, tks, rf, taur. In pbl_surface_newfrac, we
choose to set s_int to 35 for an appearing ocean point, this is
questionable. In surf_ocean, change the intent of s_int from out
to inout.

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/ocean_cpl_mod.F90 (modified) (3 diffs)

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

@@ -55 +55 @@
        radsol, snow, agesno, &
        qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-       tsurf_new, dflux_s, dflux_l, sss)
+       tsurf_new, dflux_s, dflux_l, sss, s_int)
 
 !
@@ -92 +92 @@
     REAL, DIMENSION(klon), INTENT(IN)        :: u1, v1, gustiness
     REAL, DIMENSION(klon), INTENT(IN)        :: tsurf_in
+    real, intent(in):: s_int(:) ! (knon) ocean-air interface salinity, in ppt
 
 ! In/Output arguments
@@ -205 +206 @@
 !****************************************************************************************
 
-    CALL cpl_send_ocean_fields(itime, knon, knindex, &
-         swnet, lwnet, fluxlat, fluxsens, &
-         precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, flux_u1, flux_v1, windsp,&
-         sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, tsurf_in)
-    
+    CALL cpl_send_ocean_fields(itime, knon, knindex, swnet, lwnet, fluxlat, &
+         fluxsens, precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, &
+         flux_u1, flux_v1, windsp, sens_prec_liq, sens_prec_sol, lat_prec_liq, &
+         lat_prec_sol, tsurf_in, s_int)
 
   END SUBROUTINE ocean_cpl_noice