soliton_initialization Module Reference

Detailed Description

Initial conditions for the Equatorial Rossby soliton test (Boyd).

Description of the equatorial Rossby soliton initial

conditions

Functions/Subroutines
subroutine, public	soliton_initialize_thickness (h, G, GV, US)
	Initialization of thicknesses in Equatorial Rossby soliton test. More...
subroutine, public	soliton_initialize_velocity (u, v, h, G)
	Initialization of u and v in the equatorial Rossby soliton test. More...

Variables
character(len=40)	mdl = "soliton_initialization"
	This module's name.

Function/Subroutine Documentation

soliton_initialize_thickness()

subroutine, public soliton_initialization::soliton_initialize_thickness	(	real, dimension(szi_(g),szj_(g),szk_(gv)), intent(out)	h,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US
	)

Initialization of thicknesses in Equatorial Rossby soliton test.

Parameters

[in]	g	The ocean's grid structure.
[in]	gv	The ocean's vertical grid structure.
[in]	us	A dimensional unit scaling type
[out]	h	The thickness that is being initialized [H ~> m or kg m-2].

Definition at line 32 of file soliton_initialization.F90.

  type(ocean_grid_type),   intent(in)  :: G    !< The ocean's grid structure.
  type(verticalGrid_type), intent(in)  :: GV   !< The ocean's vertical grid structure.
  type(unit_scale_type),   intent(in)  :: US   !< A dimensional unit scaling type
  real, dimension(SZI_(G),SZJ_(G),SZK_(GV)), &
                           intent(out) :: h    !< The thickness that is being initialized [H ~> m or kg m-2].
 
  integer :: i, j, k, is, ie, js, je, nz
  real    :: x, y, x0, y0
  real    :: val1, val2, val3, val4
  character(len=40) :: verticalCoordinate
 
  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
 
  call mom_mesg("soliton_initialization.F90, soliton_initialize_thickness: setting thickness")
 
  x0 = 2.0*g%len_lon/3.0
  y0 = 0.0
  val1 = 0.395
  val2 = us%m_to_Z * 0.771*(val1*val1)
 
  do j = g%jsc,g%jec ; do i = g%isc,g%iec
    do k = 1, nz
      x = g%geoLonT(i,j)-x0
      y = g%geoLatT(i,j)-y0
      val3 = exp(-val1*x)
      val4 = val2 * ( 2.0*val3 / (1.0 + (val3*val3)) )**2
      h(i,j,k) = gv%Z_to_H * (0.25*val4*(6.0*y*y + 3.0) * exp(-0.5*y*y) + g%bathyT(i,j))
    enddo
  enddo ; enddo
 

soliton_initialize_velocity()

subroutine, public soliton_initialization::soliton_initialize_velocity	(	real, dimension(szib_(g),szj_(g),szk_(g)), intent(out)	u,
		real, dimension(szi_(g),szjb_(g),szk_(g)), intent(out)	v,
		real, dimension(szi_(g),szj_(g), szk_(g)), intent(in)	h,
		type(ocean_grid_type), intent(in)	G
	)

Initialization of u and v in the equatorial Rossby soliton test.

Parameters

[in]	g	Grid structure
[out]	u	i-component of velocity [m s-1]
[out]	v	j-component of velocity [m s-1]
[in]	h	Thickness [H ~> m or kg m-2]

Definition at line 67 of file soliton_initialization.F90.

  type(ocean_grid_type),                  intent(in)     :: G  !< Grid structure
  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)), intent(out) :: u  !< i-component of velocity [m s-1]
  real, dimension(SZI_(G),SZJB_(G),SZK_(G)), intent(out) :: v  !< j-component of velocity [m s-1]
  real, dimension(SZI_(G),SZJ_(G), SZK_(G)), intent(in)  :: h  !< Thickness [H ~> m or kg m-2]
 
  real    :: x, y, x0, y0
  real    :: val1, val2, val3, val4
  integer :: i, j, k, is, ie, js, je, nz
 
  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
 
  x0 = 2.0*g%len_lon/3.0
  y0 = 0.0
  val1 = 0.395
  val2 = 0.771*(val1*val1)
 
  v(:,:,:) = 0.0
  u(:,:,:) = 0.0
 
  do j = g%jsc,g%jec ; do i = g%isc-1,g%iec+1
    do k = 1, nz
      x = 0.5*(g%geoLonT(i+1,j)+g%geoLonT(i,j))-x0
      y = 0.5*(g%geoLatT(i+1,j)+g%geoLatT(i,j))-y0
      val3 = exp(-val1*x)
      val4 = val2*((2.0*val3/(1.0+(val3*val3)))**2)
      u(i,j,k) = 0.25*val4*(6.0*y*y-9.0) * exp(-0.5*y*y)
    enddo
  enddo ; enddo
  do j = g%jsc-1,g%jec+1 ; do i = g%isc,g%iec
    do k = 1, nz
      x = 0.5*(g%geoLonT(i,j+1)+g%geoLonT(i,j))-x0
      y = 0.5*(g%geoLatT(i,j+1)+g%geoLatT(i,j))-y0
      val3 = exp(-val1*x)
      val4 = val2*((2.0*val3/(1.0+(val3*val3)))**2)
      v(i,j,k) = 2.0*val4*y*(-2.0*val1*tanh(val1*x)) * exp(-0.5*y*y)
    enddo
  enddo ; enddo