eqs_of_state::fofonoff_millard Interface Reference

Public Member Functions
subroutine	fofonoff_millard_1d (MYS, MYT, MYP, PREF, MYRHO)
subroutine	fofonoff_millard_2d (MYS, MYT, MYP, PREF, MYRHO)

Detailed Description

Definition at line 94 of file eqs_of_state.f90.

Member Function/Subroutine Documentation

fofonoff_millard_1d()

subroutine eqs_of_state::fofonoff_millard::fofonoff_millard_1d	(	real(sp), dimension(:), intent(in)	MYS,
		real(sp), dimension(:), intent(in)	MYT,
		real(sp), dimension(:), intent(in)	MYP,
		real(sp), intent(in)	PREF,
		real(sp), dimension(:), intent(inout)	MYRHO
	)

Definition at line 188 of file eqs_of_state.f90.

    
    !------------------------------------------------------------------------------|
    IMPLICIT NONE
    REAL(SP), INTENT(IN), DIMENSION(:) :: MYS,MYT,MYP
    REAL(SP), INTENT(IN) :: PREF
    REAL(SP), INTENT(INOUT), DIMENSION(:):: MYRHO
    INTEGER :: I,K,ub1,lb1
    REAL(SP)  ::PT,SVA,SIGMA
    !==============================================================================|

    ! 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 P
    if(lb1 /= lbound(myp,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myp,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!")

    DO k=lb1,ub1
       
       pt  = theta(mys(k),myt(k),myp(k),pref)
       sva = svan(mys(k),pt,myp(k),sigma)
       myrho(k) = sigma*1.e-3_sp
    END DO

    !----------------TRANSFORM TO FACE CENTER--------------------------------------

Here is the call graph for this function:

fofonoff_millard_2d()

subroutine eqs_of_state::fofonoff_millard::fofonoff_millard_2d	(	real(sp), dimension(:,:), intent(in)	MYS,
		real(sp), dimension(:,:), intent(in)	MYT,
		real(sp), dimension(:,:), intent(in)	MYP,
		real(sp), intent(in)	PREF,
		real(sp), dimension(:,:), intent(inout)	MYRHO
	)

Definition at line 134 of file eqs_of_state.f90.

    
    !------------------------------------------------------------------------------|
    IMPLICIT NONE
    REAL(SP), INTENT(IN), DIMENSION(:,:) :: MYS,MYT,MYP
    REAL(SP), INTENT(IN) :: PREF
    REAL(SP), INTENT(INOUT), DIMENSION(:,:):: MYRHO
    INTEGER :: I,K,ub1,lb1,ub2,lb2
    REAL(SP)  ::PT,SVA,SIGMA
    !==============================================================================|

    ! 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 P
    if(lb1 /= lbound(myp,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub1 /= ubound(myp,1)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")

    if(lb2 /= lbound(myp,2)) CALL fatal_error("EQS_OF_STATE:: Dimension mismatch!")
    if(ub2 /= ubound(myp,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!")

    

    DO i=lb1,ub1
       DO k=lb2,ub2  

          pt  = theta(mys(i,k),myt(i,k),myp(i,k),pref)
          sva = svan(mys(i,k),pt,myp(i,k),sigma)
          myrho(i,k) = sigma*1.e-3_sp
       END DO
    END DO

    !----------------TRANSFORM TO FACE CENTER--------------------------------------

Here is the call graph for this function:

---

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

• FVCOMg01/eqs_of_state.f90