Ignore:
Timestamp:
Jun 22, 2001, 11:51:02 AM (23 years ago)
Author:
lmdzadmin
Message:

Pb d'initialisation surle Nec
LF

Location:
LMDZ.3.3/branches/rel-LF/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/fxhyp.F

    r232 r242  
    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/branches/rel-LF/libf/dyn3d/fyhyp.F

    r232 r242  
    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.