Re: ssmi data

NCL reads the available data correctly.

I speculate the SSMI files do not contain any information
on the grid latitudes or longitudes. Hence, when you try
to plot the data on a map, it gives the errors below.

I will send you a file that should have the lat/lon information
that will allow you to plot the data.

On 9/26/12 10:25 AM, Alexander Semenov wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 26 12:05:08 2012