光学的に薄い非熱ジャイロシンクロトロン放射 （試験公開）

注意：改良についてのコメント歓迎。

放射領域物理量から放射を求めることができる。
Dulk (1985) の方法とRamaty (1969)の方法とが用意されている。 ここでは電子の分布関数が以下のようなものであると 仮定している。

$${{d NV(E)} \over {d E}} = K \left({{E}\over{E_{0}}}\right)^{-\delta}$$

$E$は電子のエネルギー（keV）・$NV(E)$はエネルギー$E$以上の電子数（particles）。 この分布はパラメータ$\delta$と$E_{0}$と$K$とによって記述されるが、$K$のかわりに $NV(E_0)=K/(\delta -1)/E_0^{\delta-1}$を与える方が物理的な意味が明確になるので そちらを実際は用いる。以下のプロシジャではDulk (1985)にならい、 $E_0 = 10\ {\rm keV}$で固定している。

(a) 物理量から放射を求める

放射領域の、電子分布べき指数$\delta$・磁場強度（G） ・視線と磁場とのなす角度（度）、さらに （エネルギー10keV以上の）非熱電子数（個）を 入力すると、放射束密度（SFU）・円偏波率が以下のように得られる。 Dulk (1985)の方法では
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc $<$CR$>$
さらに放射源の大きさ（立体角であたえる。単位はsterad）をあたえると 光学的厚さが求まる。
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau $<$CR$>$
同様にRamaty (1969)の方法では
IDL> ramaty_gysy,delta,bb,theta,nv,freq,fi,rc $<$CR$>$
IDL> ramaty_gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau $<$CR$>$

(b) 放射から物理量を求める

Dulk (1985) の近似モデルにもとづいて、放射から放射領域物理量 を求めることができる。
$\alpha$をもとめたあと（節3.7参照）、 非熱電子の分布関数のべき指数がただちに求められる。
IDL> norp_alpha,freq,fi,mvd,mvdfit,alpha_tk,alpha_tn,freqpk,fluxpk $<$CR$>$
IDL> norp_alpha2delta,alpha_tn,delta $<$CR$>$