In [1]:

    

%pylab inline
import bmi.wrapper
import matplotlib.pyplot as plt


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

model = bmi.wrapper.BMIWrapper(engine='subgrid', configfile='/Users/baart_f/models/duifp/Duifp_default.mdu')


    

In [3]:

    

model.initialize()
model.update(-1)
vars = {}
for var in ['s1', 'maxinterception']:
    vars[var] = model.get_var(var)


    

In [8]:

    

plt.imshow(vars['maxinterception'][::10,::10], vmin=0, vmax=0.1)


    

    
Out[8]:





<matplotlib.image.AxesImage at 0x119e4eb50>