Array Organisation in FASTEST

Subroutine setind

In FASTEST the Control Volumes are identified by the index ijk. These indices depend on the block and grid and are set in subroutine setind.F.

       call setind (ngr,l
#        include "listsetind.h"

  

      ngr     
:
is the actual gridlevel and

l
:
is the actual blocknumber.

Knowing these two values setind.F calculates the values for:

  

ni, nj, nk
:
The max. index in i-, j-,    k-direction.

nim,njm,nkm
:
The max. index minus ``1'' for each direction.

nij,nik,njk
:
The max. index for ij-, ik-,    jk-face.


li(i)
:
li(1)=offset for actual grid/block

 
:
li(i)= offset for jk-plane, i.e. all CVs with    constant i

lk(k)
:
lk(k)= offset for ij-plane, i.e. all CVs with    constant k

ijs, iks, jks

:
see setcon.F (2.1.3).

ijks
:
offset for grid l, block m

icst
:
first CV in grid l, block m

icen
:
last CV in grid l, block m

iis
:
=  icst ???.

iie
:
last inner CV  in grid l, block m (i.e. not    lying on block boundary)

After calling setind.F you can calculate indices ijk to identify the ith, jth, kth control volume of block m, grid l in a global array with the control volumes 1,...,nxyza of all grids and blocks (function ha(i,j,k) is defined in hahihjhk.h):
ijk=ha(i,j,k)

Subroutine setcon

       call setcon (
#        include "listsetind.h"

Calculates offsets for block m=1,...,nblo and grid l=1,...,ngr.

• ijkgit(l,m) = offset for block m, grid l
  (= first cv in actual block/grid -1)
• isbij(l,m) = offset for first ij-plane in block m, grid l
  (means k=1; i=1,...,ni; j=1,...,nj)
• isbik(l,m) = offset for first ik-plane
• isbjk(l,m) = offset for first jk-plane

Include file hahihjhk.h

In include file hahihjhk.h functions for the calculation of indices are defined.

• ha(ii,jj,kk) = lk(kk)+li(ii)+jj
  = ijk-CV in actual block and grid (according to block and grid in setind)
• hi(jj,kk) = (jj-1)*nk+kk+jks
  = first jk-plane in actual block and grid (quasi west boundary ?)
• hj(ii,kk) = (ii-1)*nk+kk+iks
  = first ik-plane
• hk(ii,jj) = (ii-1)*nk+jj+ijs
  = first ij-plane