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

Timestamp:

Feb 4, 2020, 10:36:32 PM (4 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/phys_state_var_mod.F90 (modified) (4 diffs)

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

@@ -431 +431 @@
 
       ! Ocean-atmosphere interface, subskin ocean and near-surface ocean:
+
+      REAL, ALLOCATABLE, SAVE:: s_int(:) ! interface salinity, in ppt
       
       REAL, ALLOCATABLE, SAVE:: ds_ns(:)
@@ -442 +444 @@
       ! minus foundation temperature. (Can be negative.) In K.
       
-      !$OMP THREADPRIVATE(ds_ns, dt_ns)
+      !$OMP THREADPRIVATE(s_int, ds_ns, dt_ns)
 
     CONTAINS
@@ -650 +652 @@
         du_gwd_front(:,:) = 0 !ym missing init   
       ENDIF
-      if (activate_ocean_skin >= 1) ALLOCATE(ds_ns(klon), dt_ns(klon))
+      if (activate_ocean_skin >= 1) ALLOCATE(s_int(klon), ds_ns(klon), &
+           dt_ns(klon))
 
     END SUBROUTINE phys_state_var_init
@@ -788 +791 @@
 !!! fin nrlmd le 10/04/2012
 
-      if (activate_ocean_skin >= 1) deALLOCATE(ds_ns, dt_ns)
+      if (activate_ocean_skin >= 1) deALLOCATE(s_int, ds_ns, dt_ns)
 
       is_initialized=.FALSE.