Closed arbennett closed 3 years ago
Occasionally, when using JRDN snow layering I have run into segfaults on reference to mw_exp (only when doing scalar solutions) in the line:
mw_exp
iLayerLiqFluxSnow(iLayer) = k_snow*relSaturn**mw_exp
after some debugging it turns out this loss of access to mw_exp occurs directly following the pack statement:
pack
ixLayerDesired = pack(ixLayerState, ixSnowOnlyHyd/=integerMissing)
I have not been able to determine the exact reason for this error, but have developed a workaround, which is provided in this PR.
Solves #446
Occasionally, when using JRDN snow layering I have run into segfaults on reference to
mw_exp(only when doing scalar solutions) in the line:after some debugging it turns out this loss of access to
mw_expoccurs directly following thepackstatement:I have not been able to determine the exact reason for this error, but have developed a workaround, which is provided in this PR.
Solves #446