d6/d09/scale__grid__index_8F90_source.html

 !-------------------------------------------------------------------------------
 !-------------------------------------------------------------------------------
 module scale_grid_index
   !-----------------------------------------------------------------------------
   !
   !++ used modules
   !
   use scale_precision
   use scale_stdio
   use scale_prof
   !-----------------------------------------------------------------------------
   implicit none
   private
   !-----------------------------------------------------------------------------
   !
   !++ Public procedure
   !
   public :: grid_index_setup

   !-----------------------------------------------------------------------------
   !
   !++ Public parameters & variables
   !
   integer, public, parameter :: zdir  = 1
   integer, public, parameter :: xdir  = 2
   integer, public, parameter :: ydir  = 3

 #ifdef _FIXEDINDEX_
   include "inc_index.h"
   include "inc_index_common.h"
 #else
   integer, public :: kmax   = -1
   integer, public :: imax   = -1
   integer, public :: jmax   = -1

   integer, public :: iblock = -1
   integer, public :: jblock = -1

   integer, public, parameter :: khalo  = 2
   integer, public            :: ihalo  = 2
   integer, public            :: jhalo  = 2

   integer, public :: ka
   integer, public :: ia
   integer, public :: ja

   integer, public :: ks
   integer, public :: ke
   integer, public :: is
   integer, public :: ie
   integer, public :: js
   integer, public :: je

   integer, public :: kijmax = -1
 #endif

   ! indices considering boundary
   integer, public :: imaxb
   integer, public :: jmaxb
   integer, public :: isb
   integer, public :: ieb
   integer, public :: jsb
   integer, public :: jeb
   integer, public :: ieh
   integer, public :: jeh

   integer, public :: isg
   integer, public :: ieg
   integer, public :: jsg
   integer, public :: jeg

   !-----------------------------------------------------------------------------
   !
   !++ Private procedure
   !
   !-----------------------------------------------------------------------------
   !
   !++ Private parameters & variables
   !
   !-----------------------------------------------------------------------------
 contains
   !-----------------------------------------------------------------------------
   subroutine grid_index_setup
     use scale_process, only: &
        prc_mpistop, &
        prc_myrank
     use scale_rm_process, only: &
        prc_2drank,  &
        prc_num_x,   &
        prc_num_y,   &
        prc_has_w,   &
        prc_has_e,   &
        prc_has_s,   &
        prc_has_n
     implicit none

 #ifndef _FIXEDINDEX_
     namelist / param_index / &
        kmax,   &
        imax,   &
        jmax,   &
        ihalo,  &
        jhalo,  &
        iblock, &
        jblock
 #endif

     integer :: ierr
     !---------------------------------------------------------------------------

     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[GRID_INDEX] / Categ[ATMOS-RM GRID] / Origin[SCALElib]'

 #ifdef _FIXEDINDEX_
     if( io_l ) write(io_fid_log,*) '*** No namelists.'
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** fixed index mode'
 #else
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_index,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        write(*,*) 'xxx Not appropriate names in namelist PARAM_INDEX. Check!'
        call prc_mpistop
     endif
     if( io_lnml ) write(io_fid_log,nml=param_index)

     if ( imax < ihalo ) then
        write(*,*) 'xxx number of grid size IMAX must >= IHALO! ', imax, ihalo
        call prc_mpistop
     endif
     if ( jmax < jhalo ) then
        write(*,*) 'xxx number of grid size JMAX must >= JHALO! ', jmax, jhalo
        call prc_mpistop
     endif

     ka   = kmax + khalo * 2
     ia   = imax + ihalo * 2
     ja   = jmax + jhalo * 2

     ks   = 1    + khalo
     ke   = kmax + khalo
     is   = 1    + ihalo
     ie   = imax + ihalo
     js   = 1    + jhalo
     je   = jmax + jhalo

     if( iblock == -1 ) iblock = imax
     if( jblock == -1 ) jblock = jmax

     kijmax = kmax * imax * jmax
 #endif

     !-- Block size must be divisible
     if    ( mod(imax,iblock) > 0 ) then
        write(*,*) 'xxx number of grid size IMAX must be divisible by IBLOCK! ', imax, iblock
        call prc_mpistop
     elseif( mod(jmax,jblock) > 0 ) then
        write(*,*) 'xxx number of grid size JMAX must be divisible by JBLOCK! ', jmax, jblock
        call prc_mpistop
     endif

     ! horizontal index (global domain)
     isg = ihalo + 1    + prc_2drank(prc_myrank,1) * imax
     ieg = ihalo + imax + prc_2drank(prc_myrank,1) * imax
     jsg = jhalo + 1    + prc_2drank(prc_myrank,2) * jmax
     jeg = jhalo + jmax + prc_2drank(prc_myrank,2) * jmax

     ! index considering boundary region
     imaxb = imax
     jmaxb = jmax
     isb = is
     ieb = ie
     jsb = js
     jeb = je
     ieh = ie
     jeh = je
     if ( .NOT. prc_has_w ) then
        imaxb = imaxb + ihalo
        isb = 1
     endif
     if ( .NOT. prc_has_e ) then
        imaxb = imaxb + ihalo
        ieb = ia
        ieh = ie - 1
     endif
     if ( .NOT. prc_has_s ) then
        jmaxb = jmaxb + jhalo
        jsb = 1
     endif
     if ( .NOT. prc_has_n ) then
        jmaxb = jmaxb + jhalo
        jeb = ja
        jeh = je - 1
     endif

     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** Atmosphere grid index information ***'
     if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6,A,I6)') '*** No. of Computational Grid (global)  :', &
                                                        kmax,           ' x ',                       &
                                                        imax*prc_num_x, ' x ',                       &
                                                        jmax*prc_num_y
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6,A,I6)') '*** No. of Computational Grid (local)   :', &
                                                        kmax,' x ',imax,' x ',jmax
     if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6,A,I6)') '*** No. of Grid (including HALO, local) :', &
                                                        ka," x ",ia," x ",ja
     if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6)')      '*** Global index of local grid (X)      :', &
                                                        isg," - ",ieg
     if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6)')      '*** Global index of local grid (Y)      :', &
                                                        jsg," - ",jeg

   end subroutine grid_index_setup

 end module scale_grid_index
