Pysat seasons demo 1

[rstoneback]

Demonstration of pysatSeasons seasonal averaging. Uses ICON IVM to generate an average temperature and density distribution for 'equatorial-ish' ionosphere. I have a custom function for adding log density. I don't do any spatial restrictions on the data.

[rstoneback]

Binder link https://mybinder.org/v2/gh/pysat/pysatTutorials/pysatSeasons-demo-1

[jklenzing]

Block 2 of code registers all instruments. This pops up the warning from GOLD. perhaps there should be a comment explaining that this comes from preliminary GOLD support?
Just saw this already exists in the comments

[jklenzing]

Add newline after printing "subset of data" for readability of output

# Slicing and dicing mini-examples
print('Subset of data')
print(ivm[0:5, 'Ion_Density'])
print('')
date1 = dt.datetime(2020, 1, 1)
date2 = dt.datetime(2020, 1, 1, 1, 30, 30)
print(ivm[date1:date2, 'Ion_Density'])
print('')

[jklenzing]

---

AttributeError Traceback (most recent call last)
/tmp/ipykernel_61/2148985827.py in
54 livm.data.plot(x='Magnetic_Local_Time', y='Log_Ion_Density',
55 xlim=(0,24), xticks=[],
---> 56 ax=axs[0], legend=False)
57 axs[0].set_title('ICON IVM Log Ion Density')
58 axs[0].set_ylabel('Log Ion Density (N/cc)')

AttributeError: 'AxesSubplot' object has no attribute 'rowNum'

[jklenzing]

MIGHTI wind plots should have "red line" and "green line" in titles for clarity

[jklenzing]

MIHTI Temperature block:

ValueError                                Traceback (most recent call last)
/tmp/ipykernel_61/110224154.py in <module>
     54     # Make plot
     55     a.transpose().plot(ax=axs[i], x='average_local_time', y='Tangent_Altitude',
---> 56                        cbar_kwargs={'label': 'Temperature (K)'})
     57     axs[i].set_xlim((0, 24))
     58     axs[i].set_xticks([])

ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (247,) and requested shape (247,18)

[jklenzing]

When the colorbar plots start, I'd suggest having at least one with specified max and min for the colorbar range, since this is an intro python tutorial

[rstoneback]

I get a 30 minute runtime (on Binder) for the ICON seasonal analysis. A bit long.... :(

[jklenzing]

I've just been poking buttons on binder

[rstoneback]

I get a 30 minute runtime (on Binder) for the ICON seasonal analysis. A bit long.... :(

Takes about 5 minutes on my local system. I'm currently in the middle of a big computation taking so I thought that was slow.... Gonna check on heliocloud as well.

Issue with reducing the number of days is the plots don't come out complete. Perhaps I need to a bit of text as well.

[rstoneback]

Only takes 90 seconds for 15 days on heliocloud. That works. 31 days cause ram use to exceed 4GB and the kernel restarts.

[rstoneback]

MIHTI Temperature block:

ValueError                                Traceback (most recent call last)
/tmp/ipykernel_61/110224154.py in <module>
     54     # Make plot
     55     a.transpose().plot(ax=axs[i], x='average_local_time', y='Tangent_Altitude',
---> 56                        cbar_kwargs={'label': 'Temperature (K)'})
     57     axs[i].set_xlim((0, 24))
     58     axs[i].set_xticks([])

ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (247,) and requested shape (247,18)

Thanks! This is fixed in the summer-school repo I just haven't brought it back here yet.

[rstoneback]

MIHTI Temperature block:

ValueError                                Traceback (most recent call last)
/tmp/ipykernel_61/110224154.py in <module>
     54     # Make plot
     55     a.transpose().plot(ax=axs[i], x='average_local_time', y='Tangent_Altitude',
---> 56                        cbar_kwargs={'label': 'Temperature (K)'})
     57     axs[i].set_xlim((0, 24))
     58     axs[i].set_xticks([])

ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (247,) and requested shape (247,18)

Thanks! This is fixed in the summer-school repo I just haven't brought it back here yet.

This is the result of 'Epoch' switching to 'time' in the data. Updated.

[rstoneback]

Actual merge to main pending release of pysat and such. Using review here as proxy for release to summer school.