Skip to content

Commit

Permalink
cast minimum lat/lon threshold as a real
Browse files Browse the repository at this point in the history
  • Loading branch information
Peter Bogenschutz committed Jul 17, 2024
1 parent a5aec91 commit f2a7b97
Showing 1 changed file with 2 additions and 2 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -268,9 +268,9 @@ initialize_iop_file(const util::TimeStamp& run_t0,
scorpio::read_var(iop_file,"lon",&iop_file_lon);

const Real rel_lat_err = std::fabs(iop_file_lat - m_params.get<Real>("target_latitude"))/
std::max(m_params.get<Real>("target_latitude"),0.1);
std::max(m_params.get<Real>("target_latitude"),(Real)0.1);
const Real rel_lon_err = std::fabs(std::fmod(iop_file_lon + 360.0, 360.0)-m_params.get<Real>("target_longitude"))/
std::max(m_params.get<Real>("target_longitude"),0.1);
std::max(m_params.get<Real>("target_longitude"),(Real)0.1);
EKAT_REQUIRE_MSG(rel_lat_err < std::numeric_limits<float>::epsilon(),
"Error! IOP file variable \"lat\" does not match target_latitude from IOP parameters.\n");
EKAT_REQUIRE_MSG(rel_lon_err < std::numeric_limits<float>::epsilon(),
Expand Down

0 comments on commit f2a7b97

Please sign in to comment.