erf_earth_wind Namespace Reference

Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void	grid_to_earth (amrex::Real u, amrex::Real v, amrex::Real cos_alpha, amrex::Real sin_alpha, amrex::Real &u_east, amrex::Real &v_north) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	horizontal_speed (amrex::Real u_east, amrex::Real v_north) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	meteorological_direction (amrex::Real u_east, amrex::Real v_north, amrex::Real missing_value) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void	grid_rotation_at (const amrex::Array4< const amrex::Real > &ll, int ilo, int ihi, int i, int j, amrex::Real &cos_alpha, amrex::Real &sin_alpha) noexcept

Function Documentation

grid_rotation_at()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_earth_wind::grid_rotation_at ( const amrex::Array4< const amrex::Real > &  ll, int  ilo, int  ihi, int  i, int  j, amrex::Real &  cos_alpha, amrex::Real &  sin_alpha  )

noexcept

{
    const int im = (i - 1 > ilo) ? i - 1 : ilo;
    const int ip = (i + 1 < ihi) ? i + 1 : ihi;
    const amrex::Real inv_di = (ip > im) ? amrex::Real(1.0) / amrex::Real(ip - im)
                                         : amrex::Real(0.0);
 
    // Displacement on the sphere per cell in i: east and north components, in
    // degrees of latitude (the common factor, the Earth's radius, cancels).
    // The longitude difference is taken the short way round the antimeridian.
    amrex::Real dlon = ll(ip,j,0,1) - ll(im,j,0,1);
    if (dlon >  amrex::Real(180.0)) { dlon -= amrex::Real(360.0); }
    if (dlon < -amrex::Real(180.0)) { dlon += amrex::Real(360.0); }
    const amrex::Real lat   = ll(i,j,0,0);
    const amrex::Real north = (ll(ip,j,0,0) - ll(im,j,0,0)) * inv_di;
    const amrex::Real east  = dlon * inv_di * std::cos(lat * PI / amrex::Real(180.0));
    const amrex::Real norm  = std::sqrt(east*east + north*north);
 
    if (norm > amrex::Real(0.0)) {
        cos_alpha = east  / norm;
        sin_alpha = north / norm;
    } else {
        cos_alpha = amrex::Real(1.0);
        sin_alpha = amrex::Real(0.0);
    }
}

Referenced by fill_grid_rotation(), and grid_rotation_from_latlon().

Here is the caller graph for this function:

grid_to_earth()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_earth_wind::grid_to_earth ( amrex::Real  u, amrex::Real  v, amrex::Real  cos_alpha, amrex::Real  sin_alpha, amrex::Real &  u_east, amrex::Real &  v_north  )

noexcept

{
    u_east  = u * cos_alpha - v * sin_alpha;
    v_north = u * sin_alpha + v * cos_alpha;
}

Referenced by near_surface_diagnostics::fill_wind_10m(), and ERF::sample_stations().

Here is the caller graph for this function:

horizontal_speed()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_earth_wind::horizontal_speed ( amrex::Real  u_east, amrex::Real  v_north  )

noexcept

{
    return std::sqrt(u_east * u_east + v_north * v_north);
}

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

Here is the caller graph for this function:

meteorological_direction()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_earth_wind::meteorological_direction ( amrex::Real  u_east, amrex::Real  v_north, amrex::Real  missing_value  )

noexcept

{
    if (horizontal_speed(u_east, v_north) <= amrex::Real(1.0e-12)) {
        return missing_value;
    }
    amrex::Real direction = amrex::Real(270.0)
        - std::atan2(v_north, u_east) * amrex::Real(180.0) / PI;
    while (direction < amrex::Real(0.0)) {
        direction += amrex::Real(360.0);
    }
    while (direction >= amrex::Real(360.0)) {
        direction -= amrex::Real(360.0);
    }
    return direction;
}

Referenced by ERF::sample_stations().

Here is the call graph for this function:

Here is the caller graph for this function: