Sources Esope | Cast3M

flurun
C FLURUN    SOURCE    OF166741  25/11/04    21:15:54     12349          
      SUBROUTINE FLURUN(SIGI,EPSINI,TETIN,TRINS,HPAS2,ICLFLU,NCOUR,
     . IBOU,ITYPL,BTAA,BETA,BTAEPS,BTAPAS,ICOD,DSIGK,DEPSK,SIGF,EPSINF,
     . WEP,EPSINT,PREC,ecou)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC TECOU
 
      DIMENSION SIGI(6),DSIGK(6),SIGF(6)
C
      ICOD=0
      PREBTA=2.D00
      PREZER=1.D-7
      SINI=VNMISD(SIGI,ITYPL,ecou.ALFAH,CVNMSD)
      CALL CRPLAW(VI,EPSINI,SINI,TETIN,TRINS,HPAS2,ICLFLU,NCOUR)
      DEPSK = VI*HPAS2
      IF (DEPSK.LT.PREC) GO TO 10
      BETA=BETA*0.5
      ICOD=1
      RETURN
10    CONTINUE
      EPSINF=EPSINT+DEPSK*WEP
      BTASIG=DEPSK*BTAEPS
      IF(BTASIG.GT.PREBTA) BTASIG=PREBTA
      IF(BTAA.LE.PREZER) BTASIG=BTAPAS
      ELTA=0.
      IF (SINI.NE.0.) ELTA=VI/SINI
      ELTF = ELTA*HPAS2
      DO 20 IB=1,IBOU
        DSIGK(IB) =  ecou.DSIGP0(IB)*BTAPAS
     &             + (ecou.DSIGP(IB)-ecou.DSIGP0(IB))*BTASIG
     &             - ecou.E(IB)*SIGI(IB)*ELTF
        SIGF(IB) = ecou.W2(IB)+DSIGK(IB)*WEP
20    CONTINUE
      RETURN
      END