EAMxx: hremap cannot handle fields without COL tag

[mahf708]

Details

%YAML 1.1
---
filename_prefix: new_test
horiz_remap_file: /global/cfs/cdirs/e3sm/inputdata/atm/scream/maps/map_ne4pg2_to_ne2pg2_c20240902.nc
iotype: pnetcdf
averaging_type: average
max_snapshots_per_file: 1 # only one snapshot per file
fields:
  physics_pg2:
    field_names:
    - T_mid_vert_avg_dp_weighted
    - T_mid_vert_avg_dz_weighted
    - T_mid_vert_sum_dz_weighted
output_control:
  frequency: 1
  frequency_units: nsteps
restart:
  force_new_file: true

@bartgol, you're not completely out of the woods yet ;):P

the hremap cannot handle rank-0 fields, trying T_mid_vert_sum_dz_weighted_horiz_avg in the above erros when hremapping, but is ok otherwise

0: idim>=0 && idim<m_rank
0: components/eamxx/src/share/field/field_layout.cpp:250
0: Error! Index out of bounds.

Originally posted by @mahf708 in #7507 (comment)

[bartgol]

I suspect this may be fixed by the abstract_grid.hpp mods in PR #7515 . I just merged it, so @mahf708 can you please rebase and try again and see if the error persists?

[mahf708]

Still same issue with a rebase :(

[bartgol]

Can you paste a stacktrace of some sort?

[mahf708]

Can you paste a stacktrace of some sort?

Oh, yeah, let me actually run it through gdb

[mahf708]

Does this help?

#0  0x00007ffff7e2f601 in __cxa_throw () from /usr/lib64/libstdc++.so.6
#1  0x00000000004223ff in void throw_exception<std::runtime_error>(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) ()
#2  0x0000000001aa2248 in scream::FieldLayout::reset_dim(int, int) ()
#3  0x0000000001c2c203 in scream::HorizInterpRemapperBase::create_ov_fields() ()
#4  0x0000000001c2d2f8 in scream::HorizInterpRemapperBase::registration_ends_impl() ()
#5  0x0000000001c1a9ae in scream::CoarseningRemapper::registration_ends_impl() ()
#6  0x0000000001c03f82 in scream::AbstractRemapper::registration_ends() ()
#7  0x0000000001077542 in scream::AtmosphereOutput::AtmosphereOutput(ekat::Comm const&, ekat::ParameterList const&, std::shared_ptr<scream::FieldManager const> const&, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) ()
#8  0x00000000010586fe in scream::OutputManager::setup(std::shared_ptr<scream::FieldManager> const&, std::set<std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >, std::less<std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > >, std::allocator<std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > > > const&) ()
#9  0x0000000000e03b0e in scream::control::AtmosphereDriver::initialize_output_managers() ()
#10 0x000000000058c06c in scream_init_atm ()
#11 0x000000000058888a in atm_comp_mct::atm_init_mct (eclock=..., cdata=..., x2a=..., a2x=..., nlfilename=..., _nlfilename=_nlfilename@entry=6)
    at /pscratch/sd/m/mahf708/e3sm-repo/test-pr/components/eamxx/src/mct_coupling/atm_comp_mct.F90:282
#12 0x00000000004fbf5e in component_mod::component_init_cc (eclock=..., comp=..., comp_init=0x5877d0 <atm_comp_mct::atm_init_mct>, infodata=..., nlfilename=..., seq_flds_x2c_fluxes=..., 
    seq_flds_c2x_fluxes=..., _nlfilename=6, _seq_flds_x2c_fluxes=0, _seq_flds_c2x_fluxes=0) at /pscratch/sd/m/mahf708/e3sm-repo/test-pr/driver-mct/main/component_mod.F90:258
#13 0x00000000004ec1c9 in cime_comp_mod::cime_init () at /pscratch/sd/m/mahf708/e3sm-repo/test-pr/driver-mct/main/cime_comp_mod.F90:1502
#14 0x00000000004c38ed in cime_driver () at /pscratch/sd/m/mahf708/e3sm-repo/test-pr/driver-mct/main/cime_driver.F90:122
#15 main (argc=<optimized out>, argv=<optimized out>) at /pscratch/sd/m/mahf708/e3sm-repo/test-pr/driver-mct/main/cime_driver.F90:23
#16 0x00007fffeea3e1fd in __libc_start_main () from /lib64/libc.so.6
#17 0x00000000004dd78a in _start () at ../sysdeps/x86_64/start.S:120

[bartgol]

Ok, the issue is not with the rank of the field, but rather with the fact that we are trying to remap something that doesn't have the COL dimension. I should have a fix soon.