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)
	Get Coriolis parameter.

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 )

Get Coriolis parameter.

Parameters

[in]	np	Array size
[out]	coriolis	Array storing Coriolis parameter
[in]	colioris_type	Type of Coriolis parameter [PLANE / SPHERE / NONE]
[in]	y	y-coordinate which is used when COLIORIS_type=PLANE
[in]	f0	Value of Coriolis parameter at reference latitude which is used when COLIORIS_type=PLANE
[in]	beta	Rossby parameter at reference latitude which is used when COLIORIS_type=PLANE
[in]	y0	Distance from reference latitude
[in]	lat	Latitude

Definition at line 39 of file scale_coriolis_param.F90.

 
    implicit none
 
    integer, intent(in) :: Np                 !< Array size
    real(RP), intent(out) :: coriolis(Np)     !< Array storing Coriolis parameter
    character(*), intent(in) :: COLIORIS_type !< Type of Coriolis parameter [PLANE / SPHERE / NONE]
    real(RP), intent(in), optional :: y(Np)   !< y-coordinate which is used when COLIORIS_type=PLANE
    real(RP), intent(in), optional :: f0      !< Value of Coriolis parameter at reference latitude which is used when COLIORIS_type=PLANE
    real(RP), intent(in), optional :: beta    !< Rossby parameter at reference latitude which is used when COLIORIS_type=PLANE
    real(RP), intent(in), optional :: y0      !< Distance from reference latitude
    real(RP), intent(in), optional :: lat(Np) !< Latitude
 
    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().