Differences between revisions 3 and 5 (spanning 2 versions)
Revision 3 as of 2008-09-11 10:41:58
Size: 1276
Comment:
Revision 5 as of 2008-09-23 14:51:16
Size: 445
Comment:
Deletions are marked like this. Additions are marked like this.
Line 6: Line 6:
Example of reading the ice concentration data into python: [[/ExampleRead | Example of reading the ice concentration data into python]]
Line 8: 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()
}}}
Line 36: Line 12:


= Brightness temperatures =

Get data from NSIDC
{{{wget ftp://sidads.colorado.edu/pub/DATASETS/brightness-temperatures/polar-stereo/ssmi/north/2008.tar}}}

SSM/I ASI

ftp://ftp-projects.zmaw.de/seaice/SSMI_ASI_IceConc/

Example of reading the ice concentration data into python

Animation

ftp://ftp.zmaw.de/outgoing/carmen/Arctic_corrMedian5_19920101_20080821_15000f_FlyHigher_comp5_Npole100.mpg

Brightness temperatures

Get data from NSIDC wget ftp://sidads.colorado.edu/pub/DATASETS/brightness-temperatures/polar-stereo/ssmi/north/2008.tar

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