Generate and initialize all grids at the coarsest level (level one)
subroutine initlevelone use mod_global_parameters use mod_ghostcells_update use mod_functions_connectivity, only: build_connectivity, getigrids use mod_functions_forest, only: init_forest_root use mod_amr_solution_node, only: alloc_node integer :: iigrid, igrid integer :: itimelevel levmin=1 levmax=1 call init_forest_root call getigrids call build_connectivity ! fill solution space of all root grids do iigrid=1,igridstail; igrid=igrids(iigrid); call alloc_node(igrid) ! initial_condition fetches this block's positions back itself call initial_condition(igrid) end do ! update ghost cells call getbc(global_time,0.d0,ps,iwstart,nwgc) end subroutine initlevelone