scale_grid_index::imax
integer, public imax
of computational cells: x
Definition: scale_grid_index.F90:44

scale_rm_process::prc_num_x
integer, public prc_num_x
x length of 2D processor topology
Definition: scale_rm_process.F90:43

scale_grid_index::is
integer, public is
start point of inner domain: x, local
Definition: scale_grid_index.F90:60

scale_grid_index::je
integer, public je
end point of inner domain: y, local
Definition: scale_grid_index.F90:63

scale_rm_process::prc_has_n
logical, public prc_has_n
Definition: scale_rm_process.F90:50

scale_process::prc_mpistop
subroutine, public prc_mpistop
Abort MPI.
Definition: scale_process.F90:234

scale_grid_index::jeb
integer, public jeb
Definition: scale_grid_index.F90:74

scale_grid_index::iblock
integer, public iblock
block size for cache blocking: x
Definition: scale_grid_index.F90:47

scale_stdio::io_l
logical, public io_l
output log or not? (this process)
Definition: scale_stdio.F90:59

scale_grid_index::zdir
integer, parameter, public zdir
Definition: scale_grid_index.F90:35

scale_rm_process::prc_has_e
logical, public prc_has_e
Definition: scale_rm_process.F90:51

scale_stdio
module STDIO
Definition: scale_stdio.F90:12

scale_grid_index::ydir
integer, parameter, public ydir
Definition: scale_grid_index.F90:37

scale_grid_index::ke
integer, public ke
end point of inner domain: z, local
Definition: scale_grid_index.F90:59

scale_grid_index::xdir
integer, parameter, public xdir
Definition: scale_grid_index.F90:36

scale_grid_index::grid_index_setup
subroutine, public grid_index_setup
Setup.
Definition: scale_grid_index.F90:96

scale_grid_index::imaxb
integer, public imaxb
Definition: scale_grid_index.F90:69

