PROGRAM NOD10_MAIN
C
C     SET MAXIMAL DIMENSIONS
C       MNM = MAX NUMBER OF GRIDPOINTS   (EG. 100000)
C       MNN = MAX NUMBER OF MEASUREMENTS (EG. 200)
C       MNP = MAX NUMBER OF COVARIATES   (WITHOUT ONES, EG. 5)
C       -----------
C       MNMP = MNM * MNP                   = 500000
C       MNNP = MNN * MNP                   = 500
C       MNRQ = ...                         = 6
C       MNN1 = MNN + 1                    = 101
C       MNN2 = MNN * (MNN + 1) / 2         = 5050
C       MNN3 = (MNN + 1) * (MNN + 2) / 2   = 5151
C       MNP1 = MNP + 1                    = 6
C       MNP3 = (MNP + 1) * (MNP + 2) / 2   = 21
C       MNPN = (MNP + 1) * MNN             = 600
C       MNPW = (MNP + 1) * (MNN + 1)       = 606
C
      INTEGER MNM, MNN, MNR, MNP, MNMP, MNNP, MNRQ, MNN1, MNN2, MNN3, 
     & MNP1, MNP3, MNPN, MNPW
C
C     BASIC PARAMETERS
      PARAMETER (MNM = 100000)
      PARAMETER (MNN = 200)
      PARAMETER (MNP = 5)
C
C     DERIVED PARAMETERS
      PARAMETER (MNMP = MNM * MNP)
      PARAMETER (MNNP = MNN * MNP)
      PARAMETER (MNRQ = 6)
      PARAMETER (MNN1 = MNN + 1)
      PARAMETER (MNN2 = MNN * MNN1 / 2)
      PARAMETER (MNN3 = MNN1 * (MNN + 2) / 2)
      PARAMETER (MNP1 = MNP + 1)
      PARAMETER (MNP3 = MNP1 * (MNP + 2) / 2)
      PARAMETER (MNPN = MNP1 * MNN)
      PARAMETER (MNPW = MNP1 * MNN1)
C
      INTEGER I1T(MNM)
      REAL*8 
     & VCX(MNM), VCY(MNM), VBX(MNN), VBY(MNN), MXC(MNMP), MXB(MNNP), 
     & MSGMI(MNN2), MXTSI(MNPN), MXSIX(MNP3), MSGMIF(MNN2),MXSIXF(MNP3),
     & MSGMIP(MNN3), MXTSIP(MNPW), MXSIXP(MNP3),
     & VT(MNRQ), VAA(MNN1), VBB(MNP1), VZ(MNN), VY(MNP1), 
     & VX0(MNP1), VC0(MNN1), VNP(MNP1), VNN(MNN1), 
     & F0T(MNM), FKT(MNM), F1T(MNM)
      LOGICAL VCM(MNM), VBM(MNN)
C
      CALL NOD10_SUB_MAIN(MNM,MNN,MNR,MNP,
     & MNMP,MNNP,MNRQ,MNN1,MNN2,MNN3,MNP1,MNP3,MNPN,MNPW,
     & I1T,VCX,VCY,VBX,VBY,MXC,MXB,
     & MSGMI,MXTSI,MXSIX,MSGMIF,MXSIXF,MSGMIP,MXTSIP,MXSIXP,
     & VT,VAA,VBB,VZ,VY,VX0,VC0,VNP,VNN,F0T,FKT,F1T,VCM,VBM)
C
      END