Changeset 2923 for LMDZ5

Timestamp:

Jun 30, 2017, 11:03:09 AM (7 years ago)

Author:

fhourdin

Message:

Correction sous drapeaux d'une errance sur la reevaporation de la pluie.
flag_evap_prec = 0 : pas de reevaporation

1 : comme avant
2 : nouveau ou on utilise une surface max des nuages au

dessus comme fraction de la maille dans laquelle la
pluie peut se reevaporer (alors que c'etait la fraction
de la maille juste au dessus avant).

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (10 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -103 +103 @@
   ! Options du programme:
   !
-  REAL seuil_neb ! un nuage existe vraiment au-dela
-  PARAMETER (seuil_neb=0.001)
+  REAL, SAVE :: seuil_neb=0.001 ! un nuage existe vraiment au-dela
+  !$OMP THREADPRIVATE(seuil_neb)
+
 
   INTEGER ninter ! sous-intervals pour la precipitation
   PARAMETER (ninter=5)
-  LOGICAL evap_prec ! evaporation de la pluie
-  PARAMETER (evap_prec=.TRUE.)
+  INTEGER,SAVE :: iflag_evap_prec=1 ! evaporation de la pluie
+  !$OMP THREADPRIVATE(iflag_evap_prec)
   !
   LOGICAL cpartiel ! condensation partielle
@@ -124 +125 @@
   !
   INTEGER i, k, n, kk
+  INTEGER,save::itap=0
+  !$OMP THREADPRIVATE(itap)
+
   REAL qsl, qsi
   real zct      ,zcl
@@ -154 +158 @@
   REAL t_glace_min_old
   REAL zdz(klon),zrho(klon),ztot      , zrhol(klon)
-  REAL zchau      ,zfroi      ,zfice(klon),zneb(klon)
+  REAL zchau      ,zfroi      ,zfice(klon),zneb(klon),znebprecip(klon)
   REAL zmelt, zpluie, zice
   REAL dzfice(klon)
@@ -212 +216 @@
 !CR: pour iflag_ice_thermo=2, on active que la convection
 !  ice_thermo = iflag_ice_thermo .GE. 1
+
+itap=itap+1
+znebprecip(:)=0.
+
    ice_thermo = (iflag_ice_thermo .EQ. 1).OR.(iflag_ice_thermo .GE. 3)
   zdelq=0.0
@@ -218 +226 @@
   IF (appel1er) THEN
      CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp)
+     CALL getin_p('iflag_evap_prec',iflag_evap_prec)
+     CALL getin_p('seuil_neb',seuil_neb)
      write(lunout,*)' iflag_oldbug_fisrtilp =',iflag_oldbug_fisrtilp 
      !
      WRITE(lunout,*) 'fisrtilp, ninter:', ninter
-     WRITE(lunout,*) 'fisrtilp, evap_prec:', evap_prec
+     WRITE(lunout,*) 'fisrtilp, iflag_evap_prec:', iflag_evap_prec
      WRITE(lunout,*) 'fisrtilp, cpartiel:', cpartiel
+     
      IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN
         WRITE(lunout,*) 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
@@ -387 +398 @@
      !   - zmqc: masse de precip qui doit etre thermalisee
      !
-     IF (evap_prec) THEN
+     IF (iflag_evap_prec>=1) THEN
         DO i = 1, klon
 !          S'il y a des precipitations
@@ -480 +491 @@
 !        S'il y a des precipitations
          IF (zrfl(i)+zifl(i).GT.0.) THEN
+
+         IF (iflag_evap_prec==1) THEN
+            znebprecip(i)=zneb(i)
+         ELSE
+            znebprecip(i)=MAX(zneb(i),znebprecip(i))
+         ENDIF
      
         ! Evap max pour ne pas saturer la fraction sous le nuage
-         zqev0 = MAX (0.0, (zqs(i)-zq(i))*zneb(i) )
+         zqev0 = MAX (0.0, (zqs(i)-zq(i))*znebprecip(i) )
 
          !JAM
@@ -567 +584 @@
          zrfl(i) = zrfln(i)
          zifl(i) = zifln(i)
+!        print*,'REEVAP ',itap,k,znebprecip(1),zqev0,zqev,zqevi,zrfl(1)
 
 !CR ATTENTION: deplacement de la fonte de la glace
@@ -595 +613 @@
            end if
 
+           ELSE
+              ! Si on n'a plus de pluies, on reinitialise a 0 la farcion
+              ! sous nuageuse utilisee pour la pluie.
+              znebprecip(i)=0.
            ENDIF ! (zrfl(i)+zifl(i).GT.0.)
         ENDDO
@@ -603 +625 @@
      ! Fin evaporation de la precipitation
      ! ----------------------------------------------------------------
-     ENDIF ! (evap_prec)
+     ENDIF ! (iflag_evap_prec>=1)
      !
      ! Calculer Qs et L/Cp*dQs/dT: