Changeset 3429 for LMDZ6/branches/Ocean_skin/libf/phylmd/phyetat0.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/phyetat0.F90 (modified) (4 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
+       ale_wake, ale_bl_stat, ds_ns, dt_ns
 !FC
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
@@ -29 +29 @@
   USE ocean_slab_mod, ONLY: nslay, tslab, seaice, tice, ocean_slab_init
   USE time_phylmdz_mod, ONLY: init_iteration, pdtphys, itau_phy
+  use netcdf, only: nf90_fill_real
 
   IMPLICIT none
@@ -35 +36 @@
   ! Objet: Lecture de l'etat initial pour la physique
   !======================================================================
-  include "netcdf.inc"
   include "dimsoil.h"
   include "clesphys.h"
@@ -529 +529 @@
   CALL fonte_neige_init(run_off_lic_0)
 
+  found = phyetat0_get(1, ds_ns, "ds_ns", "delta salinity near surface", 0.)
+  found = phyetat0_get(1, dt_ns, "dT_ns", "delta temperature near surface", 0.)
+
+  where (pctsrf(:, is_oce) == 0.)
+     ds_ns = nf90_fill_real
+     dt_ns = nf90_fill_real
+  end where
+
 END SUBROUTINE phyetat0