mom_eos_nemo::calculate_density_nemo Interface Reference

Detailed Description

Definition at line 41 of file MOM_EOS_NEMO.F90.

Private functions
subroutine	calculate_density_scalar_nemo (T, S, pressure, rho)
subroutine	calculate_density_array_nemo (T, S, pressure, rho, start, npts)
	This subroutine computes the in situ density of sea water (rho in units of kg/m^3) from absolute salinity (S in g/Kg), conservative temperature (T in deg C), and pressure in Pa. More...

Functions and subroutines

calculate_density_array_nemo()

subroutine mom_eos_nemo::calculate_density_nemo::calculate_density_array_nemo ( real, dimension(:), intent(in)  T, real, dimension(:), intent(in)  S, real, dimension(:), intent(in)  pressure, real, dimension(:), intent(out)  rho, integer, intent(in)  start, integer, intent(in)  npts  )

private

This subroutine computes the in situ density of sea water (rho in units of kg/m^3) from absolute salinity (S in g/Kg), conservative temperature (T in deg C), and pressure in Pa.

Parameters

[in]	t	Conservative temperature in C.
[in]	s	Absoulte salinity in g/Kg.
[in]	pressure	Pressure in Pa.
[out]	rho	In situ density in kg m-3.
[in]	start	The starting point in the arrays.
[in]	npts	The number of values to calculate.

Definition at line 213 of file MOM_EOS_NEMO.F90.

  real,    intent(in),  dimension(:) :: t     !< Conservative temperature in C.
  real,    intent(in),  dimension(:) :: s     !< Absoulte salinity in g/Kg.
  real,    intent(in),  dimension(:) :: pressure !< Pressure in Pa.
  real,    intent(out), dimension(:) :: rho   !< In situ density in kg m-3.
  integer, intent(in)                :: start !< The starting point in the arrays.
  integer, intent(in)                :: npts  !< The number of values to calculate.

! * Arguments: T - conservative temperature in C.                      *
! *  (in)      S - absoulte salinity in g/Kg.                          *
! *  (in)      pressure - pressure in Pa.                              *
! *  (out)     rho - in situ density in kg m-3.                        *
! *  (in)      start - the starting point in the arrays.               *
! *  (in)      npts - the number of values to calculate.               *

! *====================================================================*
! *  This subroutine computes the in situ density of sea water (rho in *
! *  units of kg/m^3) from absolute salinity (S in g/Kg),              *
! *  conservative temperature (T in deg C), and pressure in Pa.        *
! *====================================================================*
  real :: zp,zt , zh , zs , zr0, zn , zn0, zn1, zn2, zn3
  integer :: j

  do j=start,start+npts-1
    !Conversions
    zs = s(j) !gsw_sr_from_sp(S(j))       !Convert practical salinity to absolute salinity
    zt = t(j) !gsw_ct_from_pt(S(j),T(j))  !Convert potantial temp to conservative temp
    zp = pressure(j)* pa2db         !Convert pressure from Pascal to decibar

    !The following algorithm was provided by Roquet in a private communication.
    !It is not necessarily the algorithm used in NEMO ocean!
    zp  = zp * r1_p0 !pressure
    zt  = zt * r1_t0 !temperature
    zs  = sqrt( abs( zs + rdeltas ) * r1_s0 )   ! square root salinity
    !
    zn3 = eos013*zt   &
       &   + eos103*zs+eos003
       !
    zn2 = (eos022*zt   &
       &   + eos112*zs+eos012)*zt   &
       &   + (eos202*zs+eos102)*zs+eos002
       !
    zn1 = (((eos041*zt   &
       &   + eos131*zs+eos031)*zt   &
       &   + (eos221*zs+eos121)*zs+eos021)*zt   &
       &   + ((eos311*zs+eos211)*zs+eos111)*zs+eos011)*zt   &
       &   + (((eos401*zs+eos301)*zs+eos201)*zs+eos101)*zs+eos001
       !
    zn0 = (((((eos060*zt   &
       &   + eos150*zs+eos050)*zt   &
       &   + (eos240*zs+eos140)*zs+eos040)*zt   &
       &   + ((eos330*zs+eos230)*zs+eos130)*zs+eos030)*zt   &
       &   + (((eos420*zs+eos320)*zs+eos220)*zs+eos120)*zs+eos020)*zt   &
       &   + ((((eos510*zs+eos410)*zs+eos310)*zs+eos210)*zs+eos110)*zs+eos010)*zt   &
       &   + (((((eos600*zs+eos500)*zs+eos400)*zs+eos300)*zs+eos200)*zs+eos100)*zs+eos000
       !
    zn  = ( ( zn3 * zp + zn2 ) * zp + zn1 ) * zp + zn0
    !
    zr0 = (((((r05 * zp+r04) * zp+r03 ) * zp+r02 ) * zp+r01) * zp+r00) * zp
    !
    rho(j) =  ( zn + zr0 ) ! density

 enddo

calculate_density_scalar_nemo()

subroutine mom_eos_nemo::calculate_density_nemo::calculate_density_scalar_nemo ( real, intent(in)  T, real, intent(in)  S, real, intent(in)  pressure, real, intent(out)  rho  )

private

Definition at line 184 of file MOM_EOS_NEMO.F90.

real,    intent(in)  :: t, s, pressure
real,    intent(out) :: rho
! * Arguments: T - conservative temperature in C.                      *
! *  (in)      S - absoulte salinity in g/Kg.                          *
! *  (in)      pressure - pressure in Pa.                              *
! *  (out)     rho - in situ density in kg m-3.                        *

! *====================================================================*
! *  This subroutine computes the in situ density of sea water (rho in *
! *  units of kg/m^3) from absolute salinity (S in g/Kg), conservative *
! *  temperature (T in deg C), and pressure in Pa.                     *
! *====================================================================*

  real :: al0, p0, lambda
  integer :: j
  real, dimension(1) :: t0, s0, pressure0
  real, dimension(1) :: rho0

  t0(1) = t
  s0(1) = s
  pressure0(1) = pressure

  call calculate_density_array_nemo(t0, s0, pressure0, rho0, 1, 1)
  rho = rho0(1)

---

The documentation for this interface was generated from the following file:

• /home/adcroft/GitHub/workspace/Gaea-stats-MOM6-examples/MOM6-examples/src/MOM6/src/equation_of_state/MOM_EOS_NEMO.F90