Ignore:
Timestamp:
May 5, 2021, 12:50:37 PM (3 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:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/dyn1d/compar1d.h

    r3594 r3888  
    88      real :: nat_surf
    99      real :: tsurf
     10      real :: beta_surf
    1011      real :: rugos
    1112      real :: rugosh
     
    4546      real    :: p_nudging_u, p_nudging_v, p_nudging_w, p_nudging_t, p_nudging_qv
    4647      common/com_par1d/                                                 &
    47      & nat_surf,tsurf,rugos,rugosh,                                     &
     48     & nat_surf,tsurf,beta_surf,rugos,rugosh,                           &
    4849     & xqsol,qsurf,psurf,zsurf,albedo,time,time_ini,xlat,xlon,airefi,   &
    4950     & wtsurf,wqsurf,restart_runoff,xagesno,qsolinp,zpicinp,            &
  • LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

    r3781 r3888  
    1111       du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
    1212       falb_dir, falb_dif, &
    13        ftsol, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
     13       ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
    1414       rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, &
    1515       solsw, t_ancien, q_ancien, u_ancien, v_ancien, wake_cstar, &
     
    656656      qsol = qsolinp
    657657      qsurf = fq_sat(tsurf,psurf/100.)
     658      beta_aridity(:,:) = beta_surf
    658659      day1= day_ini
    659660      time=daytime-day
     
    795796
    796797        fder=0.
     798        snsrf(1,:)=snowmass ! masse de neige des sous surface
    797799        print *, 'snsrf', snsrf
    798         snsrf(1,:)=snowmass ! masse de neige des sous surface
    799800        qsurfsrf(1,:)=qsurf ! humidite de l'air des sous surface
    800801        fevap=0.
  • LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

    r3781 r3888  
    77       du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
    88       falb_dir, falb_dif, &
    9        ftsol, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
     9       ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
    1010       rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, &
    1111       solsw, t_ancien, q_ancien, u_ancien, v_ancien, wake_cstar, &
     
    429429      qsol = qsolinp
    430430      qsurf = fq_sat(tsurf,psurf/100.)
     431      beta_aridity(:,:) = beta_surf
    431432      day1= day_ini
    432433      time=daytime-day
Note: See TracChangeset for help on using the changeset viewer.