It was reported that the full-chemistry TM5-MP is not restartable anymore: stopping and restarting from a restart file changes results compared to a run done in one shot. First tests show that it happens with and without M7, so nothing to do with M7. Runs with only advection, convection, diffusion and wet/dry deposition are fine. So the problem is with sedimentation, boundary, chemistry or emissions.
#3 Updated by Philippe Le Sager over 5 years ago
- Status changed from In Progress to Resolved
- % Done changed from 40 to 100
Solved. The call to ss_monthly_update (sources and sinks update) had to be moved after meteo setup for a consistent vertical remapping of the NOx emissions across the two types of runs (with and without monthly restart). The call to ss_monthly_update must also be before the call to Proces_Update, because of the getDMS call down the line. This has been done, tested, and committed in r347.
#6 Updated by Philippe Le Sager over 4 years ago
I've isolated the problematic lines of code. If I comment the calls to sedimentation_calcv and to deposition_calcv, TM5 is restartable. These calls are found in:
sedimentation.F90/sedimentation_applywhich is called from