function [r,q] = map1d(xcgrid,ycgrid,qgrid)
% for 1d scatter plot, plot q vs. "rotated latitude"
% computes latitude relative to axis (x1,y1,z1) set in setprob.m

global x1 y1 z1

[xgrid,ygrid,zgrid] = mapc2m(xcgrid,ycgrid);
r = asin(xgrid(:)*x1 + ygrid(:)*y1 + zgrid(:)*z1);
q = qgrid(:);

