\documentstyle{jarticle}
\setlength{\topmargin}{-45pt}
\setlength{\oddsidemargin}{0cm}
\setlength{\evensidemargin}{-0cm}
\setlength{\textheight}{24cm}
\setlength{\textwidth}{16cm}

\newcommand{\ind}{\par\hangindent=0.5cm\hangafter=1\noindent}
\newcommand{\indd}{\par\hangindent=1.0cm\hangafter=1\hspace*{0.5cm}}
\newcommand{\idl}{\hspace*{5mm}{\tt IDL>\ }}
\newcommand{\unix}{\hspace*{5mm}{\tt unix\%\ }}
\newcommand{\crmk}{{\bf $<$CR$>$}}
\newcommand{\norp}{{\tt \$\{NORP\}}}
%\newcommand{\S}{\verb*+ +}
\newcommand{\ssw}{{\tt \$\{SSW\}}}


\begin{document}
\title{\bf 野辺山偏波計（NoRP）解析マニュアル
\\
ver. 1.0}
\author{国立天文台（野辺山）太陽電波}
\date{2000年12月21日}
\maketitle

\clearpage
\tableofcontents

\clearpage
\section{はじめに}

このマニュアルは、国立天文台（野辺山）太陽電波の
偏波計（NoRP; Nakajima et al. 1985）と
旧豊川偏波計（TYKW; Torii et al. 1979; Shibasaki et al. 1979）と
のデータ解析ためのマニュアルである。

このパッケージに関する質問や要望は以下に電子メールを送ってほしい。\\
\bigskip
\hspace*{1cm}{\tt service@solar.nro.nao.ac.jp}

