module INITIAL More...

Functions/Subroutines
subroutine, public	mktopo_setup
	Setup.

subroutine, public	mktopo (output, model_mesh, topography)
	Driver.

subroutine, public	mktopo_write (model_mesh, topography)
	Output topography data.

Variables
integer, public	mktopo_type = -1

integer, parameter, public	i_ignore = 0

integer, parameter, public	i_input_file = 1

integer, parameter, public	i_flat = 2

integer, parameter, public	i_bellshape = 3

integer, parameter, public	i_scaer = 4

integer, parameter, public	i_bellshape_global = 5

integer, parameter, public	i_scaher_global = 6

integer, parameter, public	i_barocwave_global_jw2006 = 7

Detailed Description

module INITIAL

Description: A module for preparing topography data

Author: Yuta Kawai, Team SCALE

NAMELIST

PARAM_MKTOPO

name	type	default value	comment
TOPONAME	character(len=H_SHORT)	'NONE'
OUT_BASENAME	character(len=H_LONG)	''	basename of the output file
OUT_VARNAME	character(len=H_MID)	'topo'	variable name of topo in the output file
OUT_TITLE	character(len=H_MID)	'SCALE-DG TOPOGRAPHY'	title of the output file
OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8

PARAM_MKTOPO_INPUT_FILE

name type default value comment

IN_BASENAME character(len=H_LONG)

VARNAME character(len=H_MID)
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_CX	real(RP)	0.0E0_RP	center location [m]: x-coordinate
BELL_CY	real(RP)	0.0E3_RP	center location [m]: y-coordinate
BELL_R	real(RP)	2.E3_RP	half radius [m]
BELL_HEIGHT	real(RP)	100.0_RP	height of mountain [m]
BELL_QUASI_2D	logical	.false.

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)		wavelength of mountain oscillation [m]
SCHAER_HEIGHT	real(RP)		height of mountain [m]

PARAM_MKTOPO_BELLSHAPE_GLOBAL

name	type	default value	comment
BELL_CLON	real(RP)	0.0E0_RP	center location [rad]: longitude
BELL_CLAT	real(RP)	0.0E0_RP	center location [rad]: latitude
BELL_R	real(RP)	2.E3_RP	half radius [m]
BELL_HEIGHT	real(RP)	100.0_RP	height of mountain [m]

PARAM_MKTOPO_SCHAER_GLOBAL

name	type	comment
SCHAER_CLON	real(RP)	center location [rad]: longitude
SCHAER_CLAT	real(RP)	center location [rad]: latitude
SCHAER_R	real(RP)	half-width of mountain [m]
SCHAER_LAMBDA	real(RP)	wavelength of mountain oscillation [m]
SCHAER_SHAPE_ID	integer	1: gaussian, 2: cosine bell
SCHAER_HEIGHT	real(RP)	height of mountain [m]
QUASI_2D_FLAG	logical

PARAM_MKTOPO_BAROCWAVE_GLOBAL_JW2006

name	type	default value	comment
ETA0	real(RP)	0.252_RP	The value of η at a reference level (position of the jet)
U0	real(RP)	35.E0_RP	The parameter associated with zonal jet maximum amplitude [m/s]

History Output: No history output

Function/Subroutine Documentation

◆ mktopo_setup()

subroutine, public mod_mktopo::mktopo_setup

Setup.

Definition at line 91 of file mod_mktopo.F90.

    implicit none
 
    character(len=H_SHORT) :: toponame = 'NONE'
 
    namelist / param_mktopo /  &
      toponame,                  &
      out_basename, out_varname, &
      out_title, out_dtype
 
 
    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(trim(toponame))
    case('NONE')
      mktopo_type = i_ignore
    case('INPUT_FILE')
      mktopo_type = i_input_file
    case('FLAT')
      mktopo_type = i_flat  
    case('BELLSHAPE')
      mktopo_type = i_bellshape
    case('SCHAER')
      mktopo_type = i_scaer
    case('BELLSHAPE_GLOBAL')
      mktopo_type = i_bellshape_global
    case('SCHAER_GLOBAL')
      mktopo_type = i_scaher_global
    case('BAROCWAVE_GLOBAL_JW2006')
      mktopo_type = i_barocwave_global_jw2006      
    case default
      log_error("MKTOPO_setup",*) 'Not appropriate toponame. Check!', toponame
      call prc_abort      
    end select
    
    return

References i_barocwave_global_jw2006, i_bellshape, i_bellshape_global, i_flat, i_ignore, i_input_file, i_scaer, i_scaher_global, and mktopo_type.

