Esope source files

   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 :