Changeset 1041


Ignore:
Timestamp:
Oct 30, 2008, 4:37:55 PM (16 years ago)
Author:
lmdzadmin
Message:

Correction oubli(bug) calcul eau liquide apres melange menant a des humidites specifiques negatives
JYG/NR/IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/cv3p_mixing.F

    r991 r1041  
    417417c
    418418c--     Compute liquid water static energy of mixed draughts
     419c       IF (j .GT. i) THEN
     420c        awat=elij(il,i,j)-(1.-ep(il,j))*clw(il,j)
     421c        awat=amax1(awat,0.0)
     422c       ELSE
     423c        awat = 0.
     424c       ENDIF
     425c       Hent(il,i,j) = (1.-Sigij(il,i,j))*HP(il,i)
     426c    :         + Sigij(il,i,j)*H(il,i)
     427c    :         + (LV(il,j)+(cpd-cpv)*t(il,j))*awat
     428cIM 301008 beg
     429        Hent(il,i,j) = (1.-Sigij(il,i,j))*HP(il,i)
     430     :         + Sigij(il,i,j)*H(il,i)
     431
     432        Elij(il,i,j) = Qent(il,i,j)-rs(il,j)
     433        Elij(il,i,j) = Elij(il,i,j)
     434     :    + ((h(il,j)-Hent(il,i,j))*rs(il,j)*LV(il,j)
     435     :    / ((cpd*(1.-Qent(il,i,j))+Qent(il,i,j)*cpv)
     436     :    * rrv*t(il,j)*t(il,j)))
     437        Elij(il,i,j) = Elij(il,i,j)
     438     :    / (1.+LV(il,j)*LV(il,j)*rs(il,j)
     439     :    / ((cpd*(1.-Qent(il,i,j))+Qent(il,i,j)*cpv)
     440     :    * rrv*t(il,j)*t(il,j)))
     441
     442        Elij(il,i,j) = max(elij(il,i,j),0.)
     443
     444        Elij(il,i,j) = min(elij(il,i,j),Qent(il,i,j))
     445
    419446        IF (j .GT. i) THEN
    420447         awat=elij(il,i,j)-(1.-ep(il,j))*clw(il,j)
     
    423450         awat = 0.
    424451        ENDIF
    425         Hent(il,i,j) = (1.-Sigij(il,i,j))*HP(il,i)
    426      :         + Sigij(il,i,j)*H(il,i)
    427      :         + (LV(il,j)+(cpd-cpv)*t(il,j))*awat
     452
     453c        print *,h(il,j)-hent(il,i,j),LV(il,j)*rs(il,j)/(cpd*rrv*t(il,j)*
     454c    :         t(il,j))
     455
     456        Hent(il,i,j) =  Hent(il,i,j)+(LV(il,j)+(cpd-cpv)*t(il,j))
     457     :         * awat
     458cIM 301008 end
    428459c
    429460c      print *,'mix : i,j,hent(il,i,j),sigij(il,i,j) ',
Note: See TracChangeset for help on using the changeset viewer.