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: ! ================================================ !
  ! Read in command line arguments and initialize    !
  ! ================================================ !
  
  ! ===================================== !
  ! Allocate all arrays in the modules    !
  ! ===================================== !
  subroutine init
    use state
    use constants
  
negative.code: 
    call init
  
    iter = 0
  
    do while (time .lt. tfin)
      iter = iter + 1
      if(iter .gt. itermax) then
        print*, "Exceeded maximum iterations. Exiting..."
        stop
  