Metgrid interpolator bugfixes

Source/ERF.H

@@ -1090,9 +1090,9 @@ private:
     bool        metgrid_use_below_sfc{true};
     bool        metgrid_use_sfc{true};
     bool        metgrid_retain_sfc{false};
-    amrex::Real metgrid_proximity{1000.0};
+    amrex::Real metgrid_proximity{500.0};
     int         metgrid_order{2};
-    int         metgrid_force_sfc_k{0};
+    int         metgrid_force_sfc_k{6};
 
     // 1D CDF output (for ingestion in AMR-Wind)
     static int         output_1d_column;

Source/Initialization/ERF_InitFromMetgrid.cpp

@@ -39,6 +39,13 @@ ERF::init_from_metgrid (int lev)
     if (real_width > real_set_width)
         AMREX_ALWAYS_ASSERT(real_width-real_set_width >= 3);
 
+    // Ensure a reasonable value for the order of the vertical interpolation scheme.
+    AMREX_ALWAYS_ASSERT(metgrid_order > 0 && metgrid_order <= 9);
+
+    // Odd behavior can occur if not using the surface and also omitting some near-surface levels.
+    if (metgrid_force_sfc_k > 0)
+        AMREX_ALWAYS_ASSERT(metgrid_use_sfc);
+
     // Size the SST and LANDMASK
       sst_lev[lev].resize(ntimes);
     lmask_lev[lev].resize(ntimes);