geospacelab 0.10.1__tar.gz → 0.10.2__tar.gz

{geospacelab-0.10.1/geospacelab.egg-info → geospacelab-0.10.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: geospacelab
-Version: 0.10.1
+Version: 0.10.2
 Summary: Collect, manage, and visualize geospace data.
 Home-page: https://github.com/JouleCai/geospacelab
 Author: Lei Cai
@@ -378,14 +378,9 @@ Output:
 From then on, the data source has been switched from JHUAPL to CDAWeb.
  
 ```python
-import datetime
-import matplotlib.pyplot as plt
-
-import geospacelab.visualization.mpl.geomap.geodashboards as geomap
-
 dt_fr = datetime.datetime(2015, 9, 8, 8)
 dt_to = datetime.datetime(2015, 9, 8, 23, 59)
-time1 = datetime.datetime(2015, 9, 8, 20, 21)
+time_c = datetime.datetime(2015, 9, 8, 20, 21)
 pole = 'N'
 sat_id = 'f16'
 band = 'LBHS'
@@ -396,49 +391,51 @@ dashboard = geomap.GeoDashboard(dt_fr=dt_fr, dt_to=dt_to, figure_config={'figsiz
 # If the orbit_id is specified, only one file will be downloaded. This option saves the downloading time.
 # dashboard.dock(datasource_contents=['jhuapl', 'dmsp', 'ssusi', 'edraur'], pole='N', sat_id='f17', orbit_id='46863')
 # If not specified, the data during the whole day will be downloaded.
-dashboard.dock(
-    datasource_contents=dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], 
-    pole=pole, sat_id=sat_id, orbit_id=None
-)
+ds_ssusi = dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], pole=pole, sat_id=sat_id, orbit_id=None)
 ds_s1 = dashboard.dock(
     datasource_contents=['madrigal', 'satellites', 'dmsp', 's1'],
-    dt_fr=time1 - datetime.timedelta(minutes=45),
-    dt_to=time1 + datetime.timedelta(minutes=45),
-    sat_id=sat_id
-)
+    dt_fr=time_c - datetime.timedelta(minutes=45),
+    dt_to=time_c + datetime.timedelta(minutes=45),
+    sat_id=sat_id, replace_orbit=True)
 
 dashboard.set_layout(1, 1)
 
 # Get the variables: LBHS emission intensiy, corresponding times and locations
-lbhs = dashboard.assign_variable('GRID_AUR_' + band, dataset_index=0)
-dts = dashboard.assign_variable('DATETIME', dataset_index=0).value.flatten()
-mlat = dashboard.assign_variable('GRID_MLAT', dataset_index=0).value
-mlon = dashboard.assign_variable('GRID_MLON', dataset_index=0).value
-mlt = dashboard.assign_variable(('GRID_MLT'), dataset_index=0).value
+lbhs = ds_ssusi['GRID_AUR_' + band]
+dts = ds_ssusi['DATETIME'].flatten()
+mlat = ds_ssusi['GRID_MLAT']
+mlon = ds_ssusi['GRID_MLON']
+mlt = ds_ssusi['GRID_MLT']
 
 # Search the index for the time to plot, used as an input to the following polar map
-ind_t = dashboard.datasets[0].get_time_ind(ut=time1)
+ind_t = dashboard.datasets[0].get_time_ind(ut=time_c)
+if (dts[ind_t] - time_c).total_seconds()/60 > 60:     # in minutes
+    raise ValueError("The time does not match any SSUSI data!")
 lbhs_ = lbhs.value[ind_t]
-mlat_ = mlat[ind_t]
-mlon_ = mlon[ind_t]
-mlt_ = mlt[ind_t]
+mlat_ = mlat.value[ind_t]
+mlon_ = mlon.value[ind_t]
+mlt_ = mlt.value[ind_t]
 # Add a polar map panel to the dashboard. Currently the style is the fixed MLT at mlt_c=0. See the keywords below:
-panel1 = dashboard.add_polar_map(row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM', mlt_c=0., pole=pole, ut=time1, boundary_lat=65., mirror_south=True)
+panel = dashboard.add_polar_map(
+    row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM',
+    mlt_c=0., pole=pole, ut=time_c, boundary_lat=55., mirror_south=True
+)
 
 # Some settings for plotting.
 pcolormesh_config = lbhs.visual.plot_config.pcolormesh
 # Overlay the SSUSI image in the map.
-ipm = panel1.overlay_pcolormesh(data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM',
-                                regridding=True, **pcolormesh_config)
+ipc = panel.overlay_pcolormesh(
+    data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM', 
+    regridding=False, **pcolormesh_config)
 # Add a color bar
-panel1.add_colorbar(ipm, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
+panel.add_colorbar(ipc, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
                     width=0.05, height=0.7)
 
 # Overlay the gridlines
-panel1.overlay_gridlines(lat_res=5, lon_label_separator=5)
+panel.overlay_gridlines(lat_res=5, lon_label_separator=5)
 
 # Overlay the coastlines in the AACGM coordinate
-panel1.overlay_coastlines()
+panel.overlay_coastlines()
 
 # Overlay cross-track velocity along satellite trajectory
 sc_dt = ds_s1['SC_DATETIME'].value.flatten()
@@ -448,14 +445,22 @@ sc_alt = ds_s1['SC_GEO_ALT'].value.flatten()
 sc_coords = {'lat': sc_lat, 'lon': sc_lon, 'height': sc_alt}
 
 v_H = ds_s1['v_i_H'].value.flatten()
-panel1.overlay_cross_track_vector(vector=v_H, unit_vector=1000, alpha=0.5, color='r', sc_coords=sc_coords, sc_ut=sc_dt, cs='GEOC')
+panel.overlay_cross_track_vector(
+    vector=v_H, unit_vector=1000, vector_unit='m/s', alpha=0.3, color='red',
+    sc_coords=sc_coords, sc_ut=sc_dt, cs='GEO',
+)
 # Overlay the satellite trajectory with ticks
-panel1.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEOC')
+panel.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEO')
+
+# Overlay sites
+panel.overlay_sites(
+    site_ids=['TRO', 'ESR'], coords={'lat': [69.58, 78.15], 'lon': [19.23, 16.02], 'height': 0.}, 
+    cs='GEO', marker='^', markersize=2)
 
 # Add the title and save the figure
 polestr = 'North' if pole == 'N' else 'South'
-panel1.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time1.strftime('%Y-%m-%d %H%M UT'))
-plt.savefig('DMSP_SSUSI_' + time1.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
+panel.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time_c.strftime('%Y-%m-%d %H%M UT'))
+plt.savefig('DMSP_SSUSI_' + time_c.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
 
 # show the figure
 plt.show()

{geospacelab-0.10.1 → geospacelab-0.10.2}/README.md

@@ -322,14 +322,9 @@ Output:
 From then on, the data source has been switched from JHUAPL to CDAWeb.
  
 ```python
-import datetime
-import matplotlib.pyplot as plt
-
-import geospacelab.visualization.mpl.geomap.geodashboards as geomap
-
 dt_fr = datetime.datetime(2015, 9, 8, 8)
 dt_to = datetime.datetime(2015, 9, 8, 23, 59)
-time1 = datetime.datetime(2015, 9, 8, 20, 21)
+time_c = datetime.datetime(2015, 9, 8, 20, 21)
 pole = 'N'
 sat_id = 'f16'
 band = 'LBHS'
@@ -340,49 +335,51 @@ dashboard = geomap.GeoDashboard(dt_fr=dt_fr, dt_to=dt_to, figure_config={'figsiz
 # If the orbit_id is specified, only one file will be downloaded. This option saves the downloading time.
 # dashboard.dock(datasource_contents=['jhuapl', 'dmsp', 'ssusi', 'edraur'], pole='N', sat_id='f17', orbit_id='46863')
 # If not specified, the data during the whole day will be downloaded.
-dashboard.dock(
-    datasource_contents=dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], 
-    pole=pole, sat_id=sat_id, orbit_id=None
-)
+ds_ssusi = dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], pole=pole, sat_id=sat_id, orbit_id=None)
 ds_s1 = dashboard.dock(
     datasource_contents=['madrigal', 'satellites', 'dmsp', 's1'],
-    dt_fr=time1 - datetime.timedelta(minutes=45),
-    dt_to=time1 + datetime.timedelta(minutes=45),
-    sat_id=sat_id
-)
+    dt_fr=time_c - datetime.timedelta(minutes=45),
+    dt_to=time_c + datetime.timedelta(minutes=45),
+    sat_id=sat_id, replace_orbit=True)
 
 dashboard.set_layout(1, 1)
 
 # Get the variables: LBHS emission intensiy, corresponding times and locations
