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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
!
!-----------------------------------------------------------------------
!
!Purpose:
! Fill the TabSegArc structure with information recovered in a
! surfacic file.
!
!Copyright:
! Copyright (C) 2001 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):
! G. Civario (CS-SI)
!
!-----------------------------------------------------------------------
!
module generTabSegArc
use constUtiles
use segArc
use boundCond
use SALGET_FUNS_MOD
use precision_and_kinds, only : pdb
implicit none
contains
subroutine generateTabSegArc(ipSal,sizeSA,nbNode,nbCLP,nbFlux,merg,impx)
integer,intent(inout) :: ipSal
integer,intent(in) :: nbNode,sizeSA,impx
integer,intent(out) :: nbCLP,nbFlux
integer,dimension(nbNode),intent(out) :: merg
integer, parameter :: n_datain=25, n_datare=20
integer, dimension (n_datain) :: datain
real, dimension (n_datare) :: datare
real(pdb), dimension (n_datare) :: datade
integer :: type,nber,prec,elem,i,nbMacro,fout0
integer, parameter, dimension(0:4) :: read_bc_len=(/1,1,2,3,3/)
character(len=12) :: name_geom
! internal : albedo
! vacuum surface : albedo
! specular reflexion : none
! translation : tx ty (t=translation vector)
! rotation : cx cy cos(theta) sin(theta) theta
! (c= center,theta= axis angle)
! axial symmetry : cx cy cos(theta) sin(theta) theta
! (c= center,theta= axis angle)
! central symetry : cx cy (c= center)
integer, allocatable, dimension(:) :: iflux, medium
integer,parameter :: dimTabCelluleBase = 20000
fout0=6
if(impx == 0) fout0=0
call SALGET(datain,6,ipSal,fout0,'dimensions for geometry')
if(nbNode /= datain(3)) call XABORT('g2s_generateTabSegArc: nbNode error')
if(sizeSA /= datain(4)) call XABORT('g2s_generateTabSegArc: sizeSA error')
typgeo=datain(1)
nbMacro=datain(5)
nbFlux=datain(6)
call SALGET(datain,3,ipSal,fout0,'index kndex prec')
prec=datain(3)
call SALGET(datare,1,ipSal,fout0,'eps')
call SALGET(merg,nbNode,ipSal,fout0,'flux index per node')
call SALGET(name_geom,ipSal,fout0,'names of macros')
allocate(iflux(nbFlux),stat=alloc_ok)
if (alloc_ok /= 0) call XABORT("G2S: generateTabSegArc(1) => allocation pb")
call SALGET(iflux,nbFlux,ipSal,fout0,'macro order number per flux region')
deallocate(iflux)
do elem=1,sizeSA
call SALGET(datain,3,ipSal,fout0,'integer descriptors')
type=datain(1)
tabSegArc(elem)%typ=type
if(type == 1) then
tabSegArc(elem)%noded=datain(2)
tabSegArc(elem)%nodeg=datain(3)
else
tabSegArc(elem)%nodeg=datain(2)
tabSegArc(elem)%noded=datain(3)
endif
select case (type)
case (1)
nber=4
case (2)
nber=3
case (3)
nber=5
case default
write(fout0,'(1x,''==> sal126: unknown type '',i3)')type
call xabort('g2_generateTabSegArc: unknown type')
end select
call SALGET(datade,nber,ipSal,fout0,prec,'real descriptors')
tabSegArc(elem)%x=datade(1)
tabSegArc(elem)%y=datade(2)
select case (type)
case (1)
tabSegArc(elem)%dx=datade(1)+datade(3)
tabSegArc(elem)%dy=datade(2)+datade(4)
case (2)
tabSegArc(elem)%r=datade(3)
tabSegArc(elem)%a=0.0
tabSegArc(elem)%b=0.0
case (3)
tabSegArc(elem)%r=datade(3)
tabSegArc(elem)%a=datade(4)/rad2deg
tabSegArc(elem)%b=(datade(4)+datade(5))/rad2deg
end select
enddo
call SALGET(datain,3,ipSal,fout0,'general bc data')
nbCLP=datain(2)
call SALGET(datade(1),ipSal,fout0,prec,'general albedo')
do i=1,nbCLP
call SALGET(datain,2,ipSal,fout0,'specific bc: type nber')
type=datain(1)
nber=datain(2)
SALbCDataTab(i)%SALtype=type
SALbCDataTab(i)%nber=nber
allocate(SALbCDataTab(i)%elemNb(nber),stat=alloc_ok)
if (alloc_ok /= 0) call XABORT("G2S: generateTabSegArc(3) => allocation pb")
call SALGET(SALbCDataTab(i)%elemNb,nber,ipSal,fout0,'bc elements')
! read bc motion
call SALGET(datade,read_bc_len(type),ipSal,fout0,prec,'data for specific bc condition')
select case(type)
case(0,1)
SALbCDataTab(i)%albedo=real(datade(1))
case(2)
SALbCDataTab(i)%tx=real(datade(1))
SALbCDataTab(i)%ty=real(datade(2))
case(3,4)
SALbCDataTab(i)%cx=real(datade(1))
SALbCDataTab(i)%cy=real(datade(2))
SALbCDataTab(i)%angle=real(datade(3))
end select
enddo
allocate(medium(dimTabCelluleBase))
call SALGET(medium,nbNode,ipSal,fout0,'media per node')
do i=1,sizeSA
if(tabSegArc(i)%nodeg>0) tabSegArc(i)%neutronicMixg=medium(tabSegArc(i)%nodeg)
if(tabSegArc(i)%noded>0) tabSegArc(i)%neutronicMixd=medium(tabSegArc(i)%noded)
enddo
deallocate(medium)
end subroutine generateTabSegArc
end module generTabSegArc
|
