program testmkmap
   use mkmapmodule
   implicit none
   integer,parameter        :: m=100
   integer,parameter        :: n=30
   integer,parameter        :: m2 = 40
   integer,parameter        :: n2 = 50
   real*8, dimension(m,n)   :: x,y
   real*8, dimension(n2)    :: x2,y2,xs,ys
   integer, dimension(m,n)  :: code
   integer, dimension(n2)   :: nrx,nry,iflag,nrin,covered
   real*8, dimension(4,n2)  :: w
   integer, dimension(4,n2) :: iref
   integer                  :: iprint=1
   real*8                   :: xymiss

   call mkmap(code      ,x        ,y        ,m        ,n        , &
   &                     x2        ,y2        ,n2        ,xs        ,ys        , &
   &                     nrx       ,nry       ,iflag     ,nrin      ,w         , &
   &                     iref      ,iprint    ,covered   ,xymiss)

end program testmkmap