import komod

%pylab inline

lon, lat = komod.mitbincoord(416, 480, bswap=1)

subplot(1,2,1)
imshow(lon)
subplot(1,2,2)
imshow(lat)

ndim, xdim, ydim, zdim, datatype, nrecords, timeStepNumber = komod.rmeta('AREAtave.0000002191.meta')

print('Number of dimensions = '+str(ndim))
print('Xdim = '+str(xdim))
print('Ydim = '+str(ydim))
print('Zdim = '+str(zdim))
print('Datatype - '+datatype)
print('Number of records = '+str(nrecords))
print('Time step number = '+str(timeStepNumber))

area = komod.mitbin('AREAtave.0000002191.data', 416, 480, datatype = 'float32')

pcolormesh(area[0,0,:,:])

area = komod.mitbin2('AREAtave.0000002191.data')
temp = komod.mitbin2('Ttave.0000002191.data')

figure(figsize=(10,5))
subplot(1,2,1)
pcolormesh(area[0,0,:,:])
subplot(1,2,2)
pcolormesh(temp[0,5,:,:])

komod.var_nc2d(parameters=['AREAtave'], FillValue=-999.)

!ncdump -h MIT_output_2d.nc

komod.var_nc3d(parameters=['Ttave'], FillValue=-999.)

!ncdump -h MIT_output_3d.nc