C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
C                                                                      %
C Copyright (C) 1996, The Board of Trustees of the Leland Stanford     %
C Junior University.  All rights reserved.                             %
C                                                                      %
C The programs in GSLIB are distributed in the hope that they will be  %
C useful, but WITHOUT ANY WARRANTY.  No author or distributor accepts  %
C responsibility to anyone for the consequences of using them or for   %
C whether they serve any particular purpose or work at all, unless he  %
C says so in writing.  Everyone is granted permission to copy, modify  %
C and redistribute the programs in GSLIB, but only under the condition %
C that this notice and the above copyright notice remain intact.       %
C                                                                      %
C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
c-----------------------------------------------------------------------
c
c                Variogram of Irregular Spaced 3-D Data
c                **************************************
c
c The following Parameters control static dimensioning within gamv3:
c
c   MAXDIR    maximum number of directions possible at one time
c   MAXLAG    maximum number of lags at one time
c   MXVARG    maximum number of variograms possible at one time
c
c   MXDLV     MAXDIR*MAXLAG*MXVARG - used for dimensioning
c   EPSLON    a small number to avoid dividing by zero
c
c-----------------------------------------------------------------------
c
c User Adjustable Parameters:
c
      parameter(MAXDIR =  5, MAXLAG = 200, MXVARG =  5)
c
c Fixed Parameters:
c
      parameter(MAXLG=MAXLAG+2,MXDLV=MAXDIR*MAXLG*MXVARG,
     +          EPSLON=1.0e-20,VERSION=2.901)
c
c Variable Declaration:
c
      real      azm(MAXDIR),atol(MAXDIR),bandwh(MAXDIR),dip(MAXDIR),
     +          dtol(MAXDIR),bandwd(MAXDIR),xlag,xltol,tmin,tmax
      real*8    dis(MXDLV),gam(MXDLV),hm(MXDLV),tm(MXDLV),hv(MXDLV),
     +          tv(MXDLV),np(MXDLV)
      integer   nd,nlag,ndir,nvarg,isill,test,
     +          ivtail(MXVARG),ivhead(MXVARG),ivtype(MXVARG)
      character outfl*40
c
c Common Blocks:
c
      common /datagv/ nd,nlag,xlag,xltol,ndir,nvarg,isill,
     +                azm,atol,test,bandwh,dip,dtol,bandwd,
     +                ivtail,ivhead,ivtype,tmin,tmax

      common /parmdp/ dis,gam,hm,tm,hv,tv,np

      common /chargv/ outfl
