eqs_of_state::dens2g Interface Reference

Public Member Functions
subroutine	dens2_1d (MYS, MYT, MYRHO)
subroutine	dens2_2d (MYS, MYT, MYRHO)

Detailed Description

Definition at line 99 of file eqs_of_state.f90.

Member Function/Subroutine Documentation

dens2_1d()

subroutine eqs_of_state::dens2g::dens2_1d	(	real(sp), dimension(:), intent(in)	MYS,
		real(sp), dimension(:), intent(in)	MYT,
		real(sp), dimension(:), intent(inout)	MYRHO
	)

Definition at line 335 of file eqs_of_state.f90.

    !==============================================================================|
    IMPLICIT NONE
    REAL(SP), INTENT(IN), DIMENSION(:) :: MYS,MYT
    REAL(SP), INTENT(INOUT), DIMENSION(:):: MYRHO

    INTEGER :: I,K,ub1,lb1,ub2,lb2
    !==============================================================================|

    ! SET DIMS USING MY S
    lb1 = lbound(mys,1)
    ub1 = ubound(mys,1)

     ! CHECK MY T
    if(lb1 /= lbound(myt,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myt,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    ! CHECK MY RHO
    if(lb1 /= lbound(myrho,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myrho,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    
    !
    !  CALCULATE DENSITY FROM EQUATION OF STATE
    !
    myrho =mys*mys*mys*6.76786136e-6_sp &
         & - mys*mys*4.8249614e-4_sp &
         & + mys*8.14876577e-1_sp &
         & - 0.22584586e0_sp

    myrho = myrho * &
         & ( myt*myt*myt*1.667e-8_sp &
         & - myt*myt*8.164e-7_sp &
         & + myt*1.803e-5_sp &
         & )
    
    myrho = myrho+1.0_sp &
         & - myt*myt*myt*1.0843e-6_sp &
         & + myt*myt*9.8185e-5_sp &
         & - myt*4.786e-3_sp
    
    myrho = myrho* &
         & ( mys*mys*mys*6.76786136e-6_sp &
         & - mys*mys*4.8249614e-4_sp &
         & + mys*8.14876577e-1_sp &
         & +3.895414e-2_sp &
         & )
    
    myrho = myrho &
         & - (myt-3.98_sp) *(myt-3.98_sp) * (myt+283.0_sp) / (503.57_sp*(myt+67.26_sp))
    
    !
    !  CALCULATE RHO
    !
    myrho =  myrho*1.e-3_sp

Here is the call graph for this function:

dens2_2d()

subroutine eqs_of_state::dens2g::dens2_2d	(	real(sp), dimension(:,:), intent(in)	MYS,
		real(sp), dimension(:,:), intent(in)	MYT,
		real(sp), dimension(:,:), intent(inout)	MYRHO
	)

Definition at line 266 of file eqs_of_state.f90.

    !==============================================================================|
    IMPLICIT NONE
    REAL(SP), INTENT(IN), DIMENSION(:,:) :: MYS,MYT
    REAL(SP), INTENT(INOUT), DIMENSION(:,:):: MYRHO


    INTEGER :: I,K,ub1,lb1,ub2,lb2
    !==============================================================================|

    ! SET DIMS USING MY S
    lb1 = lbound(mys,1)
    ub1 = ubound(mys,1)

    lb2 = lbound(mys,2)
    ub2 = ubound(mys,2)

    ! CHECK MY T
    if(lb1 /= lbound(myt,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myt,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    if(lb2 /= lbound(myt,2)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub2 /= ubound(myt,2)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    ! CHECK MY RHO
    if(lb1 /= lbound(myrho,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myrho,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    if(lb2 /= lbound(myrho,2)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub2 /= ubound(myrho,2)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    
    !
    !  CALCULATE DENSITY FROM EQUATION OF STATE
    !
    myrho =mys*mys*mys*6.76786136e-6_sp &
         & - mys*mys*4.8249614e-4_sp &
         & + mys*8.14876577e-1_sp &
         & - 0.22584586e0_sp

    myrho = myrho * &
         & ( myt*myt*myt*1.667e-8_sp &
         & - myt*myt*8.164e-7_sp &
         & + myt*1.803e-5_sp &
         & )
    
    myrho = myrho+1.0_sp &
         & - myt*myt*myt*1.0843e-6_sp &
         & + myt*myt*9.8185e-5_sp &
         & - myt*4.786e-3_sp
    
    myrho = myrho* &
         & ( mys*mys*mys*6.76786136e-6_sp &
         & - mys*mys*4.8249614e-4_sp &
         & + mys*8.14876577e-1_sp &
         & +3.895414e-2_sp &
         & )
    
    myrho = myrho &
         & - (myt-3.98_sp) *(myt-3.98_sp) * (myt+283.0_sp) / (503.57_sp*(myt+67.26_sp))
    
    !
    !  CALCULATE RHO
    !
    myrho =  myrho*1.e-3_sp

Here is the call graph for this function:

---

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

• FVCOMg01/eqs_of_state.f90