Grid tools

to be described More...

Namespaces
module	mod_generic_node_type
	Any node type variables.
module	mod_grid_physical_points
	Functions to manage physical points in/out the mesh.

Functions
subroutine, public	mod_grid_get_grid_volume::get_grid_volume (volumes)
	Compute the volume of each cell of a grid.

Detailed Description

Function Documentation

get_grid_volume()

subroutine, public mod_grid_get_grid_volume::get_grid_volume

(

double precision, dimension(:,:,:), intent(out)

volumes

)

Then using immersed boundaries, the volume is restricted to the inner domain. The outer domain counts as zero.

The argument volumes can be a cell field or the component of a face field. In the latter case, the routine considers "shifted cells," i. e. cells centered on face nodes. Although, on the boundary, the volume outside the physical domain and in overlaping regions is not counted. Therefore, the returned argument volumes can be used as follows:

vol = sum(volumes) !! The inner part of the local domain.
 
call global_sum(vol) !! The inner part of the global domain.
 
int = sum(field*volumes) !! The integral of `field` restricted to the
                         !! inner part of the local domain.
 
call global_sum(int) !! The integral of `field` restricted to the inner
                     !! part of the global domain.