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

Timestamp:

Jul 1, 2020, 6:57:48 PM (4 years ago)

Author:

lguez

Message:

Store delta_sst instead of sst_nff

Store as a state variable the difference between ocean-air interface
temperature and bulk SST instead of sst_nff, which can be either
interface temperature or bulk SST. This is clearer. Also, it is
analoguous to what we will do with salinity.

So replace the two dummy arguments tsurf_in and sst_nff of
procedure cpl_send_ocean_fields by a single dummy argument
delta_sst. Replace dummy argument sst_nff of procedures
ocean_cpl_noice and surf_ocean by dummy argument
delta_sst. Replace variable sst_nff of module phys_state_var_mod
by variable delta_sst. Rename local variable ysst_nff of procedure
pbl_surface to ydelta_sst. Set variable delta_sst of module
phys_state_var_mod to 0 for an appearing ocean fraction and a
missing startup field. Replace variable o_sst_nff of module
phys_output_ctrlout_mod by variable o_delta_sst.

Rename variables cpl_delta_temp and cpl_delta_temp_2D of module
cpl_mod to cpl_delta_sst and cpl_delta_sst_2D, clearer. Rename
variable ids_delta_temp of module oasis to ids_delta_sst. Change
infosend(ids_delta_temp)%name to "CODELSST".

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, sst_nff
+       ale_wake, ale_bl_stat, s_int, ds_ns, dt_ns, delta_sst
 !FC
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
@@ -544 +544 @@
         found = phyetat0_get(1, s_int, "S_int", &
              "salinity at air-ocean interface", 35.)
-        found = phyetat0_get(1, sst_nff, "T1", "bulk SST of coupled ocean", 0.)
-        if (.not. found) sst_nff = ftsol(:,is_oce)
+        found = phyetat0_get(1, delta_sst, "delta_SST", &
+             "ocean-air interface temperature minus bulk SST", 0.)
      end if
      
@@ -556 +556 @@
         ds_ns = missing_val
         dt_ns = missing_val
-        sst_nff = missing_val
+        delta_sst = missing_val
      end where
   end if