由站点经纬度找格点

matlab程序，size(A)返回矩阵维度信息。

nx=195;ny=153;sta=1551;nz=2;
output = zeros(sta,4);

for k=1:sta
    IOUT = 0;
    JOUT = 0;
    dis=1000000.0;
    for j=1:ny
        for i=1:nx
            tmp_dis=sqrt((XLAT_M(i,j)-lonlat(k,2))^2+(XLONG_M(i,j)-lonlat(k,1))^2);
            if (tmp_dis <= dis) 
                dis = tmp_dis;
                IOUT = i;
                JOUT = j;
            end
        end
    end
  output(k,1)=IOUT;
  output(k,2)=JOUT;
  output(k,3)=XLAT_M(IOUT,JOUT);
  output(k,4)=XLONG_M(IOUT,JOUT);
end