computational_domain.f90

Example showing how create different types of computational domains

#include "../src/cpp_macros.h"
 
program computational_domain
  use m_af_all
 
  implicit none
 
  integer    :: n_cell = 4 ! Boxes contain n_cell^dim cells
  type(af_t) :: tree       ! Will contain the quad/octree grid
  integer    :: grid_size(NDIM)
  real(dp)   :: domain_size(NDIM)
  logical    :: periodic(NDIM)
 
  ! Create mesh 1: two boxes along x-direction
  grid_size(:) = n_cell
  grid_size(1) = 2 * n_cell
  domain_size = 1.0_dp * grid_size
 
  call af_add_cc_variable(tree, "phi")
  call af_init(tree, n_cell, domain_size, grid_size, mem_limit_gb=0.1_dp)
  call af_write_silo(tree, "output/computational_domain_" // dimname // "_1")
  call af_destroy(tree)
 
  ! Create mesh 2: two boxes along y-direction
  grid_size(:) = n_cell
#if NDIM > 1
  grid_size(2) = 2 * n_cell
#endif
  domain_size = 1.0_dp * grid_size
 
  call af_add_cc_variable(tree, "phi")
  call af_init(tree, n_cell, domain_size, grid_size, mem_limit_gb=0.1_dp)
  call af_write_silo(tree, "output/computational_domain_" // dimname // "_2")
  call af_destroy(tree)
 
  ! Create mesh 3: Two boxes along x-direction that are fully periodic
  grid_size(:) = n_cell
  grid_size(1) = 2 * n_cell
  periodic(:)  = .true.
  domain_size  = 1.0_dp * grid_size
 
  call af_add_cc_variable(tree, "phi")
  call af_init(tree, n_cell, domain_size, grid_size, mem_limit_gb=0.1_dp)
  call af_write_silo(tree, "output/computational_domain_" // dimname // "_3")
  call af_destroy(tree)
 
end program computational_domain