anchor.diff:      end if
   
     end do
  +  !$OMP END DO
   
     return
   end subroutine enforce_BCs
  
positive.code:   twrite = frame_freq
    call cpu_time(init_stop)
  
    iter = 0
    time = 0.0_WP
    do while (time .lt. tfin)
   !   if(my_id.eq.0) 
       iter = iter + 1
  
      ! Simulate the fluid flow
      call step
      !if(my_id.eq.0)
      time = time + dt
  
      !if(my_id .eq. 0) then
        if(time .ge. twrite) then
          call cpu_time(write_start)
          twrite = twrite + frame_freq
  
negative.code:         call cpu_time(write_stop)
          write_sum = write_sum + (write_stop-write_start)
        end if
      !$OMP END MASTER
      !$OMP BARRIER
    end do
    !$OMP END PARALLEL
  
    call cpu_time(deinit_start)
    call deinit
  