function [np, gam, hm, tm, hv, tv] = vario3di (... 
          nd, x, y, z, vr, nlag, xlag, xltol, ...
          ndir, azm, atol, bandwh, dip, dtol, bandwd, nvarg, ...
          ivtail, ivhead, ivtype, nvar)
%
%      vario3di : matlab m-file
%      gamv3    : fortran mex-file (corresponds to gamv3 in GSLIB)
%      gamv3g   : fortran gateway function
%
%
%              Variogram of 3-D Irregularly Spaced Data
%              ****************************************
%
% This subroutine computes a variety of spatial continuity measures of a
% set for irregularly spaced data.  The user can specify any combination
% of direct and cross variograms using any of eight "variogram" measures
%
% INPUT VARIABLES:
%
%   nd               Number of data (no missing values)
%   x(nd)            X coordinates of the data
%   y(nd)            Y coordinates of the data
%   z(nd)            Z coordinates of the data
%   nv               The number of variables
%   vr(nd,nv)        Data values
%   tmin,tmax        Trimming limits
%   nlag             Number of lags to calculate
%   xlag             Length of the unit lag
%   xltol            Distance tolerance (if <0 then set to xlag/2)
%   ndir             Number of directions to consider
%   azm(ndir)        Azimuth angle of direction (measured positive
%                      degrees clockwise from NS).
%   atol(ndir)       Azimuth (half window) tolerances
%   bandwh           Maximum Horizontal bandwidth (i.e., the deviation 
%                      perpendicular to the defined azimuth).
%   dip(ndir)        Dip angle of direction (measured in negative
%                      degrees down from horizontal).
%   dtol(ndir)       Dip (half window) tolerances
%   bandwd           Maximum "Vertical" bandwidth (i.e., the deviation
%                      perpendicular to the defined dip).
%   nvarg            Number of variograms to compute
%   ivtail(nvarg)    Variable for the tail of each variogram
%   ivhead(nvarg)    Variable for the head of each variogram
%   ivtype(nvarg)    Type of variogram to compute:
%                      1. semivariogram
%                      2. cross-semivariogram
%                      3. covariance
%                      4. correlogram
%                      5. general relative semivariogram
%                      6. pairwise relative semivariogram
%                      7. semivariogram of logarithms
%                      8. rodogram
%                      9. madogram
%                     10. indicator semivariogram: an indicator variable
%                         is constructed in the main program.
%   nvar            Number of variables
%
% OUTPUT VARIABLES:
%
%   np()             Number of pairs
%   dis()            Distance of pairs falling into this lag
%   gam()            Semivariogram, covariance, correlogram,... value
%   hm()             Mean of the tail data
%   tm()             Mean of the head data
%   hv()             Variance of the tail data
%   tv()             Variance of the head data
%

%   INPUT VARIABLES

tmin = -1e21;
tmax = +1e21;

%   VARIOGRAMME

[np, dis, gam, hm, tm, hv, tv] = gamv3 (nd,x,y,z,vr,tmin,tmax,nlag, ...
	xlag,xltol,ndir,azm, atol,bandwh,dip,dtol,bandwd,nvarg,ivtail,ivhead,ivtype,nvar );

%   OUTPUT VARIABLES

[np, gam, hm, tm, hv, tv] = outvario(nlag+2,dis,ndir,nvarg,np,gam,hm,tm,hv,tv,ivtype);