1 from PyNGL import Nio
   2 grid = Nio.open_file('antarctic_psi_conc_antarctic_ssmi_daily_19911205_20080821_concOnly_runmedian5.nc')
   3 C=grid.variables['concentration'][0,:,:] # Day zero

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