ssmi data

Hi,

I'm trying to plot SSMI sea-ice coverage (file is in hdf format), but I get
the domain with missing values. I think NCL doesn't understand the
dimensions of a variable. Do I need to do some sort of interpolation to
aregular grid? If so, please let me know how.

pacman3 1245% ncl seaice.ncl
 Copyright (C) 1995-2011 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.0.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: ic
Type: float
Total Size: 8716288 bytes
            2179072 values
Number of Dimensions: 2
Dimensions and sizes: [fakeDim0 | 1792] x [fakeDim1 | 1216]
Coordinates:
Number Of Attributes: 3
  valid_range : ( 0, 100 )
  long_name : ASI Ice Concentration, 20120801, res: 6.25000, SSMIS (F18),
ASI Version: 5.2
  hdf_name : ASI Ice Concentration
(0) ASI Ice Concentration, 20120801, res: 6.25000, SSMIS (F18), ASI
Version: 5.2: min=nan max=nan
(0) check_for_y_lat_coord: Warning: Data either does not contain a valid
latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute equal
to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

script:

begin

 files = systemfunc("ls asi-SSMIS-n6250-201208*.hdf")
  do j = 0, dimsizes(files)-1
  fname = files(j)
  f = addfile (fname, "r")

ic=f->ASI_Ice_Concentration
printVarSummary(ic)
printMinMax(ic,0)

-- 
regards
*******************************************************
*Alex Semenov*
*
*
PhD Student - Research Assistant
International Arctic Research Center
Department of atmospheric sciences
University of Alaska Fairbanks
930 Koyukuk dr, 408c3
Fairbanks, AK, USA, 99775
work phone +19074742672

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 26 10:25:26 2012