Changeset 3720 for LMDZ6/branches

Timestamp:

Jun 13, 2020, 12:23:29 AM (4 years ago)

Author:

lguez

Message:

Output bulk sea surface salinity

Output bulk sea surface salinity when activate_ocean_skin >= 1. So
add variable sss in module phys_output_var_mod. In procedure
surf_ocean, promote sss from local variable to dummy argument.

Location:

LMDZ6/branches/Ocean_skin/libf/phylmd

Files:

• pbl_surface_mod.F90 (modified) (5 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_var_mod.F90 (modified) (2 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)
• surf_ocean_mod.F90 (modified) (3 diffs)

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

@@ -293 +293 @@
     USE ioipsl_getin_p_mod, ONLY : getin_p
     use phys_state_var_mod, only: s_int, ds_ns, dt_ns
-    use phys_output_var_mod, only: sst_nff, dter, dser, tkt, tks, taur
+    use phys_output_var_mod, only: sst_nff, dter, dser, tkt, tks, taur, sss
 #ifdef CPP_XIOS
     USE wxios, ONLY: missing_val
@@ -831 +831 @@
 
     real, DIMENSION(klon):: ysst_nff, ys_int, yds_ns, ydt_ns, ydter, ydser, &
-         ytkt, ytks, ytaur
+         ytkt, ytks, ytaur, ysss
     ! compression of sst_nff, s_int, ds_ns, dt_ns, dter, dser, tkt, tks,
-    ! taur on ocean points
+    ! taur, sss on ocean points
 
 !****************************************************************************************
@@ -1956 +1956 @@
                y_flux_u1, y_flux_v1, ysst_nff(:knon), ys_int(:knon), &
                yds_ns(:knon), ydt_ns(:knon), ydter(:knon), ydser(:knon), &
-               ytkt(:knon), ytks(:knon), ytaur(:knon))
+               ytkt(:knon), ytks(:knon), ytaur(:knon), ysss)
       IF (prt_level >=10) THEN
           print *,'arg de surf_ocean: ycdragh ',ycdragh
@@ -2570 +2570 @@
           tks = missing_val
           taur = missing_val
+          sss = missing_val
           
           s_int(ni(:knon)) = ys_int(:knon)
@@ -2580 +2581 @@
           tks(ni(:knon)) = ytks(:knon)
           taur(ni(:knon)) = ytaur(:knon)
+          sss(ni(:knon)) = ysss(:knon)
        end if
 

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

@@ -1983 +1983 @@
         "momentum flux due to rain", "Pa", '')
 
+   type(ctrl_out), save:: o_sss &
+        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'SSS', &
+        "bulk sea-surface salinity", "ppt", '')
+
 END MODULE phys_output_ctrlout_mod

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

@@ -156 +156 @@
   ! épaisseur (m) de la couche de diffusion de masse (microlayer)
   
-  REAL, SAVE, ALLOCATABLE:: taur(:)
-  ! momentum flux due to rain, in Pa
+  REAL, SAVE, ALLOCATABLE:: taur(:) ! momentum flux due to rain, in Pa
+
+  REAL, SAVE, ALLOCATABLE:: sss(:)
+  ! bulk salinity of the surface layer of the ocean, in ppt
   
-  !$OMP THREADPRIVATE(sst_nff, dter, dser, tkt, tks, taur)
+  !$OMP THREADPRIVATE(sst_nff, dter, dser, tkt, tks, taur, sss)
 
 CONTAINS
@@ -221 +223 @@
 
     if (activate_ocean_skin >= 1) allocate(sst_nff(klon), dter(klon), &
-         dser(klon), tkt(klon), tks(klon), taur(klon))
+         dser(klon), tkt(klon), tks(klon), taur(klon), sss(klon))
 
   END SUBROUTINE phys_output_var_init

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

@@ -203 +203 @@
          o_flx_co2_ocean, o_flx_co2_land, o_flx_co2_ff, o_flx_co2_bb, &
          o_sst_nff, o_s_int, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, o_tks, &
-         o_taur
+         o_taur, o_sss
 
 #ifdef CPP_StratAer
@@ -360 +360 @@
          alt_tropo, &
 !Ionela
-         ok_4xCO2atm, sst_nff, dter, dser, tkt, tks, taur
+         ok_4xCO2atm, sst_nff, dter, dser, tkt, tks, taur, sss
 
     USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
@@ -2313 +2313 @@
           CALL histwrite_phy(o_tks, tks)
           CALL histwrite_phy(o_taur, taur)
+          CALL histwrite_phy(o_sss, sss)
        end if
 

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

@@ -20 +20 @@
        z0m, z0h, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & 
        tsurf_new, dflux_s, dflux_l, lmt_bils, &
-       flux_u1, flux_v1, sst_nff, s_int, ds_ns, dt_ns, dter, dser, tkt, tks, taur)
+       flux_u1, flux_v1, sst_nff, s_int, ds_ns, dt_ns, dter, dser, tkt, tks, &
+       taur, sss)
 
     use albedo, only: alboc, alboc_cd
@@ -127 +128 @@
     ! momentum flux due to rain, in Pa
 
+    real, intent(out):: sss(:) ! (klon)
+    ! Bulk salinity of the surface layer of the ocean, in ppt. (Only
+    ! defined for subscripts 1:knon, but we have to declare it with
+    ! size klon because of the coupling machinery.)
+
     ! Local variables
     !******************************************************************************
@@ -138 +144 @@
     real rhoa(knon) ! density of moist air  (kg / m3)
     REAL sens_prec_liq(knon)
-
-    real sss(klon)
-    ! Bulk salinity of the surface layer of the ocean, in ppt. (Only
-    ! defined for subscripts 1:knon, but we have to declare it with
-    ! size klon because of the coupling machinery.)
 
     REAL t_int(knon) ! ocean-air interface temperature, in K