high_res_cube/matlab-grid-generator/map_lonlat2xyz.m

function [X,Y,Z]=map_lonlat2xyz(lon,lat)
% [X,Y,Z]=map_lonlat2xyz(lon,lat);
%
% Convert spherical coordinates [lon,lat] to 3-D coordinates [x,y,z]
%
% Assumes "lat" is positive with "z", equatorial plane
% falls at z=0  and "lon" is measured anti-clockwise (eastward)
% from x-axis (y=0) about z-axis.

% Latitude
Req=cos(lat);
Z=sin(lat);
X=Req.*cos(lon);
Y=Req.*sin(lon);
1	cnh	1.1	function [X,Y,Z]=map_lonlat2xyz(lon,lat)
2			% [X,Y,Z]=map_lonlat2xyz(lon,lat);
3			%
4			% Convert spherical coordinates [lon,lat] to 3-D coordinates [x,y,z]
5			%
6			% Assumes "lat" is positive with "z", equatorial plane
7			% falls at z=0 and "lon" is measured anti-clockwise (eastward)
8			% from x-axis (y=0) about z-axis.
9
10			% Latitude
11			Req=cos(lat);
12			Z=sin(lat);
13			X=Req.*cos(lon);
14			Y=Req.*sin(lon);