野辺山偏波計の最新の情報は国立天文台太陽電波のウェブページを
参照してほしい。URLは\\
\bigskip
\hspace*{1cm} {\tt http://solar.nro.nao.ac.jp/norp/}\\
である。

\section{解析ソフトを使うには}

ここで紹介する偏波計解析ソフトを使うには
以下の設定が必要である。

\medskip
\noindent
(1) {\bf SolarSoftware}のインストール

「SolarSoftware（SSW）」をサブパッケージ「Radio/NoRP」
を含めてインストールする。これ以外にもRadio/NoRHとYohkoh/SXTとを
含めることを強く推奨する。これらのパッケージがインストールされていない場合は
コンピュータ管理者に相談してほしい。
「SolarSoftware」の配布元は、
{\tt http://www.lmsal.com/solarsoft/}
である。このマニュアルではSSWをインストールしたトップディレクトリを
{\ssw}と定義する。またNoRPサブパッケージのトップディレクトリを
{\norp}（通常は{\ssw}{\tt /radio/norp}）と定義する。

\medskip
\noindent
(2) {\bf NoRPデータベース}のインストール

国立天文台野辺山データアーカイブ
{\tt ftp://solar.nro.nao.ac.jp/pub/norp/data}

からデータベースをコピーしてインストールする。
{\norp}{\tt /data}にインストールするのがデフォルト。

\medskip
\noindent
(3) 個人環境の設定

環境設定ファイル{\tt \char126/.cshrc}中で \\
{\tt
\hspace*{10mm} setenv SSW {\it SSWのディレクトリ}\footnote{
野辺山では {\tt /sgi1/ssw}} \\
\hspace*{10mm} setenv SSW\_INSTR ``norp'' \\
\hspace*{10mm} source \$\{SSW\}/gen/setup/setup.ssw \\
}
を適当な箇所に書き加える。ただし環境変数{\ssw}はサイトによって異なるので
注意。また他の観測機器も同時解析する場合は \\
{\tt
\hspace*{10mm} setenv SSW\_INSTR ``norp norh sxt'' \\
}
などとする。

\section{実際の解析}

\subsection{解析の開始}

この節で紹介するNoRH用IDLプロシジャはすべてSolarSoftware(SSW)に
依存している。前節の設定を行った後実際に解析する際はSSW/IDLを立ち上げる。

\unix {\bf sswidl \crmk}

\subsection{データベース}

\subsubsection{イベントリスト}

1990年から1992年5月までのイベントについては
データベースがある('evx'データ)
\footnote{1992年6月以後のイベントについては電波ヘリオグラフの
イベントリストを参照のこと}。
たとえば
1991年6月9日UTのイベントのリストを表示するには \\
\idl {\bf norp\_pr\_evx,'1991-06-09' \crmk} \\
1991年6月5日UTから6月8日UTまでであれば \\
\idl {\bf st\_day='1991-06-05' \crmk} \\
\idl {\bf ed\_day='1991-06-08' \crmk} \\
\idl {\bf norp\_pr\_evx,st\_day,ed\_day \crmk} \\
17GHzで太陽放射束密度1000SFU以上のイベントのみの場合は\\
\idl {\bf norp\_pr\_evx,st\_day,ed\_day,criterion=1000 \crmk}

\subsubsection{一日分データの平均値}

一日分データのうち静穏時のみを選びだした平均値が
データベース化されている（1987年11月以後）。プリバーストの太陽放射を得たり、
長周期変動の研究に役立ててほしい。

\noindent
旧豊川偏波計（1・2・3.75・9.4GHz）：1987年11月から1994年2月まで \\
旧豊川偏波計（3.75GHzのみ）：1994年3月から1994年4月まで \\
旧野辺山偏波計（17・35・80GHz）：1990年3月から1994年2月まで \\
野辺山偏波計（1・2・3.75・9.4・17・35・80GHz）：1994年5月から現在も稼働中 

ただし、35GHz（2400 SFU）と80GHz（9000 SFU）とは装置の設計上、
静穏太陽放射束密度は固定された値になっているので注意。

読み込みは、たとえば1999年1月1日UTから6月30日UTまでであれば \\
\idl {\bf st\_day='1999-01-01' \crmk} \\
\idl {\bf ed\_day='1999-06-30' \crmk} \\
\idl {\bf norp\_rd\_avg,st\_day,ed\_day,timavg,fiavg,fvavg \crmk} \\
表示は、1GHzの場合\\
\idl {\bf utplot,timavg,fiavg(0,*) \crmk}\\
2GHzの場合\\
\idl {\bf utplot,timavg,fiavg(1,*) \crmk}

\subsection{データのありか}

\subsubsection{キャリブレーション済データのXDR(IDL save)ファイル}

野辺山太陽電波アーカイブ(第\ref{archive}節参照)に
キャリブレーションが済んで
XDR ( IDL save )フォーマットになったファイルがいくらか置いてある。
読み込むにはIDLコマンド{\tt restore}を使う。
まずはこれらを転送して解析するのをすすめる。
定常観測については時間分解能 1 secのデータが、
重要なイベントについては時間分解能 0.1 secのデータがおいてある。 \\
\idl {\bf file=getenv('NORP\_XDR')+'/1999/08/norp19990828\_0056.xdr' \crmk}\\
\idl {\bf restore,file \crmk}\\
\idl {\bf help \crmk}\\
\hspace*{0.5cm} \begin{minipage}[t]{15cm} \small
\begin{verbatim}
DAY             STRING    = '1999-08-28'
FI              FLOAT     = Array[7, 7590]
FIAVG           FLOAT     = Array[7]
FREQ            FLOAT     = Array[7]
FV              FLOAT     = Array[7, 7590]
FVAVG           FLOAT     = Array[7]
MVD             BYTE      = Array[7, 7590]
TIM             STRUCT    = -> ANYTIM2INTS Array[7590]
\end{verbatim}
\end{minipage}\medskip\\
{\tt freq}に観測周波数がGHz単位で
（ここでは7周波数1・2・3.75・9.4・17・35・80GHz）、
{\tt tim}が観測時刻、
{\tt fi}・{\tt fv}に放射束密度I(R+L)成分・V(R--L)成分（単位SFU）が、
{\tt fiavg}・{\tt fvavg}は一日分の静穏時の平均値がはいっている。
また{\tt mvd}は、各データ点が有効・無効をあらわす配列で有効なとき1、
無効なときゼロ。


\subsubsection{生データ}

野辺山太陽電波アーカイブ(第\ref{archive}節参照)に生データの
一部が公開されている。

\subsection{生データの読み込み --- 野辺山偏波計}

野辺山偏波計データは、豊川偏波計との統合前の旧フォーマットと
統合後の新フォーマットとがあるが解析プログラムは共通である。
旧フォーマットデータには1、2、3.75、9.4 GHzのデータは含まれていないので
これらに関するオプションは無視される。


実際に読む込むには、読み込みたい日付(JST)を指定して\\
\idl {\bf day='2000-4-23'} \\
\idl {\bf norp\_rd\_dat,day,mvd,tim,fi,fv,freq \crmk} \\
または、直接ファイル名を指定して\\
\idl {\bf file='./pl921102'} \\
\idl {\bf norp\_rd\_dat,file,mvd,tim,fi,fv,freq \crmk} \\
とする。{\tt freq}に観測周波数がGHz単位で
（ここでは7周波数1・2・3.75・9.4・17・35・80GHz）、
{\tt tim}が観測時刻、
{\tt fi}・{\tt fv}に放射束密度I成分・V成分（単位SFU）がはいっている。
また{\tt mvd}は、各データ点が有効・無効をあらわす配列で有効なとき1、
無効なときゼロ。\\
%
ある時間区間だけを選びだして読み込むには \\
\idl {\bf timerange=['1992-11-2 2:00','1992-11-2 4:00'] \crmk} \\
\idl {\bf norp\_rd\_dat,day,mvd,tim,fi,fv,freq,timerange=timerange \crmk} \\
この場合キャリブレーションは必要な区間だけに行われるので読み込みが速くなる。\\
%
ある周波数だけ読み込むには以下のようにする。たとえば17GHzだけであれば \\
\idl {\bf rdfreq=[0,0,0,0,1,0,0] \crmk} \\
\idl {\bf norp\_rd\_dat,day,mvd,tim,fi,fv,freq,rdfreq=rdfreq \crmk} \\
1GHzと9.4GHzとだけであれば \\
\idl {\bf rdfreq=[1,0,0,1,0,0,0] \crmk} \\
\idl {\bf norp\_rd\_dat,day,mvd,tim,fi,fv,freq,rdfreq=rdfreq \crmk} \\
つまり{\tt rdfreq}の7つの各要素が順に1・2・3.75・9.4・17・35・80GHz
に対応しており値が1のときに読み込む。

\subsection{生データの読み込み --- 旧豊川偏波計}

旧豊川偏波計データは、 1、2、3.75、9.4 GHzのデータが含まれている。

実際に読む込むには、読み込みたい日付(JST)を指定して\\
\idl {\bf day='1992-11-2'} \\
\idl {\bf tykw\_rd\_dat,day,mvd,tim,fi,fv,freq \crmk} \\
とする。{\tt freq}に観測周波数がGHz単位で
（ここでは4周波数1・2・3.75・9.4GHz）、
{\tt tim}が観測時刻、
{\tt fi}・{\tt fv}に放射束密度I成分・V成分（単位SFU）がはいっている。
また{\tt mvd}は、各データ点が有効・無効をあらわす配列で有効なとき1、
無効なときゼロ。\\
%
直接ファイル名を指定して\\
\idl {\bf file='./ty921102.01i' \crmk} ; 1GHz I成分 \\
\idl {\bf file0pa='./ty921102.0pa' \crmk} ; キャリブレーションデータ\\
\idl {\bf filestt='./ty921102.stt' \crmk} ; ステータスデータ \\
\idl {\bf tykw\_rd\_dat,file,file0pa,filestt,mvd,tim,data \crmk} \\
とすることもできる。この場合は、各周波数のI・V成分のどちらかしか一度に
読めない。また生データのほかにキャリブレーションデータ・
ステータスデータが含まれたファイル名も同時に指定しなければならない。\\
%
ある時間区間だけを選びだして読み込むには \\
\idl {\bf timerange=['1992-11-2 2:00','1992-11-2 4:00'] \crmk} \\
\idl {\bf tykw\_rd\_dat,day,mvd,tim,fi,fv,freq,timerange=timerange \crmk} \\
この場合キャリブレーションは必要な区間だけに行われるので読み込みが速くなる。\\
%
ある周波数だけ読み込むには以下のようにする。たとえば2GHzだけであれば \\
\idl {\bf rdfreq=[0,1,0,0] \crmk} \\
\idl {\bf tykw\_rd\_dat,file,mvd,tim,fi,fv,freq,rdfreq=rdfreq \crmk} \\
1GHzと9GHzとだけであれば \\
\idl {\bf rdfreq=[1,0,0,1] \crmk} \\
\idl {\bf tykw\_rd\_dat,file,mvd,tim,fi,fv,freq,rdfreq=rdfreq \crmk} \\
つまり{\tt rdfreq}の4つの各要素が順に1、2、3.75、9.4 GHzに対応しており
値が1のときに読み込む。

\subsection{グラフ表示}

グラフ表示するには、\\
\idl {\bf mfreq=0 \crmk} \\
\idl {\bf norp\_plot,mfreq,mvd,tim,fi \crmk} \\
ここで{\tt mfreq}は、{\tt fi}（複数の周波数分を含む）に
はいっているデータのうち表示したいデータの周波数の番号である。\\
別の周波数のデータを重ねるには\\
\idl {\bf mfreq1=1 \crmk} \\
\idl {\bf norp\_plot,mfreq1,mvd,tim,fi,/over \crmk} \\
とする。直接指定して\\
\idl {\bf utplot,tim(where(mvd(mfreq,*))),fi(mfreq,where(mvd(mfreq,*))) \crmk} \\
とすることもできる。

\subsection{スペクトル解析}
\label{section:spectrum}

複数の周波数データを組み合わせることでスペクトルを求めることができる。
求まる量は$\alpha$で定義は
$$
F_{\nu} \propto {\nu}^{\alpha}
$$
である。実際は$\alpha$は全周波数帯にわたって一価ではなく、
マイクロ波帯のスペクトルは次のような関数でだいたいフィッテングできる。
\begin{equation}
F_{\nu} = \widehat{F_{\nu}}
\left(\frac{\nu}{\widehat{\nu}}\right)^{\alpha_{\rm tk}}
\left\{
1-\exp\left[
-\left(\frac{\nu}{\widehat{\nu}}\right)^{\alpha_{\rm tn}-\alpha_{\rm tk}}
\right]
\right\}
\approx
\left\{
\begin{array}{cc}
\widehat{F_{\nu}}
\left({\nu}/{\widehat{\nu}}\right)^{\alpha_{\rm tk}} 
  & \mbox{for $\nu \ll \widehat{\nu}$}\\
\widehat{F_{\nu}}
\left({\nu}/{\widehat{\nu}}\right)^{\alpha_{\rm tn}} 
  & \mbox{for $\nu \gg \widehat{\nu}$}
\end{array}
\right.
\label{equation:fitting}
\end{equation}
ここでフィッテングのパラメータは、
$\widehat{\nu}$はターンオーバー周波数・
$\widehat{F_{\nu}}$はターンオーバー放射束密度・
$\alpha_{\rm tk}$は光学的に厚い側（低周波数側）のべき指数・
$\alpha_{\rm tn}$は光学的に薄い側（高周波数側）のべき指数で
の4変数である。 
実際のフィッテングは以下のようにおこなう。\\
\idl {\bf day='2000-4-8'} \\
\idl {\bf norp\_rd\_dat,day,mvd,tim,fi,fv,freq \crmk} \\
\idl {\bf norp\_rd\_avg,day,timavg,fiavg,fvavg \crmk} \\
\idl {\bf for m=0,6 do fi(m,*)=fi(m,*)-fiavg(m) \crmk} \\
\idl {\bf norp\_alpha,freq,fi,mvd,mvdfit,alpha\_tk,alpha\_tn,freqpk,fluxpk 
\crmk} \\
フレア前静穏状態の放射束密度（一日の平均値）を
読み込んだデータから差し引いてフレア成分だけをとりだし、
それを{\tt norp\_alpha}にわたす。
{\tt mvdfit}は得られたフィッティングが有効な場合1が無効な場合
（たとえば信号が弱いなどの理由で）についてはゼロがはいっている。

\medskip
{\bf 注意：} 通常はこのような単純な手順ではフィッティングできない。
実際は(1)積分時間を長めにとってS/Nを良くする、(2)上記の関数から
明らかにはずれた点をフィッティングからはずす（理由：
低周波数は放射メカニズムが異なる、高周波数はS/Nが良くない）
などの手続が必要である。手順は節 \ref{section:spectrum_method}を参照のこと。


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

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

放射領域物理量から放射を求めることができる。\\
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}$で固定している。 \\

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

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

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

Dulk (1985) の近似モデルにもとづいて、放射から放射領域物理量
を求めることができる。\\
$\alpha$をもとめたあと（節\ref{section:spectrum}参照）、
非熱電子の分布関数のべき指数がただちに求められる。\\
\idl {\bf norp\_alpha,freq,fi,mvd,mvdfit,alpha\_tk,alpha\_tn,freqpk,fluxpk 
\crmk} \\
\idl {\bf norp\_alpha2delta,alpha\_tn,delta \crmk} 

\subsubsection{光学的に薄い熱制動放射 （試験公開）}

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

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

Dulk (1985) の近似モデルにもとづいて、放射領域物理量から放射
を求めることができる。\\
放射領域の電子温度（K）・視線方向磁場強度（G）
・ボリュームエミッションメジャ（cm$^{-3}$）を
入力すると、放射束密度（SFU）・円偏波率が以下のように得られる。 \\
\idl {\bf dulk\_frfr,te,b\_loc,vem,freq,fi,rc \crmk} \\
さらに放射源のサイズ（sterad）をあたえると光学的厚みも得られる。\\
\idl {\bf dulk\_frfr,te,b\_loc,vem,freq,fi,rc,omega,tau \crmk}

%\bigskip
%\noindent
%(b) 放射から物理量を求める
%
%Dulk (1985) の近似モデルにもとづいて、放射から放射領域物理量
%を求めることができる。\\
%17GHz放射束密度のI成分・V成分をあたえると、視線方向磁場と$VEM/\sqrt{T_e}$
%（$VEM$はボリュームエミッションメジャ・$T_e$は電子温度）とが求まる。 \\
%\idl {\bf norp\_frfr\_inv,freq,fi,fv,b\_los,vem\_over\_sqrtte \crmk} \\
%電子温度かボリュームエミッションメジャを仮定するともう一方が求まる。\\
%\idl {\bf norp\_frfr\_inv,freq,fi,fv,b\_los,tein,vem,/assumete \crmk} \\
%\idl {\bf norp\_frfr\_inv,freq,fi,fv,b\_los,vemin,te,/assumevem \crmk} \\
%さらに放射源の大きさ（立体角であたえる。単位はsterad）をあたえると
%光学的厚さが求まる。\\
%\idl {\bf norp\_frfr\_inv,freq,fi,fv,b\_los,tein,vem,/assumete,omega,tau 
%\crmk} \\
%\idl {\bf norp\_frfr\_inv,freq,fi,fv,b\_los,vemin,te,/assumevem,omega,tau \crmk}

\subsection{具体的なデータ解析の例}

この節では解析の例を記す。

\subsubsection{イベント選択}

重要なイベントについては、偏波計のウェブサイトに紹介されている。
まずはそこをみて概要をつかむのがのぞましい。URLは \\
\hspace*{1cm}{\tt http://solar.nro.nao.ac.jp/norp/}

\subsubsection{光度曲線}

イベントを選びだしたら、光度曲線をみてみよう。
重要なイベントについてはキャリブレーションの済んだデータが
XDR（IDL save）フォーマットで保存してあるのでそれを読む。\\
\idl {\bf file=getenv('NORP\_XDR')+'/1999/08/norp19990828\_0056.xdr' \crmk}\\
\idl {\bf restore,file \crmk} \\
グラフ表示するには、\\
\idl {\bf mfreq=0 \crmk} \\
\idl {\bf norp\_plot,mfreq,mvd,tim,fi \crmk} \\
ここで{\tt mfreq}は、表示したいデータの周波数の{\tt freq}の中での順番である。\\
別の周波数のデータを重ねるには\\
\idl {\bf mfreq1=1 \crmk} \\
\idl {\bf norp\_plot,mfreq1,mvd,tim,fi,/over \crmk} \\
とする。

\subsubsection{スペクトル解析}
\label{section:spectrum_method}

スペクトルを計算する際は積分時間を長くとった方がよい。300点（30秒間）分の
データを積分するには\\
\idl {\bf mint=300 \crmk} \\
\idl {\bf norp\_mkint,mint,mvd,tim,fi,fv,mvdav,timav,fiav,fvav \crmk} \\
次にフレア前静穏電波放射束を差し引く。\\
\idl {\bf for m=0,6 do fiav(m,*)=fiav(m,*)-fiavg(m) \crmk} \\
えられた結果をスペクトルフィッティングする。\\
\idl {\bf norp\_alpha,freq,fiav,mvdav,mvdfit,alpha\_tk,alpha\_tn,freqpk,fluxpk 
\crmk} \\
放射がジャイロシンクロトロン放射であると仮定すると非熱電子のべき$\delta$
がただちにもとまる。\\
\idl {\bf norp\_alpha2delta,alpha\_tn,delta \crmk}\\

\medskip
また、フィッティング関数（方程式\ref{equation:fitting}）から
明らかにはずれた観測点があるとフィッティングがうまくいかないことがある。
理由は、低周波数では放射メカニズムが異なる（ジャイロシンクロトロンではなく
プラズマ放射）のと、高周波数はS/Nが良くない（特に80GHzはよくない）ためである。
そのようなときはその観測点をフィッティングから以下のようにしてはずす。\\
\idl {\bf mvfreq=[2,3,4,5,6] \crmk} \\
\idl {\bf norp\_alpha,freq(mvfreq),fiav(mvfreq,*),mvdav(mvfreq,*) \$ \crmk}\\
\idl {\bf ,mvdfit,alpha\_tk,alpha\_tn,freqpk,fluxpk \crmk} \\

\medskip
フィッティングした結果を確かめるために、フィッティング関数と観測点とを
重ね合わせてみる。たとえば{\tt n}ステップ目のデータを見るには \\
\idl {\bf norp\_funcp,alpha\_tk(n),alpha\_tn(n),freqpk(n)
,fluxpk(n),freqfit,fifit \crmk} \\
\idl {\bf plot,freq,fiav(*,n),/xlog,/ylog,psym=5 \crmk} \\
\idl {\bf oplot,freqfit,fifit \crmk}

\subsection{その他}

\subsubsection{生データを直接よみこむ}

生データを直接よみこむ方法は、野辺山データは\\
\idl {\bf norp\_rd\_rdt,file,rdata,date,version \crmk} \\
豊川データは\\
\idl {\bf tykw\_rd\_rdt,file,rdata \crmk} \\
データに対応する時刻（JST）を得るには\\
\idl {\bf timej=norp\_rdt2timej(rdata) \crmk}


\appendix

\section{データアーカイブ匿名FTPサイト}
\label{archive}

観測データベースファイルは日々蓄積されていくので
適宜更新する必要がある。国立天文台野辺山太陽電波では
生データを含むデータベースを匿名FTPサーバで公開している。
NoRPデータのURLは

\bigskip
{\tt ftp://solar.nro.nao.ac.jp/pub/norp}
\bigskip

\noindent
である。
また大量のデータ転送はネットワーク負荷を軽減するために
できるだけCD-ROMによる郵送をお願いしたい。
電子メールでアドレス

\bigskip
{\tt service@solar.nro.nao.ac.jp}

\bigskip
\noindent
に報せてもらえればただちに作成して送る。

\section*{謝辞}

以下の方々に感謝します。

SolarSoftWareは、
太陽観測データ解析用のIDLパッケージで、
NASAのさまざまなプロジェクトのもとで
支援を受けています。
IDLはResearch Systems Inc.が製造・販売しているデータ解析ソフトウェアです。

\section*{参考文献}

\noindent
Dulk, G. A., 1985, ARAA, {\bf 23}, 169 \\
%
Nakajima et al. "The Radiometer and Polarimeters at 80, 35, and 17 GHz
for Solar Observations at Nobeyama", PASJ, 37, 163 (1985) \\
%
Ramaty, 1969, ApJ, {\bf 158}, 753 \\
%
Shibasaki et al. "Solar Radio Data Acquisition and Communication System
(SORDACS) of Toyokawa Observatory",
Proc. of the Res. Inst. of Atmospherics,
Nagoya Univ., 26, 117 (1979) \\
%
Torii et al. "Full-Automatic Radiopolarimeters for Solar Patrol at
Microwave Frequencies", Proc. of the Res. Inst. of Atmospherics,
Nagoya Univ., 26, 129 (1979) \\

\end{document}
