function msaselmap

pos=[0,140,35786000]; % MTSAT-1R PRN129

lat=20:2:50; lon=120:2:160;
for i=1:length(lon)
    for j=1:length(lat)
        azel=satazel(geodtoecef(pos),geodtoecef([lat(j),lon(i),0]));
        el(j,i)=azel(2)*180/pi;
    end
end
figure('color','w')
gmt('mmap','proj','miller','cent',[136,37],'scale',12,'pos',[0,0,1,0.95]);
gmt('mcoast','ccolor',[.5,.5,.5]);
[x,y]=gmt('lltoxy',lon,lat);
[c,h]=contour(x,y,el,0:90,'r'); clabel(c,h);
title('MTSAT Elavation Angle (deg) (MTSAT-1R PRN129 140E)');
