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 [ ]: