'''plot end-member fractions'''
from pylab import *
import gemma
import gemma.plotting as gp
import sys, os
import datetime as dt
sys.path.append(os.path.abspath('pyplots'))

g = gemma.gemma('gemma_output.bin', version=34)
g.read_data(100)
g.solutes.sort()

hiatuses = [[dt.datetime(2011,11,6),dt.datetime(2011,11,9)],[dt.datetime(2012,1,31),dt.datetime(2012,2,20)],
            [dt.datetime(2012,5,16),dt.datetime(2012,5,28)],[dt.datetime(2012,8,6),dt.datetime(2012,8,23)]]

fig, ax = gp.endmembers_time((5,50,95), result=g, fracload='load', hiatuses=hiatuses,
                             separate=True, returnfig=True, figsize=(7,10))

show()