From:	"Jason Roberts" <>
To:	"'Regis GALLON'" <>
Cc:	<>
Subject:	RE: [mget-help] help with Netcdf file converstion
Date:	Tue, 22 May 2012 11:18:12 -0400

Dear Régis,

I am sorry about this problem. It is most likely a bug in MGET relating to
the fact that the surface_solar_irradiance variable has three dimensions.
The MGET tool is designed to convert variables that have only two
dimensions. But because the time dimension of surface_solar_irradiance has a
length of 1, the tool is supposed to accept it as a special case and perform
the conversion of the single time slice. I suspect the bug relates to
handling that special case.

Could you please email me an example netCDF file so I can try to convert it
myself? I do not have access to that data. I will debug the problem and send
you a patch.

Best regards,

Jason

-----Original Message-----
From: Regis GALLON 
[mailto:]
 
Sent: Tuesday, May 22, 2012 10:10 AM
To: 

Subject: [mget-help] help with Netcdf file converstion

Dear,

I try the MGET tool and I have got some difficulties to use "the Netcdf to
ArcGIS raster tool"
When I run the tool I get the error "not yet implemented for numpy.ma
<http://numpy.ma/> arrays"


My netcdf header is :
dimensions:
    time = 1 ;
    depth = 1 ;
    latitude = 1500 ;
    longitude = 1451 ;
variables:
    int time(time) ;
        time:long_name = "time" ;
        time:units = "hours since 1900-1-1 0:0:0" ;
    float depth(depth) ;
        depth:long_name = "depth" ;
        depth:units = "meters" ;
        depth:positive = "up" ;
    float latitude(latitude) ;
        latitude:long_name = "latitude" ;
        latitude:units = "degrees_N" ;
    float longitude(longitude) ;
        longitude:long_name = "longitude" ;
        longitude:units = "degrees_E" ;
    short surface_solar_irradiance(time, latitude, longitude) ;
        surface_solar_irradiance:long_name = "surface solar irradiance" ;
        surface_solar_irradiance:units = "watt/meter2" ;
        surface_solar_irradiance:scale_factor = 0.1 ;
        surface_solar_irradiance:add_offset = 0.f ;
        surface_solar_irradiance:valid_min = 0s ;
        surface_solar_irradiance:valid_max = 400s ;
        surface_solar_irradiance:missing_value = 32767s ;
        surface_solar_irradiance:_FillValue = 32767s ;
    short quality_index(time, latitude, longitude) ;
        quality_index:long_name = "quality index" ;
        quality_index:units = "n/a" ;
        quality_index:scale_factor = 1.f ;
        quality_index:add_offset = 0.f ;
        quality_index:valid_min = -32767s ;
        quality_index:valid_max = 32767s ;
        quality_index:missing_value = 32767s ;
        quality_index:_FillValue = 32767s ;

// global attributes:
        :CONVENTIONS = "COARDS" ;
        :long_name = "Merged Atlantic Product / 24-hourly surface solar
irradiance" ;
        :short_name = "MAP-SSI-24" ;
        :producer_agency = "O&SI SAF" ;
        :netcdf_version_id = "3.4" ;
        :product_version = "1.0" ;
        :time_resolution = "24-hourly mean" ;
        :spatial_resolution = "0.1 degree" ;
        :platform_id = "GOES/MSG/NOAA" ;
        :south_latitude = -60. ;
        :north_latitude = 89.9 ;
        :west_longitude = -100. ;
        :east_longitude = 45. ;
}
 

To convert to Raster I use the following parameters :

variable name : surface_solar_irradiance
x coordinate -100
y coordinate -60
cell size 0.1
NODATA value 32767

Is it correct?

Anything you can offer up would be warmly received.

Sincerely

Régis GALLON

< date >
< thread >

Archives powered by MHonArc.