Ignore:
Timestamp:
Jun 8, 2004, 6:47:45 PM (20 years ago)
Author:
lmdzadmin
Message:

Petit bug sur le signe de la valeur seuil du gradient tair-tsol
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/stdlevvar.F90

    r509 r537  
    9898        ustar(i) = sqrt(cdram(i) * speed(i) * speed(i))
    9999        zdte(i) = tpot(i) - ts1(i)
     100!       PRINT*,'AVANT i,zdte',i,zdte(i)
    100101!IM cf FH : on prend le max : pour eviter le plantage sur SUN
    101         zdte(i) = max(zdte(i),1.e-10)
     102!IM BUG BUG BUG       zdte(i) = max(zdte(i),1.e-10)
     103        zdte(i) = sign(max(abs(zdte(i)),1.e-10),zdte(i))
     104!       PRINT*,'APRES i,zdte',i,zdte(i)
     105!
    102106        zdq(i) = max(q1(i),0.0) - max(qsurf(i),0.0)
    103107!
Note: See TracChangeset for help on using the changeset viewer.