1 : C FLURUN SOURCE CHAT 05/01/13 00:04:56 5004 2 : SUBROUTINE FLURUN(SIGI,EPSINI,TETIN,TRINS,HPAS2,ICLFLU,NCOURB, 3 : . IBOU,ITYP,BTAA,BETA,BTAEPS,BTAPAS,ICOD,DSIGK,DEPSK,SIGF,EPSINF, 4 : . WEP,EPSINT,PREC,ecou) 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8(A-H,O-Z) 7 : DIMENSION SIGI(6),DSIGK(6),SIGF(6) 8 : C 9 : SEGMENT ECOU 10 : *** COMMON/ECOU/TEST,ALFAH, 11 : REAL*8 TEST, ALFAH, 12 : 1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6), 13 : 2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),SIG5(6),SIGINT(6), 14 : 1 DALPHA(6),DSIGO(6),E(12),XINV(3), 15 : 2 SIPLAD(6),DSIGP0(6),TET,TETI 16 : ENDSEGMENT 17 : COMMON/ECOU/TEST,ALFAH, 18 : 1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6), 19 : 2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),SIG5(6),SIGINT(6), 20 : 1 DALPHA(6),DSIGO(6),E(12),XINV(3), 21 : 2 SIPLAD(6),DSIGP0(6),TET,TETI 22 : C 23 : ICOD=0 24 : PREBTA=2.D00 25 : PREZER=1.D-7 26 : SINI=VNMISD(SIGI,ITYP,ALFAH,CVNMSD) 26 : CALL CRPLAW(VI,EPSINI,SINI,TETIN,TRINS,HPAS2,ICLFLU,NCOURB) 27 : DEPSK=VI*HPAS2 28 : IF(DEPSK.LT.PREC) GO TO 10 29 : BETA=BETA*0.5 30 : ICOD=1 31 : RETURN 32 : 10 EPSINF=EPSINT+DEPSK*WEP 33 : BTASIG=DEPSK*BTAEPS 34 : IF(BTASIG.GT.PREBTA) BTASIG=PREBTA 35 : IF(BTAA.LE.PREZER) BTASIG=BTAPAS 36 : ELTA=0. 37 : IF(SINI.NE.0.) ELTA=VI/SINI 38 : ELTF=ELTA*HPAS2 39 : DO 20 IB=1,IBOU 40 : DSIGK(IB)=DSIGP0(IB)*BTAPAS+(DSIGP(IB)-DSIGP0(IB))*BTASIG 41 : . -E(IB)*SIGI(IB)*ELTF 42 : SIGF(IB)=SIGINT(IB)+DSIGK(IB)*WEP 43 : 20 CONTINUE 44 : RETURN 45 : END 46 : 47 : 48 :
© Cast3M 2003 - All rights reserved.
Disclaimer