;; Slide #6 ; Set the time range to 9-18 UT on Sep. 21, 2019 timespan,'2019-09-21/09:00:00',9,/hour ; Load Arase l2 def orbit data erg_load_orb ; Plot Arase orbit in the SM coord. tplot,'erg_orb_l2_pos_sm' ;; Slide #7 ; Set the time range to 9-18 UT on Sep. 21, 2019 timespan,'2019-09-21/09:00:00',9,/hour ; Load Arase l2 def orbit data erg_load_orb ; Plot Arase orbit on the SM X-Y plane tplotxy,'erg_orb_l2_pos_sm' ;; Slide #8 ; Plot Arase orbit on the X-Y plane of - 6 Re < x,y < 6 Re tplotxy,'erg_orb_l2_pos_sm',versus='xy',xrange = [-6,6],yrange = [-6,6] ; The same but on the X-Z plane tplotxy,'erg_orb_l2_pos_sm',versus='xz',xrange = [-6,6],yrange = [-6,6] ;; Slide #10 ; Initialize the map2d environment and Set the coordinate system to geographical coordinate map2d_init & map2d_coord, 'geo' ; Generate an empty window window, 0, xsize=800, ysize=640 & erase ; Draw the grid (10 deg in lat & 15 deg in log) on the window and Draw the coast lines map2d_set & overlay_map_coast ;; Slide #11 ; get l2 data for ionospheric footprints on the northern hemisphere erg_load_orb get_data, 'erg_orb_l2_pos_iono_north', data=data_foot ; Store latitude and longitude separately. store_data,'IGRF_ifoot_glat',data= {x:data_foot.x,y:reform(data_foot.y[*,0])} store_data,'IGRF_ifoot_glog',data= {x:data_foot.x,y:reform(data_foot.y[*,1])} ; Overplot Arase footprint on the map overlay_map_sc_ifoot, 'IGRF_ifoot_glat', 'IGRF_ifoot_glog' ;; Slide #12 ; get l3 OP77Q orbit data for ionospheric footprints erg_load_orb_l3,model = 'op' & get_data, 'erg_orb_l3_pos_iono_north_op', data=data_foot_op store_data, 'OP77Q_ifoot_glat',data= {x:data_foot_op.x, y:reform(data_foot_op.y[*,0])} store_data, 'OP77Q_ifoot_glog',data= {x:data_foot_op.x, y:reform(data_foot_op.y[*,1])} ; get l3 T89 orbit data for ionospheric footprints erg_load_orb_l3,model = 't89' & get_data, 'erg_orb_l3_pos_iono_north_t89', data=data_foot_t89 store_data, 'T89_ifoot_glat',data= {x:data_foot_t89.x, y:reform(data_foot_t89.y[*,0])} store_data, 'T89_ifoot_glog',data= {x:data_foot_t89.x,y:reform(data_foot_t89.y[*,1])} ; get l3 TS04 orbit data for ionospheric footprints erg_load_orb_l3,model = 'ts04' & get_data, 'erg_orb_l3_pos_iono_north_TS04', data=data_foot_ts04 store_data, 'TS04_ifoot_glat',data= {x:data_foot_ts04.x, y:reform(data_foot_ts04.y[*,0])} store_data, 'TS04_ifoot_glog',data= {x:data_foot_ts04.x,y:reform(data_foot_ts04.y[*,1])} ;; Slide #13 ; Generate a new map window window, 0, xsize=800, ysize=640 & erase map2d_set,glonc = 280,glatc=60,scale = 1e7 & overlay_map_coast ; Overplot Arase footprints using different models overlay_map_sc_ifoot, 'IGRF_ifoot_glat', 'IGRF_ifoot_glog' overlay_map_sc_ifoot, 'OP77Q_ifoot_glat', 'OP77Q_ifoot_glog',trace_color=2 overlay_map_sc_ifoot, 'T89_ifoot_glat', 'T89_ifoot_glog',trace_color=4 overlay_map_sc_ifoot, 'TS04_ifoot_glat', 'TS04_ifoot_glog',trace_color=6