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

Timestamp:

Jan 2, 2019, 7:44:20 PM (5 years ago)

Author:

lguez

Message:

Create subdirectory Ocean_skin in libf. For now, Ocean_skin is under
control of git, not subversion.

Add variable eps_w to common YOMCST.

For now, continue to read ocean skin parameters from a namelist in
config_ocean_skin.

The parameterisation is called from procedure surf_ocean.

Add two prognostic variables for the parameterisation: dt_ns and
ds_ns. Add eight diagnostic variables: t_int, s_int, dter, dser, tkt,
tks, rf, taur. The ten variables are only defined on ocean surface,
elsewhere they are set to nf90_fill_real. In pbl_surface, we can
initialize the eight diagnostic variables to nf90_fill_real before the
loop on sub-surfaces, but we need to keep the old values of dt_ns and
ds_ns as input of the parameterisation so we set dt_ns and ds_ns to
nf90_fill_real after the call to surf_ocean. Define ten corresponding
compressed variables in pbl_surface. Define ten corresponding NetCDF
output variables in phys_output_ctrlout_mod.

In procedure pbl_surface_newfrac, for an appearing ocean sub-surface,
dt_ns and ds_ns are set to 0. In phyetat0, also set initial dt_ns and
ds_ns to 0 if they are not in start file.

In procedure surf_ocean, for now, we use a constant specific latent
heat of vaporization, as elsewhere in LMDZ, and a constant bulk
salinity.

File:

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

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

@@ -196 +196 @@
 ! Tropopause
          o_p_tropopause, o_z_tropopause, o_t_tropopause,  &
-         o_col_O3_strato, o_col_O3_tropo               ! Added ThL
-
+         o_col_O3_strato, o_col_O3_tropo, &               ! Added ThL
+         o_t_int, o_s_int, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, o_tks, &
+         o_rf, o_taur
 
 #ifdef CPP_StratAer
@@ -240 +241 @@
          ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
          rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, &
-         vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90
+         vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90, ds_ns, dt_ns
 
     USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
@@ -348 +349 @@
          alt_tropo, &
 !Ionela
-         ok_4xCO2atm
+         ok_4xCO2atm, t_int, s_int, dter, dser, tkt, tks, rf, taur
 
     USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
@@ -2110 +2111 @@
        ENDIF   !(iflag_phytrac==1)
 
+       CALL histwrite_phy(o_t_int, t_int)
+       CALL histwrite_phy(o_s_int, s_int)
+       CALL histwrite_phy(o_ds_ns, ds_ns)
+       CALL histwrite_phy(o_dt_ns, dt_ns)
+       CALL histwrite_phy(o_dter, dter)
+       CALL histwrite_phy(o_dser, dser)
+       CALL histwrite_phy(o_tkt, tkt)
+       CALL histwrite_phy(o_tks, tks)
+       CALL histwrite_phy(o_rf, rf)
+       CALL histwrite_phy(o_taur, taur)
 
        IF (.NOT.vars_defined) THEN