mod_amr_grid.t

Go to the documentation of this file.

module mod_amr_grid
 
  implicit none
  private
 
  public :: settree
  public :: resettree
  public :: resettree_convert
 
contains
 
  !> Build up AMR
  subroutine settree
    use mod_global_parameters
    use mod_fix_conserve
    use mod_coarsen_refine
    use mod_errest, only: errest 
 
    ! create and initialize grids on all levels > 1. On entry, all
    ! level=1 grids have been formed and initialized. This subroutine
    ! creates and initializes new level grids
 
    integer :: igrid, iigrid, levnew
 
    ! when only one level allowed, there is nothing to do anymore
    if (refine_max_level == 1) return
 
    do levnew=2,refine_max_level
 
       call errest
 
       call amr_coarsen_refine
 
       if (.not.reset_grid) then
         ! if no finer level grids created: exit
         if (levmax/=levnew) exit
       end if
    end do
 
    ! set up boundary flux conservation arrays
    if(levmax>levmin) call allocatebflux
 
  end subroutine settree
 
  !> reset AMR and (de)allocate boundary flux storage at level changes
  subroutine resettree
    use mod_global_parameters
    use mod_fix_conserve
    use mod_amr_fct
    use mod_coarsen_refine
    use mod_errest, only: errest 
    use mod_particles
    use mod_eos, only: prepare_eos_w_fields
 
    if(levmax>levmin) call deallocatebflux
    if(stagger_grid) call deallocatebfaces
 
    ! deallocate old grids for particles
    if(use_particles) then
      call finish_gridvars()
    end if
  
    call errest
 
    call amr_coarsen_refine
 
    ! set up boundary flux conservation arrays
    if(levmax>levmin) call allocatebflux
 
    ! generate new grids for particles
    if(use_particles) then
      call init_gridvars()
    end if
 
    ! after getbc in amr_coarsen_refine, eos variables of ghost cells are zero for new blocks
    ! update the additional w state eos variables including ghost cells
    call prepare_eos_w_fields()
 
  end subroutine resettree
 
  !> Force the tree to desired level(s) from level_io(_min/_max)
  !> used for conversion to vtk output
  subroutine resettree_convert
    use mod_global_parameters
    use mod_ghostcells_update
    use mod_coarsen_refine
    use mod_errest, only: forcedrefine_grid_io 
 
    integer  :: igrid,iigrid, my_levmin, my_levmax
 
    if(level_io > 0) then
      my_levmin = level_io
      my_levmax = level_io
    else
      my_levmin = max(1,level_io_min)
      my_levmax = min(refine_max_level,level_io_max)
    end if
 
    do while(levmin<my_levmin.or.levmax>my_levmax)
     call getbc(global_time,0.d0,ps,iwstart,nwgc)
     do iigrid=1,igridstail; igrid=igrids(iigrid);
        call forcedrefine_grid_io(igrid,ps(igrid)%w)
     end do
  
     call amr_coarsen_refine
    end do
 
  end subroutine resettree_convert
  
end module mod_amr_grid