Télécharger crois.eso

Retour à la liste

Numérotation des lignes :

  1. C CROIS SOURCE STRU 08/12/22 21:15:07 6246
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE CROIS (TT,FI,FIT, PCROI, VEPSCR)
  5. C --------------------------------------------------------------------------
  6. C Calcul des vitesses de deformation de croissance VEPSCR(IJ)
  7. C en fonction :
  8. C - de la temperature TT,
  9. C - du flux de neutrons rapides FI,
  10. C - de la fluence FIT.
  11. C --------------------------------------------------------------------------
  12. IMPLICIT INTEGER (I-N)
  13. IMPLICIT REAL*8 (A-H, O-Z)
  14. PARAMETER ( NVC=2 )
  15. DIMENSION PCROI(1:*)
  16. DIMENSION VCR(1:3,1:NVC),C(1:3),CFIT(1:3)
  17. DIMENSION VEPSCR(1:*)
  18. FITPL = FIT+PCROI(1)
  19. NPCROI = 1
  20. DO IVC = 1,NVC
  21. NPCROI = NPCROI+1
  22. KCI = NINT(PCROI(NPCROI))
  23. IF (KCI.EQ.0) THEN
  24. DO IJ =1,3
  25. VCR(IJ,IVC) = 0.
  26. END DO
  27. ELSE
  28. IF (KCI.EQ.-1) THEN
  29. DO IJ =1,3
  30. NPCROI = NPCROI+1
  31. C(IJ) = PCROI(NPCROI)
  32. END DO
  33. ELSE IF (KCI.EQ.1) THEN
  34. DO IJ =1,3
  35. C(IJ) =
  36. & PCROI(NPCROI+1)+(PCROI(NPCROI+2)+PCROI(NPCROI+3)*TT)*TT
  37. NPCROI = NPCROI+3
  38. END DO
  39. ELSE IF (KCI.EQ.2) THEN
  40. COEF = EXP( PCROI(NPCROI+4)/TT
  41. & +(PCROI(NPCROI+5)+(PCROI(NPCROI+6)+PCROI(NPCROI+7)*TT)*TT)*TT
  42. & )
  43. DO IJ =1,3
  44. C(IJ) = PCROI(NPCROI+IJ)*COEF
  45. END DO
  46. NPCROI = NPCROI+7
  47. END IF
  48. COEF = (PCROI(NPCROI+4)+PCROI(NPCROI+5)*TT)*TT
  49. DO IJ =1,3
  50. CFIT(IJ) = PCROI(NPCROI+IJ)+COEF
  51. END DO
  52. NPCROI = NPCROI+6
  53. DO IJ =1,3
  54. VCR(IJ,IVC) = C(IJ)*FITPL**CFIT(IJ)*FI**PCROI(NPCROI)
  55. END DO
  56. END IF
  57. END DO
  58. IF (KCI.EQ.0) THEN
  59. DO IJ = 1,3
  60. VEPSCR(IJ) = VCR(IJ,1)
  61. END DO
  62. ELSE
  63. FITTC = PCROI(NPCROI+2)+PCROI(NPCROI+3)*TT
  64. FT = 1./(1.+EXP(PCROI(NPCROI+1)*(FIT/FITTC-1.)))
  65. DO IJ = 1,3
  66. VEPSCR(IJ) = VCR(IJ,1)*FT+VCR(IJ,2)*(1.-FT)
  67. END DO
  68. END IF
  69. DO IJ = 4,6
  70. VEPSCR(IJ) = 0.
  71. END DO
  72. RETURN
  73. END
  74.  
  75.  

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