scale_ocean_grid_icoA_index.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_ocean_grid_icoa_index
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_grid_icoa_index_setup
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  integer, public :: okmax = -1 ! # of computational cells: z for ocean
  integer, public :: oimax = -1 ! # of computational cells: x for ocean
  integer, public :: ojmax = -1 ! # of computational cells: y for ocean
  !$acc declare create(OKMAX,OIMAX,OJMAX)
 
  integer, public :: oka       ! # of total grids: z for ocean, local
  integer, public :: oks       ! start point of inner domain: z for ocean, local
  integer, public :: oke       ! end   point of inner domain: z for ocean, local
  !$acc declare create(OKA,OKS,OKE)
 
  integer, public :: oia       ! # of total grids: x for ocean, local
  integer, public :: ois       ! start point of inner domain: x for ocean, local
  integer, public :: oie       ! end   point of inner domain: x for ocean, local
  !$acc declare create(OIA,OIS,OIE)
 
  integer, public :: oja       ! # of total grids: y for ocean, local
  integer, public :: ojs       ! start point of inner domain: y for ocean, local
  integer, public :: oje       ! end   point of inner domain: y for ocean, local
  !$acc declare create(OJA,OJS,OJE)
 
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine ocean_grid_icoa_index_setup
    use scale_prc, only: &
         prc_abort
    use scale_atmos_grid_icoa_index, only: &
       imax,       &
       ia, is, ie, &
       jmax,       &
       ja, js, je
    implicit none
 
    namelist / param_ocean_grid_icoa_index / &
       okmax
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("OCEAN_GRID_ICOA_INDEX_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_ocean_grid_icoa_index,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("OCEAN_GRID_ICOA_INDEX_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("OCEAN_GRID_ICOA_INDEX_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_GRID_ICOA_INDEX. Check!'
       call prc_abort
    endif
    log_nml(param_ocean_grid_icoa_index)
 
    if ( okmax < 1 ) then
       log_error("OCEAN_GRID_ICOA_INDEX_setup",*) 'OKMAX must be >= 1 ', okmax
       call prc_abort
    end if
 
    oka  = okmax
    oks  = 1
    oke  = okmax
 
    log_newline
    log_info("OCEAN_GRID_ICOA_INDEX_setup",*) 'Ocean grid index information '
    log_info_cont('(1x,A,I6,A,I6,A,I6)') 'z-axis levels :', okmax
 
 
    ! at this moment horizontal grid is same as that in atmosphere
    oimax = imax
    oia = ia
    ois = is
    oie = ie
 
    ojmax = jmax
    oja = ja
    ojs = js
    oje = je
 
    !$acc update device(OKMAX,OIMAX,OJMAX)
    !$acc update device(OKA,OKS,OKE)
    !$acc update device(OIA,OIS,OIE)
    !$acc update device(OJA,OJS,OJE)
    return
  end subroutine ocean_grid_icoa_index_setup
 
end module scale_ocean_grid_icoa_index