!jcw
      integer :: ixr, iyr, ix, iy, ncomp, numdata

      parameter(ixr=300,iyr=300,ix=21,iy=1000,ncomp=1,numdata=10000)