[gradsusr] Data not being displayed

Syed Ahsan Ali ahsanshah01 at gmail.com
Sun Jul 27 03:07:20 EDT 2014


I was trying to open a netcdf file  and there was error SDF file has no
discernible X coordinate. I wrote a grads description file the data file as
below

DSET ^GlobSnow_SWE_L3B_monthly_201301_v2.0.nc
TITLE GlobalSWE
UNDEF 99999.0
XDEF 1637 LINEAR -180 0.22
YDEF 200 LINEAR 35 0.22
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR jan2013 1mo
VARS 2
SWE_avg 1 99
SWE_max 1 99
ENDVARS



The file is opened but data is not being displayed rather some abnormal
values like -8e+38  are shown.

File is present on dropbox link
https://www.dropbox.com/s/fp6rft5g26hej49/GlobSnow_SWE_L3B_monthly_201301_v2.0.nc
. output from ncdump is attached. Please guide me.

Regards
Ahsan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140727/feb44de4/attachment.html 
-------------- next part --------------
[pmdtest at r720 ~]$ ncdump -h GlobSnow_SWE_L3B_monthly_201301_v2.0.nc 
netcdf GlobSnow_SWE_L3B_monthly_201301_v2.0 {
dimensions:
        x = 721 ;
        y = 721 ;
variables:
        char lambert_azimuthal_equal_area ;
                lambert_azimuthal_equal_area:grid_mapping_name = "lambert_azimuthal_equal_area" ;
                lambert_azimuthal_equal_area:false_easting = 0. ;
                lambert_azimuthal_equal_area:false_northing = 0. ;
                lambert_azimuthal_equal_area:latitude_of_projection_origin = 90. ;
                lambert_azimuthal_equal_area:longitude_of_projection_origin = 0. ;
                lambert_azimuthal_equal_area:longitude_of_prime_meridian = 0. ;
                lambert_azimuthal_equal_area:semi_major_axis = 6371228. ;
                lambert_azimuthal_equal_area:inverse_flattening = 0. ;
                lambert_azimuthal_equal_area:spatial_ref = "PROJCS[\"NSIDC EASE-Grid North\",GEOGCS[\"Unspecified datum based upon the International 1924 Authalic Sphere\",DATUM[\"Not_specified_based_on_International_1924_Authalic_Sphere\",SPHEROID[\"International 1924 Authalic Sphere\",6371228,0,AUTHORITY[\"EPSG\",\"7057\"]],AUTHORITY[\"EPSG\",\"6053\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4053\"]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_center\",90],PARAMETER[\"longitude_of_center\",0],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"3408\"]]" ;
                lambert_azimuthal_equal_area:GeoTransform = "-9036842.762 25067.525 0 9036842.763000002 0 -25067.525 " ;
        double x(x) ;
                x:standard_name = "projection_x_coordinate" ;
                x:long_name = "x coordinate of projection" ;
                x:units = "m" ;
        double y(y) ;
                y:standard_name = "projection_y_coordinate" ;
                y:long_name = "y coordinate of projection" ;
                y:units = "m" ;
        float lat(y, x) ;
                lat:long_name = "WGS84 latitude coordinates, center of pixel" ;
                lat:_FillValue = -99999.f ;
                lat:grid_mapping = "lambert_azimuthal_equal_area" ;
                lat:units = "degrees" ;
        float lon(y, x) ;
                lon:long_name = "WGS84 longitude coordinates, center of pixel" ;
                lon:_FillValue = -99999.f ;
                lon:grid_mapping = "lambert_azimuthal_equal_area" ;
                lon:units = "degrees" ;
        float SWE_avg(y, x) ;
                SWE_avg:long_name = "Monthly mean Snow Water Equivalent" ;
                SWE_avg:_FillValue = -99999.f ;
                SWE_avg:grid_mapping = "lambert_azimuthal_equal_area" ;
                SWE_avg:units = "mm" ;
        float SWE_max(y, x) ;
                SWE_max:long_name = "Monthly maximum Snow Water Equivalent" ;
                SWE_max:_FillValue = -99999.f ;
                SWE_max:grid_mapping = "lambert_azimuthal_equal_area" ;
                SWE_max:units = "mm" ;

// global attributes:
                :GDAL_AREA_OR_POINT = "Area" ;
                :Conventions = "CF-1.5" ;
                :GDAL = "GDAL 1.10.1, released 2013/08/26" ;
                :history = "Sat Jan 04 16:29:16 2014: GDAL CreateCopy( D:\\GlobSnow_SWE_processing\\SWE_mat_to_geotiff_nc\\temp\\GlobSnow_SWE_L3B_monthly_201301_v2.0.nc, ... )" ;
                :title = "GLOBSNOW FPS SWE" ;
                :sensor = "SSM/I" ;
                :Data\ Date = "2013-01" ;
                :Processing\ Date = "2014-01-04" ;
                :Coordinate\ System = "Equal-Area Scalable Earth Grid (EASE-Grid) - Northern Hemisphere" ;
                :Latitude\ Range = "35N - 85N" ;
                :Longitude\ Range = "180W - 180E" ;
                :Spatial\ Resolution = "25 X 25 sq.km" ;
                :Processing\ Software\ Name = "FMI assimilation algorithm (Pulliainen 2006)" ;
                :Processing\ Software\ Version = "v 2.0" ;
                :Processing\ Organisation = "Finnish Meteorological Institute" ;
                :Weather\ Station\ Data\ Date = "" ;
                :Land\ Mask = "GLC-2000 derived land classification mask" ;
                :Landmask\ Version = "v 2.0" ;
                :Mountain\ Mask = "ETOPO-5 derived mountain mask" ;
                :Mountain\ Mask\ Version = "v 2.0" ;
                :Forest\ Mask = "GLC-2000 derived forest mask" ;
                :project_id = "ESA GlobSnow" ;
                :contact = "Kari Luojus <Kari.Luojus at fmi.fi>" ;


More information about the gradsusr mailing list