diff options
Diffstat (limited to 'Dragon/src/MCGFFIR.f')
| -rw-r--r-- | Dragon/src/MCGFFIR.f | 131 |
1 files changed, 131 insertions, 0 deletions
diff --git a/Dragon/src/MCGFFIR.f b/Dragon/src/MCGFFIR.f new file mode 100644 index 0000000..91ad708 --- /dev/null +++ b/Dragon/src/MCGFFIR.f @@ -0,0 +1,131 @@ +*DECK MCGFFIR + SUBROUTINE MCGFFIR(SUBSCH,K,KPN,M,N,H,NOM,NZON,XST,S,NREG,KEYFLX, + 1 KEYCUR,F,B,W,OMEGA2,IDIR,NSOUT,XSI) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Solution of transport equation on a track +* ray-tracing (isotropic scattering case,'source term isolation' on). +* +*Copyright: +* Copyright (C) 2002 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version +* +*Author(s): I. Suslov and R. Le Tellier +* +*Parameters: input +* SUBSCH track coefficients calculation subroutine. +* K total number of volumes for which specific values +* of the neutron flux and reactions rates are required. +* KPN total number of unknowns in vectors F. +* M number of material mixtures. +* N number of elements in the current track. +* H vector containing the lenght of the different segments of this +* track. +* NOM vector containing the region number of the different segments +* of this track. +* NZON index-number of the mixture type assigned to each volume. +* XST macroscopic total cross section. +* S total source vector. +* NREG number of volumes. +* KEYFLX position of flux elements in PHI vector. +* KEYCUR position of current elements in PHI vector. +* W weight associated with this track. +* OMEGA2 square x, y and z-component of the direction +* Omega for 2D geometry. +* IDIR direction of fundamental current for TIBERE with MoC +* (=0,1,2,3). +* NSOUT number of outer surfaces. +* XSI x,y and z component of the shape parameter for TIBERE. +* +*Parameters: input/output +* F vector containing the zonal scalar flux. +* +*Parameters: scratch +* B undefined. +* +*----------------------------------------------------------------------- +* + IMPLICIT NONE +*---- +* SUBROUTINE ARGUMENTS +*---- + INTEGER K,KPN,M,N,NOM(N),NZON(K),NREG,KEYFLX(NREG,1), + 1 KEYCUR(K-NREG),IDIR,NSOUT + REAL XST(0:M) + DOUBLE PRECISION W,H(N),S(KPN),F(KPN),B(N),OMEGA2(3) + DOUBLE PRECISION XSI(NSOUT) + EXTERNAL SUBSCH +*--- +* LOCAL VARIABLES +*--- + DOUBLE PRECISION F0,RP,RM,WW,SSS + INTEGER I,J,NOMI,IND,IND1,INDN,NOMJ,INDC +* + WW=DBLE(W) +*---- +* Calculation coefficients for this track. +*---- +* MCGSCA: Step-Characteristics Scheme with Tabulated Exponentials +* MCGDDF: Diamond-Differencing Scheme +* MCGSCE: Step-Characteristics Scheme with Exact Exponentials + CALL SUBSCH(N,K,M,NOM,NZON,H,XST,B) +*---- +* Summation along the track in both directions +*---- +* incoming flux in + direction + IND1=KEYCUR(-NOM(1)) + RP=S(IND1) +* incoming flux in - direction + INDN=KEYCUR(-NOM(N)) + RM=S(INDN) + IF(IDIR.GT.0) THEN +* incoming flux in + direction + RP=RP*OMEGA2(IDIR)/XSI(IND1-NREG) +* incoming flux in - direction + RM=RM*OMEGA2(IDIR)/XSI(INDN-NREG) + ENDIF + DO I=2,N-1 +* + direction + NOMI=NOM(I) + SSS=S(NOMI) + IF(IDIR.GT.0) THEN + SSS=SSS*OMEGA2(IDIR) + ENDIF + F0=B(I)*RP + RP=RP+B(I)*(SSS-DBLE(XST(NZON(NOMI)))*RP) + IND=KEYFLX(NOMI,1) + F(IND)=F(IND)+F0*WW + INDC=KPN/2+IND +* Calculate Xi, Yi and Zi for TIBERE + IF(IDIR.GT.0) THEN + F(INDC)=F(INDC)+F0*WW*OMEGA2(IDIR) + ENDIF +* - direction + J=N+1-I + NOMJ=NOM(J) + SSS=S(NOMJ) + IF(IDIR.GT.0) THEN + SSS=SSS*OMEGA2(IDIR) + ENDIF + F0=B(J)*RM + RM=RM+B(J)*(SSS-DBLE(XST(NZON(NOMJ)))*RM) + IND=KEYFLX(NOMJ,1) + F(IND)=F(IND)+F0*WW + INDC=KPN/2+IND +* Calculate Xi, Yi and Zi for TIBERE + IF(IDIR.GT.0) THEN + F(INDC)=F(INDC)+F0*WW*OMEGA2(IDIR) + ENDIF + ENDDO +* outgoing flux in + direction + F(INDN)=F(INDN)+RP*WW +* outgoing flux in - direction + F(IND1)=F(IND1)+RM*WW +* + RETURN + END |