Quick view of netCDF file (param/lat/lon type)

nc-quick_view

import gdal
import netCDF4
import sys
import matplotlib.pyplot as plt


def plot(file_path,  title=""):
    gdal_data = gdal.Open(file_path)
    arr = gdal_data.ReadAsArray()
    #print "shapes", arr.shape, arr.dtype
    plt.imshow(arr)
    plt.colorbar()
    plt.title(title)
    plt.show()

var_name = 'airtemp'
file_path_out = '/path/to/%s/daily/' % var_name
root_name = 'ecmwf_2mtemperature_fulldaily_0.125deg_-22S16E-35S33E_'
year = 2001
day = 181

file_name = root_name + str(year) + '_avg_%03d.nc' % day
data_path = file_path_out + file_name
plot(data_path, file_name)

file_name = root_name + str(year) + '_avg_%03d_studyarea.nc' % day
data_path = file_path_out + file_name
plot(data_path, file_name)