ERF_Tagging.cpp File Reference

#include <ERF.H>
#include <ERF_Derive.H>

Include dependency graph for ERF_Tagging.cpp:

Functions
void	tag_on_distance_from_eye (const Geometry &cgeom, TagBoxArray *tags, const Real eye_x, const Real eye_y, const Real rad_tag)
	Tag cells based on 2D distance from the storm eye. More...

Function Documentation

tag_on_distance_from_eye()

void tag_on_distance_from_eye	(	const Geometry &	cgeom,
		TagBoxArray *	tags,
		const Real	eye_x,
		const Real	eye_y,
		const Real	rad_tag
	)

Tag cells based on 2D distance from the storm eye.

Parameters

[in]	cgeom	Geometry defining the domain.
[out]	tags	Array of tagged cells for refinement.
[in]	eye_x	X-coordinate of the eye center.
[in]	eye_y	Y-coordinate of the eye center.
[in]	rad_tag	Radius within which cells are tagged.

Tag cells based on distance from the hurricane eye.

Parameters

	cgeom	Geometry used for coordinates.
[out]	tags	TagBoxArray to be filled.
	eye_x	x-coordinate of the eye.
	eye_y	y-coordinate of the eye.
	rad_tag	Radius within which to tag.

{
    const auto dx      = cgeom.CellSizeArray();
    const auto prob_lo = cgeom.ProbLoArray();
 
    for (MFIter mfi(*tags); mfi.isValid(); ++mfi) {
        TagBox& tag = (*tags)[mfi];
        auto tag_arr = tag.array();  // Get device-accessible array
 
        const Box& tile_box = mfi.tilebox(); // The box for this tile
 
        ParallelFor(tile_box, [=] AMREX_GPU_DEVICE(int i, int j, int k) {
            // Compute cell center coordinates
            Real x = prob_lo[0] + (i + myhalf) * dx[0];
            Real y = prob_lo[1] + (j + myhalf) * dx[1];
 
            Real dist = std::sqrt((x - eye_x)*(x - eye_x) + (y - eye_y)*(y - eye_y));
 
            if (dist < rad_tag) {
                tag_arr(i,j,k) = TagBox::SET;
            }
        });
    }
}

Referenced by ERF::ErrorEst(), and ERF::HurricaneTracker().

Here is the call graph for this function:

Here is the caller graph for this function: