Sources Esope | Cast3M

flurun
C FLURUN    SOURCE    CHAT      05/01/13    00:04:56     5004
      SUBROUTINE FLURUN(SIGI,EPSINI,TETIN,TRINS,HPAS2,ICLFLU,NCOURB,
     . IBOU,ITYP,BTAA,BETA,BTAEPS,BTAPAS,ICOD,DSIGK,DEPSK,SIGF,EPSINF,
     . WEP,EPSINT,PREC,ecou)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION SIGI(6),DSIGK(6),SIGF(6)
C
        SEGMENT ECOU
***     COMMON/ECOU/TEST,ALFAH,
       REAL*8 TEST, ALFAH,
     1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6),
     2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),SIG5(6),SIGINT(6),
     1 DALPHA(6),DSIGO(6),E(12),XINV(3),
     2 SIPLAD(6),DSIGP0(6),TET,TETI
      ENDSEGMENT
      COMMON/ECOU/TEST,ALFAH,
     1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6),
     2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),SIG5(6),SIGINT(6),
     1 DALPHA(6),DSIGO(6),E(12),XINV(3),
     2 SIPLAD(6),DSIGP0(6),TET,TETI
C
      ICOD=0
      PREBTA=2.D00
      PREZER=1.D-7
      SINI=VNMISD(SIGI,ITYP,ALFAH,CVNMSD)
      CALL CRPLAW(VI,EPSINI,SINI,TETIN,TRINS,HPAS2,ICLFLU,NCOURB)
      DEPSK=VI*HPAS2
      IF(DEPSK.LT.PREC) GO TO 10
      BETA=BETA*0.5
      ICOD=1
      RETURN
10    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)=DSIGP0(IB)*BTAPAS+(DSIGP(IB)-DSIGP0(IB))*BTASIG
     . -E(IB)*SIGI(IB)*ELTF
      SIGF(IB)=SIGINT(IB)+DSIGK(IB)*WEP
20    CONTINUE
      RETURN
      END