source: trunk/src/scripts_Laura/read_netcdf_AMSUA_test.py @ 26

#!/usr/bin/env python
# -*- coding: utf-8 -*-
import string
import numpy as np
import matplotlib.pyplot as plt
from pylab import *
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.basemap import shiftgrid, cm
from netCDF4 import Dataset


#####################
## read daily data ##
#####################
month = array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
len_month = np.array([31, 28, 31, 30, 31, 30, 31])

## JANUARY ##
imo = 0
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_JAN = ncfile.variables['day'][:]
emis_JAN = ncfile.variables['emissivity'][:]
ncfile.close()

## FEBRUARY ##
imo = 1
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_FEB = ncfile.variables['day'][:]
emis_FEB = ncfile.variables['emissivity'][:]
ncfile.close()

## MARCH ##
imo = 2
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_MAR = ncfile.variables['day'][:]
emis_MAR = ncfile.variables['emissivity'][:]
ncfile.close()

## APRIL ##
imo = 3
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_APR = ncfile.variables['day'][:]
emis_APR = ncfile.variables['emissivity'][:]
ncfile.close()

## MAY ##
imo = 4
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_MAY = ncfile.variables['day'][:]
emis_MAY = ncfile.variables['emissivity'][:]
ncfile.close()

## JUNE ##
imo = 5
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_JUN = ncfile.variables['day'][:]
emis_JUN = ncfile.variables['emissivity'][:]
ncfile.close()

## JULY ##
imo = 6
ncfile = Dataset('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] + '_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
jrs_JUL = ncfile.variables['day'][:]
emis_JUL = ncfile.variables['emissivity'][:]
ncfile.close()





#######################
## read monthly data ##
#######################

ncfile = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_JAN-MAR-MAY-JUN_ANTARC.nc', 'r', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
mon = ncfile.variables['month'][:]
emis_moy1 = ncfile.variables['emissivity'][:]
ncfile.close()


ncfile = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_FEB-APR-JUL_ANTARC.nc', 'w', format = 'NETCDF4')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
chan = ncfile.variables['channels'][:]
mon = ncfile.variables['month'][:]
emis_moy2 = ncfile.variables['emissivity'][:]
ncfile.close()