|
Size: 1276
Comment:
|
← Revision 6 as of 2009-01-27 16:51:59 ⇥
Size: 555
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| #acl AdminGroup:admin,read,write,delete,revert EditorGroup:read,write All:read #format wiki #language en |
|
| Line 6: | Line 10: |
| Example of reading the ice concentration data into python: | [[/ExampleRead | Example of reading the ice concentration data into python]] |
| Line 8: | Line 12: |
| {{{#!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 16: |
= 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
Brightness temperatures
Get data from NSIDC wget ftp://sidads.colorado.edu/pub/DATASETS/brightness-temperatures/polar-stereo/ssmi/north/2008.tar