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