Ignore:
Timestamp:
Jul 22, 2017, 6:24:00 PM (7 years ago)
Author:
fhourdin
Message:

Nouvelle option pour la reevaporation des pluies.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/fisrtilp.F90

    r2956 r2962  
    519519         ! Formulation en racine du flux de precip
    520520         ! (Klemp & Wilhelmson, 1978; Sundqvist, 1988)
     521         IF (iflag_evap_prec==3) THEN
     522         zqevt = znebprecip(i)*coef_eva*(1.0-zq(i)/qsl) &
     523              *SQRT(zrfl(i)/max(1.e-4,znebprecip(i))) &
     524              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
     525         ELSE
    521526         zqevt = 1.*coef_eva*(1.0-zq(i)/qsl)*SQRT(zrfl(i)) &
    522527              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
     528         ENDIF
     529
     530
    523531         zqevt = MAX(0.0,MIN(zqevt,zrfl(i))) &
    524532              *RG*dtime/(paprs(i,k)-paprs(i,k+1))
     
    532540         ! Calcul de la sublimation du flux de precip solide herite
    533541         !   d'au-dessus
     542         IF (iflag_evap_prec==3) THEN
     543         zqevti = znebprecip(i)*coef_eva*(1.0-zq(i)/qsi) &
     544              *SQRT(zifl(i)/max(1.e-4,znebprecip(i))) &
     545              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
     546         ELSE
    534547         zqevti = 1.*coef_eva*(1.0-zq(i)/qsi)*SQRT(zifl(i)) &
    535548              *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
     549         ENDIF
    536550         zqevti = MAX(0.0,MIN(zqevti,zifl(i))) &
    537551              *RG*dtime/(paprs(i,k)-paprs(i,k+1))   
Note: See TracChangeset for help on using the changeset viewer.