anchor.diff:  
   ! output value of PI and time taken
   
  -  write(*,'(A,1F12.10,A)'),"Obtained value of PI: ", mypi
  -  write(*,'(A,1F12.5,A)'),"Time taken:           ",(stop-start), " seconds"
  +  write(*,'(A,1F12.10,A)') "Obtained value of PI: ", mypi
  +  write(*,'(A,1F12.5,A)') "Time taken:           ",(stop-start), " seconds"
   
   end program pi_openmp
   
  
positive.code:   implicit none
  
    double precision  :: step, x, sum, mypi, start, stop, omp_get_wtime
    integer(kind=8)   :: num_steps, i
    character(len=32) :: arg
    integer, external :: omp_get_num_threads 
    integer, external :: omp_get_max_threads
  
negative.code: 
  ! Output start message
  
    write(*,'(A)') "Calculating PI using:"
    write(*,'(A,1I16,A)') "                  ",num_steps, " slices"
    write(*,'(A,1I16,A)') "                  ",omp_get_max_threads()," OpenMP threads"
  
  ! Initialise time counter and sum: set step size
  
  