dump_results Subroutine

  subroutine dump_results(this)
    type(pf_results_t), intent(inout) :: this
    integer :: i, j, k
    
    open(unit=this%p_index, file=this%fname_r, form='formatted')
    do k = 1, this%nlevels
       do j = 1, this%nblocks
          do i = 1 , this%niters
             write(this%p_index, '(I4, I4, I4, e21.14)') i,j,k,this%residuals(i, j, k)
          end do
       end do
    enddo
    close(this%p_index)

  end subroutine dump_results