#!/opt/canopy-1.3.0/Canopy_64bit/User/bin/python


# calcul de la surface du modele peru12 a partir du meshmask

import matplotlib.pyplot as plt
from netCDF4 import Dataset as netcdf

nc=netcdf('mesh_mask_peru12.L75_from_trop12.nc','r')
mask=nc.variables['tmask'][0,0,:,:]
e1t=nc.variables['e1t'][0,:,:]
e2t=nc.variables['e2t'][0,:,:]
nc.close()

surf=e1t*e2t
msurf=surf*mask
area=sum(sum(msurf[:,:]))
print "calcul de la surface du domaine peru12 a partir du mesh_mask = %4.1f Petametres" % (area/1e+12)
plt.pcolor(msurf)
plt.colorbar()
plt.show()