scale_ocean_grid_icoa Module Reference

module ocean / grid / icosahedralA More...

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

Variables
real(rp), dimension(:), allocatable, public	ocean_grid_icoa_cz
	center coordinate [m]: z, local=global More...
real(rp), dimension(:), allocatable, public	ocean_grid_icoa_fz
	face coordinate [m]: z, local=global More...
real(rp), dimension(:), allocatable, public	ocean_grid_icoa_cdz
	z-length of control volume [m] More...

Detailed Description

module ocean / grid / icosahedralA

Description

Grid module for icosahedral coordinate for ocean

Author

Team SCALE

NAMELIST

• PARAM_OCEAN_GRID_ICOA

  name	type	default value	comment
  OCEAN_GRID_ICOA_IN_BASENAME	character(len=H_LONG)	''
  OCEAN_GRID_ICOA_IN_AGGREGATE	logical
  FZ	real(RP), dimension(FZ_MAX)		face coordinate without surface (=0 m)

  
  

History Output

No history output

Function/Subroutine Documentation

ocean_grid_icoa_setup()

subroutine, public scale_ocean_grid_icoa::ocean_grid_icoa_setup

Setup.

Definition at line 60 of file scale_ocean_grid_icoA.F90.

    use scale_prc, only: &
       prc_abort
    use scale_file, only: &
       file_aggregate
    implicit none
 
    namelist / param_ocean_grid_icoa / &
       ocean_grid_icoa_in_basename,  &
       ocean_grid_icoa_in_aggregate, &
       fz
 
    integer :: ierr
    integer :: k
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("OCEAN_GRID_ICOA_setup",*) 'Setup'
 
    if ( okmax < 1 ) then
       log_info("OCEAN_GRID_ICOA_setup",*) 'Skip because OKMAX < 1'
       return
    end if
 
    fz(:) = 0.0_rp
 
    ocean_grid_icoa_in_aggregate = file_aggregate
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_ocean_grid_icoa,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("OCEAN_GRID_ICOA_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("OCEAN_GRID_ICOA_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_GRID_ICOA. Check!'
       call prc_abort
    endif
    log_nml(param_ocean_grid_icoa)
 
    allocate( ocean_grid_icoa_cz(oks  :oke) )
    allocate( ocean_grid_icoa_fz(oks-1:oke) )
    allocate( ocean_grid_icoa_cdz(oks  :oke) )
    !$acc enter data create(OCEAN_GRID_ICOA_CZ,OCEAN_GRID_ICOA_FZ,OCEAN_GRID_ICOA_CDZ)
 
    log_newline
    log_info("OCEAN_GRID_ICOA_setup",*) 'Ocean grid information '
 
    if ( ocean_grid_icoa_in_basename /= '' ) then
       call ocean_grid_icoa_read
    else
       log_info("OCEAN_GRID_ICOA_setup",*) 'Not found input grid file. Grid position is calculated.'
 
       call ocean_grid_icoa_generate
    endif
 
    if ( oke == oks ) then
       log_newline
       log_info("OCEAN_GRID_ICOA_setup",*) 'Single layer. ODZ = ', ocean_grid_icoa_cdz(1)
    else
       log_newline
       log_info("OCEAN_GRID_ICOA_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_icoa_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_icoa_cz(k),'        ',ocean_grid_icoa_cdz(k),'     | '
       log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',ocean_grid_icoa_fz(k),'       |',k,' | '
       enddo
       k = oke
       log_info_cont('(1x,A,I4,F8.3,A,F8.3,A)') '|',k,ocean_grid_icoa_cz(k),'        ',ocean_grid_icoa_cdz(k),'     | '
       log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',ocean_grid_icoa_fz(k),'        ',k,' | layer of no motion'
       log_info_cont('(1x,A)')                  '|=================================|'
    endif
 
    !$acc update device(OCEAN_GRID_ICOA_CZ,OCEAN_GRID_ICOA_FZ,OCEAN_GRID_ICOA_CDZ)
    return

References scale_file::file_aggregate, scale_file::file_open(), scale_io::io_fid_conf, ocean_grid_icoa_cdz, ocean_grid_icoa_cz, ocean_grid_icoa_fz, scale_ocean_grid_icoa_index::oka, scale_ocean_grid_icoa_index::oke, scale_ocean_grid_icoa_index::okmax, scale_ocean_grid_icoa_index::oks, scale_prc::prc_abort(), and scale_prc::prc_myrank.

Here is the call graph for this function:

Variable Documentation

ocean_grid_icoa_cz

real(rp), dimension (:), allocatable, public scale_ocean_grid_icoa::ocean_grid_icoa_cz

center coordinate [m]: z, local=global

Definition at line 35 of file scale_ocean_grid_icoA.F90.

  real(RP), public, allocatable :: OCEAN_GRID_ICOA_CZ (:)

Referenced by ocean_grid_icoa_setup().

ocean_grid_icoa_fz

real(rp), dimension (:), allocatable, public scale_ocean_grid_icoa::ocean_grid_icoa_fz

face coordinate [m]: z, local=global

Definition at line 36 of file scale_ocean_grid_icoA.F90.

  real(RP), public, allocatable :: OCEAN_GRID_ICOA_FZ (:)

Referenced by ocean_grid_icoa_setup().

ocean_grid_icoa_cdz

real(rp), dimension(:), allocatable, public scale_ocean_grid_icoa::ocean_grid_icoa_cdz

z-length of control volume [m]

Definition at line 37 of file scale_ocean_grid_icoA.F90.

  real(RP), public, allocatable :: OCEAN_GRID_ICOA_CDZ(:)

Referenced by ocean_grid_icoa_setup().