-lbhs = dashboard.assign_variable('GRID_AUR_' + band, dataset_index=0)
-dts = dashboard.assign_variable('DATETIME', dataset_index=0).value.flatten()
-mlat = dashboard.assign_variable('GRID_MLAT', dataset_index=0).value
-mlon = dashboard.assign_variable('GRID_MLON', dataset_index=0).value
-mlt = dashboard.assign_variable(('GRID_MLT'), dataset_index=0).value
+lbhs = ds_ssusi['GRID_AUR_' + band]
+dts = ds_ssusi['DATETIME'].flatten()
+mlat = ds_ssusi['GRID_MLAT']
+mlon = ds_ssusi['GRID_MLON']
+mlt = ds_ssusi['GRID_MLT']
 
 # Search the index for the time to plot, used as an input to the following polar map
-ind_t = dashboard.datasets[0].get_time_ind(ut=time1)
+ind_t = dashboard.datasets[0].get_time_ind(ut=time_c)
+if (dts[ind_t] - time_c).total_seconds()/60 > 60:     # in minutes
+    raise ValueError("The time does not match any SSUSI data!")
 lbhs_ = lbhs.value[ind_t]
-mlat_ = mlat[ind_t]
-mlon_ = mlon[ind_t]
-mlt_ = mlt[ind_t]
+mlat_ = mlat.value[ind_t]
+mlon_ = mlon.value[ind_t]
+mlt_ = mlt.value[ind_t]
 # Add a polar map panel to the dashboard. Currently the style is the fixed MLT at mlt_c=0. See the keywords below:
