Changeset 23 for LMDZ.3.3


Ignore:
Timestamp:
Jan 18, 2000, 5:37:55 PM (25 years ago)
Author:
lmdz
Message:

modif de physiq et fistrlp.F pour recuperer les flux d'eau precipitante (cf Olivier)

  1. Hauglustaine
Location:
LMDZ.3.3/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/phylmd/fisrtilp.F

    r2 r23  
    11      SUBROUTINE fisrtilp(dtime,paprs,pplay,t,q,
    2      s                   d_t, d_q, d_ql, rneb, radliq, rain, snow)
     2     s                   d_t, d_q, d_ql, rneb, radliq, rain, snow,
     3     s                   prfl, psfl)
    34c
    45      IMPLICIT none
    … …  
    2829      REAL rain(klon) ! pluies (mm/s)
    2930      REAL snow(klon) ! neige (mm/s)
     31      REAL prfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s)
     32      REAL psfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s)
    3033c
    3134c Options du programme:
    … …  
    103106c Initialiser les sorties:
    104107c
     108      DO k = 1, klev+1
     109      DO i = 1, klon
     110         prfl(i,k) = 0.0
     111         psfl(i,k) = 0.0
     112      ENDDO
     113      ENDDO
    105114      DO k = 1, klev
    106115      DO i = 1, klon
    … …  
    268277         zrfl(i) = zrfl(i)+ MAX(zcond(i)-zoliq(i),0.0)
    269278     .                    * (paprs(i,k)-paprs(i,k+1))/(RG*dtime)
    270       ENDIF
    271       ENDDO
     279         IF (zt(i).LT.RTT) THEN
     280           psfl(i,k)=zrfl(i)
     281         ELSE
     282           prfl(i,k)=zrfl(i)
     283         ENDIF
     284      ENDIF
     285      ENDDO
     286
    272287c
    273288c Calculer les tendances de q et de t:

  • LMDZ.3.3/trunk/libf/phylmd/physiq.F

    r15 r23  
    276276      REAL yv1(klon)            ! vents dans la premiere couche V
    277277      LOGICAL offline           ! Controle du stockage ds "physique"
    278       PARAMETER (offline=.TRUE.)
     278      PARAMETER (offline=.FALSE.)
    279279      REAL pfrac_impa(klon,klev)! Produits des coefs lessivage impaction
    280280      save pfrac_impa
    … …  
    432432      INTEGER kcbot(klon), kctop(klon), kdtop(klon)
    433433      REAL pmflxr(klon,klev+1), pmflxs(klon,klev+1)
     434      REAL prfl(klon,klev+1), psfl(klon,klev+1)
    434435c
    435436      INTEGER ibas_con(klon), itop_con(klon)
    … …  
    15531554     .           rain_lsc, snow_lsc,
    15541555     .           pfrac_impa, pfrac_nucl, pfrac_1nucl,
    1555      .           frac_impa, frac_nucl )
     1556     .           frac_impa, frac_nucl,
     1557     .           prfl, psfl)
    15561558      DO k = 1, klev
    15571559      DO i = 1, klon
Note: See TracChangeset for help on using the changeset viewer.