[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