Télécharger la plaquette Cast3M


fr en

Télécharger crois.eso

Retour à la liste

Numérotation des lignes : 

crois

  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