import komod
import Nio
from IPython.display import Image
import numpy as np

!wget http://data.nodc.noaa.gov/thredds/fileServer/woa/WOA09/NetCDFdata/temperature_annual_5deg.nc

ff =Nio.open_file('temperature_annual_5deg.nc')

temp = ff.variables['t_mn'][:]
lat  = ff.variables['lat'][:]
lon  = ff.variables['lon'][:]
lev  = ff.variables['depth'][:]

temp.shape

komod.arctpl(lon, lat, temp[0,0,:,:], region = 'Global')

!convert -scale 50% -rotate -90 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'Global')

!convert -scale 50% -rotate -90 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'Global', raster_fill=True)

!convert -scale 50% -rotate -90 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'Global', raster_fill=True, 
             colormap_name='posneg_1', start_color=4, end_color=-2)

!convert -scale 50% -rotate -90 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'Global', 
             minLon=100, maxLon=290 , minLat=-80 , maxLat=75)

!convert -scale 60% -rotate -0 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'NAtlanticOcean')

!convert -scale 60% -rotate -0 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0,:,:], 
             datamin=-2, datamax=30, datastep=1,
             vtitle='T, deg. C', region = 'NAtlanticOcean', add_cyclic=True)

!convert -scale 60% -rotate -0 output.ps output.png
Image(filename='output.png')

komod.arctpl(lon, lat, temp[0,0:5,:,:], region = 'Global')

Image(filename='earthpy_mitplot1.png')

komod.arctpl(lon, lat, temp[0,0:5,:,:], region = 'Global',
             levon=True, llevel=lev[0:5])

Image(filename='earthpy_mitplot2.png')

komod.arctpltnc('temperature_annual_5deg.nc', 't_mn', 
                lon='lon', lat='lat', llevel="depth", 
                region = 'Global', add_cyclic=True, levon=True)

Image(filename='earthpy_mitplot3.png')