Monitoring the convergence of the iterative solution method along a vertical cut at r = 10 au as a function of the vertical hydrogen nuclei column density for the model with medium strong mixing, α_mix = 10⁻³. The upper plots show the changing gas and dust temperature as a function of iteration number. The lower plots shows the concentrations of a few selected oxygen-bearing gas and ice species. The left row shows the first few iterations and the right row shows the convergence towards iteration 1000, which is considered as the final solution. The different full, dashed and dotted lines correspond to the different numbers of iterations carried out as indicated.