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

Initialisations pour eviter des divisions par 0
LF

File:
1 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.
Note: See TracChangeset for help on using the changeset viewer.