function [lon2,lat2]=caploc(lonlat,rad,N,method)
%  [lon2,lat2]=CAPLOC([lon1 lat1],rad,N,method)
%
% Calculates the location of a circle on the sphere
%
% INPUT:
%
% [lon1 lat1]      Longitude, latitude of center(s), in degrees
% rad              Radius, in degrees
% N                Number of points [default: 100]
% method           1 Cartesian [default]
%                  2 Mollweide projection
%
% OUPUT:
%
% [lon2,lat2]      The requested points, in degrees
%
% See also: LATITUDE, LONGITUDE
%
% Last modified by fjsimons-at-alum.mit.edu, Feb 20th, 2004

defval('N',100)
angl=linspace(0,360,N);
defval('method',1)

delta=rad*fralmanac('DegDis','Earth')/1000;

[lon2,lat2]=grcazim(lonlat,delta,angl);

if method==2
   [lon2,lat2]=mollweide(lon2*pi/180,lat2*pi/180);
end