Changeset 2962 for LMDZ5

Timestamp:

Jul 22, 2017, 6:24:00 PM (7 years ago)

Author:

fhourdin

Message:

Nouvelle option pour la reevaporation des pluies.

File:

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

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -519 +519 @@
          ! Formulation en racine du flux de precip
          ! (Klemp & Wilhelmson, 1978; Sundqvist, 1988)
+         IF (iflag_evap_prec==3) THEN
+         zqevt = znebprecip(i)*coef_eva*(1.0-zq(i)/qsl) &
+              *SQRT(zrfl(i)/max(1.e-4,znebprecip(i))) &
+              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
+         ELSE
          zqevt = 1.*coef_eva*(1.0-zq(i)/qsl)*SQRT(zrfl(i)) &
               *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
+         ENDIF
+
+
          zqevt = MAX(0.0,MIN(zqevt,zrfl(i))) &
               *RG*dtime/(paprs(i,k)-paprs(i,k+1)) 
@@ -532 +540 @@
          ! Calcul de la sublimation du flux de precip solide herite
          !   d'au-dessus
+         IF (iflag_evap_prec==3) THEN
+         zqevti = znebprecip(i)*coef_eva*(1.0-zq(i)/qsi) &
+              *SQRT(zifl(i)/max(1.e-4,znebprecip(i))) &
+              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
+         ELSE
          zqevti = 1.*coef_eva*(1.0-zq(i)/qsi)*SQRT(zifl(i)) &
               *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
+         ENDIF
          zqevti = MAX(0.0,MIN(zqevti,zifl(i))) &
               *RG*dtime/(paprs(i,k)-paprs(i,k+1))