Changeset 2778 for LMDZ5


Ignore:
Timestamp:
Jan 26, 2017, 3:49:44 PM (8 years ago)
Author:
musat
Message:

add max to be consistent with the rest

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/cdrag.F90

    r2772 r2778  
    172172! the negative qsurf and q1 are set to zero (as in coefcdrag)
    173173     ztsolv = tsurf(i) * (1.0+RETV*max(qsurf(i),0.0)) ! negative qsurf set to zero
    174      ztvd = (t1(i)+zgeop1(i)/RCPD/(1.+RVTMP2*q1(i))) &
     174     ztvd = (t1(i)+zgeop1(i)/RCPD/(1.+RVTMP2*max(q1(i),0.0))) &
    175175          *(1.+RETV*max(q1(i),0.0)) ! negative q1 set to zero
    176176     zri(i) = zgeop1(i)*(ztvd-ztsolv)/(zdu2*ztvd)
Note: See TracChangeset for help on using the changeset viewer.