Changeset 2236 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Mar 17, 2015, 12:04:12 PM (9 years ago)

Author:

fhourdin

Message:

Changement du drapeau iflag_cldth en iflag_cld_th
et du défaut de iflag_cld_cv pour le choix entre schema
de nuages convectifs

lognormal = 0
bigaussien = 1

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (11 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -639 +639 @@
   !$OMP THREADPRIVATE(fact_cldcon,facttemps)
 
-  integer iflag_cldth
-  save iflag_cldth
-  !$OMP THREADPRIVATE(iflag_cldth)
+  integer iflag_cld_th
+  save iflag_cld_th
+  !$OMP THREADPRIVATE(iflag_cld_th)
   logical ptconv(klon,klev)
   !IM cf. AM 081204 BEG
@@ -919 +919 @@
           solarlong0,seuil_inversion, &
           fact_cldcon, facttemps,ok_newmicro,iflag_radia, &
-          iflag_cldth,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
+          iflag_cld_th,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
           ok_ade, ok_aie, ok_cdnc, aerosol_couple,  &
           flag_aerosol, flag_aerosol_strat, new_aod, &
@@ -1027 +1027 @@
      print*,'CYCLE_DIURNE', cycle_diurne
      !
-     IF (iflag_con.EQ.2.AND.iflag_cldth.GT.-1) THEN
-        abort_message = 'Tiedtke needs iflag_cldth=-2 or -1'
+     IF (iflag_con.EQ.2.AND.iflag_cld_th.GT.-1) THEN
+        abort_message = 'Tiedtke needs iflag_cld_th=-2 or -1'
         CALL abort_gcm (modname,abort_message,1)
      ENDIF
@@ -1143 +1143 @@
                 ,alp_bl_prescr, ale_bl_prescr)
            ! 11/09/06 rajout initialisation ALE et ALP du wake et PBL(YU)
-           !        print*,'apres ini_wake iflag_cldth=', iflag_cldth
+           !        print*,'apres ini_wake iflag_cld_th=', iflag_cld_th
         endif
 
@@ -2298 +2298 @@
      !   calcul des proprietes des nuages convectifs
      clwcon0(:,:)=fact_cldcon*clwcon0(:,:)
-     IF (iflag_cld_cv <= 1) THEN
+     IF (iflag_cld_cv == 0) THEN
      call clouds_gno &
           (klon,klev,q_seri,zqsat,clwcon0,ptconv,ratqsc,rnebcon0)
@@ -2539 +2539 @@
   END IF
 
-  !      print*,'apres callwake iflag_cldth=', iflag_cldth
+  !      print*,'apres callwake iflag_cld_th=', iflag_cld_th
   !
   !===================================================================
@@ -2860 +2860 @@
   ! water distribution
   CALL  calcratqs(klon,klev,prt_level,lunout,        &
-       iflag_ratqs,iflag_con,iflag_cldth,pdtphys,  &
+       iflag_ratqs,iflag_con,iflag_cld_th,pdtphys,  &
        ratqsbas,ratqshaut,tau_ratqs,fact_cldcon,   &
        ptconv,ptconvth,clwcon0th, rnebcon0th,     &
@@ -2882 +2882 @@
        frac_impa, frac_nucl, beta_prec_fisrt, &
        prfl, psfl, rhcl,  &
-       zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cldth, &
+       zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, &
        iflag_ice_thermo)
   !
@@ -2938 +2938 @@
   !
   !IM cf FH
-  !     IF (iflag_cldth.eq.-1) THEN ! seulement pour Tiedtke
-  IF (iflag_cldth.le.-1) THEN ! seulement pour Tiedtke
+  !     IF (iflag_cld_th.eq.-1) THEN ! seulement pour Tiedtke
+  IF (iflag_cld_th.le.-1) THEN ! seulement pour Tiedtke
      snow_tiedtke=0.
      !     print*,'avant calcul de la pseudo precip '
-     !     print*,'iflag_cldth',iflag_cldth
-     if (iflag_cldth.eq.-1) then
+     !     print*,'iflag_cld_th',iflag_cld_th
+     if (iflag_cld_th.eq.-1) then
         rain_tiedtke=rain_con
      else
@@ -2976 +2976 @@
      ENDDO
 
-  ELSE IF (iflag_cldth.ge.3) THEN
+  ELSE IF (iflag_cld_th.ge.3) THEN
      !  On prend pour les nuages convectifs le max du calcul de la
      !  convection et du calcul du pas de temps precedent diminue d'un facteur
@@ -3066 +3066 @@
      !   On prend la somme des fractions nuageuses et des contenus en eau
 
-     if (iflag_cldth>=5) then
+     if (iflag_cld_th>=5) then
 
         do k=1,klev