calc_gcc_seg_pa Module Function

pure elemental module function calc_gcc_seg_pa(self, state) result(suction)

Arguments

Type IntentOptional Attributes Name
class(type_gcc_segregation_pa), intent(in) :: self
type(type_state), intent(in) :: state

Return Value real(kind=real64)


Called by

proc~~calc_gcc_seg_pa~~CalledByGraph proc~calc_gcc_seg_pa calc_gcc_seg_pa interface~calc_gcc_seg_pa type_gcc_segregation_pa%calc_gcc_seg_pa interface~calc_gcc_seg_pa->proc~calc_gcc_seg_pa

Source Code

    module pure elemental function calc_gcc_seg_pa(self, state) result(suction)
        implicit none
        class(type_gcc_segregation_pa), intent(in) :: self
        type(type_state), intent(in) :: state
        real(real64) :: suction

        if (state%temperature <= self%Tf) then
            suction = (state%density_ice / state%density_water - 1.0d0) * state%pressure &
                      - self%Lf * state%density_ice * log((state%temperature + self%TtoK) / (self%Tf + self%TtoK))
        else
            suction = 0.0d0
        end if

    end function calc_gcc_seg_pa