Changeset 3719 for LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Jun 12, 2020, 9:31:02 PM (4 years ago)

Author:

lguez

Message:

Change name and meaning of t_int to sst_nff

Rename variable t_int of module phys_output_var_mod to
sst_nff. Rename local variable yt_int of pbl_surface to
ysst_nff. Rename type o_t_int of module phys_output_ctrlout_mod
to o_sst_nff. In procedure phys_output_write, write sst_nff if
activate_ocean_skin >= 1, while t_int was written only for
activate_ocean_skin == 1.

In procedure surf_ocean, rename dummy argument t_int to
sst_nff. Create local variable t_int. Set sst_nff to t_int or
tsurf_new depending on activate_ocean_skin.

Motivation: we want to output the bulk SST if activate_ocean_skin

2. Before this revision, the bulk SST in variable tsurf_new in

procedure surf_ocean was overwritten with interface temperature and
lost. We already had two output variables for SST: t_int and
ftsol, containing the same values if activate_ocean_skin

2. Instead of creating a third variable for bulk SST, change the

name and meaning of t_int to sst_nff, the SST not used to compute
surface fluxes. That is either the interface temperature or the bulk
SST depending on activate_ocean_skin. The SST used to compute
surface fluxes is always ftsol, meaning either the interface
temperature or the bulk SST depending on activate_ocean_skin.

File:

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

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

@@ -202 +202 @@
 !--interactive CO2
          o_flx_co2_ocean, o_flx_co2_land, o_flx_co2_ff, o_flx_co2_bb, &
-         o_t_int, o_s_int, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, o_tks, &
+         o_sst_nff, o_s_int, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, o_tks, &
          o_taur
 
@@ -360 +360 @@
          alt_tropo, &
 !Ionela
-         ok_4xCO2atm, t_int, dter, dser, tkt, tks, taur
+         ok_4xCO2atm, sst_nff, dter, dser, tkt, tks, taur
 
     USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
@@ -2304 +2304 @@
 
        if (activate_ocean_skin >= 1) then
-          if (activate_ocean_skin == 1) CALL histwrite_phy(o_t_int, t_int)
-          ! (no need to duplicate with netCDF variable "tsol_oce" if
-          ! activate_ocean_skin == 2)
-          
+          CALL histwrite_phy(o_sst_nff, sst_nff)
           CALL histwrite_phy(o_s_int, s_int)
           CALL histwrite_phy(o_ds_ns, ds_ns)