So we probably need to check if internal_heat is associated to avoid the segfault that occurs if DO_GEOTHERMAL = False, and either fail with a descriptive warning or assume geothermal heating is 0 (currently we assume it is 0 anyways in the regional models).
This one has been around for a while. COBALT uses the internal (geothermal) heat argument to calculate an iron flux:
https://github.com/NOAA-CEFI-Regional-Ocean-Modeling/ocean_BGC/blob/1bc8466335573dc62cf89eb69ec982339b0413b2/generic_tracers/generic_COBALT.F90#L4326
However if you trace it back far enough, internal heat is a null pointer by default that only gets allocated if
DO_GEOTHERMAL = Truein the MOM_input: https://github.com/NOAA-CEFI-Regional-Ocean-Modeling/MOM6/blob/3624434308f2863985353091cd02ee274cd18b7b/src/core/MOM.F90#L2802So we probably need to check if internal_heat is associated to avoid the segfault that occurs if DO_GEOTHERMAL = False, and either fail with a descriptive warning or assume geothermal heating is 0 (currently we assume it is 0 anyways in the regional models).