EISCATデータ講習用　コマンドリスト（2012/02/22 版）

------------------------------------------------------
（Lesson1）EISCATサマリープロットの表示
------------------------------------------------------

（1-1）初期化：
IDL> thm_init

（1-2）プロットする日にちの入力：
THEMIS> timespan, '2012-01-22',1

（1-3）データの読み込み（まずはTromso UHF データ）：
THEMIS> iug_load_eiscat, site='tro_uhf'

（1-4）ロードデータ名の確認：
THEMIS> tplot_names

（1-5）データのプロット：
THEMIS> tplot,['eiscat_trouhf_beat_all_ne','eiscat_trouhf_beat_all_te','eiscat_trouhf_beat_all_ti','eiscat_trouhf_beat_all_vi']

（1-6）軸の変更（Y軸の高度やカラースケールの変更）：
THEMIS> ylim, 'eiscat_trouhf_beat_*', 90, 600
THEMIS> zlim, 'eiscat_trouhf_beat_vi', -400, 400
THEMIS> tplot

------------------------------------------------------
（Lesson2）関連する複数のプロットを同時表示
------------------------------------------------------

（2-1）データの読み込み（すべてのEISCATデータ）：
THEMIS > iug_load_eiscat
（注：現状では、
THEMIS > iug_load_eiscat, site='tro_uhf'
THEMIS > iug_load_eiscat, site='esr_42m'
のように、個別に読み込む方が処理が早い状況。）

（2-2）ロードデータ名の確認：
THEMIS > tplot_names

（2-3）電子密度データのプロット（ここではESR42m、トロムソUHFレーダーデータ）：
THEMIS> tplot, 'eiscat_*_all_ne'
(tplot, ['eiscat_esr42m_ipy0_all_ne','eiscat_trouhf_beat_all_ne']と同じ)

（2-4）IMF Bz データのプロットを追加
THEMIS> omni_hro_load
THEMIS> tplot_names
THEMIS> tplot,['eiscat_*_all_ne','OMNI_HRO_1min_BZ_GSM']

（2-5）光学データのケオグラムを追加
THEMIS> iug_load_ask_nipr, site='tro'
THEMIS> iug_load_ask_nipr, site='lyr'
THEMIS> tplot_names
THEMIS> ylim, 'iug_ask_*', -75, 75
THEMIS> tplot,['eiscat_*_all_ne','OMNI_HRO_1min_BZ_GSM','iug_ask_*']
THEMIS> zlim, 'iug_ask_lyr', 0, 150
THEMIS> zlim, 'iug_ask_tro', 0, 200
THEMIS> tplot

------------------------------------------------------
（Lesson3）詳細なデータ比較1（@ロングイアビン 06-08 UT）
------------------------------------------------------

（3-1）データの読み込み（ESR42mデータとロングイアビン光学データ）：
THEMIS> timespan, '2012-01-22',1
THEMIS> iug_load_eiscat, site='esr_42m'
THEMIS> iug_load_ask_nipr, site='lyr'
THEMIS> omni_hro_load

（3-2）ロードデータ名の確認とプロット：
THEMIS > tplot_names
THEMIS> tplot,['eiscat_esr42m_ipy0_all_ne','eiscat_esr42m_ipy0_all_te','iug_ask_lyr','OMNI_HRO_1min_BZ_GSM']

（3-3）時刻や軸、カラースケールの指定
THEMIS> tlimit, '2012-01-22/06:00:00', '2012-01-22/08:00:00'
THEMIS> ylim, 'eiscat_esr42m_ipy0_all_*', 90, 500
THEMIS> ylim, 'iug_ask_lyr', -75, 75
THEMIS> zlim, 'iug_ask_lyr', 0, 100
THEMIS> tplot

------------------------------------------------------
（Lesson4）詳細なデータ比較2（@トロムソ 20-23 UT）
------------------------------------------------------

（4-1）データの読み込み（トロムソUHFデータと光学データ）：
THEMIS> timespan, '2012-01-22',1
THEMIS> iug_load_eiscat, site='tro_uhf'
THEMIS> iug_load_ask_nipr, site='tro'
THEMIS> omni_hro_load

（4-2）ロードデータ名の確認とプロット：
THEMIS > tplot_names
THEMIS> tplot,['eiscat_trouhf_beat_all_ne','eiscat_trouhf_beat_all_te','eiscat_trouhf_beat_all_ti','iug_ask_tro','OMNI_HRO_1min_BZ_GSM']

（4-3）時刻や軸、カラースケールの指定
THEMIS> tlimit, '2012-01-22/20:00:00', '2012-01-22/23:00:00'
THEMIS> ylim, 'eiscat_trouhf_beat_all_*', 90, 700
THEMIS> ylim, 'iug_ask_tro', -75, 75
THEMIS> zlim, 'iug_ask_tro', 0, 200
THEMIS> tplot

------------------------------------------------------
（Lesson5）ある時刻の高度分布を表示
------------------------------------------------------

（5-1）特定の時刻のためのインデックスを見つける：
THEMIS > inds=nn('eiscat_trouhf_beat_all_ne','2012-01-22/20:56:00')
THEMIS > print, inds

（5-2）tplotの変数をIDLの変数に変換する：
THEMIS > get_data, 'eiscat_trouhf_beat_all_alt', data=alt1
THEMIS > get_data, 'eiscat_trouhf_beat_all_ne', data=ne1
THEMIS > help, ne1, /struct

（5-3）新しいウィンドウを開き、プロットを表示する：
THEMIS > window,1
THEMIS > plot, ne1.y(inds, *), alt1.y(inds, *), /xlog, xtitle='Ne', ytitle='Altitude'

（5-4）その他の方法として、マウスで表示したい時刻を選択：
THEMIS > window,2,xsize=256,ysize=256
THEMIS > ctime,/cut
（その後、右クリックでコマンドプロンプトに戻る。）

------------------------------------------------------
（Lesson6）その他
------------------------------------------------------
（6-1）複数日にまたがる連続観測データの表示例：
THEMIS> timespan, '2012-01-21',3
THEMIS> iug_load_eiscat , site='tro_uhf'
THEMIS> iug_load_ask_nipr, site='tro'
THEMIS> tplot,['eiscat_trouhf_beat_all_ne','iug_ask_tro','OMNI_HRO_1min_BZ_GSM']

------------------------------------------------------