function [zepo,Kg]=zpsd(th,ku,ky)
%ZPSD	See TH2ZP

% *** Set up default values ***
if nargin<3,ky=[];end
if nargin<2,ku=[];end
[zepo,Kg] = th2zp(th,ku,ky);
