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

Timestamp:

Feb 3, 2020, 2:27:46 PM (4 years ago)

Author:

lguez

Message:

If the ocean skin parameterization is working (passively or actively,
activate_ocean_skin >= 1) and we are coupled to the ocean then
receive bulk salinity of the surface layer of the ocean from the ocean
and feed it to procedure bulk_flux instead of the constant
value 35. If the ocean skin parameterization is working actively
(activate_ocean_skin == 2) and we are coupled to the ocean then send
ocean-air interface temperature to the ocean. We can only send
interface temperature from the previous time-step since communication
with the ocean is before the call to bulk_flux. In module cpl_mod,
define cpl_t_int with rank 1: no dimension for cpl_index because
t_int is only defined over ocean. New dummy argument sss of
procedures cpl_receive_ocean_fields and ocean_cpl_noice. New dummy
argument t_int of cpl_send_ocean_fields. In procedure
surf_ocean, rename local variable s1 to sss and give it the size
klon, which is required by the coupling machinery.

File:

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

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

@@ -57 +57 @@
     REAL, DIMENSION(klon), INTENT(IN)        :: rmu0  
     REAL, DIMENSION(klon), INTENT(IN)        :: fder
-    REAL, DIMENSION(klon), INTENT(IN)        :: tsurf_in
+    REAL, INTENT(IN):: tsurf_in(klon) ! defined only for subscripts 1:knon
     REAL, DIMENSION(klon), INTENT(IN)        :: p1lay,z1lay ! pression (Pa) et altitude (m) du premier niveau
     REAL, DIMENSION(klon), INTENT(IN)        :: cdragh
@@ -102 +102 @@
     REAL, DIMENSION(klon), INTENT(OUT)       :: flux_u1, flux_v1
 
-    REAL, intent(out):: t_int(:) ! (knon) interface temperature, in K
-    real, intent(out):: s_int(:) ! (knon) interface salinity, in ppt
+    REAL, intent(out):: t_int(:) ! (knon) ocean-air interface temperature, in K
+    real, intent(out):: s_int(:) ! (knon) ocean-air interface salinity, in ppt
 
     REAL, intent(out):: dter(:) ! (knon)
@@ -111 +111 @@
 
     REAL, intent(out):: dser(:) ! (knon)
-    ! Temperature variation in the diffusive microlayer, that is
-    ! subskin temperature minus ocean-air interface temperature. In K.
+    ! Salinity variation in the diffusive microlayer, that is
+    ! ocean-air interface salinity minus subskin salinity. In ppt.
 
     REAL, intent(out):: tkt(:) ! (knon)
@@ -140 +140 @@
     real xlv(knon) ! chaleur latente d'évaporation (J / kg)
     real precip_tot(knon) ! rain + snow
-    real S1(knon) ! salinity at depth_1, in ppt
+
+    real sss(klon)
+    ! Bulk salinity of the surface layer of the ocean, in ppt. (Only
+    ! defined for subscripts 1:knon, but we have to declare it with
+    ! size klon because of the coupling machinery.)
 
     ! End definition
@@ -186 +190 @@
             radsol, snow, agesno, &
             qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-            tsurf_new, dflux_s, dflux_l)
+            tsurf_new, dflux_s, dflux_l, sss)
 
     CASE('slab')
@@ -319 +323 @@
        rf =  sens_heat_rain(precip_tot, temp_air(:knon), spechum(:knon), rhoa, &
             xlv, tsurf_in(:knon), ps(:knon))
-       s1 = 35.
+       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 = s1, &
+            u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = sss(:knon), &
             rain = precip_tot, hf = - fluxsens(:knon), hlb = - fluxlat(:knon), &
             rnl = - lwnet(:knon), &