function [map, lat, lon ] = fdmap(data, LonBeg, LonEnd)

Lnum = 144;
Znum = 72;

lat = [-Znum/2:Znum/2-1]*360/Znum+360/Lnum/2;
lon = [-Lnum/2:Lnum/2-1]*360/Lnum+360/Lnum/2;

map = zeros(Znum,Lnum);
Icells = [3,9,16,22,28,34,40,46,52,58,64,69,75,80,85,90,95,100,104,108, ...
            112,116,120,123,126,129,132,134,136,138,140,141,142,143,144, ...
            144,144,144,143,142,141,140,138,136,134,132,129,126,123,120, ...
            116,112,108,104,100,95,90,85,80,75,69,64,58,52,46,40,34,28,22,16, ...
            9,3];
        

ibox = 1;
for( lat = 1:Znum )
    for( lon = 1:Icells(lat) )
        lb = LonBeg(ibox)+1;
        le = LonEnd(ibox)+1;
        for( ilon = lb:le )
            % shifts map from LON=0-360 to LON= -180,+180
            lonsq = ilon;
            lonsq = lonsq + Lnum/2;
            if( lonsq > Lnum )
                lonsq = lonsq - Lnum;
            end
            map(lat,lonsq) = data(ibox);
        end
        ibox = ibox + 1;
    end
end