Setting myself notes so I have code that I can easily pick up when I come back in the fall.
Finished up a run of the code when controlling when the flux limitation in substrate is implemented to only occur. It seems to not be the proper fix so I think I’m going to revert to the original if statement spanning over all of domain 1. Flux is now going the wrong way (from particle into substrate).
Since I don’t want to make any huge changes to the model before I leave for the summer, I simply did a quick check run to make sure everything functions properly for the moment.
Have my current summary of what I’ve done in a powerpoint now for my reference later – if I don’t have time during finals to come in, then this may be the last update until the fall.