Changeset 509 for LMDZ.3.3

Timestamp:

Mar 24, 2004, 10:21:49 AM (21 years ago)

Author:

lmdzadmin

Message:

On limite le dt a 1e-10 pour eviter le plantage sur SUN
FH/IM

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/stdlevvar.F90 (modified) (8 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/stdlevvar.F90

@@ -1 +1 @@
       SUBROUTINE stdlevvar(klon, knon, nsrf, zxli, &
-                           u1, v1, t1, q1, z1, &
-                           ts1, qsurf, rugos, psol, pat1, &
-                           t_2m, q_2m, u_10m) 
+ &                         u1, v1, t1, q1, z1, &
+ &                         ts1, qsurf, rugos, psol, pat1, &
+ &                         t_2m, q_2m, u_10m) 
       IMPLICIT NONE
 !-------------------------------------------------------------------------
@@ -87 +87 @@
       okri=.FALSE.
       CALL coefcdrag(klon, knon, nsrf, zxli, &
-                     speed, t1, q1, z1, psol, &
-                     ts1, qsurf, rugos, okri, ri1,  &         
-                     cdram, cdrah, cdran, zri1, pref)            
+ &                   speed, t1, q1, z1, psol, &
+ &                   ts1, qsurf, rugos, okri, ri1,  &         
+ &                   cdram, cdrah, cdran, zri1, pref)            
 !
 !---------Star variables----------------------------------------------------
@@ -98 +98 @@
         ustar(i) = sqrt(cdram(i) * speed(i) * speed(i))
         zdte(i) = tpot(i) - ts1(i)
+!IM cf FH : on prend le max : pour eviter le plantage sur SUN
+        zdte(i) = max(zdte(i),1.e-10)
         zdq(i) = max(q1(i),0.0) - max(qsurf(i),0.0)
 !
@@ -103 +105 @@
         qstar(i) = (cdrah(i) * zdq(i) * speed(i))/ustar(i)
         lmon(i) = (ustar(i) * ustar(i) * tpot(i))/ &
-                  (RKAR * RG * testar(i))
+ &                (RKAR * RG * testar(i))
       ENDDO
 !
@@ -109 +111 @@
       zref = 2.0
       CALL screenp(klon, knon, nsrf, speed, tpot, q1, &
-                   ts1, qsurf, rugos, lmon, &
-                   ustar, testar, qstar, zref, &
-                   delu, delte, delq)
+ &                 ts1, qsurf, rugos, lmon, &
+ &                 ustar, testar, qstar, zref, &
+ &                 delu, delte, delq)
 !
       DO i = 1, knon
@@ -129 +131 @@
         okri=.TRUE.
         CALL screenc(klon, knon, nsrf, zxli, &
-                     u_zref, temp, q_zref, zref, &
-                     ts1, qsurf, rugos, psol, &           
-                     ustar, testar, qstar, okri, ri1, &
-                     pref, delu, delte, delq) 
+ &                   u_zref, temp, q_zref, zref, &
+ &                   ts1, qsurf, rugos, psol, &           
+ &                   ustar, testar, qstar, okri, ri1, &
+ &                   pref, delu, delte, delq) 
 !
         DO i = 1, knon
@@ -211 +213 @@
       zref = 10.0
       CALL screenp(klon, knon, nsrf, speed, tpot, q1, &
-                   ts1, qsurf, rugos, lmon, &
-                   ustar, testar, qstar, zref, &
-                   delu, delte, delq)
+ &                 ts1, qsurf, rugos, lmon, &
+ &                 ustar, testar, qstar, zref, &
+ &                 delu, delte, delq)
 !
       DO i = 1, knon
@@ -231 +233 @@
         okri=.TRUE.
         CALL screenc(klon, knon, nsrf, zxli, &
-                     u_zref, temp, q_zref, zref, &
-                     ts1, qsurf, rugos, psol, &
-                     ustar, testar, qstar, okri, ri1, &
-                     pref, delu, delte, delq)
+ &                   u_zref, temp, q_zref, zref, &
+ &                   ts1, qsurf, rugos, psol, &
+ &                   ustar, testar, qstar, okri, ri1, &
+ &                   pref, delu, delte, delq)
 !
         DO i = 1, knon