Changeset 2534

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:

• calcratqs.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/calcratqs.F90

@@ -1 +1 @@
 SUBROUTINE calcratqs(klon,klev,prt_level,lunout,       &
            iflag_ratqs,iflag_con,iflag_cld_th,pdtphys, &
-           ratqsbas,ratqshaut,tau_ratqs,fact_cldcon,   &
+           ratqsbas,ratqshaut,ratqsp0,ratqsdp, &
+           tau_ratqs,fact_cldcon,   &
            ptconv,ptconvth,clwcon0th, rnebcon0th,      &
            paprs,pplay,q_seri,zqsat,fm_therm,          &
@@ -21 +22 @@
 integer,intent(in) :: iflag_con,iflag_cld_th,iflag_ratqs
 real,intent(in) :: pdtphys,ratqsbas,ratqshaut,fact_cldcon,tau_ratqs
+real,intent(in) :: ratqsp0, ratqsdp
 real, dimension(klon,klev+1),intent(in) :: paprs
 real, dimension(klon,klev),intent(in) :: pplay,q_seri,zqsat,fm_therm
@@ -125 +127 @@
          do k=1,klev
            ratqss(:,k)=ratqsbas+0.5*(ratqshaut-ratqsbas) &
-           *( tanh( (50000.-pplay(:,k))/20000.) + 1.)
+!          *( tanh( (50000.-pplay(:,k))/20000.) + 1.)
+           *( tanh( (ratqsp0-pplay(:,k))/ratqsdp) + 1.)
          enddo
 

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -679 +679 @@
     save ratqsbas,ratqshaut,tau_ratqs
     !$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
+    REAL, SAVE :: ratqsp0=50000., ratqsdp=20000.
+    !$OMP THREADPRIVATE(ratqsp0, ratqsdp)
 
     ! Parametres lies au nouveau schema de nuages (SB, PDF)
@@ -1051 +1053 @@
        CALL getin_p('random_notrig_max',random_notrig_max)
        CALL getin_p('ok_adjwk',ok_adjwk)
+       CALL getin_p('ratqsp0',ratqsp0)
+       CALL getin_p('ratqsdp',ratqsdp)
     ENDIF
 
@@ -2851 +2855 @@
     CALL  calcratqs(klon,klev,prt_level,lunout,        &
          iflag_ratqs,iflag_con,iflag_cld_th,pdtphys,  &
-         ratqsbas,ratqshaut,tau_ratqs,fact_cldcon,   &
+         ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
+         tau_ratqs,fact_cldcon,   &
          ptconv,ptconvth,clwcon0th, rnebcon0th,     &
          paprs,pplay,q_seri,zqsat,fm_therm, &