SSM/I ASI

Animation

IfmWiki: SeaiceGroup/SSMI (last edited 2009-01-27 16:51:59 by LarsKaleschke)

-  ⇤ ← Revision 3 as of 2008-09-11 10:41:58 → 
  Size: 1276
  Editor: XiangshanTianKunze
  Comment:
+   ← Revision 4 as of 2008-09-11 10:45:09 → ⇥
  Size: 281
  Editor: LarsKaleschke
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 6:
-Example of reading the ice concentration data into python:
+[[/ExampleRead | Example of reading the ice concentration data into python]]
 Line 8:
-{{{#!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()
}}}