scale_ocean_grid_cartesC_index.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_ocean_grid_cartesc_index
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_grid_cartesc_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

  integer, public :: oka   = -1 ! # 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

  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

  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

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine ocean_grid_cartesc_index_setup
    use scale_prc, only: &
         prc_abort
    use scale_atmos_grid_cartesc_index, only: &
       imax,       &
       ia, is, ie, &
       jmax,       &
       ja, js, je
    implicit none

    namelist / param_ocean_grid_cartesc_index / &
       okmax

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

    log_newline
    log_info("OCEAN_GRID_CARTESC_INDEX_setup",*) 'Setup'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_ocean_grid_cartesc_index,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("OCEAN_GRID_CARTESC_INDEX_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("OCEAN_GRID_CARTESC_INDEX_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_GRID_CARTESC_INDEX. Check!'
       call prc_abort
    endif
    log_nml(param_ocean_grid_cartesc_index)

    if ( okmax < 1 ) then
       log_error("OCEAN_GRID_CARTESC_INDEX_setup",*) 'OKMAX must be >= 1 ', okmax
       call prc_abort
    end if

    oka  = okmax
    oks  = 1
    oke  = okmax

    log_newline
    log_info("OCEAN_GRID_CARTESC_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

    return
  end subroutine ocean_grid_cartesc_index_setup

end module scale_ocean_grid_cartesc_index