C MSHEA3    SOURCE    CHAT      05/01/13    01:53:01     5004
      SUBROUTINE MSHEA3(BETA,EXPN, EPLUS,EMOIN, APLUS,AMOIN)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
C======================================================================
C       calcul de APLUS et AMOIN
C======================================================================
C=======================================================================
C
C LISTE D'ECHANGE
C ---------------
C
C BETA         = parametre beta
C EXPN         = exposant n
C EPLUS        = energie (sens positif)
C EMOIN        =         (sens négatif)
C
C APLUS        = attenuation cyclique (sens positif)
C AMOIN        =                      (sens négatif)
C
C=======================================================================
      REAL*8 BETA,EXPN, EPLUS,EMOIN, APLUS,AMOIN
      PARAMETER (XZER=0.D0,UN=1.D0)
*
      IF(EXPN.EQ.XZER)THEN
        APLUS=TANH(BETA*EPLUS)
        AMOIN=TANH(BETA*EMOIN)
      ELSE
        APLUS=((BETA*EPLUS)**EXPN)/(UN+((BETA*EPLUS)**EXPN))
        AMOIN=((BETA*EMOIN)**EXPN)/(UN+((BETA*EMOIN)**EXPN))
      ENDIF
*
      RETURN
      END
C+2004 ND