CUI操作の紹介
2015年8月19日
第2回「太陽地球環境データ解析に基づく超高層大気の空間・時間変動の解明」
データ解析セッション田中良昌（極地研）
1

SPEDASソースコードのダウンロード
THEMISウェブサイトから、SPEDASのソースコードをダウンロードする。
1. THEMIS satellite software ウェブサイトにアクセス
http://themis.ssl.berkeley.edu/software.shtml
2. “Download” セクションまでスクロールダウン

3. Download TDAS 9.0 + SPEDAS 1.00
source (13 MB) をクリックし、ソースコードをダウンロードする。

2

SPEDASのインストール
Windowsの例（Macでもほぼ同様）：

メニューからウィンドウ→設定 を選択。

 ダウンロード＆解凍したSPEDASフォルダをIDLパスに挿入する。
 左のチェックボックスにチェックを入れる。

3

UDAS最新版の利⽤⽅法
IUGONETのウェブサイトから、UDASをダウンロード＆インストールすることにより、最新のUDASのコマンドが利用できるようになる。
THEMISウェブサイト
http://themis.ssl.berkeley.edu/index.shtml

SPEDAS
UDASのほとんどのロードルーチンが使⽤できる。

IUGONETウェブサイト

http://www.iugonet.org/

インストール

ユーザのPC

UDAS
最新のUDASのロードルーチンを使えるようになる。
4

UDAS最新版のインストール
Windowsの例（Macでもほぼ同様）：

メニューからウィンドウ→設定 を選択。

 ダウンロード＆解凍したUDASフォルダをIDLパスに挿入する。
 既にインストールしたSPEDASよりも、UDASが上に来るようにする（UDASが先に読まれるようにする）。

SPEDASを利用したデータ読込描画概要
インターネット上のデータサーバ

ユーザーのPC
自動ダウンロードdata

データはtplot
変数として読み込まれるSSL, Berkeley,THEMIS, GBO
CDAWeb, OMNI,ACE, Wind, etc.

ダウンロードディレクト
リは自動的に作成される.

SPEDAS/TDAS

data
data

SPEDAS-CUI ツールを使うと、3つの基本コマンドで簡単にデータ読み込み・描画をすることが出来ます1. 時間幅を設定
2. ロードプロシージャを実行
3. プロットプロシージャを実行

timespan, ‘yyyy-mm-dd’
iug_load_***tplot, +++6

UDASに含まれるロードコマンド
現在、UDAS s1.00.1（for SPEDAS v1.00）を公開中
No.

Instrument Type

Load routines1

太陽磁場活動望遠鏡（SMART）データ

iug_load_smart

2

太陽VHF/UHF電波スペクトルデータ

iug_load_iprt

3

HF帯太陽・木星電波広帯域スペクトルデータ

iug_load_hf_tohokuu

4

自動気象観測装置データ

iug_load_aws_rish

5

境界層レーダーデータ

iug_load_blr_rish

6

Lバンド下部対流圏レーダーデータ

iug_load_ltr_rish

7

赤道大気レーダーデータ

iug_load_ear

8

MUレーダーデータ

iug_load_mu

9

流星レーダーデータ

iug_load_meteor_rish10

MFレーダーデータ

iug_load_mf_rish

11

ウィンドプロファイラーレーダーデータ

iug_load_wpr_rish

12

イオノゾンデ（信楽）データ

iug_load_ionosonde_rish

13

ラジオゾンデデータ

iug_load_radiosonde_rish

 (*)は、ERG-SCで開発されたロードコマンドのエイリアス。

7

UDASに含まれるロードコマンド
No.

Instrument Type

Load routines

14

SuperDARNレーダーデータ(*)

iug_load_sdfit (*)

15

EISCATレーダーデータ

iug_load_eiscat

16

EISCATレーダーデータ（イオン速度/電場ベクトル）

iug_load_eiscat_vief

17

イメージングリオメータデータ

iug_load_irio_nipr

18

LF帯標準電波観測データ

iug_load_lfrto

19

アジアVLF観測ネットワーク (AVON/VLF-B) データ

iug_load_avon_vlfb

20

超高層大気イメージングシステム (OMTI) データ

iug_load_camera_omti_asi (*)

21

全天イメージャデータ

iug_load_asi_nipr

22

全天イメージャケオグラムデータ

iug_load_ask_nipr

23

地磁気指数（AE, Dst, ASY/SYM）、WDC地磁気観測所
データ

iug_load_gmag_wdc

24

昭和基地、アイスランド地磁気データ、南極無人磁力計ネットワーク

iug_load_gmag_nipr

25

210°地磁気観測網データ(*)

iug_load_gmag_mm210 (*)

26

MAGDAS地磁気1秒値データ

iug_load_gmag_magdas_1sec (*)

27

STEL誘導磁力計観測網データ(*)

iug_load_gmag_stel_induction (*)

28

昭和基地、アイスランド誘導磁力計データ

iug_load_gmag_nipr_induction

29

九大GCMシミュレーションデータ

Iug_load_kyushugcm8

HTMLコマンドリファレンス①
http://themis.ssl.berkeley.edu/software.shtml

HTML documentationをクリック

9

HTMLコマンドリファレンス②
コマンドの使い方の説明が表示される。

ここをクリック

プログラムの中身（IDL言語）が表示される。

10

マニュアル、講習会資料
http://www.iugonet.org/software/documentation.html

逆引きユーザマニュアル（日本語）UDASロードプロシジャのリスト。
使い方動画へのリンク。
過去の講習会資料（PDF, 日本語）

11

ユーザーズガイドなど（日本語）
SPEDAS/UDAS逆引きユーザーマニュアル
SPEDAS/TDAS情報
https://github.com/spedas-j