Referenced by mod_dg_prep::dg_prep().

◆ mktopo()

subroutine, public mod_mktopo::mktopo	(	logical, intent(out)	output,
		class(atmosmesh), intent(in), target	model_mesh,
		class(meshtopography), intent(inout)	topography )

Driver.

Definition at line 147 of file mod_mktopo.F90.

    use scale_model_var_manager, only: modelvarmanager
    implicit none
 
    logical, intent(out) :: output
    class(AtmosMesh), target, intent(in) :: model_mesh
    class(MeshTopography), intent(inout) :: topography
 
    integer :: n
    integer :: ke
    class(LocalMesh3D), pointer :: lcmesh3D
    class(MeshBase3D), pointer :: mesh
    class(LocalMesh2D), pointer :: lcmesh2D
    class(MeshBase2D), pointer :: mesh2D
    !---------------------------------------------------------------------------
 
    mesh => model_mesh%ptr_mesh
 
    if ( mktopo_type == i_ignore ) then
      log_newline
      log_progress(*) 'skip  making topography data'
      output = .false.
    else
      log_newline
      log_progress(*) 'start making topography data'
 
      call prof_rapstart('_MkTOPO_main', 3)   
      
      call mesh%GetMesh2D( mesh2d )
 
      select case( mktopo_type )
      case ( i_input_file )
        call mktopo_input_file( mesh2d, topography%topo )
      case ( i_flat )
        call mktopo_flat( mesh2d, topography%topo )
      case ( i_bellshape )
        call mktopo_bellshape( mesh2d, topography%topo )
      case ( i_scaer )
        call mktopo_schaer_type_mountain( mesh2d, topography%topo )
      case ( i_bellshape_global )
        call mktopo_bellshape_global( mesh2d, topography%topo )
      case ( i_scaher_global )
        call mktopo_schaer_type_global( mesh2d, topography%topo )        
      case ( i_barocwave_global_jw2006 )
        call mktopo_barocwave_global_jw2006( mesh2d, topography%topo )
      end select
 
      call prof_rapend  ('_MkTOPO_main', 3)
      log_progress(*) 'end making topography data'
 
      output = .true.
    end if
 
    return

References i_barocwave_global_jw2006, i_bellshape, i_bellshape_global, i_flat, i_ignore, i_input_file, i_scaer, i_scaher_global, and mktopo_type.

Referenced by mod_dg_prep::dg_prep().

◆ mktopo_write()

subroutine, public mod_mktopo::mktopo_write	(	class(atmosmesh), intent(in), target	model_mesh,
		class(meshtopography), intent(inout)	topography )

Output topography data.

Definition at line 206 of file mod_mktopo.F90.

 
    use scale_time, only: &
      nowdate   => time_nowdate,   &
      nowsubsec => time_nowsubsec, &
      nowdaysec => time_nowdaysec
    use scale_prc, only: prc_myrank
 
    use scale_file_base_meshfield, only: file_base_meshfield
    use scale_mesh_rectdom2d, only: meshrectdom2d
    use scale_mesh_cubedspheredom2d, only: meshcubedspheredom2d
    use scale_mesh_base2d, only: &
      mftype2d_xy => meshbase2d_dimtypeid_xy
 
    implicit none
    class(AtmosMesh), target, intent(in) :: model_mesh
    class(MeshTopography), intent(inout) :: topography
 
    type(FILE_base_meshfield) :: file
    class(MeshBase2D), pointer :: mesh
 
    logical :: file_existed
    integer, parameter :: vid_topo = 1
    !--------------------------------------------------
   
    nullify( mesh )
    call model_mesh%ptr_mesh%GetMesh2D( mesh )
    select type( mesh )
    class is (meshcubedspheredom2d)
      call file%Init( 1, meshcubedsphere2d=mesh )
    class is (meshrectdom2d)      
      call file%Init( 1, mesh2d=mesh )
    end select
       
    call file%Create( out_basename, out_title, out_dtype, &
                      file_existed,                       &
                      myrank=prc_myrank                   )
 
    if ( .not. file_existed ) then
      call file%Put_GlobalAttribute_time( nowdate, nowsubsec )
    end if
 
    call file%Def_Var( topography%topo, "TOPOGRAPHY", vid_topo, mftype2d_xy, 'XY' )
    call file%End_def()
    call file%Write_var2D( vid_topo, topography%topo, nowdaysec, nowdaysec )
    call file%Close()
    call file%Final()
 
    return