anchor.diff:      call step
       time = time + dt
   
  +    if(mod(itermax,10) .eq. 0) then
  +      print*,"Iteration: ",iter," Time: ",time
  +      print*,"Max U: ", maxval(vx(:))," Max V: ",maxval(vy(:))," Max W: ",maxval(vz(:))
  +      print*,"  "
  +    end if
  +
     end do
   
     call deinit
  
positive.code: ! ===================================== !
  program sph_run
    use state
    implicit none
  
    time = 0.0_WP
  
negative.code: ! ================================================ !
  ! Read in command line arguments and initialize    !
  ! ================================================ !
  subroutine init
    use state
    use constants
  