In [25]:
%pylab notebook
import glob
import pandas as pd
import matplotlib.pyplot as plt


Populating the interactive namespace from numpy and matplotlib
WARNING: pylab import has clobbered these variables: ['all']
`%matplotlib` prevents importing * from pylab and numpy

In [27]:
%cd /Users/brodzik/projects/MODICE/data/MASCONS
#%ls


/Users/brodzik/projects/MODICE/data/MASCONS

In [28]:
def read_modice_mascon_area(filename):
    df = pd.read_csv(filename, delim_whitespace=True)
    df.set_index('Year', drop=True, inplace=True, verify_integrity=True)
    return(df)

In [29]:
region_names = ['NHem',
                'Arctic',
                'Alaskag', 
                'Altaigl',
                'Baffing',
                'Ellesme',
                'FrnJLnd',
                'Greenland',
                'HghMtnA',
                'Him+Kar',
                'Iceland',
                'NW_Amer',
                'Nov_Zem',
                'Nthasia',
                'Pam+Kun',
                'Scandin',
                'Sev_Zem',
                'Sib+Kam',
                'Svalbar',
                'Tianshn',
                'Tib+Qil']

In [30]:
all = pd.DataFrame()
for region_name in region_names:
    for nstrikes in [1, 2, 3]:
        list = glob.glob("MODICE.v0.4.%dstrikes.MOD44W.%s.area.txt" % (nstrikes, region_name))
        filename = list[0]
        print("Next: %s" % filename)
        df = read_modice_mascon_area(filename)
        all["%s_%dstrike_MODICE_area_km^2" % (region_name, nstrikes)] = df['MODICE_area_km^2']
        all["%s_%dstrike_MODICE_NS_km^2" % (region_name, nstrikes)] = df['MODICE_NS_km^2']


Next: MODICE.v0.4.1strikes.MOD44W.NHem.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.NHem.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.NHem.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Arctic.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Arctic.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Arctic.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Alaskag.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Alaskag.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Alaskag.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Altaigl.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Altaigl.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Altaigl.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Baffing.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Baffing.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Baffing.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Ellesme.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Ellesme.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Ellesme.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.FrnJLnd.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.FrnJLnd.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.FrnJLnd.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Greenland.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Greenland.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Greenland.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.HghMtnA.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.HghMtnA.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.HghMtnA.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Him+Kar.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Him+Kar.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Him+Kar.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Iceland.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Iceland.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Iceland.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.NW_Amer.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.NW_Amer.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.NW_Amer.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Nov_Zem.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Nov_Zem.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Nov_Zem.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Nthasia.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Nthasia.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Nthasia.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Pam+Kun.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Pam+Kun.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Pam+Kun.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Scandin.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Scandin.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Scandin.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Sev_Zem.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Sev_Zem.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Sev_Zem.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Sib+Kam.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Sib+Kam.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Sib+Kam.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Svalbar.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Svalbar.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Svalbar.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Tianshn.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Tianshn.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Tianshn.area.txt
Next: MODICE.v0.4.1strikes.MOD44W.Tib+Qil.area.txt
Next: MODICE.v0.4.2strikes.MOD44W.Tib+Qil.area.txt
Next: MODICE.v0.4.3strikes.MOD44W.Tib+Qil.area.txt

In [16]:
all


