Changeset 3429 for LMDZ6/branches/Ocean_skin/libf/phylmd/pbl_surface_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/pbl_surface_mod.F90 (modified) (8 diffs)

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

@@ -290 +290 @@
     USE print_control_mod,  ONLY : prt_level,lunout
     USE ioipsl_getin_p_mod, ONLY : getin_p
+    use phys_state_var_mod, only: ds_ns, dt_ns
+    use phys_output_var_mod, only: t_int, s_int, dter, dser, tkt, tks, rf, taur
+    use netcdf, only: nf90_fill_real
 
     IMPLICIT NONE
@@ -821 +824 @@
     ! Martin
 
+    real, DIMENSION(klon):: yt_int, ys_int, yds_ns, ydt_ns, ydter, ydser, &
+         ytkt, ytks, yrf, ytaur
+    ! compression of t_int, s_int, ds_ns, dt_ns, dter, dser, tkt, tks, rf,
+    ! taur on ocean points
+
 !****************************************************************************************
 ! End of declarations
@@ -954 +962 @@
 !!! jyg le 10/02/2012
     rh2m_x(:) = 0. ; qsat2m_x(:) = 0. ; rh2m_w(:) = 0. ; qsat2m_w(:) = 0.
-!!!
+
+    t_int = nf90_fill_real
+    s_int = nf90_fill_real
+    dter = nf90_fill_real
+    dser = nf90_fill_real
+    tkt = nf90_fill_real
+    tks = nf90_fill_real
+    rf = nf90_fill_real
+    taur = nf90_fill_real
 
 ! 2b) Initialization of all local variables that will be compressed later
@@ -1414 +1430 @@
           ENDDO
        ENDIF
+
+       if (nsrf == is_oce) then
+          yds_ns(:knon) = ds_ns(ni(:knon))
+          ydt_ns(:knon) = dt_ns(ni(:knon))
+       end if
        
 !****************************************************************************************
@@ -1934 +1955 @@
                yz0m, yz0h, SFRWL,yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
                ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, &
-               y_flux_u1, y_flux_v1)
+               y_flux_u1, y_flux_v1, yt_int(:knon), ys_int(:knon), &
+               yds_ns(:knon), ydt_ns(:knon), ydter(:knon), ydser(:knon), &
+               ytkt(:knon), ytks(:knon), yrf(:knon), ytaur(:knon))
       IF (prt_level >=10) THEN
           print *,'arg de surf_ocean: ycdragh ',ycdragh
@@ -2538 +2561 @@
        ENDIF
 
+       if (nsrf == is_oce) then
+          dt_ns = nf90_fill_real
+          ds_ns = nf90_fill_real
+          ds_ns(ni(:knon)) = yds_ns(:knon)
+          dt_ns(ni(:knon)) = ydt_ns(:knon)
+
+          t_int(ni(:knon)) = yt_int(:knon)
+          s_int(ni(:knon)) = ys_int(:knon)
+          dter(ni(:knon)) = ydter(:knon)
+          dser(ni(:knon)) = ydser(:knon)
+          tkt(ni(:knon)) = ytkt(:knon)
+          tks(ni(:knon)) = ytks(:knon)
+          rf(ni(:knon)) = yrf(:knon)
+          taur(ni(:knon)) = ytaur(:knon)
+       end if
+
 !****************************************************************************************
 ! 14) Calculate the temperature and relative humidity at 2m and the wind at 10m 
@@ -3141 +3180 @@
 
 !albedo SB >>>
-SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &
-     evap, z0m, z0h, agesno,                                  &
-     tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke)  
-!albedo SB <<<
+  SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &
+       evap, z0m, z0h, agesno,                                  &
+       tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke)  
+    !albedo SB <<<
     ! Give default values where new fraction has appread
 
     USE indice_sol_mod
+    use phys_state_var_mod, only: ds_ns, dt_ns
 
     INCLUDE "dimsoil.h"
@@ -3234 +3274 @@
                       alb_dif(i,k,nsrf) = 0.06
                    ENDDO
+                   ds_ns(i) = 0.
+                   dt_ns(i) = 0.
                 ELSE IF (nsrf.EQ.is_sic) THEN
                    tsurf(i,nsrf) = 271.35