C     IF YOU CHANGE ANY OF THESE PARAMETERS, YOU WILL ALSO NEED
c     TO CHANGE THE DIMENSIONS OF THE ARRAYS IN THE MATR SUBROUTINE
      PARAMETER (NPN1=100, NPNG1=300, NPNG2=2*NPNG1, NPN2=2*NPN1,  
     &           NPL=NPN2+1, NPN3=NPN1+1,  
     &           NPN4=80, NPN5=2*NPN4, NPN6=NPN4+1, NPL1=NPN5+1)