-
IDL/SPEDAS で解析 map の作り方2024
2024/9/12 の講習会の後、堀さんが作ってくれたマップの作り方例
次の内容をコピペして densitymap.pro とか名前を付けて保存
<![CDATA[
;;;ここから -------------------------------------
pro densitymap
;; 2024/9/12 の講習会 tutorial で教えてもらった手順を打ち込んでみた
dens =''
xsm =''
ysm =''
for i=0, 200 do begin
timespan, time_double('2017-04-01') + i *86400L
del_data , 'erg_*'
erg_load_pwe_hfa, level='l3'
erg_load_orb
if tnames('erg_pwe_hfa_l3_1min_ne_mgf') eq '' or tnames('erg_orb_l2_pos_sm') eq '' then continue
get_data, 'erg_pwe_hfa_l3_1min_ne_mgf',data=d_ne
tinterpol_mxn, 'erg_orb_l2_pos_sm','erg_pwe_hfa_l3_1min_ne_mgf'
get_data, 'erg_orb_l2_pos_sm_interp', data=d_pos
append_array, dens, d_ne.y
append_array, xsm, d_pos.y[*,0]
append_array, ysm, d_pos.y[*,1]
endfor
bin2d, xsm,ysm, dens, xrange=[-6,6],yrange=[-6,6],binsize=[0.2,0.2],xc=xc,yc=yc, ave=ave, binhist=num
;; bin2dについてに説明は
;; https://github.com/spedas-j/member_contrib/wiki/spedas_useful_command#bin2d
plotxyz, xc,yc,ave, /zlog, xtitle='Xsm [Re]', ytitle='Ysm [Re]', ztitle ='Ne [/cc]'
end
;;; ここまで --------------------------
]]>
できたdensitymap.pto を
自分の
c:/ユーザー/n/IDLWorkspace/Default の下におき
- IDLを立ち上げて
IDL> erg_init
コンパイルして実行すると、