Esope source files

   1 : C MSHEA3    SOURCE    CHAT      05/01/13    01:53:01     5004           
   2 :       SUBROUTINE MSHEA3(BETA,EXPN, EPLUS,EMOIN, APLUS,AMOIN)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : C
   6 : C======================================================================
   7 : C       calcul de APLUS et AMOIN
   8 : C======================================================================
   9 : C=======================================================================
  10 : C
  11 : C LISTE D'ECHANGE
  12 : C ---------------
  13 : C
  14 : C BETA         = parametre beta
  15 : C EXPN         = exposant n
  16 : C EPLUS        = energie (sens positif)
  17 : C EMOIN        =         (sens négatif)
  18 : C
  19 : C APLUS        = attenuation cyclique (sens positif)
  20 : C AMOIN        =                      (sens négatif)
  21 : C
  22 : C=======================================================================
  23 :       REAL*8 BETA,EXPN, EPLUS,EMOIN, APLUS,AMOIN
  24 :       PARAMETER (XZER=0.D0,UN=1.D0)
  25 : *
  26 :       IF(EXPN.EQ.XZER)THEN
  27 :         APLUS=TANH(BETA*EPLUS)
  28 :         AMOIN=TANH(BETA*EMOIN)
  29 :       ELSE
  30 :         APLUS=((BETA*EPLUS)**EXPN)/(UN+((BETA*EPLUS)**EXPN))
  31 :         AMOIN=((BETA*EMOIN)**EXPN)/(UN+((BETA*EMOIN)**EXPN))
  32 :       ENDIF
  33 : *
  34 :       RETURN
  35 :       END
  36 : C+2004 ND
  37 :  
  38 :