diff options
Diffstat (limited to 'Trivac/src/FLDBH1.f')
| -rwxr-xr-x | Trivac/src/FLDBH1.f | 55 |
1 files changed, 55 insertions, 0 deletions
diff --git a/Trivac/src/FLDBH1.f b/Trivac/src/FLDBH1.f new file mode 100755 index 0000000..aea864f --- /dev/null +++ b/Trivac/src/FLDBH1.f @@ -0,0 +1,55 @@ +*DECK FLDBH1 + SUBROUTINE FLDBH1 (NEL,NUN,LL4,EVECT,VOL,IDL,KN,QFR) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Calculation of the averaged flux with a mesh centered finite +* difference method in hexagonal geometry with triangular +* mesh-splitting. +* +*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): A. Hebert +* +*Parameters: input +* NEL number of hexagons. +* NUN number of unknowns per energy group. +* LL4 order of the system matrices. +* EVECT variational coefficients of the flux. The information is +* contained in position EVECT(1) to EVECT(LL4). +* VOL volume of each element. +* IDL position of the average flux component associated with each +* volume. +* KN element-ordered unknown list. +* QFR element-ordered information. +* +*Parameters: output +* EVECT averaged fluxes. The information is contained in positions +* EVECT(IDL(I)). +* +*----------------------------------------------------------------------- +* +*---- +* SUBROUTINE ARGUMENTS +*---- + INTEGER NEL,NUN,LL4,IDL(NEL),KN(4*LL4) + REAL EVECT(NUN),VOL(NEL),QFR(4*LL4) +* + NSURF=3 + DO 10 K=1,NEL + IF(IDL(K).NE.0) EVECT(IDL(K))=0.0 + 10 CONTINUE + NUM1=0 + DO 20 IND1=1,LL4 + K=KN(NUM1+NSURF+1) + EVECT(IDL(K))=EVECT(IDL(K))+QFR(NUM1+NSURF+1)*EVECT(IND1)/VOL(K) + NUM1=NUM1+NSURF+1 + 20 CONTINUE + RETURN + END |