Previously additional microphysics fields that are to be read from the forcing (rain, snow and graupel) were not updated along the boundaries, leaving them at their initial values. this commit should fix this. if mean_fields=True the fields are now also averaged correctly along the boundary.
As discussed here https://github.com/NCAR/icar/pull/61#issuecomment-552944452 I've built a separate fix for this. I will remove the handling of the ice field from the N_from_forcing branch to avoid merging conflicts later down the line.
Previously additional microphysics fields that are to be read from the forcing (rain, snow and graupel) were not updated along the boundaries, leaving them at their initial values. this commit should fix this. if mean_fields=True the fields are now also averaged correctly along the boundary.
As discussed here https://github.com/NCAR/icar/pull/61#issuecomment-552944452 I've built a separate fix for this. I will remove the handling of the ice field from the N_from_forcing branch to avoid merging conflicts later down the line.