Changeset 3888 for LMDZ6/trunk/libf/phylmd/dyn1d

Timestamp:

May 5, 2021, 12:50:37 PM (4 years ago)

Author:

jyg

Message:

New provisional version of the splitting of the
diffusive boundary layer into inwake and offwake
PBLs. The splitting of the diffuse BL should NOT
be activated yet for general purpose simulations.

The splitting is activated by:
mod(iflag_pbl_split,10)=1 for the option with
fixed surface temperature and
mod(iflag_pbl_split,10)=2 for the option with
coupled surface temperature.

iflag_pbl_split=0 ==> no splittingat all.
iflag_pbl_split=10 ==> splitting of thermals.
iflag_pbl_split=11 ==> splitting of thermals and
of vertical diffusion (fixed surf. temp.).
iflag_pbl_split=12 ==> splitting of thermals and
of vertical diffusion (coupled surf. temp.).

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• compar1d.h (modified) (2 diffs)
• old_lmdz1d.F90 (modified) (3 diffs)
• scm.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/compar1d.h

@@ -8 +8 @@
       real :: nat_surf
       real :: tsurf
+      real :: beta_surf
       real :: rugos
       real :: rugosh
@@ -45 +46 @@
       real    :: p_nudging_u, p_nudging_v, p_nudging_w, p_nudging_t, p_nudging_qv
       common/com_par1d/                                                 &
-     & nat_surf,tsurf,rugos,rugosh,                                     &
+     & nat_surf,tsurf,beta_surf,rugos,rugosh,                           &
      & xqsol,qsurf,psurf,zsurf,albedo,time,time_ini,xlat,xlon,airefi,   &
      & wtsurf,wqsurf,restart_runoff,xagesno,qsolinp,zpicinp,            &

LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

@@ -11 +11 @@
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
        falb_dir, falb_dif, &
-       ftsol, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
+       ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
        rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, &
        solsw, t_ancien, q_ancien, u_ancien, v_ancien, wake_cstar, &
@@ -656 +656 @@
       qsol = qsolinp
       qsurf = fq_sat(tsurf,psurf/100.)
+      beta_aridity(:,:) = beta_surf
       day1= day_ini
       time=daytime-day
@@ -795 +796 @@
 
         fder=0.
+        snsrf(1,:)=snowmass ! masse de neige des sous surface
         print *, 'snsrf', snsrf
-        snsrf(1,:)=snowmass ! masse de neige des sous surface
         qsurfsrf(1,:)=qsurf ! humidite de l'air des sous surface
         fevap=0.

LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

@@ -7 +7 @@
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
        falb_dir, falb_dif, &
-       ftsol, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
+       ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
        rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, &
        solsw, t_ancien, q_ancien, u_ancien, v_ancien, wake_cstar, &
@@ -429 +429 @@
       qsol = qsolinp
       qsurf = fq_sat(tsurf,psurf/100.)
+      beta_aridity(:,:) = beta_surf
       day1= day_ini
       time=daytime-day