Initialize grid variables for particles
subroutine init_gridvars() use mod_global_parameters use mod_timing integer :: igrid, iigrid tpartc_grid_0=MPI_WTIME() do iigrid=1,igridstail; igrid=igrids(iigrid); allocate(gridvars(igrid)%w(ixGlo1:ixGhi1,ixGlo2:ixGhi2,ixGlo3:ixGhi3,& 1:ngridvars),gridvars(igrid)%wold(ixGlo1:ixGhi1,ixGlo2:ixGhi2,& ixGlo3:ixGhi3,1:ngridvars)) gridvars(igrid)%w = 0.0d0 gridvars(igrid)%wold = 0.0d0 end do call particles_fill_gridvars() if (associated(particles_define_additional_gridvars)) then call particles_fill_additional_gridvars() end if tpartc_grid = tpartc_grid + (MPI_WTIME()-tpartc_grid_0) end subroutine init_gridvars