mesh_cart_t – PARALLAX

type, public :: mesh_cart_t

Datatype for a mesh_2d

Components

Type	Visibility	Attributes		Name
real(kind=FP),	public		::	spacing_c	Grid spacing on the actual level (coarse) NOTE: This variable is set public only temporarily and must not be changed from outside of this module. In normal circumstances use the routine get_spacing_c. (It needs to be accessed directly in helmholtz_operators, to circumvent the known issue of OMP scaling with intel compiler)
integer,	public		::	n_points	Number of total mesh points NOTE: This variable is set public only temporarily and must not be changed from outside of this module. In normal circumstances use the routine get_n_points. (It needs to be accessed directly in helmholtz_operators, to circumvent the known issue of OMP scaling with intel compiler)
integer,	public,	allocatable, dimension(:,:,:)	::	index_neighbor	Index of neighbor points
integer,	public,	allocatable, dimension(:)	::	inner_indices	Indices of inner points
integer,	public,	allocatable, dimension(:)	::	boundary_indices	Indices of boundary points
integer,	public,	allocatable, dimension(:)	::	ghost_indices	Indices of ghost points
integer,	public,	allocatable, dimension(:)	::	redblack_indices	Index list of redblack points, eventually needed for red-black gauss seidel smoother First nred entries are indices of red points Following nblack entries are indices of black points Final nboundary+nghost entries are indices of boundary and ghost points

Type-Bound Procedures

procedure, public :: get_phi

interface

