#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <fstream>
#include <iomanip>
#include <sstream>
#include <utility>
#include <vector>
#include <AMReX_ParmParse.H>
#include <AMReX_Utility.H>
#include <ERF.H>
#include <ERF_StationSampler.H>
#include <ERF_Plotfile2DCatalog.H>
#include <ERF_Plotfile2DSampledLevel.H>
#include <ERF_Plotfile2DUtils.H>
#include <ERF_Plotfile2DWaterPath.H>
#include <ERF_Constants.H>
#include <ERF_EpochTime.H>
#include <ERF_LatLonMap.H>
#include <ERF_TerrainSurfaceSlab.H>
#include <ERF_Utils.H>

Include dependency graph for ERF_StationSampler.cpp:

Functions
static Real	terrain_at_station (const Array4< const Real > &nd, int klo, Real x, Real y, const GpuArray< Real, AMREX_SPACEDIM > &problo, const GpuArray< Real, AMREX_SPACEDIM > &dx, int i0, int i1, int j0, int j1)

Function Documentation

◆ terrain_at_station()

static Real terrain_at_station	(	const Array4< const Real > &	nd,
		int	klo,
		Real	x,
		Real	y,
		const GpuArray< Real, AMREX_SPACEDIM > &	problo,
		const GpuArray< Real, AMREX_SPACEDIM > &	dx,
		int	i0,
		int	i1,
		int	j0,
		int	j1
	)

static

 {
     const Real fx = (x - problo[0]) / dx[0];
     const Real fy = (y - problo[1]) / dx[1];
     const int  in = std::clamp(static_cast<int>(std::floor(fx)), i0, i1);
     const int  jn = std::clamp(static_cast<int>(std::floor(fy)), j0, j1);
     const Real tx = std::clamp(fx - Real(in), Real(0.0), Real(1.0));
     const Real ty = std::clamp(fy - Real(jn), Real(0.0), Real(1.0));
     return (Real(1.0)-ty) * ( (Real(1.0)-tx)*nd(in,jn  ,klo) + tx*nd(in+1,jn  ,klo) )
          +            ty  * ( (Real(1.0)-tx)*nd(in,jn+1,klo) + tx*nd(in+1,jn+1,klo) );
 }

Referenced by ERF::resolve_station_stencils(), and ERF::sample_stations().

Here is the caller graph for this function:

Functions

Function Documentation

◆ terrain_at_station()