Esope source files

   1 : C PARDEC    SOURCE    CHAT      05/01/13    02:09:08     5004           
   2 : C PARDEC     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE PARDEC(ST,SIE,EPSC,EQSTR,EDC,IPLA,DEFR0,RB,ALPHA,EPO1,
   4 :      & EX,EMAX,ICAL)
   5 : C
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 : C
   9 : C             ***********************************************
  10 : C             **  CALCUL :  EPSC , EQSTR , EDC             **
  11 : C             ***********************************************
  12 : C
  13 :       FPT=ALPHA*RB
  14 :       DEFRL = RB/EX-EPO1
  15 :       DEFR0 = ST-SIE/EX
  16 : C           
  17 :       IF(DEFR0.GT.DEFRL) THEN
  17 :          CALL EQUIV1(DEFR0,EPSC,EQSTR,EX,RB,EPO1,ICAL)
  18 :       ELSE
  19 :          A3  = RB/(EPO1-EMAX)
  20 :          B1  = -1.D0*EX*DEFR0
  21 :          B3  = A3*EMAX
  22 :          EPSC=(B3-B1)/(EX-A3)
  23 :          EQSTR=EX*EPSC+B1
  24 :          IPLA=2
  25 :       ENDIF 
  26 :       EDC=(DABS(EQSTR)+RB)/(DABS(EPSC)+(RB/EX))
  27 : C
  28 :       IF(EDC.GT.EX.OR.EDC.LT.1.D-5) EDC=EX
  29 :       DEFR0 = EPSC-EQSTR/EDC
  30 :       EQSTR=DABS(EQSTR)
  31 :       IF(DABS(DEFR0).LT.1.D-8) DEFR0=0.D0
  32 : C     
  33 :       RETURN
  34 :       END
  35 :           
  36 :           
  37 :  
  38 :