Ｅ偏光でA1対称な局在モードの電場分布をFDTD法で算出するプログラム。

　基本的にはsa1.fと同じプログラムであり，双極子輻射された電磁場について電磁
エネルギーの代わりに電場分布を出力する。sa1.fで求めた双極子輻射スペクトルの
ピーク周波数をOM1とOM2に入力する。

　プログラム中のパラメータの意味は以下の通りです。

OM1　：双極子輻射スペクトルを計算する規格化周波数の下限。正の実数。
OM2　：双極子輻射スペクトルを計算する規格化周波数の上限。正の実数。
NPOINT　：双極子輻射スペクトルを計算する周波数の個数。正の整数。
NX　：格子定数当たりのデータ点の数で正の整数。
LX　：x軸方向の単位格子の数で正の整数。
NT　：１周期当たりの時間の分割数で正の整数。NTはNXの２倍程度にとって下さい。
　　　　これより小さくすると，計算が不安定になり発散します。そのときは，さらに
　　　　NTを大きくして下さい。
NC　：双極子輻射を計算する時間の長さで正の整数。単位は双極子の振動の周期。
A　：格子定数。
R　：円柱の半径。
QEPS1　：円柱の誘電率（=屈折率の自乗）で実数。（複素数も可ですが計算が不安定）
QEPS2　：背景の誘電率で実数。（複素数も可ですが計算が不安定）
QEPSD　：欠陥部の誘電率で実数。（複素数も可ですが計算が不安定）
X0　：振動双極子のx座標で実数。
Y0　：振動双極子のy座標で実数。0≦Y0≦X0＜LX×Aであることが必要
EMU　：振動双極子の振幅で実数。通常1.0とします。

　計算結果はOPEN文で指定したファイルに保存されます。ファイルに保存されるデータの
並びは以下の通りです。

１行目〜７行目：　計算に用いたパラメータ
８行目：　規格化周波数と放射エネルギー
９行目以降：　x座標，y座標（単位は格子定数），電場のz成分（最大値を１に規格化）

　sa1p.datはsa1p.fを実行した結果。これを図示したものがsa1p.jpgです。