gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Add new global variables to the configuration.

Browse Source 
CONFIGURATION

 - add two new variables which determined the half and the double of
   the ghost cell size (nh and nd, respectively);
 - add a variable to determine the half size of the block (ih, jh, and
   kh);
This commit is contained in:
Grzegorz Kowal 2010-09-18 11:50:24 +02:00
parent d54b57e1eb
commit db478af9e4
1 changed files with 13 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/config.F90
View File

@ -34,14 +34,14 @@ module config
! block dimensions, number of ghost zones ! block dimensions, number of ghost zones
! !
integer(kind=4), save :: ncells = 8, nghost = 4, ngrids = 16 & integer(kind=4), save :: ncells = 8, nghost = 4, ngrids = 16 &
, igrids = 16, jgrids = 16, kgrids = 16 , igrids = 16, jgrids = 16, kgrids = 16
! derived dimensional variables ! derived dimensional variables
! !
integer(kind=4), save :: in, im, ib, ie, ibl, ibu, iel, ieu, ng & integer(kind=4), save :: in, im, ib, ie, ibl, ibu, iel, ieu, ih, ng &
, jn, jm, jb, je, jbl, jbu, jel, jeu & , jn, jm, jb, je, jbl, jbu, jel, jeu, jh, nd &
, kn, km, kb, ke, kbl, kbu, kel, keu , kn, km, kb, ke, kbl, kbu, kel, keu, kh, nh
! the dimensions of the lowest level ! the dimensions of the lowest level
! !
@ -304,9 +304,12 @@ module config
#endif /* */ #endif /* */
ng = nghost ng = nghost
nd = ng * 2
nh = ng / 2
in = ncells in = ncells
im = in + 2 * ng im = in + nd
ih = im / 2
ib = ng + 1 ib = ng + 1
ie = ng + in ie = ng + in
ibl = ib - 1 ibl = ib - 1
@ -315,7 +318,8 @@ module config
ieu = ie + 1 ieu = ie + 1
jn = ncells jn = ncells
jm = jn + 2 * ng jm = jn + nd
jh = jm / 2
jb = ng + 1 jb = ng + 1
je = ng + jn je = ng + jn
jbl = jb - 1 jbl = jb - 1
@ -326,6 +330,7 @@ module config
#if NDIMS == 2 #if NDIMS == 2
kn = 1 kn = 1
km = 1 km = 1
kh = 1
kb = 1 kb = 1
ke = 1 ke = 1
kbl = 1 kbl = 1
@ -335,7 +340,8 @@ module config
#endif /* NDIMS == 2 */ #endif /* NDIMS == 2 */
#if NDIMS == 3 #if NDIMS == 3
kn = ncells kn = ncells
km = kn + 2 * ng km = kn + nd
kh = km / 2
kb = ng + 1 kb = ng + 1
ke = ng + kn ke = ng + kn
kbl = kb - 1 kbl = kb - 1