Jupyter notebook atms391geodata/Week 12/Homework 12.ipynb
ATMS 391: Geophysical Data Analysis
Homework 12
(1) Recreate Figure 1 from Horel and Wallace (1981) using the NCEP reanalysis data.
(2) Create normalized index values (mean=0, sd=1) of Sea Level Pressure between Tahiti and Darwin similar to Fig. 4. Invert the index (multiply by -1, as in the paper). Plot the time series.
(3) Create normalized index values of 200 hPa height at the stations listed in Fig 6. Plot the time series.
(4) Compute the PNA and WP indices as on p. 823 of the paper. Plot the time series.
(5) Compute standardized rainfall values at the locations in Fig. 3. Plot the time series.
(6) Reproduce Fig 8 with the top 4 and bottom 4 values of PNA index (that you find in part 4), except contour the values you have rather than plot points as in the paper.
(7) Recreate Figure 9 using your data.
---------------------------------------------------------------------------
IndexError Traceback (most recent call last)
<ipython-input-34-90e6cb8a4272> in <module>()
----> 1 outmap=spatial_corr(fanningprecip_standard.values,h700['hgt'].values)
2 from mpl_toolkits.basemap import Basemap
3 import numpy as np
4 import matplotlib.pyplot as plt
5
<ipython-input-31-76b51efb0c17> in spatial_corr(index, mapvals)
17 df_bis = pd.DataFrame({'maptser': maptser,
18 'index': np.squeeze(index)})
---> 19 outmap[i,j]=df_bis.corr().values[1,0]
20
21 return outmap
IndexError: index 1 is out of bounds for axis 0 with size 1