Changes

Summary

Commented out the k=1 boundary conditions in diffusion_zm_lhs and diffusion_zt_lhs, since these are not currently fed into the final LHS matrix. I included comments to describe why this was done. Results are b4b. (details)

Commit db7584133b87452bedc0538fed00f339dca49ae3 by benjamin.andrew.stephens

Commented out the k=1 boundary conditions in diffusion_zm_lhs and diffusion_zt_lhs, since these are not currently fed into the final LHS matrix.  I included comments to describe why this was done.  Results are b4b.

See https://github.com/larson-group/clubb/issues/920#issuecomment-942735433.

diffusion.F90 (diff)