本研究は一部日本学術振興会科学研究費補助金と 高度情報科学技術研究機構(RIST) 技術振興調整費の援助を受けた. プログラムの開発とテスト計算には一部 国立天文台の共同研究ならびに国立天文台データ解析計算センターの 共同利用計算機を利用させて頂いた.