coord_zlike Module Reference

Detailed Description

Regrid columns for a z-like coordinate (z-star, z-level)

Data Types
type	zlike_cs
	Control structure containing required parameters for a z-like coordinate. More...

Functions/Subroutines
subroutine, public	init_coord_zlike (CS, nk, coordinateResolution)
	Initialise a zlike_CS with pointers to parameters. More...
subroutine, public	end_coord_zlike (CS)
	Deallocates the zlike control structure. More...
subroutine, public	set_zlike_params (CS, min_thickness)
	Set parameters in the zlike structure. More...
subroutine, public	build_zstar_column (CS, depth, total_thickness, zInterface, z_rigid_top, eta_orig, zScale)
	Builds a z* coordinate with a minimum thickness. More...

Function/Subroutine Documentation

build_zstar_column()

subroutine, public coord_zlike::build_zstar_column	(	type(zlike_cs), intent(in)	CS,
		real, intent(in)	depth,
		real, intent(in)	total_thickness,
		real, dimension(cs%nk+1), intent(inout)	zInterface,
		real, intent(in), optional	z_rigid_top,
		real, intent(in), optional	eta_orig,
		real, intent(in), optional	zScale
	)

Builds a z* coordinate with a minimum thickness.

Parameters

[in]	cs	Coordinate control structure
[in]	depth	Depth of ocean bottom (positive in the output units)
[in]	total_thickness	Column thickness (positive in the same units as depth)
[in,out]	zinterface	Absolute positions of interfaces
[in]	z_rigid_top	The height of a rigid top (negative in the same units as depth)
[in]	eta_orig	The actual original height of the top in the same units as depth
[in]	zscale	Scaling factor from the target coordinate resolution in Z to desired units for zInterface, perhaps Z_to_H

Definition at line 65 of file coord_zlike.F90.

  type(zlike_CS),           intent(in)    :: CS !< Coordinate control structure
  real,                     intent(in)    :: depth !< Depth of ocean bottom (positive in the output units)
  real,                     intent(in)    :: total_thickness !< Column thickness (positive in the same units as depth)
  real, dimension(CS%nk+1), intent(inout) :: zInterface !< Absolute positions of interfaces
  real, optional,           intent(in)    :: z_rigid_top !< The height of a rigid top (negative in the
                                                         !! same units as depth)
  real, optional,           intent(in)    :: eta_orig !< The actual original height of the top in the
                                                      !! same units as depth
  real, optional,           intent(in)    :: zScale !< Scaling factor from the target coordinate resolution
                                                    !! in Z to desired units for zInterface, perhaps Z_to_H
  ! Local variables
  real :: eta, stretching, dh, min_thickness, z0_top, z_star, z_scale
  integer :: k
  logical :: new_zstar_def
 
  z_scale = 1.0 ; if (present(zscale)) z_scale = zscale
 
  new_zstar_def = .false.
  min_thickness = min( cs%min_thickness, total_thickness/real(cs%nk) )
  z0_top = 0.
  if (present(z_rigid_top)) then
    z0_top = z_rigid_top
    new_zstar_def = .true.
  endif
 
  ! Position of free-surface (or the rigid top, for which eta ~ z0_top)
  eta = total_thickness - depth
  if (present(eta_orig)) eta = eta_orig
 
  ! Conventional z* coordinate:
  !   z* = (z-eta) / stretching   where stretching = (H+eta)/H
  !   z = eta + stretching * z*
  ! The above gives z*(z=eta) = 0, z*(z=-H) = -H.
  ! With a rigid top boundary at eta = z0_top then
  !   z* = z0 + (z-eta) / stretching   where stretching = (H+eta)/(H+z0)
  !   z = eta + stretching * (z*-z0) * stretching
  stretching = total_thickness / ( depth + z0_top )
 
  if (new_zstar_def) then
    ! z_star is the notional z* coordinate in absence of upper/lower topography
    z_star = 0. ! z*=0 at the free-surface
    zinterface(1) = eta ! The actual position of the top of the column
    do k = 2,cs%nk
      z_star = z_star - cs%coordinateResolution(k-1)*z_scale
      ! This ensures that z is below a rigid upper surface (ice shelf bottom)
      zinterface(k) = min( eta + stretching * ( z_star - z0_top ), z0_top )
      ! This ensures that the layer in inflated
      zinterface(k) = min( zinterface(k), zinterface(k-1) - min_thickness )
      ! This ensures that z is above or at the topography
      zinterface(k) = max( zinterface(k), -depth + real(cs%nk+1-k) * min_thickness )
    enddo
    zinterface(cs%nk+1) = -depth
 
  else
    ! Integrate down from the top for a notional new grid, ignoring topography
    ! The starting position is offset by z0_top which, if z0_top<0, will place
    ! interfaces above the rigid boundary.
    zinterface(1) = eta
    do k = 1,cs%nk
      dh = stretching * cs%coordinateResolution(k)*z_scale ! Notional grid spacing
      zinterface(k+1) = zinterface(k) - dh
    enddo
 
    ! Integrating up from the bottom adjusting interface position to accommodate
    ! inflating layers without disturbing the interface above
    zinterface(cs%nk+1) = -depth
    do k = cs%nk,1,-1
      if ( zinterface(k) < (zinterface(k+1) + min_thickness) ) then
        zinterface(k) = zinterface(k+1) + min_thickness
      endif
    enddo
  endif
 

