anchor.diff:  			call calc_OBJPOT_obj
   		end if
   		if (enableTrap) then
  -			trVelx = (TVXDASH/2.0d0)*(tanh((-6.0d0*(DBLE(TIME)-(3.0d0*TTM/4.0d0))/(TTM/4.0d0))+3.0d0)+1.0d0)
  -			trVely = (TVYDASH/2.0d0)*(tanh((-6.0d0*(DBLE(TIME)-(3.0d0*TTM/4.0d0))/(TTM/4.0d0))+3.0d0)+1.0d0)
  +			trVelx = (TVXDASH/2.0d0)*(tanh((-6.0d0*(TIME-(3.0d0*TTM/4.0d0))/(TTM/4.0d0))+3.0d0)+1.0d0)
  +			trVely = (TVYDASH/2.0d0)*(tanh((-6.0d0*(TIME-(3.0d0*TTM/4.0d0))/(TTM/4.0d0))+3.0d0)+1.0d0)
   			TXDASH  = TXDASH  + (trVelx*DBLE(DT))
   			TYDASH  = TYDASH  + (trVely*DBLE(DT))  
   			do i = -NX/2,NX/2
  
positive.code: 		end if
  	end if
  	!Section 2 - trap moving at terminal velocity
  	if(TIME .gt. (TTM/4.0d0) .and. TIME .lt. (TTM/2.0d0)) then
  		if (enablePot) then
  			call calc_OBJPOT_obj
  		end if
  
negative.code: 			call calc_OBJPOT_obj
  		end if
  		if (enableTrap) then
  			trVelx = (TVXDASH/2.0d0)*(tanh((6.0d0*TIME/(TTM*4.0d0))-3.0d0)+1.0d0)
  			trVely = (TVYDASH/2.0d0)*(tanh((6.0d0*TIME/(TTM*4.0d0))-3.0d0)+1.0d0)
  			TXDASH  = TXDASH  + (trVelx*DBLE(DT))
  			TYDASH  = TYDASH  + (trVely*DBLE(DT))  
  			do i = -NX/2,NX/2
  