      parameter(mnx=801 ,mny=801,mlsour=600,mnmat=20)
      parameter(mngx=801, mngy=801)
      parameter(ntstep=10000,mnrec=150,mtsnap=10)
