Numérotation des lignes :

C CROIS     SOURCE    STRU      08/12/22    21:15:07     6246c     Sous-Programme du module MISTRAL0C       --------------------------------------------------------------------------        SUBROUTINE CROIS (TT,FI,FIT, PCROI, VEPSCR)C       --------------------------------------------------------------------------C       Calcul des vitesses de deformation de croissance VEPSCR(IJ)C       en fonction :C       - de la temperature TT,C       - du flux de neutrons rapides FI,C       - de la fluence FIT.C       --------------------------------------------------------------------------      IMPLICIT INTEGER (I-N)      IMPLICIT REAL*8 (A-H, O-Z)        PARAMETER ( NVC=2  )      DIMENSION PCROI(1:*)        DIMENSION VCR(1:3,1:NVC),C(1:3),CFIT(1:3)        DIMENSION VEPSCR(1:*)        FITPL = FIT+PCROI(1)      NPCROI = 1        DO IVC = 1,NVC       NPCROI = NPCROI+1       KCI = NINT(PCROI(NPCROI))       IF (KCI.EQ.0) THEN        DO IJ =1,3         VCR(IJ,IVC) = 0.        END DO       ELSE          IF (KCI.EQ.-1) THEN         DO IJ =1,3          NPCROI = NPCROI+1          C(IJ) = PCROI(NPCROI)         END DO        ELSE IF (KCI.EQ.1) THEN         DO IJ =1,3          C(IJ) =     &     PCROI(NPCROI+1)+(PCROI(NPCROI+2)+PCROI(NPCROI+3)*TT)*TT          NPCROI = NPCROI+3         END DO        ELSE IF (KCI.EQ.2) THEN         COEF = EXP( PCROI(NPCROI+4)/TT     &    +(PCROI(NPCROI+5)+(PCROI(NPCROI+6)+PCROI(NPCROI+7)*TT)*TT)*TT     &    )         DO IJ =1,3          C(IJ) = PCROI(NPCROI+IJ)*COEF         END DO         NPCROI = NPCROI+7        END IF        COEF = (PCROI(NPCROI+4)+PCROI(NPCROI+5)*TT)*TT        DO IJ =1,3         CFIT(IJ) = PCROI(NPCROI+IJ)+COEF        END DO        NPCROI = NPCROI+6        DO IJ =1,3         VCR(IJ,IVC) = C(IJ)*FITPL**CFIT(IJ)*FI**PCROI(NPCROI)        END DO       END IF        END DO      IF (KCI.EQ.0) THEN       DO IJ = 1,3        VEPSCR(IJ) = VCR(IJ,1)       END DO      ELSE       FITTC = PCROI(NPCROI+2)+PCROI(NPCROI+3)*TT       FT = 1./(1.+EXP(PCROI(NPCROI+1)*(FIT/FITTC-1.)))       DO IJ = 1,3        VEPSCR(IJ) = VCR(IJ,1)*FT+VCR(IJ,2)*(1.-FT)       END DO      END IF        DO IJ = 4,6       VEPSCR(IJ) = 0.        END DO        RETURN        END

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