mg_copy_from_tree Subroutine

public subroutine mg_copy_from_tree(iw_from, iw_to, state_to)

Uses

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

Copy a variable from the multigrid tree

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: iw_from

Variable to use as right-hand side
integer, intent(in) :: iw_to

Copy to this variable
type(state), intent(inout), optional, target :: state_to(max_blocks)

If present, use this as the output state

Source Code

  subroutine mg_copy_from_tree(iw_from, iw_to, state_to)
    use mod_global_parameters
    use mod_forest
    integer, intent(in)      :: iw_from
    !< Variable to use as right-hand side
    integer, intent(in)      :: iw_to
    !< Copy to this variable
    !> If present, use this as the output state
    type(state), intent(inout), optional, target :: state_to(max_blocks)
    integer                  :: iigrid, igrid, id
    integer                  :: nc, lvl
    type(tree_node), pointer :: pnode

    if (.not. mg%is_allocated) error stop "mg_copy_from_tree: tree not allocated yet"

    do iigrid = 1, igridstail
       igrid =  igrids(iigrid);
       pnode => igrid_to_node(igrid, mype)%node
       id    =  pnode%id
       lvl   =  mg%boxes(id)%lvl
       nc    =  mg%box_size_lvl(lvl)

       if (present(state_to)) then
          state_to(igrid)%w(ixMlo1:ixMhi1,ixMlo2:ixMhi2,ixMlo3:ixMhi3,&
              iw_to) = mg%boxes(id)%cc(1:nc,1:nc,1:nc, iw_from)
       else
          ps(igrid)%w(ixMlo1:ixMhi1,ixMlo2:ixMhi2,ixMlo3:ixMhi3,&
              iw_to) = mg%boxes(id)%cc(1:nc,1:nc,1:nc, iw_from)
       end if
    end do
  end subroutine mg_copy_from_tree