Fortran module from Python

from mod_arstats import ar_stats as arstats
import numpy as np
import pylab as pl


a = np.ones((2800,2800)) * np.arange(2800)
s = 2
b = np.ones((a.shape[0]/s, a.shape[1]/s))

arstats.wm = a
arstats.scf = s
arstats.new_wm = b

arstats.mean()

pl.imshow(arstats.new_wm)
pl.show()

exit