scale_coriolis_param Module Reference

module common / Coriolis parameter More...

Functions/Subroutines
subroutine, public	get_coriolis_parameter (coriolis, colioris_type, np, y, f0, beta, y0, lat)

Detailed Description

module common / Coriolis parameter

Description

Setup coriolis parameter (for regional model)

Reference

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

get_coriolis_parameter()

subroutine, public scale_coriolis_param::get_coriolis_parameter	(	real(rp), dimension(np), intent(out)	coriolis,
		character(*), intent(in)	colioris_type,
		integer, intent(in)	np,
		real(rp), dimension(np), intent(in), optional	y,
		real(rp), intent(in), optional	f0,
		real(rp), intent(in), optional	beta,
		real(rp), intent(in), optional	y0,
		real(rp), dimension(np), intent(in), optional	lat )

Definition at line 38 of file scale_coriolis_param.F90.

 
    implicit none
 
    integer, intent(in) :: Np
    real(RP), intent(out) :: coriolis(Np)
    character(*), intent(in) :: COLIORIS_type
    real(RP), intent(in), optional :: y(Np)
    real(RP), intent(in), optional :: f0
    real(RP), intent(in), optional :: beta
    real(RP), intent(in), optional :: y0
    real(RP), intent(in), optional :: lat(Np)
 
    integer :: i
    !------------------------------------------
 
    if ( trim(colioris_type) == 'PLANE' ) then
      if ( ( .not. present(f0) ) .or. ( .not. present(beta) )    &
            .or. ( .not. present(y) ) .or. ( .not. present(y0) ) ) then
        log_error('get_coriolis_parameter',*) 'If COLIORIS_type is set to PLANE, f0 and beta must be passed. Check!'
        call prc_abort
      end if
      !$omp parallel do
      do i=1, np
        coriolis(i) = f0 + beta * ( y(i) - y0 )
      end do
    else if ( trim(colioris_type) == 'SPHERE' ) then
      if ( .not. present(lat) ) then
        log_error('get_coriolis_parameter',*) 'If COLIORIS_type is set to SPHERE, lat must be passed. Check!'
        call prc_abort
      end if
      !$omp parallel do
      do i=1, np
        coriolis(i) =  2.0_rp * ohm * sin( lat(i) )
      end do
    else if ( trim(colioris_type) == 'NONE' ) then
      !$omp parallel do
      do i=1, np
        coriolis(i) =  0.0_rp
      end do      
    else
      log_error('get_coriolis_parameter',*) 'Unexpected COLIORIS_type is specified. Check! COLIORIS_type=', colioris_type
      call prc_abort
    end if  
 
    return

Referenced by mod_atmos_dyn::atmosdyn_setup().