Gets the current value of beta_prime (or equivalent var) Really just gets the variable that we can use to control the pressure gradient.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(out) | :: | bp_out |
subroutine get_beta_prime(bp_out)
use geometry, only: bishop, p_prime_input, invLp_input, beta_prime_input, alpha_input, dp_mult
implicit none
real, intent(out) :: bp_out
select case (bishop)
case(2)
bp_out = p_prime_input
case(3)
bp_out = invLp_input
case(4)
bp_out = beta_prime_input
case(5)
bp_out = alpha_input
case(6)
bp_out = beta_prime_input
case(7)
bp_out = dp_mult
case(8)
bp_out = beta_prime_input
case(9)
bp_out = beta_prime_input
case default
bp_out = 0.0
print*,"ERROR: For bishop not in 2-9 cannot modify the pressure gradient used"
end select
end subroutine get_beta_prime