I had a problem with the gas mixing ratio which changed U etc in the 5th timestep (one after UKCA and the same as radiation).
| alf,cs(1,1,1) | ||
| ↓ | EG_BICGSTAB | |
| exner_prime(1,1,1) | ||
| ↓ | EG_ADD_PRESSURE | |
| thetav(1,1,1) (aka thetav_np1) | ||
| ↓ | EG_SISL_INIT_UVW | |
| theta_wet(1,1,1) | ||
| ↓ | EG_SISL_INIT_UVW | |
| work1(1,1,1) | ||
| ↓ | EG_SISL_INIT_UVW | |
| R_u_a(1,1,1) | I think this is U in Cartesian coordinates | |
| ↓ | EG_ADD_PRESSURE | |
| u_np1(1,1,1) | U at next timestep | |
| ↓ | ATM_STEP_4A | Standard update of timestep |
| u(1,1,1) |
