Esope source files

   1 : C EQUIV1    SOURCE    CHAT      05/01/12    23:44:44     5004           
   2 : C EQUIV1     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE EQUIV1(DEFR0,EPSEQC,EQSTRC,EX,RB,EPO1,ICAL)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 : C            ********************************************
   9 : C            **  CALCUL  APPROCHE DE EPSEQ1 ET EQSTR1 ***
  10 : C            ** COURBE DE COMPRESSION  INSA-HYPERBOLE ***
  11 : C            ********************************************
  12 : C
  13 :       IF(ICAL.EQ.0) THEN
  14 :         NBIT=200
  15 :         PREC=1.D-8
  16 :         EPEQ0=DABS(DEFR0)
  17 :         EPEQ1=EPO1
  18 :         DO 10 I=1,NBIT
  19 :           EPEQ=EPEQ0+(EPEQ1-EPEQ0)/2.D0
  20 :           SEQ1=EX*(EPEQ-DABS(DEFR0))
  21 :           SEQ2=EX*EPEQ/(1.D0+(EX/RB-2.D0/EPO1)*EPEQ+(EPEQ/EPO1)**2)
  22 :           CONT=DABS((SEQ1-SEQ2)/SEQ1)
  23 :           IF(CONT.LT.PREC) GOTO 20 
  24 :           IF(SEQ1.GT.SEQ2) THEN
  25 :             EPEQ1=EPEQ
  26 :           ELSE
  27 :             EPEQ0=EPEQ
  28 :           ENDIF
  29 :  10     CONTINUE
  30 :  20     CONTINUE
  31 :       ELSE
  32 :         EPEQ=EPO1
  33 :         SEQ2=RB
  34 :       ENDIF
  35 : C      
  36 :       IF(EPEQ.GT.EPO1) EPEQ=EPO1
  37 :       IF(SEQ2.GT.RB) SEQ2=RB
  38 :       EPSEQC=-1.D0*EPEQ
  39 :       EQSTRC=-1.D0*SEQ2
  40 :       DEFR0 = EPSEQC-EQSTRC/EX
  41 : C      
  42 :       RETURN
  43 :       END
  44 :  
  45 :