public pure module function get_phi(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value real(kind=fp)

procedure, public :: set_phi

interface

public module subroutine set_phi(self, new_phi)

Sets phi value of mesh to a new value

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(inout)			::	self	Instance of the type
real(kind=FP),	intent(in)			::	new_phi	New phi value

procedure, public :: get_lvl

interface

public pure module function get_lvl(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_lvst

interface

public pure module function get_lvst(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_n_points

interface

public pure module function get_n_points(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_spacing_f

interface

public pure module function get_spacing_f(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value real(kind=fp)

procedure, public :: get_spacing_c

interface

public pure module function get_spacing_c(self)

Returns mesh distance of level, spacing_c = spacing_f*lvst

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value real(kind=fp)

procedure, public :: get_size_neighbor

interface

public pure module function get_size_neighbor(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_cart_i

interface

public pure module function get_cart_i(self, ind)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value integer

procedure, public :: get_cart_j

interface

public pure module function get_cart_j(self, ind)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value integer

procedure, public :: get_x

interface

public pure module function get_x(self, ind)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value real(kind=fp)

procedure, public :: get_y

interface

public pure module function get_y(self, ind)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value real(kind=fp)

procedure, public :: get_n_points_inner

interface

public pure module function get_n_points_inner(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_n_points_boundary

interface

public pure module function get_n_points_boundary(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_n_points_ghost

interface

public pure module function get_n_points_ghost(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_xmin

interface

public pure module function get_xmin(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value real(kind=fp)

procedure, public :: get_ymin

interface

public pure module function get_ymin(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value real(kind=fp)

procedure, public :: get_n_points_red

interface

public pure module function get_n_points_red(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: get_n_points_black

interface

public pure module function get_n_points_black(self)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

Return Value integer

procedure, public :: is_inner_point

interface

public pure module function is_inner_point(self, ind)

Indicates if point is inner grid point

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value logical

procedure, public :: is_boundary_point

interface

public pure module function is_boundary_point(self, ind)

Indicates if point is boundary point

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value logical

procedure, public :: is_ghost_point

interface

public pure module function is_ghost_point(self, ind)

Indicates if point is ghost point

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value logical

procedure, public :: get_district

interface

public pure module function get_district(self, ind)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
integer,	intent(in)			::	ind	Mesh index

Return Value integer

procedure, public :: cart_to_mesh_index

interface

public module function cart_to_mesh_index(self, i, j)

Given Cartesian indizes (i,j) returns mesh index if no point could be found 0 is returned

Arguments

Type	Intent		Name
class(mesh_cart_t),	intent(in)	::	self	Instance of the type
integer,	intent(in)	::	i	Horizontal Cartesian index
integer,	intent(in)	::	j	Vertical Cartesian index

Return Value integer

procedure, public :: build_shaped_cart_arrays

interface

public module subroutine build_shaped_cart_arrays(self)

Builds a 2D array where shaped_cart_l(i,j) = l gives the index where cart_i(l) == i and cart_j(l) == j, or zero if no point is found. Used to accelerate cart_to_mesh_index

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(inout)			::	self

procedure, public :: deallocate_shaped_cart_arrays

interface

public module subroutine deallocate_shaped_cart_arrays(self)

Deallocates the shaped_cart_l array to reduce the memory requirement for the grid

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(inout)			::	self

procedure, public :: get_index_neighbor

interface

public pure module function get_index_neighbor(self, ishift, jshift, ind)

Arguments

Type	Intent		Name
class(mesh_cart_t),	intent(in)	::	self	Instance of the type
integer,	intent(in)	::	ishift	Shift in horizontal direction w.r.t. point ind
integer,	intent(in)	::	jshift	Shift in vertical direction w.r.t. point ind
integer,	intent(in)	::	ind	Mesh index

Return Value integer

procedure, public :: determine_interpolation_stencil

interface

public module subroutine determine_interpolation_stencil(self, xp, yp, intorder, ns_max_nearest, ns_stencil, intorder_actual, ind_interpol_stencil)

Finds ns_stencil x ns_stencil indices of interpolation stencil around (xp, yp), if possible of order intorder

Arguments

Type	Intent	Attributes		Name
class(mesh_cart_t),	intent(in)		::	self	Instance of the type
real(kind=FP),	intent(in)		::	xp	x-coordinate of point
real(kind=FP),	intent(in)		::	yp	y-coordinate of point
integer,	intent(in)		::	intorder	Requested stencil order
integer,	intent(in)		::	ns_max_nearest	Size of area that is searched for nearest point
integer,	intent(out)		::	ns_stencil	Size of determined stencil
integer,	intent(out)		::	intorder_actual	Order of determined stencil
integer,	intent(out),	allocatable, dimension(:,:)	::	ind_interpol_stencil	Indices of determined stencil

procedure, public :: get_surrounding_indices

interface

public module subroutine get_surrounding_indices(self, xp, yp, ns, ind_sur, ns_complete, ind_nearest)

Finds ns x ns indices of points surrounding (xp, yp)

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type
real(kind=FP),	intent(in)			::	xp	x-coordinate of point
real(kind=FP),	intent(in)			::	yp	y-coordinate of point
integer,	intent(in)			::	ns	Dimension of surrounding local grid
integer,	intent(out),	optional,	dimension(ns, ns)	::	ind_sur	Indices of surrounding grid points
integer,	intent(out),	optional		::	ns_complete	Dimension of surrounding local grid, that is complete (w/o missing points)
integer,	intent(out),	optional		::	ind_nearest	Index of nearest point on surrounding grid

procedure, public :: normal_to_boundary

interface

public module subroutine normal_to_boundary(self, ind, nbx, nby)

Returns the discrete normal vector to the boundary

Arguments

Type	Intent		Name
class(mesh_cart_t),	intent(in)	::	self	Instance of the type
integer,	intent(in)	::	ind	Mesh index, must be index of a boundary point
integer,	intent(out)	::	nbx	x-component of discrete normal vector
integer,	intent(out)	::	nby	y-component of discrete normal vector

procedure, public :: display => display_mesh_cart

interface

public module subroutine display_mesh_cart(self)

Displays information for mesh

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in)			::	self	Instance of the type

procedure, public :: write_netcdf => write_netcdf_mesh_cart

interface

public module subroutine write_netcdf_mesh_cart(self, fgid)

Writes the mesh to netcdf

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(inout)			::	self	Instance of the type
integer,	intent(in)			::	fgid	File or group id number of existing Netcdf4 file

procedure, public :: read_netcdf => read_netcdf_mesh_cart

interface

public module subroutine read_netcdf_mesh_cart(self, fgid)

Reads mesh from netcdf

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(inout)			::	self	Instance of the type
integer,	intent(in)			::	fgid	File or group id number of existing Netcdf4 file

procedure, public :: init => init_mesh_cart

interface

public module subroutine init_mesh_cart(self, equi, phi, lvl, spacing_f, size_neighbor, size_ghost_layer, extend_beyond_wall, dbgout)

Builds complete mesh

Arguments

Type	Intent	Optional		Name
class(mesh_cart_t),	intent(out)		::	self	Instance of the type
class(equilibrium_t),	intent(inout)		::	equi	Equilibrium
real(kind=FP),	intent(in)		::	phi	Toroidal angle
integer,	intent(in)		::	lvl	Coarseness level
real(kind=FP),	intent(inout)		::	spacing_f	Grid spacing on finest level
integer,	intent(in)		::	size_neighbor	Number of neighbor points determined in each direction
integer,	intent(in)		::	size_ghost_layer	Depth of ghost point layer
logical,	intent(in)		::	extend_beyond_wall	If true, mesh is extended beyond the wall
integer,	intent(in),	optional	::	dbgout	Debug output level

procedure, public :: reorder

interface

public module subroutine reorder(self, grid_reorder_block, dbgout)

Reorders mesh indices for cache optimisation

Arguments

Type	Intent	Optional		Name
class(mesh_cart_t),	intent(inout)		::	self	Instance of the type
integer,	intent(in)		::	grid_reorder_block	Size of basic unshuffled grid square
integer,	intent(in),	optional	::	dbgout	Debug output level

procedure, public :: exchange_mesh_mpi

interface

public module subroutine exchange_mesh_mpi(self, comm, step, mesh_received)

MPI communication of mesh Sends mesh to rank-step and receives from rank+step

Arguments

Type	Intent		Name
class(mesh_cart_t),	intent(in)	::	self	Instance of the type
integer,	intent(in)	::	comm	MPI communicator
integer,	intent(in)	::	step	Step size of MPI communication
type(mesh_cart_t),	intent(out)	::	mesh_received	Mesh received from rank+step

procedure, public :: expose_data

interface

public module subroutine expose_data(self, data_object)

Exposes mesh data through mesh_cart_data_t object

Arguments

Type	Intent	Optional	Attributes		Name
class(mesh_cart_t),	intent(in),		target	::	self	Instance of the type
type(mesh_cart_data_t),	intent(out)			::	data_object	Destination mesh cart data object. Warning: the purpose of this method is to allow access to the data of the instance.

mesh_cart_t Derived Type

Contents

Variables

Type-Bound Procedures

type, public :: mesh_cart_t

Components

Type-Bound Procedures

procedure, public :: get_phi

interface

public pure module function get_phi(self)

Arguments

Return Value real(kind=fp)

procedure, public :: set_phi

interface

public module subroutine set_phi(self, new_phi)

Arguments

procedure, public :: get_lvl

interface

public pure module function get_lvl(self)

Arguments

Return Value integer

procedure, public :: get_lvst

interface

public pure module function get_lvst(self)

Arguments

Return Value integer

procedure, public :: get_n_points

interface

public pure module function get_n_points(self)

Arguments

Return Value integer

procedure, public :: get_spacing_f

interface

public pure module function get_spacing_f(self)

Arguments

Return Value real(kind=fp)

procedure, public :: get_spacing_c

interface

public pure module function get_spacing_c(self)

Arguments

Return Value real(kind=fp)

procedure, public :: get_size_neighbor

interface

public pure module function get_size_neighbor(self)

Arguments

Return Value integer

procedure, public :: get_cart_i

interface

public pure module function get_cart_i(self, ind)

Arguments

Return Value integer

procedure, public :: get_cart_j

interface

public pure module function get_cart_j(self, ind)

Arguments

Return Value integer

procedure, public :: get_x

interface

public pure module function get_x(self, ind)

Arguments

Return Value real(kind=fp)

procedure, public :: get_y

interface

public pure module function get_y(self, ind)

Arguments

Return Value real(kind=fp)

procedure, public :: get_n_points_inner

interface

public pure module function get_n_points_inner(self)

Arguments

Return Value integer

procedure, public :: get_n_points_boundary

interface

public pure module function get_n_points_boundary(self)

Arguments

Return Value integer

procedure, public :: get_n_points_ghost

interface

public pure module function get_n_points_ghost(self)

Arguments