func dsgrid(m1,m2,n,nwin,ngwin)
 ;11/29/88 r. shine
 ;1/10/94  mods for nxg != nyg
 ;just returns the gridmatch result
 
 nx=dimen(m1,0)
 ny=dimen(m1,1)
 if !narg lt 5 then ngw=fix(2.*(nx>ny)/n)  else ngw=ngwin
 if !narg lt 4 then nw=fix(1.25*ngw)       else nw=nwin
 nxg=n
 ;nyg=n
 nyg=rfix(float(n)*ny/nx)
 ;ty,'nxg, nyg =',nxg,nyg
 ;note the wx calculation should be FP, 3/18/89
 wx=float(nx)/nxg
 wy=float(ny)/nyg
 ;this is supposed to compute the grid the same way as Stu's programs
 gx=indgen(fltarr(nxg,nyg),0)*wx+wx/2.-1
 gy=indgen(fltarr(nxg,nyg),1)*wy+wy/2.-1
 dx=nw
 dy=nw
 gwid=ngw
 ;ty,'gx,gy,dx,dy,gwid =', gx,gy,dx,dy,gwid
 return,gridmatch(m1,m2,gx,gy,dx,dy,gwid)
 endfunc
 ;============================================================================