diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Dragon/src/MCGDDFT.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/src/MCGDDFT.f')
| -rw-r--r-- | Dragon/src/MCGDDFT.f | 59 |
1 files changed, 59 insertions, 0 deletions
diff --git a/Dragon/src/MCGDDFT.f b/Dragon/src/MCGDDFT.f new file mode 100644 index 0000000..a545f84 --- /dev/null +++ b/Dragon/src/MCGDDFT.f @@ -0,0 +1,59 @@ +*DECK MCGDDFT + SUBROUTINE MCGDDFT(N,K,M,NOM,NZON,H,XST,B) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Calculate coefficients of a track for the characteristics integration. +* Diamond-Differencing scheme without fix-up. with +* 'MOCC/MCI' integration strategy. +* +*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): R. Le Tellier +* +*Parameters: input +* N number of elements in the current track. +* K total number of volumes for which specific values +* of the neutron flux and reactions rates are required. +* M number of material mixtures. +* NOM vector containing the region number of the different segments +* of this track. +* NZON index-number of the mixture type assigned to each volume. +* H vector containing the lenght of the different segments of this +* track. +* XST macroscopic total cross section. +* and step characteristics (SC). +* +*Parameters: output +* B DD1 coefficients. +* +*----------------------------------------------------------------------- +* + IMPLICIT NONE +*---- +* SUBROUTINE ARGUMENTS +*---- + INTEGER N,K,M,NOM(N),NZON(K) + REAL XST(0:M) + DOUBLE PRECISION H(N),B(N) +*--- +* LOCAL VARIABLES +*--- + INTEGER I,NOMI,NZI + DOUBLE PRECISION TAUD +* + DO I=2,N-1 + NOMI=NOM(I) + NZI=NZON(NOMI) + TAUD=H(I)*XST(NZI) + B(I)=2.D0*TAUD/(2.D0+TAUD) + ENDDO +* + RETURN + END |