%      test the conversion to x,y from lat lon
global SECS COSROT SINROT ONCE
global UTMXOFF UTMYOFF XOFF YOFF ROTANGLE IZONE ISPHERE
global AXIS BXIS AK0 RADSEC

lat = [-66  -69  -67  -68  -67 ];
lon = [-67  -77  -65  -71  -69 ];

%     correct grid  values 
cx = [571 044  525 288 428];
cy = [112 240 -022 125 122];

[x,y]=ll2grid(lat,lon);


for i=1:5

   disp(['location: (' num2str(lat(i)) ', ' num2str(lon(i)) ') is grid ' num2str(x(i)) ...
    '.' num2str(y(i))  ]);
   disp(['      correct is ' num2str(cx(i)) '.' num2str(cy(i))]);
   disp('    ');

end



 