Referenced by mom_regridding::build_zstar_grid(), and mom_diag_remap::diag_remap_update().

Here is the caller graph for this function:

end_coord_zlike()

subroutine, public coord_zlike::end_coord_zlike

(

type(zlike_cs), pointer

CS

)

Deallocates the zlike control structure.

Parameters

cs	Coordinate control structure

Definition at line 44 of file coord_zlike.F90.

  type(zlike_CS), pointer :: CS !< Coordinate control structure
 
  ! Nothing to do
  if (.not. associated(cs)) return
  deallocate(cs%coordinateResolution)
  deallocate(cs)

Referenced by mom_regridding::end_regridding().

Here is the caller graph for this function:

init_coord_zlike()

subroutine, public coord_zlike::init_coord_zlike	(	type(zlike_cs), pointer	CS,
		integer, intent(in)	nk,
		real, dimension(:), intent(in)	coordinateResolution
	)

Initialise a zlike_CS with pointers to parameters.

Parameters

	cs	Unassociated pointer to hold the control structure
[in]	nk	Number of levels in the grid
[in]	coordinateresolution	Target coordinate resolution [Z ~> m]

Definition at line 30 of file coord_zlike.F90.

  type(zlike_CS),     pointer    :: CS !< Unassociated pointer to hold the control structure
  integer,            intent(in) :: nk !< Number of levels in the grid
  real, dimension(:), intent(in) :: coordinateResolution !< Target coordinate resolution [Z ~> m]
 
  if (associated(cs)) call mom_error(fatal, "init_coord_zlike: CS already associated!")
  allocate(cs)
  allocate(cs%coordinateResolution(nk))
 
  cs%nk                   = nk
  cs%coordinateResolution = coordinateresolution

References mom_error_handler::mom_error().

Referenced by mom_regridding::initcoord().

Here is the call graph for this function:

Here is the caller graph for this function:

set_zlike_params()

subroutine, public coord_zlike::set_zlike_params	(	type(zlike_cs), pointer	CS,
		real, intent(in), optional	min_thickness
	)

Set parameters in the zlike structure.

Parameters

	cs	Coordinate control structure
[in]	min_thickness	Minimum allowed thickness [H ~> m or kg m-2]

Definition at line 54 of file coord_zlike.F90.

  type(zlike_CS), pointer    :: CS !< Coordinate control structure
  real, optional, intent(in) :: min_thickness !< Minimum allowed thickness [H ~> m or kg m-2]
 
  if (.not. associated(cs)) call mom_error(fatal, "set_zlike_params: CS not associated")
 
  if (present(min_thickness)) cs%min_thickness = min_thickness

References mom_error_handler::mom_error().

Referenced by mom_regridding::set_regrid_params().

Here is the call graph for this function:

Here is the caller graph for this function: