Changeset 251 for LMDZ.3.3/trunk


Ignore:
Timestamp:
Jul 4, 2001, 12:46:06 PM (23 years ago)
Author:
lmdz
Message:

Initialisations pour eviter des divisions par 0
LF

Location:
LMDZ.3.3/trunk/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F

    r216 r251  
    7979
    8080
    81        IF( 200.* fb .LT. - fa )   THEN
    82          fhyp ( i) = - 1.
    83        ELSEIF( 200. * fb .LT. fa ) THEN
    84          fhyp ( i) =   1.
    85        ELSE
    86          fhyp(i) =  TANH ( fa/fb )
    87        ENDIF
     81       if ( xtild(i) .ne. 0. .and. xtild(i) .ne. pi) then
     82         IF( 200.* fb .LT. - fa )   THEN
     83           fhyp ( i) = - 1.
     84         ELSEIF( 200. * fb .LT. fa ) THEN
     85           fhyp ( i) =   1.
     86         ELSE
     87           fhyp(i) =  TANH ( fa/fb )
     88         ENDIF
     89       endif
    8890
    8991       IF ( xtild(i).EQ. 0. )  fhyp(i) =  1.

  • LMDZ.3.3/trunk/libf/dyn3d/fyhyp.F

    r212 r251  
    1 c
     1C
    22c $Header$
    33c
    … …  
    100100       y0max    = 2.*y0*heavyy0  + pis2
    101101
     102       fa = 999.999
     103       fb = 999.999
     104
    102105       DO i = 0, nmax2
    103106        IF( yt(i).LT.y0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.