scale_grid_index::jsg
integer, public jsg
start point of the inner domain: y, global
Definition: scale_grid_index.F90:80

scale_grid_index::jmaxb
integer, public jmaxb
Definition: scale_grid_index.F90:70

scale_rm_process::prc_has_s
logical, public prc_has_s
Definition: scale_rm_process.F90:52

scale_grid_index::ieb
integer, public ieb
Definition: scale_grid_index.F90:72

scale_grid_index::ieg
integer, public ieg
end point of the inner domain: x, global
Definition: scale_grid_index.F90:79

scale_rm_process::prc_num_y
integer, public prc_num_y
y length of 2D processor topology
Definition: scale_rm_process.F90:44

scale_grid_index::isg
integer, public isg
start point of the inner domain: x, global
Definition: scale_grid_index.F90:78

scale_grid_index
module grid index
Definition: scale_grid_index.F90:14

scale_grid_index::ia
integer, public ia
of x whole cells (local, with HALO)
Definition: scale_grid_index.F90:55

scale_grid_index::ka
integer, public ka
of z whole cells (local, with HALO)
Definition: scale_grid_index.F90:54

scale_grid_index::jblock
integer, public jblock
block size for cache blocking: y
Definition: scale_grid_index.F90:48

scale_grid_index::kmax
integer, public kmax
of computational cells: z
Definition: scale_grid_index.F90:43

scale_grid_index::jhalo
integer, public jhalo
of halo cells: y
Definition: scale_grid_index.F90:52

scale_grid_index::js
integer, public js
start point of inner domain: y, local
Definition: scale_grid_index.F90:62

scale_process
module PROCESS
Definition: scale_process.F90:14

scale_grid_index::ks
integer, public ks
start point of inner domain: z, local
Definition: scale_grid_index.F90:58

scale_process::prc_myrank
integer, public prc_myrank
process num in local communicator
Definition: scale_process.F90:78

scale_grid_index::khalo
integer, parameter, public khalo
of halo cells: z
Definition: scale_grid_index.F90:50

scale_grid_index::jeh
integer, public jeh
end point of inner domain: y, local (half level)
Definition: scale_grid_index.F90:76

scale_grid_index::ieh
integer, public ieh
end point of inner domain: x, local (half level)
Definition: scale_grid_index.F90:75

scale_grid_index::kijmax
integer, public kijmax
of computational cells: z*x*y
Definition: scale_grid_index.F90:65

scale_rm_process
module RM PROCESS
Definition: scale_rm_process.F90:10

scale_prof
module profiler
Definition: scale_prof.F90:10

scale_grid_index::ie
integer, public ie
end point of inner domain: x, local
Definition: scale_grid_index.F90:61

scale_grid_index::jeg
integer, public jeg
end point of the inner domain: y, global
Definition: scale_grid_index.F90:81

scale_stdio::io_lnml
logical, public io_lnml
output log or not? (for namelist, this process)
Definition: scale_stdio.F90:60

scale_precision
module PRECISION
Definition: scale_precision.F90:13

scale_rm_process::prc_2drank
integer, dimension(:,:), allocatable, public prc_2drank
node index in 2D topology
Definition: scale_rm_process.F90:46

scale_grid_index::isb
integer, public isb
Definition: scale_grid_index.F90:71

scale_stdio::io_fid_conf
integer, public io_fid_conf
Config file ID.
Definition: scale_stdio.F90:55

scale_stdio::io_fid_log
integer, public io_fid_log
Log file ID.
Definition: scale_stdio.F90:56

scale_grid_index::jsb
integer, public jsb
Definition: scale_grid_index.F90:73

scale_rm_process::prc_has_w
logical, public prc_has_w
Definition: scale_rm_process.F90:49

scale_grid_index::jmax
integer, public jmax
of computational cells: y
Definition: scale_grid_index.F90:45

scale_grid_index::ihalo
integer, public ihalo
of halo cells: x
Definition: scale_grid_index.F90:51

scale_grid_index::ja
integer, public ja
of y whole cells (local, with HALO)
Definition: scale_grid_index.F90:56