anchor.diff:          
   
       write(*,*) "Starting Solver"
  -    call system_clock(startTime, clock_rate, clock_max)
  +    call system_clock(COUNT_RATE=clock_rate, COUNT_MAX=clock_max)
  +    call system_clock(startTime)
       call solveOrder2(tEnd,nOuts,tDel,n,row,col,M,C,yLen,xDel,yDel,&
          c0,m0,y,nu,gam,alpha,beta,k,delta,ndiag,e1,e2,nit,esoln,dSelect,&
          fSelect,gSelect,avgIters,maxIters,avgNit,maxNit)
  -    call system_clock(endTime, clock_rate, clock_max)
  -    
  +    call system_clock(endTime)
  +   
  +    ! Convert times into seconds
  +    elapsedTime = real(endTime - startTime)/real(clock_rate)
  +
  +
   !    call solveOrder1(tEnd,nOuts,tDel,n,row,col,M,C,xLen,yLen,xDel,yDel&
   !       ,c0,m0,nu,gam,alpha,beta,k,delta,ndiag,e1,e2,nit,dSelect,fSelect)
     
  
positive.code:     write(*,*) "-------------------------------------------------------------"
      write(*,*) "Statsitcs:"
      write(*,*) "-------------------------------------------------------------"
      write(*,*) "Time to compute = ", endTime-startTime
      write(*,*) ""
      write(*,*) "Avg Iters for iterating betn. soln. =", avgIters
      write(*,*) "Max Iters for iterating betn. soln. =", maxIters
      write(*,*) "Avg Iters for linear solver =", avgNit
      write(*,*) "Max Iters for linear solver =", maxNit  
      write(*,*) "Writing statistics to file"
      call reportStats(avgIters,maxIters,avgNit,maxNit,endTime-startTime)
      
      write(*,*) "Execution completed"
    
  
negative.code:     real :: avgIters, maxIters
      real :: avgNit, maxNit
      integer :: startTime, endTime, clock_rate, clock_max
      real :: elapsedTime
      
      write(*,*) "Enter parameter file name: ex. 'parameter.txt' "
      write(*,'(A)', advance="no") "    "
  