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

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 68 of file mod_mktopo.F90.

References i_bellshape, i_flat, i_ignore, i_schaer, scale_io::io_fid_conf, mktopo_type, and scale_prc::prc_abort().

Referenced by mod_rm_prep::rm_prep().

     implicit none
 
     character(len=H_SHORT) :: mktopo_name = 'NONE'
 
     namelist / param_mktopo / &
        mktopo_name
 
     integer :: ierr
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("MKTOPO_setup",*) 'Setup'
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mktopo,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("MKTOPO_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("MKTOPO_setup",*) 'Not appropriate names in namelist PARAM_MKTOPO. Check!'
        call prc_abort
     endif
     log_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
        log_error("MKTOPO_setup",*) 'Unsupported TYPE:', trim(mktopo_name)
        call prc_abort
     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 112 of file mod_mktopo.F90.

References scale_const::const_pi, i_bellshape, i_flat, i_ignore, i_schaer, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::js, mktopo_type, scale_prc::prc_abort(), scale_topography::topo_write(), and scale_topography::topo_zsfc.

Referenced by mod_rm_prep::rm_prep().

     use scale_topography, only: &
        topo_write
     implicit none
     !---------------------------------------------------------------------------
 
     if ( mktopo_type == i_ignore ) then
        log_newline
        log_info("MKTOPO",*) 'SKIP  MAKING TOPOGRAPHY DATA'
     else
        log_newline
        log_info("MKTOP",*) '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
           log_error("MKTOPO",*) 'Unsupported TYPE:', mktopo_type
           call prc_abort
        endselect
 
        log_info("MKTOPO",*) 'END   MAKING TOPOGRAPHY DATA'
 
        ! output topography file
        call topo_write
     endif
 
     return