eqs_of_state Module Reference

Data Types
interface	dens2g
interface	fofonoff_millard
interface	jacket_mcdougall

Functions/Subroutines
subroutine	dens1
subroutine	dens2
subroutine	dens3

Function/Subroutine Documentation

dens1()

subroutine eqs_of_state::dens1

(

)

Definition at line 229 of file eqs_of_state.f90.

    !------------------------------------------------------------------------------|
    IMPLICIT NONE
    INTEGER :: K
    REAL(SP), PARAMETER ::PR = 0.0_sp
    REAL(SP), DIMENSION(0:MT,1:KB) :: RZU
    !==============================================================================|
    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "START: DENS1"

    ! The thickness of the water column
    ! Is not the depth realtive to z=0
    DO k=1,kbm1
       rzu(:,k) = -grav_n*1.025_sp*(zz(:,k)*d(:))*0.1_sp
    END DO

    CALL fofonoff_millard_2d(s1,t1,rzu,pr,rho1)
    rho1(:,kb)=0.0_sp
    rho1(0,:)=0.0_sp


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

    CALL n2e3d(rho1,rho)

    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "END: DENS1"
    RETURN

Here is the call graph for this function:

Here is the caller graph for this function:

dens2()

subroutine eqs_of_state::dens2

(

)

Definition at line 393 of file eqs_of_state.f90.

    !==============================================================================|
    IMPLICIT NONE
    !==============================================================================|
    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "START: DENS2"

    !
    !  CALCULATE DENSITY FROM EQUATION OF STATE
    !
    CALL dens2_2d(s1,t1,rho1)
    rho1(:,kb)=0.0_sp
    rho1(0,:)=0.0_sp


    !
    !  AVERAGE FROM NODES TO FACE CENTERS
    !
    CALL n2e3d(rho1,rho)

    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "END: DENS2"
    RETURN

Here is the call graph for this function:

Here is the caller graph for this function:

dens3()

subroutine eqs_of_state::dens3

(

)

Definition at line 688 of file eqs_of_state.f90.

!==============================================================================|
    IMPLICIT NONE
    REAL(SP), DIMENSION(0:MT,KB) :: MYP
    INTEGER :: K
    
    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "START: DENS3"

    DO k=1,kbm1
       myp(:,k) = -grav_n(:)*1.025_sp*zz(:,k)*d(:) *0.01_sp
    END DO

    CALL jacket_mcdougall_2d(s1,t1,myp,rho1)
    rho1(:,kb)=0.0_sp
    rho1(0,:)=0.0_sp

    !
    !  AVERAGE FROM NODES TO FACE CENTERS
    !
    CALL n2e3d(rho1,rho)

    IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "END: DENS3"
    RETURN

Here is the call graph for this function:

Here is the caller graph for this function: