set_geometry_index_ranges Subroutine

private subroutine set_geometry_index_ranges()

Uses

  • proc~~set_geometry_index_ranges~~UsesGraph proc~set_geometry_index_ranges set_geometry_index_ranges module~mod_global_parameters mod_global_parameters proc~set_geometry_index_ranges->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

Set the index range over which the cell metrics are known. With an even number of ghost cells this is just the ghosted block; with an odd number (or a staggered grid) ghost-cell prolongation needs one extra layer.

Arguments

None

Calls

proc~~set_geometry_index_ranges~~CallsGraph proc~set_geometry_index_ranges set_geometry_index_ranges proc~mpistop mpistop proc~set_geometry_index_ranges->proc~mpistop mpi_abort mpi_abort proc~mpistop->mpi_abort

Called by

proc~~set_geometry_index_ranges~~CalledByGraph proc~set_geometry_index_ranges set_geometry_index_ranges proc~initialize_vars initialize_vars proc~initialize_vars->proc~set_geometry_index_ranges proc~initialize_amrvac initialize_amrvac proc~initialize_amrvac->proc~initialize_vars proc~main main proc~main->proc~initialize_amrvac program~agile agile program~agile->proc~main

Source Code

  subroutine set_geometry_index_ranges()
    use mod_global_parameters

    integer :: icase

    icase = mod(nghostcells,2)
    if (stagger_grid) icase = 1
    select case (icase)
    case (0)
       ixGextmin1=ixGlo1;ixGextmin2=ixGlo2;ixGextmin3=ixGlo3
       ixGextmax1=ixGhi1;ixGextmax2=ixGhi2;ixGextmax3=ixGhi3;
    case (1)
       ixGextmin1=ixGlo1-1;ixGextmin2=ixGlo2-1;ixGextmin3=ixGlo3-1
       ixGextmax1=ixGhi1+1;ixGextmax2=ixGhi2+1;ixGextmax3=ixGhi3+1;
    case default
       call mpistop("no such case")
    end select

  end subroutine set_geometry_index_ranges