-panel1 = dashboard.add_polar_map(row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM', mlt_c=0., pole=pole, ut=time1, boundary_lat=65., mirror_south=True)
+panel = dashboard.add_polar_map(
+    row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM',
+    mlt_c=0., pole=pole, ut=time_c, boundary_lat=55., mirror_south=True
+)
 
 # Some settings for plotting.
 pcolormesh_config = lbhs.visual.plot_config.pcolormesh
 # Overlay the SSUSI image in the map.
-ipm = panel1.overlay_pcolormesh(data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM',
-                                regridding=True, **pcolormesh_config)
+ipc = panel.overlay_pcolormesh(
+    data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM', 
+    regridding=False, **pcolormesh_config)
 # Add a color bar
-panel1.add_colorbar(ipm, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
+panel.add_colorbar(ipc, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
                     width=0.05, height=0.7)
 
 # Overlay the gridlines
-panel1.overlay_gridlines(lat_res=5, lon_label_separator=5)
+panel.overlay_gridlines(lat_res=5, lon_label_separator=5)
 
 # Overlay the coastlines in the AACGM coordinate
-panel1.overlay_coastlines()
+panel.overlay_coastlines()
 
 # Overlay cross-track velocity along satellite trajectory
 sc_dt = ds_s1['SC_DATETIME'].value.flatten()
@@ -392,14 +389,22 @@ sc_alt = ds_s1['SC_GEO_ALT'].value.flatten()
 sc_coords = {'lat': sc_lat, 'lon': sc_lon, 'height': sc_alt}
 
 v_H = ds_s1['v_i_H'].value.flatten()
-panel1.overlay_cross_track_vector(vector=v_H, unit_vector=1000, alpha=0.5, color='r', sc_coords=sc_coords, sc_ut=sc_dt, cs='GEOC')
+panel.overlay_cross_track_vector(
+    vector=v_H, unit_vector=1000, vector_unit='m/s', alpha=0.3, color='red',
+    sc_coords=sc_coords, sc_ut=sc_dt, cs='GEO',
+)
 # Overlay the satellite trajectory with ticks
-panel1.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEOC')
+panel.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEO')
+
+# Overlay sites
+panel.overlay_sites(
+    site_ids=['TRO', 'ESR'], coords={'lat': [69.58, 78.15], 'lon': [19.23, 16.02], 'height': 0.}, 
+    cs='GEO', marker='^', markersize=2)
 
 # Add the title and save the figure
 polestr = 'North' if pole == 'N' else 'South'
-panel1.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time1.strftime('%Y-%m-%d %H%M UT'))
-plt.savefig('DMSP_SSUSI_' + time1.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
+panel.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time_c.strftime('%Y-%m-%d %H%M UT'))
+plt.savefig('DMSP_SSUSI_' + time_c.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
 
 # show the figure
 plt.show()

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/__init__.py

@@ -6,7 +6,7 @@ __author__ = "Lei Cai"
 __copyright__ = "Copyright 2021, GeospaceLAB"
 __credits__ = ["Lei Cai"]
 __license__ = "BSD-3-Clause License"
-__version__ = "0.10.1"
+__version__ = "0.10.2"
 __maintainer__ = "Lei Cai"
 __email__ = "lei.cai@oulu.fi"
 __status__ = "Developing"

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/isr/eiscat/loader.py

@@ -133,13 +133,14 @@ class Loader:
                             print("Note: the number of range gates doesn't match nrec!")
                         var = var.reshape(num_row, num_col)
                         if self.gate_num is None:
-                            self.gate_num = num_col
+                            num_gates=self.gate_num = num_col
                         else:
-                            var_array = np.empty((num_row, num_gates))
-                            var_array[::] = np.nan
-                            for i in range(num_row):
-                                var_array[i, 0:num_col] = var[i, :]
-                            var = var_array
+                            num_gates = self.gate_num
+                        var_array = np.empty((num_row, num_gates))
+                        var_array[::] = np.nan
+                        for i in range(num_row):
+                            var_array[i, 0:num_col] = var[i, :]
+                        var = var_array
                     elif nrec_group == 'par1d':
                         if self.gate_num is None:
                             self.gate_num = num_gates = int(np.max(nrec))

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/isr/pfisr/fitted/__init__.py

@@ -46,7 +46,7 @@ default_variable_names = [
     'DATETIME', 'AZ', 'EL', 'PULSE_LENGTH',
     'P_Tx', 'n_e', 'n_e_err', 'T_i', 'T_i_err', 'T_e', 'T_e_err',
     'v_i_los', 'v_i_los_err', 'comp_mix', 'comp_mix_err',
-    'HEIGHT', 'RANGE', 'CGM_LAT', 'CGM_LON'
+    'HEIGHT', 'RANGE', 'CGM_LAT', 'CGM_LON', 'BEAM_ID', 'CHISQ'
 ]
 
 # default_data_search_recursive = True
@@ -55,7 +55,7 @@ default_attrs_required = []
 
 pulse_code_dict = {
     'alternating code': 'AC',
-    'long pulse': 'PL',
+    'long pulse': 'LP',
 }
 
 
@@ -138,8 +138,8 @@ class Dataset(datahub.DatasetSourced):
         self.beam_az = load_obj.beam_az
         self.beam_el = load_obj.beam_el
 
-        if self.add_APEX or self.add_AACGM:
-            self.calc_lat_lon()
+        # if self.add_APEX or self.add_AACGM:
+        self.calc_lat_lon()
 
         if self['HEIGHT'].value is None:
             self['HEIGHT'] = self['GEO_ALT']
@@ -353,7 +353,7 @@ class Dataset(datahub.DatasetSourced):
             pulse_code=self.pulse_code,
             dry_run=dry_run,
             data_file_root_dir=self.data_root_dir,
-            include_exp_ids=self.exp_ids,
+            include_exp_ids=self.experiment_ids,
             include_exp_name_patterns=include_exp_name_patterns)
         return download_obj.data_file_paths
 

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/isr/pfisr/fitted/downloader.py

@@ -15,7 +15,7 @@ EXCLUDE_FILE_TYPE_PATTERNS = [['from power'], ['velocity'], ['uncorrected']]
 
 pulse_code_dict = {
     'alternating code': 'AC',
-    'long pulse': 'PL',
+    'long pulse': 'LP',
 }
 
 

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/isr/pfisr/fitted/loader.py

@@ -37,6 +37,8 @@ var_name_dict = {
     'CGM_LAT': 'cgm_lat',
     'CGM_LON': 'cgm_long',
     'HEIGHT': 'gdalt',
+    'BEAM_ID': 'beamid',
+    'CHISQ': 'chisq',
 }
 
 

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/__init__.py

@@ -205,8 +205,8 @@ class Dataset(datahub.DatasetSourced):
         sin_glon_1 = np.sin(glon_1 * factor)
         sin_glon_2 = np.sin(glon_2 * factor)
         cos_glon_2 = np.cos(glon_2 * factor)
-        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
-        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
+        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
+        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
         sin_glon_2_i = itpf_sin(sectime_1)
         sin_glon_2_i = np.where(sin_glon_2_i > 1., 1., sin_glon_2_i)
         sin_glon_2_i = np.where(sin_glon_2_i < -1., -1., sin_glon_2_i)

{geospacelab-0.10.1 → geospacelab-0.10.2}/geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/__init__.py

@@ -177,8 +177,8 @@ class Dataset(datahub.DatasetSourced):
         sin_glon_1 = np.sin(glon_1 * factor)
         sin_glon_2 = np.sin(glon_2 * factor)
         cos_glon_2 = np.cos(glon_2 * factor)
-        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
-        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
+        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
+        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
         sin_glon_2_i = itpf_sin(sectime_1)
         cos_glon_2_i = itpf_cos(sectime_1)
         rad = np.sign(sin_glon_2_i) * (np.pi / 2 - np.arcsin(cos_glon_2_i))

{geospacelab-0.10.1 → geospacelab-0.10.2/geospacelab.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: geospacelab
-Version: 0.10.1
+Version: 0.10.2
 Summary: Collect, manage, and visualize geospace data.
 Home-page: https://github.com/JouleCai/geospacelab
 Author: Lei Cai
@@ -378,14 +378,9 @@ Output:
 From then on, the data source has been switched from JHUAPL to CDAWeb.
  
 ```python
-import datetime
-import matplotlib.pyplot as plt
-
-import geospacelab.visualization.mpl.geomap.geodashboards as geomap
-
 dt_fr = datetime.datetime(2015, 9, 8, 8)
 dt_to = datetime.datetime(2015, 9, 8, 23, 59)
-time1 = datetime.datetime(2015, 9, 8, 20, 21)
+time_c = datetime.datetime(2015, 9, 8, 20, 21)
 pole = 'N'
 sat_id = 'f16'
 band = 'LBHS'
@@ -396,49 +391,51 @@ dashboard = geomap.GeoDashboard(dt_fr=dt_fr, dt_to=dt_to, figure_config={'figsiz
 # If the orbit_id is specified, only one file will be downloaded. This option saves the downloading time.
 # dashboard.dock(datasource_contents=['jhuapl', 'dmsp', 'ssusi', 'edraur'], pole='N', sat_id='f17', orbit_id='46863')
 # If not specified, the data during the whole day will be downloaded.
-dashboard.dock(
-    datasource_contents=dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], 
-    pole=pole, sat_id=sat_id, orbit_id=None
-)
+ds_ssusi = dashboard.dock(datasource_contents=['cdaweb', 'dmsp', 'ssusi', 'edr_aur'], pole=pole, sat_id=sat_id, orbit_id=None)
 ds_s1 = dashboard.dock(
     datasource_contents=['madrigal', 'satellites', 'dmsp', 's1'],
-    dt_fr=time1 - datetime.timedelta(minutes=45),
-    dt_to=time1 + datetime.timedelta(minutes=45),
-    sat_id=sat_id
-)
+    dt_fr=time_c - datetime.timedelta(minutes=45),
+    dt_to=time_c + datetime.timedelta(minutes=45),
+    sat_id=sat_id, replace_orbit=True)
 
 dashboard.set_layout(1, 1)
 
 # Get the variables: LBHS emission intensiy, corresponding times and locations
-lbhs = dashboard.assign_variable('GRID_AUR_' + band, dataset_index=0)
-dts = dashboard.assign_variable('DATETIME', dataset_index=0).value.flatten()
-mlat = dashboard.assign_variable('GRID_MLAT', dataset_index=0).value
-mlon = dashboard.assign_variable('GRID_MLON', dataset_index=0).value
-mlt = dashboard.assign_variable(('GRID_MLT'), dataset_index=0).value
+lbhs = ds_ssusi['GRID_AUR_' + band]
+dts = ds_ssusi['DATETIME'].flatten()
+mlat = ds_ssusi['GRID_MLAT']
+mlon = ds_ssusi['GRID_MLON']
+mlt = ds_ssusi['GRID_MLT']
 
 # Search the index for the time to plot, used as an input to the following polar map
-ind_t = dashboard.datasets[0].get_time_ind(ut=time1)
+ind_t = dashboard.datasets[0].get_time_ind(ut=time_c)
+if (dts[ind_t] - time_c).total_seconds()/60 > 60:     # in minutes
+    raise ValueError("The time does not match any SSUSI data!")
 lbhs_ = lbhs.value[ind_t]
-mlat_ = mlat[ind_t]
-mlon_ = mlon[ind_t]
-mlt_ = mlt[ind_t]
+mlat_ = mlat.value[ind_t]
+mlon_ = mlon.value[ind_t]
+mlt_ = mlt.value[ind_t]
 # Add a polar map panel to the dashboard. Currently the style is the fixed MLT at mlt_c=0. See the keywords below:
-panel1 = dashboard.add_polar_map(row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM', mlt_c=0., pole=pole, ut=time1, boundary_lat=65., mirror_south=True)
+panel = dashboard.add_polar_map(
+    row_ind=0, col_ind=0, style='mlt-fixed', cs='AACGM',
+    mlt_c=0., pole=pole, ut=time_c, boundary_lat=55., mirror_south=True
+)
 
 # Some settings for plotting.
 pcolormesh_config = lbhs.visual.plot_config.pcolormesh
 # Overlay the SSUSI image in the map.
-ipm = panel1.overlay_pcolormesh(data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM',
-                                regridding=True, **pcolormesh_config)
+ipc = panel.overlay_pcolormesh(
+    data=lbhs_, coords={'lat': mlat_, 'lon': mlon_, 'mlt': mlt_}, cs='AACGM', 
+    regridding=False, **pcolormesh_config)
 # Add a color bar
-panel1.add_colorbar(ipm, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
+panel.add_colorbar(ipc, c_label=band + " (R)", c_scale=pcolormesh_config['c_scale'], left=1.1, bottom=0.1,
                     width=0.05, height=0.7)
 
 # Overlay the gridlines
-panel1.overlay_gridlines(lat_res=5, lon_label_separator=5)
+panel.overlay_gridlines(lat_res=5, lon_label_separator=5)
 
 # Overlay the coastlines in the AACGM coordinate
-panel1.overlay_coastlines()
+panel.overlay_coastlines()
 
 # Overlay cross-track velocity along satellite trajectory
 sc_dt = ds_s1['SC_DATETIME'].value.flatten()
@@ -448,14 +445,22 @@ sc_alt = ds_s1['SC_GEO_ALT'].value.flatten()
 sc_coords = {'lat': sc_lat, 'lon': sc_lon, 'height': sc_alt}
 
 v_H = ds_s1['v_i_H'].value.flatten()
-panel1.overlay_cross_track_vector(vector=v_H, unit_vector=1000, alpha=0.5, color='r', sc_coords=sc_coords, sc_ut=sc_dt, cs='GEOC')
+panel.overlay_cross_track_vector(
+    vector=v_H, unit_vector=1000, vector_unit='m/s', alpha=0.3, color='red',
+    sc_coords=sc_coords, sc_ut=sc_dt, cs='GEO',
+)
 # Overlay the satellite trajectory with ticks
-panel1.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEOC')
+panel.overlay_sc_trajectory(sc_ut=sc_dt, sc_coords=sc_coords, cs='GEO')
+
+# Overlay sites
+panel.overlay_sites(
+    site_ids=['TRO', 'ESR'], coords={'lat': [69.58, 78.15], 'lon': [19.23, 16.02], 'height': 0.}, 
+    cs='GEO', marker='^', markersize=2)
 
 # Add the title and save the figure
 polestr = 'North' if pole == 'N' else 'South'
-panel1.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time1.strftime('%Y-%m-%d %H%M UT'))
-plt.savefig('DMSP_SSUSI_' + time1.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
+panel.add_title(title='DMSP/SSUSI, ' + band + ', ' + sat_id.upper() + ', ' + polestr + ', ' + time_c.strftime('%Y-%m-%d %H%M UT'))
+plt.savefig('DMSP_SSUSI_' + time_c.strftime('%Y%m%d-%H%M') + '_' + band + '_' + sat_id.upper() + '_' + pole, dpi=300)
 
 # show the figure
 plt.show()