Segmentation fault when interpolating regional initial conditions

[daniloceano]

Dear all,

I’m running into a blocking error while generating regional initial conditions for Regional MPAS.

What I’m doing

• MPAS‑A version: 8.2.2
• Global variable‑resolution mesh 200 km (outer) → 10 km (inner) created with JIGSAW.
• Used "create_region" to clip the high‑resolution patch over my domain; static fields created successfully.
• Input data: hourly ERA5 GRIB (pressure‑level + single‑level) downloaded via CDS‑API. WPS handled link_grib.csh and ungrib.exe with Vtable.ECMWF.
• When I launch
  mpiexec ‑n <nproc> ./init_atmosphere init_atmosphere.nl
  to interpolate regional ICs, the run dies with segmentation fault

The full error log is on "nohup.out.txt"

Any help would be greatly appreciated.

Thanks,

Danilo

 

[Ming Chen]

Your namelist.init_atmosphere looks fine to me.

The error message "Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation" indicates fundamental problems in your case.

Without looking at your datafile, I am not sure whether the 200 km → 10 km mesh you created is correct or not.

As a good starting point, can you run a test case of regional MPAS using a variable-resolution mesh we provide in MPAS home page?

This test case would serve as the basis that can at least demonstrate your approach is correct. If this test case is done successfully, we will know for sure that you follow the correct steps to run regional MPAS. Then we can further explore what is wrong with your failed case.