C dim.p2
      INTEGER MAXSTP              !array size for density
      PARAMETER (MAXSTP=1000)