1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
*DECK FMAC04
SUBROUTINE FMAC04(NGPRT,NGP,NPART,NK,H)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Reorganize boundary information from FMAC-M array.
*
*Copyright:
* Copyright (C) 2020 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): C. Bienvenue
*
*Parameters: input
* NGPRT number of energy groups per particle.
* NGP number of energy groups for all particles.
* NPART number of particles in FMAC-M file.
* NK number of mixtures
* H boundary informations array from FMAC-M file.
*
*Parameters: output
* H boundary informations array organized by particle.
*
*-----------------------------------------------------------------------
*
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER NGPRT(NPART),NGP,NPART
REAL H(NK,NGP+NPART),HTEMP(NK,NGP+NPART)
*
NGT1=0
NGT2=0
DO I=1,NPART
DO J=1,NGPRT(I)+1
DO K=1,NK
IF(J.LE.NGPRT(I)) THEN
HTEMP(k,NGT2+J)=H(k,NGT1+J)
ELSE
HTEMP(k,NGT2+J)=H(k,NGP+NPART+1-I)
ENDIF
ENDDO
ENDDO
NGT1=NGT1+NGPRT(I)
NGT2=NGT2+NGPRT(I)+1
ENDDO
H=HTEMP
RETURN
END
|
