Changeset 251 for LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F

Timestamp:

Jul 4, 2001, 12:46:06 PM (23 years ago)

Author:

lmdz

Message:

Initialisations pour eviter des divisions par 0
LF

File:

• LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F (modified) (1 diff)

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

@@ -79 +79 @@
 
 
-       IF( 200.* fb .LT. - fa )   THEN
-         fhyp ( i) = - 1.
-       ELSEIF( 200. * fb .LT. fa ) THEN
-         fhyp ( i) =   1.
-       ELSE
-         fhyp(i) =  TANH ( fa/fb )
-       ENDIF
+       if ( xtild(i) .ne. 0. .and. xtild(i) .ne. pi) then
+         IF( 200.* fb .LT. - fa )   THEN
+           fhyp ( i) = - 1.
+         ELSEIF( 200. * fb .LT. fa ) THEN
+           fhyp ( i) =   1.
+         ELSE
+           fhyp(i) =  TANH ( fa/fb )
+         ENDIF
+       endif
 
        IF ( xtild(i).EQ. 0. )  fhyp(i) =  1.