pressure_lower Function

  1. COESA.f90
  2. COESA_module
  3. pressure_lower

private pure function pressure_lower(H)

Arguments

TypeIntentOptionalAttributesName
real(kind=wp), intent(in) :: H

Return Value real(kind=wp)


Calls

proc~~pressure_lower~~CallsGraph proc~pressure_lower pressure_lower proc~find find proc~pressure_lower->proc~find

Called by

proc~~pressure_lower~~CalledByGraph proc~pressure_lower pressure_lower proc~coesa_atmosphere COESA_atmosphere proc~coesa_atmosphere->proc~pressure_lower proc~coesa_density COESA_density proc~coesa_density->proc~pressure_lower

Contents

Source Code

pressure_lower

Source Code

pure real(wp) function pressure_lower(H)
    real(wp),intent(in) :: H
    integer :: i
    i = find(H, Hb)
    if (Lmb(i) == 0) then
        pressure_lower = Pb(i) * exp(-g0 * M0 * (H - Hb(i)) / (Rstar * Tmb(i)))
    else
        pressure_lower = Pb(i) * (Tmb(i) / (Tmb(i) + Lmb(i) * (H - Hb(i)))) ** (g0 * M0 / (Rstar * Lmb(i)))
    end if
end function pressure_lower