Changeset 165 for LMDZ.3.3/trunk


Ignore:
Timestamp:
Nov 17, 2000, 2:34:28 PM (24 years ago)
Author:
lmdz
Message:

Petite correction "physique" dans la correction precedente. O.Boucher & J.Quaas
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/phylmd/conflx.F

    r164 r165  
    10431043cjq The variable maxpdmfdp(klon) has been introduced by Olivier Boucher
    10441044cjq 14/11/00 to fix the problem with the negative precipitation.     
    1045       REAL pdmfdp(klon,klev), maxpdmfdp(klon)
     1045      REAL pdmfdp(klon,klev), maxpdmfdp(klon,klev)
    10461046      REAL prfl(klon), psfl(klon)
    10471047      REAL pmflxr(klon,klev+1), pmflxs(klon,klev+1)
    10481048      INTEGER  kcbot(klon), kctop(klon), ktype(klon)
    10491049      LOGICAL  ldland(klon), ldcum(klon)
    1050       INTEGER k, i
     1050      INTEGER k, kp, i
    10511051      REAL zcons1, zcons2, zcucov, ztmelp2
    10521052      REAL pdtime, zdp, zzp, zfac, zsnmlt, zrfl, zrnew
     
    11731173c        panache descendant pour eliminer la negativite
    11741174         IF ((pmflxr(i,k+1)+pmflxs(i,k+1)).LT.0.0) THEN
    1175 cjq            write (*,*) 'total precip negative. rain ',pmflxr(i,k+1),
    1176 cjq     .                  'snow ',pmflxs(i,k+1)
    11771175            pdmfdp(i,k) = -pmflxr(i,k)-pmflxs(i,k)-pdmfup(i,k)
    11781176            pmflxr(i,k+1) = 0.0
     
    11921190cjq as well as the evaporation.
    11931191cjq     
     1192      DO k = 1, klev
    11941193      DO i = 1, klon
    1195          maxpdmfdp(i)=0.0
    1196       ENDDO
    1197       DO k = ktopm2, klev
    1198          DO i = 1, klon
    1199          if (k.GE.kcbot(i)) then
    1200          maxpdmfdp(i)=maxpdmfdp(i)+pdmfdp(i,k)
    1201          endif
    1202          ENDDO
     1194         maxpdmfdp(i,k)=0.0
     1195      ENDDO
     1196      ENDDO
     1197      DO k = 1, klev
     1198       DO kp = k, klev
     1199        DO i = 1, klon
     1200         maxpdmfdp(i,k)=maxpdmfdp(i,k)+pdmfdp(i,kp)
     1201        ENDDO
     1202       ENDDO
    12031203      ENDDO
    12041204cjq End of initialization
     
    12211221cjq be evaporated (surely the evaporation can't be positive):           
    12221222            zdrfl=MAX(zdrfl,
    1223      .            MIN(-pmflxr(i,k)-pmflxs(i,k)-maxpdmfdp(i),0.0))
     1223     .            MIN(-pmflxr(i,k)-pmflxs(i,k)-maxpdmfdp(i,k),0.0))
    12241224cjq End of insertion
    12251225c           
     
    12421242            pdmfdp(i,k) = 0.0
    12431243            pdpmel(i,k) = 0.0
    1244          ENDIF
     1244         ENDIF         
     1245         if (pmflxr(i,k) + pmflxs(i,k).lt.-1.e-26)
     1246     .    write(*,*) 'precip. < 1e-16 ',pmflxr(i,k) + pmflxs(i,k)
    12451247      ENDIF
    12461248      ENDDO
Note: See TracChangeset for help on using the changeset viewer.