Changeset 509


Ignore:
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:
1 edited

Legend:

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

    r472 r509  
    11      SUBROUTINE stdlevvar(klon, knon, nsrf, zxli, &
    2                            u1, v1, t1, q1, z1, &
    3                            ts1, qsurf, rugos, psol, pat1, &
    4                            t_2m, q_2m, u_10m)
     2 &                         u1, v1, t1, q1, z1, &
     3 &                         ts1, qsurf, rugos, psol, pat1, &
     4 &                         t_2m, q_2m, u_10m)
    55      IMPLICIT NONE
    66!-------------------------------------------------------------------------
     
    8787      okri=.FALSE.
    8888      CALL coefcdrag(klon, knon, nsrf, zxli, &
    89                      speed, t1, q1, z1, psol, &
    90                      ts1, qsurf, rugos, okri, ri1,  &         
    91                      cdram, cdrah, cdran, zri1, pref)           
     89 &                   speed, t1, q1, z1, psol, &
     90 &                   ts1, qsurf, rugos, okri, ri1,  &         
     91 &                   cdram, cdrah, cdran, zri1, pref)           
    9292!
    9393!---------Star variables----------------------------------------------------
     
    9898        ustar(i) = sqrt(cdram(i) * speed(i) * speed(i))
    9999        zdte(i) = tpot(i) - ts1(i)
     100!IM cf FH : on prend le max : pour eviter le plantage sur SUN
     101        zdte(i) = max(zdte(i),1.e-10)
    100102        zdq(i) = max(q1(i),0.0) - max(qsurf(i),0.0)
    101103!
     
    103105        qstar(i) = (cdrah(i) * zdq(i) * speed(i))/ustar(i)
    104106        lmon(i) = (ustar(i) * ustar(i) * tpot(i))/ &
    105                   (RKAR * RG * testar(i))
     107 &                (RKAR * RG * testar(i))
    106108      ENDDO
    107109!
     
    109111      zref = 2.0
    110112      CALL screenp(klon, knon, nsrf, speed, tpot, q1, &
    111                    ts1, qsurf, rugos, lmon, &
    112                    ustar, testar, qstar, zref, &
    113                    delu, delte, delq)
     113 &                 ts1, qsurf, rugos, lmon, &
     114 &                 ustar, testar, qstar, zref, &
     115 &                 delu, delte, delq)
    114116!
    115117      DO i = 1, knon
     
    129131        okri=.TRUE.
    130132        CALL screenc(klon, knon, nsrf, zxli, &
    131                      u_zref, temp, q_zref, zref, &
    132                      ts1, qsurf, rugos, psol, &           
    133                      ustar, testar, qstar, okri, ri1, &
    134                      pref, delu, delte, delq)
     133 &                   u_zref, temp, q_zref, zref, &
     134 &                   ts1, qsurf, rugos, psol, &           
     135 &                   ustar, testar, qstar, okri, ri1, &
     136 &                   pref, delu, delte, delq)
    135137!
    136138        DO i = 1, knon
     
    211213      zref = 10.0
    212214      CALL screenp(klon, knon, nsrf, speed, tpot, q1, &
    213                    ts1, qsurf, rugos, lmon, &
    214                    ustar, testar, qstar, zref, &
    215                    delu, delte, delq)
     215 &                 ts1, qsurf, rugos, lmon, &
     216 &                 ustar, testar, qstar, zref, &
     217 &                 delu, delte, delq)
    216218!
    217219      DO i = 1, knon
     
    231233        okri=.TRUE.
    232234        CALL screenc(klon, knon, nsrf, zxli, &
    233                      u_zref, temp, q_zref, zref, &
    234                      ts1, qsurf, rugos, psol, &
    235                      ustar, testar, qstar, okri, ri1, &
    236                      pref, delu, delte, delq)
     235 &                   u_zref, temp, q_zref, zref, &
     236 &                   ts1, qsurf, rugos, psol, &
     237 &                   ustar, testar, qstar, okri, ri1, &
     238 &                   pref, delu, delte, delq)
    237239!
    238240        DO i = 1, knon
Note: See TracChangeset for help on using the changeset viewer.