setpar Subroutine

   impure subroutine setpar()
   !! This routine sets the basic problem parameters which are passed to the various routines
   !! via the module [[web_par]].
      integer :: i, j

      ax = one
      ay = one

      aa = one
      ee = 1e4_rk
      gg = 0.5e-6_rk
      bb = one
      dprey = one
      dpred = 0.05_rk
      alpha = 5e1_rk
      beta = 3e2_rk

      dx = ax/(mx - 1)
      dy = ay/(my - 1)

      do j = 1, np
         do i = 1, np
            acoef(np + i, j) = ee
            acoef(i, np + j) = -gg
         end do
         acoef(j, j) = -aa
         acoef(np + j, np + j) = -aa
         bcoef(j) = bb
         bcoef(np + j) = -bb
         diff(j) = dprey
         diff(np + j) = dpred
      end do

      do i = 1, ns
         cox(i) = diff(i)/dx**2
         coy(i) = diff(i)/dy**2
      end do

   end subroutine setpar