逆引きできる SPEDASに関する情報が日本語で記載されている。
 SPEDASユーザが便利なプログラム等が置かれている。
12

動画によるUDASの使い方の解説
http://www.iugonet.org/software/loadprocedures.html
からリンクしている。

Crib Sheet（サンプルスクリプト）•Crib Sheetとは、 用途毎にSPEDASのコマンドの例が書かれているスクリプト。
iug_crib_aws_rish.proの例• コピー＆ペーストして実行するか、;+
THEMIS> .run Crib_Sheet_名;PROCEDURE: IUG_CRIB_AWS_RISH.PRO;で実行する。
•

spedas_1_00/idl/ディレクトリ以下の次のような場所に置かれている。
general/examples
general/missions/rbsp/efw/examples
projects/erg/examples
projects/iugonet/examples
projects/themis/examples/advanced
projects/themis/examples/basic
projects/themis/examples/deprecated;Initializes system variables for themis:;=====================================
thm_init;Specify timespan:;=================
timespan,'1994-05-01',7,/day;Load zonal, meridional and vertical winds atShigaraki in timespan:;=====================================
iug_load_aws_rish, site = 'sgk';Plot time-height distribution of zonal and meridional
winds:;=====================================
tplot,['iug_aws_sgk_uwnd','iug_aws_sgk_vwnd']

例１：太陽風パラメータ（OMNI）のプロット
1. IDL> thm_init （初期化）
2. THEMIS> timespan, '2015-3-15', 4 （日時設定）
3. THEMIS> omni_hro_load （OMNIデータをロード）
4. THEMIS> calc, "'OMNI_HRO_1min_Vxn'=-1*'OMNI_HRO_1min_Vx‘”（簡単な計算。太陽風速度Vxが太陽方向なので、-1をかけて地球方向にする）
5. THEMIS> tplot, ['OMNI_HRO_1min_Vxn', 'OMNI_HRO_1min_proton_density', $'OMNI_HRO_1min_BZ_GSM']

15

例２：地磁気データのプロット
1. THEMIS> kyodo_load_dst（Dst指数をロード）
2. THEMIS> kyodo_load_ae（AE指数をロード）
3. THEMIS> iug_load_gmag_nipr, site='syo‘（昭和基地の地磁気をロード）
4. THEMIS> tplot, $['OMNI_HRO_1min_Vxn’, $'OMNI_HRO_1min_proton_density', $'OMNI_HRO_1min_BZ_GSM', $'kyoto_dst', $'kyoto_ae', $'nipr_mag_syo_1sec']（プロット）
5. THEMIS> tlimit, '2015-3-17/00', $'2015-3-18/12‘（時間軸を変更）16

例３：SuperDARNレーダーデータのプロット
1. IDL> thm_init （初期化）
2. THEMIS> timespan, '2015-3-17' （日時設定）
3. THEMIS> iug_load_sdfit, site=‘hok’ （SuperDARN HOKデータをロード）
4. THEMIS> tplot, ['sd_hok_pwr_1', 'sd_hok_vlos_1', 'sd_hok_spec_width_1']（SDのエコーパワー、ドップラー速度、スペクトル幅をプロット）
5. THEMIS> zlim, 'sd_hok_vlos_1', -400, 400（ドップラー速度のカラースケールを変更）
6. THEMIS> tplot

17

例４：EISCATレーダーデータのプロット
1. IDL> thm_init （初期化）
2. THEMIS> timespan,‘2013-7-9’, 3（日時設定）
3. THEMIS> iug_load_eiscat（EISCATレーダーの全てのデータをロード）
4. THEMIS> tplot, $['eiscat_trovhf_ne', $'eiscat_trovhf_te', $'eiscat_trovhf_ti', $'eiscat_trovhf_vi']（EISCATレーダーの基本物理データ（電子密度、電子温度、イオン温度、イオン速度）をプロット）18

例５：ケオグラムのプロット
1. IDL> thm_init （初期化）
2. THEMIS> timespan,‘2012-01-22’, 1（日時設定）
3. THEMIS> iug_load_ask_nipr, site=‘hus’（Husafellの全天イメージャのケオグラムをロード）
4. THEMIS> tplot, ['nipr_ask_hus_0000_ns’, ‘nipr_ask_hus_0000_ew']（ケオグラムをプロット。tplot, [1, 2] でも同じ。）
6. THEMIS> tlimit, '2012-01-22/20', '2012-01-22/22‘（時刻の範囲を絞る。）

19

例６：全天イメージャ（ASI）画像のロード
1. THEMIS> timespan,‘2012-01-22/20:30’, /min, 30 （時刻範囲を変更）
2. THEMIS> iug_load_asi_nipr, site=‘hus‘ （HusafellのASIデータをロード。）
3. THEMIS> iug_load_asi_nipr, site=‘tro‘ （TromsoのASIデータをロード。）
4. THEMIS> window, 1, xsize=600, ysize=600 & erase（新しいウィンドウを開く）
5. THEMIS> map2d_time, '2012-1-22/20:50‘（2次元データをプロットする時刻を設定）
6. THEMIS> map2d_set, glatc=65., glonc=0., $scale=30e+6, /label（地理座標を表示する。中心座標を緯度65°、経度0°に指定。）
7. THEMIS> overlay_map_asi_nipr, $['nipr_asi_hus_0000', 'nipr_asi_tro_0000']（ASIデータを地図上にプロット）
8. THEMIS> overlay_map_coast（海岸線をプロット）20

例７：THEMIS GBOの全天イメージャ画像群のプロット
1. IDL> thm_init （初期化）
2. THEMIS> thm_asi_create_mosaic, '2008-02-10/07:13:30', /verbose,exclude='kian'21