redi_absolute_position Function

private pure function redi_absolute_position(nk, Pint, Karr, frac) result(p)

Absolute (pressure) position of neutral surface ks in a column (MOM6 absolute_position): Pint(K) + frac*(Pint(K+1)-Pint(K)).

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nk
real(kind=wp), intent(in) :: Pint(nk+1)

interface pressures
integer, intent(in) :: Karr

layer index for this surface
real(kind=wp), intent(in) :: frac

fractional position within the layer

Return Value real(kind=wp)


Called by

proc~~redi_absolute_position~~CalledByGraph proc~redi_absolute_position redi_absolute_position proc~redi_neutral_positions_continuous redi_neutral_positions_continuous proc~redi_neutral_positions_continuous->proc~redi_absolute_position proc~redi_face_coeffs redi_face_coeffs proc~redi_face_coeffs->proc~redi_neutral_positions_continuous proc~redi_calc_coeffs_x redi_calc_coeffs_x proc~redi_calc_coeffs_x->proc~redi_face_coeffs proc~redi_calc_coeffs_y redi_calc_coeffs_y proc~redi_calc_coeffs_y->proc~redi_face_coeffs proc~redi_calc_coeffs redi_calc_coeffs proc~redi_calc_coeffs->proc~redi_calc_coeffs_x proc~redi_calc_coeffs->proc~redi_calc_coeffs_y proc~ocean_dyn_step_split ocean_dyn_step_split proc~ocean_dyn_step_split->proc~redi_calc_coeffs

Source Code

   pure function redi_absolute_position(nk, Pint, Karr, frac) result(p)
      !$acc routine seq
      integer, intent(in) :: nk
      real(wp), intent(in) :: Pint(nk + 1)  !! interface pressures
      integer, intent(in) :: Karr         !! layer index for this surface
      real(wp), intent(in) :: frac        !! fractional position within the layer
      real(wp) :: p

      p = Pint(Karr) + frac*(Pint(Karr + 1) - Pint(Karr))
   end function redi_absolute_position