ClimaAtmos.jl is a library for building atmospheric circulation models that is designed from the outset to leverage data assimilation and machine learning tools. We welcome contributions!
1785 set the iteration to 3 and tolerance to 0.003 without an initial guess. With the previous temperature as an initial guess the simulation breaks, see this build for example. It seems temperature can change by 4K in one timestep, which is not physical. One idea is to check the previous q_tot to see if the change is from precipitation.
We should use 1 iteration for the thermodynamic state, and use the temperature at the previous time step as the initial condition for the iterations.