Télécharger flurun.eso

Retour à la liste

Numérotation des lignes :

  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)
  27. CALL CRPLAW(VI,EPSINI,SINI,TETIN,TRINS,HPAS2,ICLFLU,NCOURB)
  28. DEPSK=VI*HPAS2
  29. IF(DEPSK.LT.PREC) GO TO 10
  30. BETA=BETA*0.5
  31. ICOD=1
  32. RETURN
  33. 10 EPSINF=EPSINT+DEPSK*WEP
  34. BTASIG=DEPSK*BTAEPS
  35. IF(BTASIG.GT.PREBTA) BTASIG=PREBTA
  36. IF(BTAA.LE.PREZER) BTASIG=BTAPAS
  37. ELTA=0.
  38. IF(SINI.NE.0.) ELTA=VI/SINI
  39. ELTF=ELTA*HPAS2
  40. DO 20 IB=1,IBOU
  41. DSIGK(IB)=DSIGP0(IB)*BTAPAS+(DSIGP(IB)-DSIGP0(IB))*BTASIG
  42. . -E(IB)*SIGI(IB)*ELTF
  43. SIGF(IB)=SIGINT(IB)+DSIGK(IB)*WEP
  44. 20 CONTINUE
  45. RETURN
  46. END
  47.  
  48.  
  49.  
  50.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales