mod_copytopo Module Reference

module Copy topography More...

Functions/Subroutines
subroutine, public	copytopo (topo_cd)
	Setup and Main. More...

Variables
integer, public	cnvtopo_type = -1
integer, parameter, public	i_ignore = 0
integer, parameter, public	i_gtopo30 = 1
integer, parameter, public	i_dem50m = 2
integer, parameter, public	i_gmted2010 = 3

Detailed Description

module Copy topography

Description

subroutines for preparing topography data (copy from parent domain)

Author

Team SCALE

NAMELIST

• PARAM_COPYTOPO

  name	type	default value	comment
  COPYTOPO_IN_BASENAME	character(len=H_LONG)	''
  COPYTOPO_TRANSITION_DX	real(RP)	-1.0_RP	thickness of transition region [m]: x
  COPYTOPO_TRANSITION_DY	real(RP)	-1.0_RP	thickness of transition region [m]: y
  COPYTOPO_TRANSFACT	real(RP)	-1.0_RP	stretch factor of transition region
  COPYTOPO_FRACX	real(RP)	1.0_RP	fraction of transition region (x) (0-1)
  COPYTOPO_FRACY	real(RP)	1.0_RP	fraction of transition region (y) (0-1)
  COPYTOPO_TAUX	real(RP)	1.0_RP	maximum value for mixing tau (x) [s]
  COPYTOPO_TAUY	real(RP)	1.0_RP	maximum value for mixing tau (y) [s]
  COPYTOPO_ENTIRE_REGION	logical	.false.	copy parent topo over an entire region
  COPYTOPO_LINEAR_H	logical	.true.	linear or non-linear profile of relax region
  COPYTOPO_EXP_H	real(RP)	2.0_RP	factor of non-linear profile of relax region

  
  

History Output

No history output

Function/Subroutine Documentation

copytopo()

subroutine, public mod_copytopo::copytopo

(

real(rp), dimension(:,:), intent(inout)

topo_cd

)

Setup and Main.

Definition at line 84 of file mod_copytopo.f90.

References scale_grid::buffer_dx, scale_grid::buffer_dy, scale_grid::bufffact, scale_const::const_d2r, scale_const::const_eps, scale_const::const_radius, scale_grid::dx, scale_grid::dy, scale_grid::grid_cbfxg, scale_grid::grid_cbfyg, scale_grid::grid_cxg, scale_grid::grid_cyg, scale_grid::grid_fxg, scale_grid::grid_fyg, scale_grid_index::ia, scale_grid_index::iag, scale_grid_index::ihalo, scale_grid_index::imax, scale_interpolation_nest::intrpnest_domain_compatibility(), scale_interpolation_nest::intrpnest_interp_2d, scale_interpolation_nest::intrpnest_interp_fact_latlon(), 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::ja, scale_grid_index::jag, scale_grid_index::jhalo, scale_grid_index::jmax, scale_grid_nest::nest_domain_shape(), scale_grid_nest::nest_interp_level, scale_grid_nest::nest_tile_id, scale_grid_nest::nest_tile_num_x, scale_grid_nest::nest_tile_num_y, scale_grid_nest::parent_imax, scale_grid_nest::parent_jmax, scale_rm_process::prc_2drank, scale_process::prc_mpistop(), scale_process::prc_myrank, scale_rm_process::prc_num_x, scale_rm_process::prc_num_y, scale_grid_real::real_lat, and scale_grid_real::real_lon.

Referenced by mod_cnvtopo::cnvtopo().

    use scale_process, only: &
       prc_mpistop
    use scale_grid, only: &
       dx,         &
       dy,         &
       buffer_dx,  &
       buffer_dy,  &
       bufffact
    use scale_grid_nest, only: &
       nest_interp_level
    implicit none

    real(RP), intent(inout) :: topo_cd(:,:)

    namelist / param_copytopo / &
       copytopo_in_basename,    &
       copytopo_transition_dx,  &
       copytopo_transition_dy,  &
       copytopo_transfact,      &
       copytopo_fracx,          &
       copytopo_fracy,          &
       copytopo_taux,           &
       copytopo_tauy,           &
       copytopo_entire_region,  &
       copytopo_linear_h,       &
       copytopo_exp_h

    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ Module[COPYTOPO]/Categ[COPYTOPO]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_copytopo,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_COPYTOPO. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_copytopo)

    if ( copytopo_transition_dx < 0.0_rp ) then
       copytopo_transition_dx = buffer_dx
    end if
    if ( copytopo_transition_dy < 0.0_rp ) then
       copytopo_transition_dy = buffer_dy
    end if
    if ( copytopo_transfact < 0.0_rp ) then
       copytopo_transfact = bufffact
    end if

    allocate( ctrx(ia) )
    allocate( ctry(ja) )
    allocate( copytopo_alpha(ia,ja) )
    allocate( topo_pd(ia,ja) )
    copytopo_alpha(:,:) = 0.0_rp

    itp_nh = int( nest_interp_level )

    ! copy topography from parent domain to transition region
    call copytopo_transgrid

    call copytopo_setalpha

    call copytopo_input_data( topo_pd ) ! (out)

    call copytopo_mix_data( topo_cd, &  ! (inout)
                            topo_pd  )  ! (in)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

cnvtopo_type

integer, public mod_copytopo::cnvtopo_type = -1

Definition at line 36 of file mod_copytopo.f90.

  integer, public            :: CNVTOPO_TYPE = -1

i_ignore

integer, parameter, public mod_copytopo::i_ignore = 0

Definition at line 38 of file mod_copytopo.f90.

  integer, public, parameter :: I_IGNORE     =  0

i_gtopo30

integer, parameter, public mod_copytopo::i_gtopo30 = 1

Definition at line 39 of file mod_copytopo.f90.

  integer, public, parameter :: I_GTOPO30    =  1

i_dem50m

integer, parameter, public mod_copytopo::i_dem50m = 2

Definition at line 40 of file mod_copytopo.f90.

  integer, public, parameter :: I_DEM50M     =  2

i_gmted2010

integer, parameter, public mod_copytopo::i_gmted2010 = 3

Definition at line 41 of file mod_copytopo.f90.

  integer, public, parameter :: I_GMTED2010  =  3