anchor.diff:  		if((vortexKillX .and. abs(xloc) > vortexKillXDist) .or.&
   		   (vortexKillY .and. abs(yloc) > vortexKillYDist)) then
   			write(6,*) "Killing positive vortex at: ", xloc, ", ", yloc
  -			call insert_vortex(xloc,yloc,-1.0d0)
  +			call insert_vortex(xloc,yloc,-1)
   			call set_homg_sq(nint(xloc/DSPACE),nint(yloc/DSPACE),8)
   		end if
   	end do
  
positive.code: 	end do
  end subroutine
  
  subroutine insert_rand_vortex(n,r,circ)
  	use params
  	implicit none
  	integer :: i,n
  	double precision ::circ,r
  	do i = 1,n
  		call insert_vortex(dble(r*RAND() - r/2.0),dble(r*RAND() - r/2.0),circ) 
  	end do
  end
  
  subroutine impose_rand_vortex(n,r,circ,seed)
  	use params
  	implicit none
  	integer :: i,j,k,n,seed
  	double precision ::circ,rv,tv,xloc,yloc,rs,xx,yy,r
  	call srand(seed+nint(circ)+1)
  
  	do k=1,n
  		rv = rand()*r
  
negative.code: 		if((vortexKillX .and. abs(xloc) > vortexKillXDist) .or.&
  		   (vortexKillY .and. abs(yloc) > vortexKillYDist)) then
  			write(6,*) "Killing negative vortex at: ", xloc, ", ", yloc
  			call insert_vortex(xloc,yloc,1)
  			call set_homg_sq(nint(xloc/DSPACE),nint(yloc/DSPACE),8)
  		end if
  	end do
  