{{{#!python
from PyNGL import Nio
grid = Nio.open_file('antarctic_psi_conc_antarctic_ssmi_daily_19911205_20080821_concOnly_runmedian5.nc')
C=grid.variables['concentration'][0,:,:] # Day zero
}}}
{{{#!python
from PyNGL import Nio
from scipy import *
from pylab import *
import os,sys,time
arg=sys.argv[1]  #input date
year0,month0,day0,hour0=int(arg[0:4]),int(arg[4:6]),int(arg[6:8]),int(arg[8:10])
t0=805812.0 #the hour of 1991.12.5 00:00 since 1900-01-01 00:00
lt0=time.mktime((1991,12,5,0,0,0,0,0,0))
lt2=time.mktime((year0,month0,day0,0,0,0,0,0,0))
t2=t0+(lt2-lt0)/3600.0  #the hour of input date since 1900-01-01 00:00
#grid = Nio.open_file('antarctic_psi_conc_antarctic_ssmi_daily_19911205_20080821_concOnly_runmedian5.nc')
grid=Nio.open_file('arctic_psi_conc_daily_19911205_20080821_runmedian5_Npole100.nc')
t=grid.variables['time'][:] #time series of the data
index=abs(t-t2).argmin() #find the index of input date
C=grid.variables['concentration'][index,:,:] 
imshow(C)
show()
}}}