Changeset 3719 for LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_var_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_var_mod.F90 (modified) (3 diffs)

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

@@ -135 +135 @@
   ! Ocean-atmosphere interface, subskin ocean and near-surface ocean:
   
-  REAL, ALLOCATABLE, SAVE:: t_int(:) ! interface temperature, in K
+  REAL, ALLOCATABLE, SAVE:: sst_nff(:)
+  ! SST not used to compute surface fluxes, in K. If
+  ! activate_ocean_skin == 0 then it is not allocated; if
+  ! activate_ocean_skin == 1 then it is the ocean-air interface
+  ! temperature; if activate_ocean_skin == 2 then it is the bulk SST.
 
   REAL, ALLOCATABLE, SAVE:: dter(:)
@@ -155 +159 @@
   ! momentum flux due to rain, in Pa
   
-  !$OMP THREADPRIVATE(t_int, dter, dser, tkt, tks, taur)
+  !$OMP THREADPRIVATE(sst_nff, dter, dser, tkt, tks, taur)
 
 CONTAINS
@@ -216 +220 @@
     IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
 
-    if (activate_ocean_skin >= 1) allocate(t_int(klon), dter(klon), &
+    if (activate_ocean_skin >= 1) allocate(sst_nff(klon), dter(klon), &
          dser(klon), tkt(klon), tks(klon), taur(klon))