Compute exospheric temperature from F10.7 values using the Jacchia-Roberts formula
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(jacchia_roberts_type), | intent(in) | :: | me | |||
| real(kind=dp), | intent(in) | :: | f107_daily |
Daily F10.7 value |
||
| real(kind=dp), | intent(in) | :: | f107a_avg |
81-day average F10.7a |
||
| real(kind=dp), | intent(in) | :: | tkp |
Kp index |
Exospheric temperature (K)
pure function parms_to_temp(me, f107_daily, f107a_avg, tkp) result(xtemp) class(jacchia_roberts_type), intent(in) :: me real(dp), intent(in) :: f107_daily !! Daily F10.7 value real(dp), intent(in) :: f107a_avg !! 81-day average F10.7a real(dp), intent(in) :: tkp !! Kp index real(dp) :: xtemp !! Exospheric temperature (K) xtemp = 379.0_dp + 3.24_dp * f107a_avg + 1.3_dp * (f107_daily - f107a_avg) end function parms_to_temp