Out[16]:
NHem_1strike_MODICE_area_km^2 NHem_1strike_MODICE_NS_km^2 NHem_2strike_MODICE_area_km^2 NHem_2strike_MODICE_NS_km^2 NHem_3strike_MODICE_area_km^2 NHem_3strike_MODICE_NS_km^2 Arctic_1strike_MODICE_area_km^2 Arctic_1strike_MODICE_NS_km^2 Arctic_2strike_MODICE_area_km^2 Arctic_2strike_MODICE_NS_km^2 ... Tianshn_2strike_MODICE_area_km^2 Tianshn_2strike_MODICE_NS_km^2 Tianshn_3strike_MODICE_area_km^2 Tianshn_3strike_MODICE_NS_km^2 Tib+Qil_1strike_MODICE_area_km^2 Tib+Qil_1strike_MODICE_NS_km^2 Tib+Qil_2strike_MODICE_area_km^2 Tib+Qil_2strike_MODICE_NS_km^2 Tib+Qil_3strike_MODICE_area_km^2 Tib+Qil_3strike_MODICE_NS_km^2
Year
2000 2326273.2 6782.4 2437080.0 7052.8 2545934.5 7174.5 2219722.9 6753.4 2301854.9 7022.3 ... 15558.7 0.0 18209.9 0.0 5972.9 0.0 6958.2 0.0 7666.3 0.0
2001 2271053.6 6737.7 2353716.9 6981.3 2424138.9 7111.4 2179881.4 6709.4 2238214.2 6950.2 ... 12116.0 0.0 14201.0 0.0 5403.4 0.0 6243.1 0.0 6778.1 0.0
2002 2257186.8 6517.9 2343400.4 6748.9 2421538.1 6848.7 2155740.0 6490.0 2214405.4 6720.5 ... 11697.0 0.0 13390.8 0.0 6001.9 0.0 6992.3 0.0 7766.8 0.0
2003 2242640.7 6471.3 2324957.1 6785.1 2395010.8 6997.9 2135469.8 6442.6 2193148.8 6756.0 ... 12122.6 0.0 13794.6 0.0 6469.0 0.0 7229.9 0.0 7738.4 0.0
2004 2266576.4 6605.5 2371329.9 6839.9 2474045.7 6932.4 2155598.5 6577.6 2230702.5 6811.5 ... 11784.5 0.0 13543.5 0.0 6039.2 0.0 6759.2 0.0 7251.4 0.0
2005 2235608.9 6491.5 2314377.5 6828.7 2384378.5 6956.0 2129632.6 6463.2 2182380.0 6799.7 ... 12604.5 0.0 14679.4 0.0 6577.1 0.0 7468.4 0.0 8125.5 0.0
2006 2260105.5 6643.7 2346194.2 6871.0 2419323.9 6981.1 2172806.9 6616.0 2236814.0 6842.2 ... 10021.3 0.0 11581.1 0.0 5380.6 0.0 6291.9 0.0 6861.8 0.0
2007 2238029.0 6384.4 2322117.4 6614.3 2392410.8 6754.4 2147029.4 6356.7 2206923.2 6585.9 ... 9749.4 0.0 11177.7 0.0 5952.1 0.0 6815.4 0.0 7442.2 0.0
2008 2250219.9 6358.6 2333971.9 6577.1 2402240.1 6678.5 2158051.7 6330.9 2216167.7 6549.5 ... 9919.4 0.0 11712.6 0.0 6783.4 0.0 7750.3 0.0 8459.7 0.0
2009 2235194.0 6372.8 2313571.9 6579.9 2375536.3 6665.8 2133774.2 6344.5 2183825.1 6549.9 ... 16410.2 0.0 18984.2 0.0 6382.7 0.0 7288.7 0.0 7973.9 0.0
2010 2178647.1 6373.2 2270225.2 6608.9 2329746.0 6720.7 2082241.7 6345.5 2141383.8 6581.2 ... 11639.0 0.0 13670.5 0.0 4762.0 0.0 5806.9 0.0 6348.1 0.0
2011 2189851.0 6314.2 2259022.8 6517.9 2311196.9 6610.4 2106700.6 6286.5 2153282.6 6490.2 ... 11043.5 0.0 12737.4 0.0 5808.7 0.0 6636.4 0.0 7169.6 0.0
2012 2223326.6 6432.7 2301822.3 6640.5 2361511.6 6750.2 2120971.3 6405.0 2173795.6 6612.8 ... 9435.3 0.0 11020.8 0.0 6180.2 0.0 6947.6 0.0 7487.7 0.0
2013 2215347.1 6382.4 2291580.5 6617.9 2355050.4 6717.5 2133591.5 6350.0 2191479.2 6585.5 ... 9758.0 0.0 11125.5 0.0 4922.8 0.0 5656.0 0.0 6104.7 0.0
2014 2232360.5 885.5 2315498.9 923.7 2385634.7 936.1 2140724.6 855.0 2201786.0 893.0 ... 9038.4 0.0 10551.3 0.0 5225.7 0.0 6039.6 0.0 6571.8 0.0
2015 2216014.4 880.7 2282289.4 919.6 2334466.5 929.7 2129312.1 850.7 2176819.7 889.5 ... 8518.3 0.0 10417.6 0.0 5426.8 0.0 6178.1 0.0 6712.2 0.0

16 rows × 126 columns


In [17]:
%pwd
#all.to_csv('mascon_areas_by_year.csv')


Out[17]:
u'/Users/brodzik/projects/MODICE/data/MASCONS'

In [18]:
all.to_csv('mascon_areas_by_year.v02.csv')

In [33]:
HMA1 = all['NHem_1strike_MODICE_area_km^2'] - all['Arctic_1strike_MODICE_area_km^2']
HMA2 = all['NHem_2strike_MODICE_area_km^2'] - all['Arctic_2strike_MODICE_area_km^2']
HMA3 = all['NHem_3strike_MODICE_area_km^2'] - all['Arctic_3strike_MODICE_area_km^2']

In [48]:
fig, ax = plt.subplots(2,1, figsize=(9,6))
HMA1.plot(ax=ax[0], label='1strike')
HMA2.plot(ax=ax[0], label='2strike')
HMA3.plot(ax=ax[0], label='3strike')
ax[0].set_title("MODICE in High Mountain Asia (12 tiles)")
ax[0].set_ylabel("Area ($km^2$)")
ax[0].legend()
all['Arctic_1strike_MODICE_area_km^2'].plot(ax=ax[1], label='1strike')
all['Arctic_2strike_MODICE_area_km^2'].plot(ax=ax[1], label='2strike')
all['Arctic_3strike_MODICE_area_km^2'].plot(ax=ax[1], label='3strike')
ax[1].set_title("All Arctic MODICE (46 tiles)")
ax[1].set_ylabel("Area ($km^2$)")
ax[1].legend()
fig.tight_layout()
fig.savefig("MODICE.2000-2015.png")



In [45]:



Out[45]:
<matplotlib.legend.Legend at 0x11ec8fd10>

In [ ]: