Changeset 539

Timestamp:

Jun 9, 2004, 3:44:33 PM (21 years ago)

Author:

lmdzadmin

Message:

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

File:

• LMDZ4/trunk/libf/phylmd/stdlevvar.F90 (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/stdlevvar.F90

@@ -101 +101 @@
         ustar(i) = sqrt(cdram(i) * speed(i) * speed(i))
         zdte(i) = tpot(i) - ts1(i)
+!       PRINT*,'AVANT i,zdte',i,zdte(i)
 !IM cf FH : on prend le max : pour eviter le plantage sur SUN
-        zdte(i) = max(zdte(i),1.e-10)
+!IM BUG BUG BUG       zdte(i) = max(zdte(i),1.e-10)
+        zdte(i) = sign(max(abs(zdte(i)),1.e-10),zdte(i))
+!       PRINT*,'APRES i,zdte',i,zdte(i)
+!
         zdq(i) = max(q1(i),0.0) - max(qsurf(i),0.0)
 !