We approximate the infitine sums, which include terms like j_n(b), as finite sums from n=-Nlarge to Nlarge. Since j_n(b) is pretty small when b = n/2, we just make sure that n is large enough for all b (which is related to vperp) values
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in) | :: | vperpmax |
max value to be computed in normalized velocity space |
||
real, | intent(in) | :: | tau_s |
T_ref/T_s|_parallel |
||
real, | intent(in) | :: | mu_s |
m_ref/m_s |
||
real, | intent(in) | :: | aleph_r |
T_perp/T_parallel_R |