#737: SCHISM Coastal Formulation

[PhilMiller]

Resolves #737

Additions

• CoastalFormulation abstract base class
• Support classes for relevant DataProvider instantiations that CoastalFormulation will need
• SCHISM-based implementation of CoastalFormulation

Changes

• parallel_utils.h gains the MPI size/rank global variables, and has its contents split from the header to a source file
• Bmi_Fortran_Adapter gains a constructor that implements the prototype MPI communicator setup protocol

Testing

Screenshots

Notes

Todos

Checklist

• PR has an informative and human-readable title
• Changes are limited to a single goal (no scope creep)
• Code can be automatically merged (no conflicts)
• Code follows project standards (link if applicable)
• Passes all existing automated tests
• Any change in functionality is tested
• New functions are documented (with a description, list of inputs, and expected output)
• Placeholder code is flagged / future todos are captured in comments
• Project documentation has been updated (including the "Unreleased" section of the CHANGELOG)
• Reviewers requested with the Reviewers tool ➡️

Testing checklist (automated report can be put here)

Target Environment support

• Linux
• macOS

[PhilMiller]

This is totally not ready for review, but it's now a demonstration of code that drives SCHISM through the NGen Fortran BMI adapter, initializing, setting variables, updating, and finalizing.

I think for our deliverable, we'll be able to use NetCDF files and the associated provider rather than the forcing engine, but I'll have to see. Justin, in the meanwhile, this should give a decent sense of what I'd like for an unstructured mesh version of the ForcingsEngineDataProvider. Don't look at anything I've written on that front, especially, as set in stone. It's all meant to be just enough to fill out all aspects of the design.

Note that there's still a shortfall in the SCHISM BMI code that limits this to running on just 1 MPI rank for now, since in partitioned form, some arrays are split while others aren't.

[PhilMiller]

I've found that I need to disable Fortran array bounds checking in the BMI implementation to work with this on GNU Fortran, because the array arguments that come in have. reported bounds of 0.

This maybe arose new in gfortran-14. I haven't gone back to older versions to check.