module ocean / grid / cartesianC More...

Functions/Subroutines
subroutine, public	ocean_grid_cartesc_setup
	Setup. More...

subroutine, public	ocean_grid_cartesc_finalize
	Finalize. More...

Variables
real(rp), dimension(:), allocatable, public	ocean_grid_cartesc_cz
	center coordinate [m]: z, local=global More...

real(rp), dimension(:), allocatable, public	ocean_grid_cartesc_fz
	face coordinate [m]: z, local=global More...

real(rp), dimension(:), allocatable, public	ocean_grid_cartesc_cdz
	z-length of control volume [m] More...

Detailed Description

module ocean / grid / cartesianC

Description: Grid module for cartesian coordinate for ocean

Author: Team SCALE

NAMELIST

PARAM_OCEAN_GRID_CARTESC

name type default value comment

OCEAN_GRID_CARTESC_IN_BASENAME character(len=H_LONG) ''

OCEAN_GRID_CARTESC_IN_AGGREGATE logical

ODZ real(RP), dimension(ODZ_MAX) layer thickness for input

History Output: No history output

Function/Subroutine Documentation

◆ ocean_grid_cartesc_setup()

subroutine, public scale_ocean_grid_cartesc::ocean_grid_cartesc_setup

Setup.

Definition at line 62 of file scale_ocean_grid_cartesC.F90.

     use scale_prc, only: &
        prc_abort
     use scale_file, only: &
        file_aggregate
     implicit none
  
     namelist / param_ocean_grid_cartesc / &
        ocean_grid_cartesc_in_basename,  &
        ocean_grid_cartesc_in_aggregate, &
        odz
  
     integer :: ierr
     integer :: k
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("OCEAN_GRID_CARTESC_setup",*) 'Setup'
  
     if ( okmax < 1 ) then
        log_info("OCEAN_GRID_CARTESC_setup",*) 'Skip because OKMAX < 1'
        return
     end if
  
     odz(:) = 10.0_rp
  
     ocean_grid_cartesc_in_aggregate = file_aggregate
  
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_ocean_grid_cartesc,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("OCEAN_GRID_CARTESC_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("OCEAN_GRID_CARTESC_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_GRID_CARTESC. Check!'
        call prc_abort
     endif
     log_nml(param_ocean_grid_cartesc)
  
     allocate( ocean_grid_cartesc_cz(oks  :oke) )
     allocate( ocean_grid_cartesc_fz(oks-1:oke) )
     allocate( ocean_grid_cartesc_cdz(oks  :oke) )
  
     log_newline
     log_info("OCEAN_GRID_CARTESC_setup",*) 'Ocean grid information '
  
     if ( ocean_grid_cartesc_in_basename /= '' ) then
        call ocean_grid_cartesc_read
     else
        log_info("OCEAN_GRID_CARTESC_setup",*) 'Not found input grid file. Grid position is calculated.'
  
        call ocean_grid_cartesc_generate
     endif
     !$acc enter data create(OCEAN_GRID_CARTESC_CZ,OCEAN_GRID_CARTESC_FZ,OCEAN_GRID_CARTESC_CDZ)
  
     if ( oke == oks ) then
        log_newline
        log_info("OCEAN_GRID_CARTESC_setup",*) 'Single layer. ODZ = ', ocean_grid_cartesc_cdz(1)
     else
        log_newline
        log_info("OCEAN_GRID_CARTESC_setup",'(1x,A)') 'Vertical Coordinate'
        log_info_cont('(1x,A)')                  '|   k       z      zh      dz   k |'
        log_info_cont('(1x,A)')                  '|         [m]     [m]     [m]     |'
        k = oks-1
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',ocean_grid_cartesc_fz(k),'        ',k,' | Atmosphere interface'
        do k = oks, oke-1
        log_info_cont('(1x,A,I4,F8.3,A,F8.3,A)') '|',k,ocean_grid_cartesc_cz(k),'        ',ocean_grid_cartesc_cdz(k),'     | '
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',ocean_grid_cartesc_fz(k),'       |',k,' | '
        enddo
        k = oke
        log_info_cont('(1x,A,I4,F8.3,A,F8.3,A)') '|',k,ocean_grid_cartesc_cz(k),'        ',ocean_grid_cartesc_cdz(k),'     | '
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',ocean_grid_cartesc_fz(k),'        ',k,' | layer of no motion'
        log_info_cont('(1x,A)')                  '|=================================|'
     endif
  
     !$acc update device(OCEAN_GRID_CARTESC_CZ,OCEAN_GRID_CARTESC_FZ,OCEAN_GRID_CARTESC_CDZ)
  
     return

References scale_file::file_aggregate, scale_io::io_fid_conf, ocean_grid_cartesc_cdz, ocean_grid_cartesc_cz, ocean_grid_cartesc_fz, scale_ocean_grid_cartesc_index::oke, scale_ocean_grid_cartesc_index::okmax, scale_ocean_grid_cartesc_index::oks, and scale_prc::prc_abort().

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ocean_grid_cartesc_finalize()

subroutine, public scale_ocean_grid_cartesc::ocean_grid_cartesc_finalize

Finalize.

Definition at line 144 of file scale_ocean_grid_cartesC.F90.

     implicit none
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("OCEAN_GRID_CARTESC_finalize",*) 'Finalize'
  
     !$acc exit data delete(OCEAN_GRID_CARTESC_CZ,OCEAN_GRID_CARTESC_FZ,OCEAN_GRID_CARTESC_CDZ)
     deallocate( ocean_grid_cartesc_cz  )
     deallocate( ocean_grid_cartesc_fz  )
     deallocate( ocean_grid_cartesc_cdz )
  
     return