mod_mktopo_util.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module mod_mktopo_util
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prc
 
  use scale_const, only: &
    pi => const_pi
 
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: mktopoutil_barocwave_global_jw2006_calc_topo
 
contains
 
  !-----------------------------------------------------------------------------
!OCL SERIAL
  subroutine mktopoutil_barocwave_global_jw2006_calc_topo( topo, &
    U0, ETA0, lat, Np )
 
    use scale_const, only: &
      ohm => const_ohm,      &
      grav => const_grav,    &
      rplanet => const_radius
 
    implicit none
 
    integer, intent(in) :: np
    real(rp), intent(out) :: topo(np)
    real(rp), intent(in) :: u0
    real(rp), intent(in) :: eta0
    real(rp), intent(in) :: lat(np)
 
    real(rp) :: sin_lat(np)
    real(rp) :: sin_lat_pow_6(np)
    real(rp) :: cos_lat(np)
 
    real(rp ) :: tmp
    !-------------------------------------------
 
    sin_lat(:) = sin(lat(:))
    sin_lat_pow_6(:) = sin_lat(:)**6
    cos_lat(:) = cos(lat(:))
 
    tmp = cos( (1.0_rp - eta0) * 0.5_rp * pi )
    tmp = u0 * tmp * sqrt(tmp)
 
    ! Calc horizontal variation of geopotential height
    topo(:) = tmp * & 
      (   tmp * ( - 2.0_rp * sin_lat_pow_6(:) * ( cos_lat(:)**2 + 1.0_rp / 3.0_rp ) + 10.0_rp / 63.0_rp )          &
        + rplanet * ohm * ( 8.0_rp / 5.0_rp * cos_lat(:)**3 * ( sin_lat(:)**2 + 2.0_rp / 3.0_rp ) - 0.25_rp * pi ) &
      ) / grav
 
    return
  end subroutine mktopoutil_barocwave_global_jw2006_calc_topo
 
end module mod_mktopo_util