It seems like int_snow will be inconsistent when interpolating from initial conditions with different snow layers
- Dominant language
- Fortran
- Stars
- 352
- Forks
- 361
- Avg merge
- 2d 21h
- Merged PRs (30d)
- 7
Description
If you have an initial conditions file that was generated with (say) h2osno_max = 10000 and 12 snow layers, and you use this in a run with a smaller h2osno_max and/or fewer snow layers, it seems like the int_snow value on the initial conditions file will be inconsistent with the new run. In particular, you could have a large value of int_snow that is much larger than would be allowed by the new values of h2osno_max and the number of snow layers.
This is not just a theoretical issue: For example, the NWP compset uses fewer snow layers with an initial conditions file that had more snow layers, so I think this issue could occur there.
I wonder if we should somehow detect if int_snow is too large given the current initial conditions and parameter settings, and if so, reset it?
I'm unsure whether (a) this is really a problem, (b) if so, how big of a problem it is, and (c) what would be the best way to resolve this. @swensosc, do you have any sense here?
Contributor guide
Assessment
This issue has not been assessed yet.