Changeset 3767 for LMDZ6/branches/Ocean_skin/libf/phylmd/phyetat0.F90

Timestamp:

Jul 19, 2020, 5:25:06 PM (4 years ago)

Author:

lguez

Message:

Store delta_sal instead of s_int

Revision analoguous to revision [3744], for salinity. Store as a state
variable and send to the ocean delta_sal, the difference between
ocean-air interface salinity and bulk salinity, instead of s_int,
the interface salinity.

So replace dummy argument s_int of procedure cpl_send_ocean_fields
by dummy argument delta_sal. Replace dummy argument s_int of
procedures ocean_cpl_noice and surf_ocean by dummy argument
delta_sal. Replace variable s_int of module phys_state_var_mod
by variable delta_sal. Rename local variable ys_int of procedure
pbl_surface to ydelta_sal. Set variable delta_sal of module
phys_state_var_mod to 0 for an appearing ocean fraction and a
missing startup field. Replace variable o_s_int of module
phys_output_ctrlout_mod by variable o_delta_sal.

Rename variables cpl_s_int and cpl_s_int_2D of module cpl_mod to
cpl_delta_sal and cpl_delta_sal_2D. Rename variable ids_s_int of
module oasis to ids_delta_sal. Change infosend(ids_delta_sal)%name
to "CODELSSS".

File:

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

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

@@ -19 +19 @@
        wake_s, wake_dens, zgam, zmax0, zmea, zpic, zsig, &
        zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m, treedrg, &
-       ale_wake, ale_bl_stat, s_int, ds_ns, dt_ns, delta_sst
+       ale_wake, ale_bl_stat, ds_ns, dt_ns, delta_sst, delta_sal
 !FC
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
@@ -525 +525 @@
   if (activate_ocean_skin >= 1) then
      if (activate_ocean_skin == 2 .and. type_ocean == 'couple') then
-        found = phyetat0_get(1, s_int, "S_int", &
-             "salinity at air-ocean interface", 35.)
+        found = phyetat0_get(1, delta_sal, "delta_sal", &
+             "ocean-air interface salinity minus bulk salinity", 0.)
         found = phyetat0_get(1, delta_sst, "delta_SST", &
              "ocean-air interface temperature minus bulk SST", 0.)
@@ -536 +536 @@
 
      where (pctsrf(:, is_oce) == 0.)
-        s_int = missing_val
         ds_ns = missing_val
         dt_ns = missing_val
         delta_sst = missing_val
+        delta_sal = missing_val
      end where
   end if