anchor.diff:         write(17, "(A, I0)") 'nProcX = ', nProcX
          write(17, "(A, I0)") 'nProcY = ', nProcY
       end if
  +    call clean_stop(0, .FALSE.)
       stop 1
     end if
   
  
positive.code:   end if
    return
  end subroutine write_output
  
negative.code:       hypre_grid)
  
    ! Finalize MPI
    call clean_stop(nTimeSteps, .TRUE.)
  
    stop
  end program aronnax
  