Global & Local Temperature

This example uses historical data since 1880 on average global surface temperature changes from NASA'sGISS Surface Temperature Analysis (GISTEMP) (original file) That site has lots of other data sets in csv format, too.

Step 1: Click on the "Save Icon". If you have a problem and need to completely restart, then click "File" and "Revert to Checkpoint" to reset.

Click the "play" icon or press shift+enter to execute each cell.

Pre-Questions

Answer questions 1 & 2 from your coding booklet.


In [ ]:
# import the software packages needed
import pandas as pd
import numpy as np
%matplotlib inline
import matplotlib as mpl
import matplotlib.pyplot as plt
inline_rc = dict(mpl.rcParams)
import datetime as dt
import matplotlib.dates as mdates

Importing a data set

This next cell may take a little while to run if it's grabbing a pretty big data set. The cell label to the left will look like "In [*]" while it's still thinking and "In [2]" when it's finished.


In [ ]:
# Combined land and ocean temperature averages (LOTI: Land Ocean Temperature Index)
data1 = pd.read_csv('http://github.com/scpsscience/data/raw/master/LOTI.csv', header=1).replace(to_replace="***", value=np.NaN)
data_LOTI = data1.apply(lambda x: pd.to_numeric(x, errors='ignore'))

# Only land temperature averages
data2 = pd.read_csv('http://github.com/scpsscience/data/raw/master/LAND.csv', header=1).replace(to_replace="***", value=np.NaN)
data_LAND = data2.apply(lambda x: pd.to_numeric(x, errors='ignore'))

# Temperature Data for Orlando, Florida
data3 = pd.read_csv('http://github.com/scpsscience/data/raw/master/Orlando_Temps.csv')

dates = ['1/1/1975','1/2/1975','01/3/1975']
x = [dt.datetime.strptime(d,'%m/%d/%Y').date() for d in dates]
y = "range(len(x))"

PART 1: Local Temperature Changes & Climate

This section will be used for questions 3-5.

We can view the first few rows of data from the file we just imported.


In [ ]:
# The .head(n) command displays the first n rows of the local temperature data file.
data3.head(4)

In [ ]:
data3.columns

The code below will display a graph displaying the average temperatures at Orlando International Airport. You can double click on the graph to make it bigger.


In [ ]:
ax = data3.plot('date', 'avg', title="Average Temperature at Orlando International Airport (1975 - Present)", legend=False, color='b')
ax.set(xlabel="Date", ylabel="Average Temperature (Fahrenheit)")


# Get current size
fig_size = plt.rcParams["figure.figsize"]

 
# Set figure width to 20 and height to 6.5
fig_size[0] = 100
fig_size[1] = 6.5
plt.rcParams["figure.figsize"] = fig_size

plt.show()

The code below can show you the Highest and Lowest Temperatures Recorded at Orlando International Airport


In [ ]:
#This section of code will show you the lowest temperatures at Orlando International Airport. 
data3.sort_values('low', ascending=True).head(5)

#This section of code will show you the Highest temperatures at Orlando International Airport. 
#Add a # above and remove the one below to see the highest temperatures.

#data3.sort_values('high', ascending=False).head(5)

PART 2: Global Temperature Changes & Climate

This section will be used for questions 6-8.


In [ ]:
# The .head(n) command displays the first n rows of the file.
data_LAND.head(5)

In [ ]:
x1 = data_LOTI.Year
y1 = data_LOTI.JanDec

# plt.plot() makes a line graph, by default
fig = plt.figure(figsize=(10, 5))
plt.plot(x1, y1)
plt.title('Average land and ocean temperature readings')
plt.xlabel('Year')
plt.ylabel('Percent temp change')

In [ ]:
x2 = data_LAND.Year
y2 = data_LAND.JanDec

# plt.plot() makes a line graph, by default
fig = plt.figure(figsize=(10, 5))
plt.plot(x2, y2)
plt.title('Land temperature readings')
plt.xlabel('Year')
plt.ylabel('Percent temp change')

In [ ]:
# Wow, this needs a title and axis labels!
fig = plt.figure(figsize=(10, 5))
plt.plot(x1, y1, label="Land and Ocean")
plt.plot(x2, y2, label="Land only")


plt.legend()
plt.show()

Part 3: Unstructured Coding

Modify the codes for Part 1 above to answer questions 9-12. The graph in Part 1 will be useful to view precipitaion. You can modify the data charts to view the precipation totals as well!

DONE?!?!

Try editing some code and re-running the cell to see the effects. For more information on formatting the markdown text in a cell like this one, go to Help > Markdown > Basic Writing and Formatting Text.

To save your work: go to File > Save and Checkpoint. That only saves your edits as long as you're working inthis notebook. All is lost after 10 minutes of inactivity.

To really save your work: go to File > Download as > iPython notebook (or save as pdf if you just want to show someone a snapshot of what your code and output look like.

References