Changeset 2534 for LMDZ5


Ignore:
Timestamp:
Jun 2, 2016, 4:24:31 PM (8 years ago)
Author:
musat
Message:

Add ratqsp0 and ratqsdp flags to control ratqs profile via
physiq.def for iflag_ratqs=4.
Default values are as before: ratqsp0=50000., ratqsdp=20000.
FH/IM

Location:
LMDZ5/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/calcratqs.F90

    r2236 r2534  
    11SUBROUTINE calcratqs(klon,klev,prt_level,lunout,       &
    22           iflag_ratqs,iflag_con,iflag_cld_th,pdtphys, &
    3            ratqsbas,ratqshaut,tau_ratqs,fact_cldcon,   &
     3           ratqsbas,ratqshaut,ratqsp0,ratqsdp, &
     4           tau_ratqs,fact_cldcon,   &
    45           ptconv,ptconvth,clwcon0th, rnebcon0th,      &
    56           paprs,pplay,q_seri,zqsat,fm_therm,          &
    … …  
    2122integer,intent(in) :: iflag_con,iflag_cld_th,iflag_ratqs
    2223real,intent(in) :: pdtphys,ratqsbas,ratqshaut,fact_cldcon,tau_ratqs
     24real,intent(in) :: ratqsp0, ratqsdp
    2325real, dimension(klon,klev+1),intent(in) :: paprs
    2426real, dimension(klon,klev),intent(in) :: pplay,q_seri,zqsat,fm_therm
    … …  
    125127         do k=1,klev
    126128           ratqss(:,k)=ratqsbas+0.5*(ratqshaut-ratqsbas) &
    127            *( tanh( (50000.-pplay(:,k))/20000.) + 1.)
     129!          *( tanh( (50000.-pplay(:,k))/20000.) + 1.)
     130           *( tanh( (ratqsp0-pplay(:,k))/ratqsdp) + 1.)
    128131         enddo
    129132

  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2530 r2534  
    679679    save ratqsbas,ratqshaut,tau_ratqs
    680680    !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
     681    REAL, SAVE :: ratqsp0=50000., ratqsdp=20000.
     682    !$OMP THREADPRIVATE(ratqsp0, ratqsdp)
    681683
    682684    ! Parametres lies au nouveau schema de nuages (SB, PDF)
    … …  
    10511053       CALL getin_p('random_notrig_max',random_notrig_max)
    10521054       CALL getin_p('ok_adjwk',ok_adjwk)
     1055       CALL getin_p('ratqsp0',ratqsp0)
     1056       CALL getin_p('ratqsdp',ratqsdp)
    10531057    ENDIF
    10541058
    … …  
    28512855    CALL  calcratqs(klon,klev,prt_level,lunout,        &
    28522856         iflag_ratqs,iflag_con,iflag_cld_th,pdtphys,  &
    2853          ratqsbas,ratqshaut,tau_ratqs,fact_cldcon,   &
     2857         ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
     2858         tau_ratqs,fact_cldcon,   &
    28542859         ptconv,ptconvth,clwcon0th, rnebcon0th,     &
    28552860         paprs,pplay,q_seri,zqsat,fm_therm, &
Note: See TracChangeset for help on using the changeset viewer.