module TOPOGRAPHY More...

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

subroutine, public	topo_fillhalo (Zsfc, FILL_BND)
	HALO Communication. More...

subroutine, public	topo_write
	Write topography. More...

Variables
logical, public	topo_exist = .false.
	topography exists? More...

real(rp), dimension(:,:), allocatable, public	topo_zsfc
	absolute ground height [m] More...

Detailed Description

module TOPOGRAPHY

Description: Topography module

Author: Team SCALE

NAMELIST

PARAM_TOPO

name	type	default value	comment
TOPO_IN_BASENAME	character(len=H_LONG)	''	basename of the input file
TOPO_IN_AGGREGATE	logical		> switch to use aggregated file
TOPO_IN_CHECK_COORDINATES	logical	.false.	> switch for check of coordinates
TOPO_OUT_BASENAME	character(len=H_LONG)	''	basename of the output file
TOPO_OUT_AGGREGATE	logical		> switch to use aggregated file
TOPO_OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8

History Output: No history output

Function/Subroutine Documentation

◆ topo_setup()

subroutine, public scale_topography::topo_setup ( )

Setup.

Definition at line 62 of file scale_topography.F90.

References scale_file::file_aggregate, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_prc::prc_abort(), and topo_zsfc.

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

     use scale_file, only: &
        file_aggregate
     use scale_prc, only: &
        prc_abort
     implicit none
 
     namelist / param_topo / &
        topo_in_basename,          &
        topo_in_aggregate,         &
        topo_in_check_coordinates, &
        topo_out_basename,         &
        topo_out_aggregate,        &
        topo_out_dtype
 
     integer :: ierr
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("TOPO_setup",*) 'Setup'
 
     topo_in_aggregate  = file_aggregate
     topo_out_aggregate = file_aggregate
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_topo,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("TOPO_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("TOPO_setup",*) 'Not appropriate names in namelist PARAM_TOPO. Check!'
        call prc_abort
     endif
     log_nml(param_topo)
 
     allocate( topo_zsfc(ia,ja) )
     topo_zsfc(:,:) = 0.0_rp
 
     ! read from file
     call topo_read
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ topo_fillhalo()

subroutine, public scale_topography::topo_fillhalo	(	real(rp), dimension(ia,ja), intent(inout), optional	Zsfc,
		logical, intent(in), optional	FILL_BND
	)

HALO Communication.

Definition at line 108 of file scale_topography.F90.

References scale_file_cartesc::file_cartesc_close(), scale_file_cartesc::file_cartesc_flush(), scale_file_cartesc::file_cartesc_open(), scale_prc::prc_abort(), topo_exist, and topo_zsfc.

Referenced by mod_cnvtopo::cnvtopo(), and topo_write().

     use scale_comm_cartesc, only: &
        comm_vars8, &
        comm_wait
     implicit none
 
     real(RP), intent(inout), optional :: zsfc(ia,ja)
     logical,  intent(in),    optional :: fill_bnd
 
     logical :: fill_bnd_
     !---------------------------------------------------------------------------
 
     fill_bnd_ = .false.
     if ( present(fill_bnd) ) fill_bnd_ = fill_bnd
 
     if ( present(zsfc) ) then
        call comm_vars8( zsfc(:,:), 1 )
        call comm_wait ( zsfc(:,:), 1, fill_bnd_ )
     else
        call comm_vars8( topo_zsfc(:,:), 1 )
        call comm_wait ( topo_zsfc(:,:), 1, fill_bnd_ )
     end if
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ topo_write()

subroutine, public scale_topography::topo_write ( )

Write topography.

Definition at line 181 of file scale_topography.F90.

References topo_fillhalo(), and topo_zsfc.

Referenced by mod_cnvtopo::cnvtopo(), and mod_mktopo::mktopo().

     use scale_file_cartesc, only: &
        file_cartesc_write
     implicit none
     !---------------------------------------------------------------------------
 
     if ( topo_out_basename /= '' ) then
 
        log_newline
        log_info("TOPO_write",*) 'Output topography file '
 
        call topo_fillhalo( fill_bnd=.false. )
 
        call file_cartesc_write( topo_zsfc(:,:), topo_out_basename, topo_out_title, & ! [IN]
                                 'TOPO', 'Topography', 'm', 'XY',   topo_out_dtype, & ! [IN]
                                 standard_name="surface_altitude",                  & ! [IN]
                                 haszcoord=.false., aggregate=topo_out_aggregate    ) ! [IN]
 
     endif
 
     return