function [data,lat,lon] = fdread(fdname)

fid = fopen(fdname, 'r', 'ieee-be');
if(fid < 0)
    error(['Could not open file ' fdname]);
end

Inull = -1000;
GRID = 6596;
data = fread(fid, GRID, 'float32');

data(find(data == Inull)) = NaN;

lon = [-72:71]*360/144+360/288;
lon = [-72:71]*360/144+360/288;