fid = fopen('lf_dbz','r','b');
test = fread(fid,'float');
lf = reshape(test,240,240);

fid = fopen('lf_coords','r','b');
test = fread(fid,'float');
coords = reshape(test,240,2);

[lat,lon] = meshgrid(coords(:,1),coords(:,2));

%xl = (lon - slon)*111*cos(slat*d2r);   %km
%yl = (lat - slat)*111;                 %km

figure(5)
v = 20:3:52;

contourf(lon.',lat.',lf.',v)     %the ".'" rotates the array 180 degrees about the x-axis
colorbar