errest Subroutine

public subroutine errest()

Uses

  • proc~~errest~~UsesGraph proc~errest errest module~mod_forest mod_forest proc~errest->module~mod_forest module~mod_global_parameters mod_global_parameters proc~errest->module~mod_global_parameters module~mod_basic_types mod_basic_types module~mod_global_parameters->module~mod_basic_types module~mod_connectivity mod_connectivity module~mod_global_parameters->module~mod_connectivity module~mod_constants mod_constants module~mod_global_parameters->module~mod_constants module~mod_physicaldata mod_physicaldata module~mod_global_parameters->module~mod_physicaldata module~mod_variables mod_variables module~mod_global_parameters->module~mod_variables mpi mpi module~mod_global_parameters->mpi module~mod_variables->module~mod_basic_types

Do all local error estimation which determines (de)refinement

Arguments

None

Calls

proc~~errest~~CallsGraph proc~errest errest proc~forcedrefine_grid forcedrefine_grid proc~errest->proc~forcedrefine_grid proc~lohner_grid lohner_grid proc~errest->proc~lohner_grid proc~mpistop mpistop proc~errest->proc~mpistop usr_refine_grid usr_refine_grid proc~forcedrefine_grid->usr_refine_grid mpi_abort mpi_abort proc~mpistop->mpi_abort

Called by

proc~~errest~~CalledByGraph proc~errest errest proc~resettree resettree proc~resettree->proc~errest proc~settree settree proc~settree->proc~errest proc~main main proc~main->proc~settree proc~timeintegration timeintegration proc~main->proc~timeintegration proc~timeintegration->proc~resettree program~agile agile program~agile->proc~main

Source Code

  subroutine errest
    use mod_forest, only: refine, coarsen
    use mod_global_parameters

    if (igridstail==0) return

    select case (refine_criterion)
    case (1)
       ! all refinement solely based on user routine usr_refine_grid
    case (3)
       ! Error estimation is based on Lohner's scheme
       call lohner_grid()
    case default
       call mpistop("Unknown error estimator")
    end select

    if ( refine_usr ) then
       call forcedrefine_grid()
    end if

    

  end subroutine errest