module MKTOPO More...

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

subroutine, public	mktopo
	Driver. More...

Variables
integer, public	mktopo_type = -1

integer, parameter, public	i_ignore = 0

integer, parameter, public	i_flat = 1

integer, parameter, public	i_bellshape = 2

integer, parameter, public	i_schaer = 3

Detailed Description

module MKTOPO

Description: subroutines for preparing topography data (ideal case)

Author: Team SCALE

History

2012-12-26 (H.Yashiro) [new]

NAMELIST

PARAM_MKTOPO

name type default value comment

MKTOPO_NAME character(len=H_SHORT) 'NONE'
PARAM_MKTOPO_FLAT

name type default value comment

FLAT_HEIGHT real(RP) 100.0_RP height of mountain [m]

PARAM_MKTOPO_BELLSHAPE

name	type	default value	comment
BELL_EACHNODE	logical	.false.	Arrange mountain at each node? [kg/kg]
BELL_CX	real(RP)	2.E3_RP	center location [m]: x
BELL_CY	real(RP)	2.E3_RP	center location [m]: y
BELL_RX	real(RP)	2.E3_RP	bubble radius [m]: x
BELL_RY	real(RP)	2.E3_RP	bubble radius [m]: y
BELL_HEIGHT	real(RP)	100.0_RP	height of mountain [m]

PARAM_MKTOPO_SCHAER

name	type	default value	comment
SCHAER_CX	real(RP)	25.E3_RP	center location [m]: x
SCHAER_RX	real(RP)	5.E3_RP	bubble radius [m]: x
SCHAER_LAMBDA	real(RP)	4.E3_RP	wavelength of wavelike perturbation [m]: x
SCHAER_HEIGHT	real(RP)	250.0_RP	height of mountain [m]
SCHAER_SWAPXY	logical	.false.

History Output: No history output

Function/Subroutine Documentation

◆ mktopo_setup()

subroutine, public mod_mktopo::mktopo_setup ( )

Setup.

Definition at line 70 of file mod_mktopo.f90.

References i_bellshape, i_flat, i_ignore, i_schaer, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, mktopo_type, and scale_process::prc_mpistop().

Referenced by mod_rm_prep::scalerm_prep().

     implicit none
 
     character(len=H_SHORT) :: MKTOPO_name = 'NONE'
 
     namelist / param_mktopo / &
        mktopo_name
 
     integer :: ierr
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[make topo] / Categ[preprocess] / Origin[SCALE-RM]'
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mktopo,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_MKTOPO. Check!'
        call prc_mpistop
     endif
     if( io_nml ) write(io_fid_nml,nml=param_mktopo)
 
     select case(mktopo_name)
     case('NONE')
        mktopo_type = i_ignore
     case('FLAT')
        mktopo_type = i_flat
     case('BELLSHAPE')
        mktopo_type = i_bellshape
     case('SCHAER')
        mktopo_type = i_schaer
     case default
        write(*,*) 'xxx Unsupported TYPE:', trim(mktopo_name)
        call prc_mpistop
     endselect
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ mktopo()

subroutine, public mod_mktopo::mktopo ( )

Driver.

Definition at line 114 of file mod_mktopo.f90.

References scale_const::const_pi, i_bellshape, i_flat, i_ignore, i_schaer, scale_grid_index::ia, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, scale_grid_index::is, scale_grid_index::ja, scale_grid_index::js, mktopo_type, scale_process::prc_mpistop(), scale_topography::topo_write(), and scale_topography::topo_zsfc.

Referenced by mod_rm_prep::scalerm_prep().

     use scale_topography, only: &
        topo_write
     implicit none
     !---------------------------------------------------------------------------
 
     if ( mktopo_type == i_ignore ) then
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '++++++ SKIP  MAKING TOPOGRAPHY DATA ++++++'
     else
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '++++++ START MAKING TOPOGRAPHY DATA ++++++'
 
        select case(mktopo_type)
        case(i_flat)
           call mktopo_flat
 
        case(i_bellshape)
           call mktopo_bellshape
 
        case(i_schaer)
           call mktopo_schaer
 
        case default
           write(*,*) 'xxx Unsupported TYPE:', mktopo_type
           call prc_mpistop
        endselect
 
        if( io_l ) write(io_fid_log,*) '++++++ END   MAKING TOPOGRAPHY DATA ++++++'
 
        ! output topography file
        call topo_write
     endif
 
     return