I changed the tile IDs for the none plant tiles and it changed results. It turned out that because the tile IDs, which are used as vertical levels, are changed then the data in `ainitial' file are not read in - because the tile IDs/vertical levels don't match. This included stash item 00,384 the `SNOW LAYER TEMPERATURE ON TILES (K)', which writes the field tsnowlayer_p in the reconfiguration.
| tsnowlayer_p(27,9,1) | ||
| ↓ | ATMOS_PHYSICS2 | |
| tsnow/tsnow_surft(l=27,n=9,1) | ||
| ↓ | SF_EXCH (JULES) | |
| ts1_elev/tsurf(l=27) | ||
| ↓ | SF_FLUX (JULES) | |
| surf_ht_flux (with l=27) | ||
| ↓ | SF_FLUX (JULES) | |
| dtstar(27) (with i=3 & j=5) | ||
| ↓ | SF_FLUX (JULES) | |
| ftl(3,5,1) | ||
| ↓ | BDY_IMPL4 | |
| dtl(3,5,1) | ||
| ↓ | BDY_IMPL4 | |
| tl(3,5,1) | ||
| ↓ | BDY_IMPL4 | |
| T_latest(3,5,1) | ||
| ↓ | ATMOS_PHYSICS2 | |
| theta_star(3,5,1) | ||
| ↓ | EG_R_S | |
| r_theta(3,5,0:1)/s_thetav(3,5,0:1) | ||
| ↓ | EG_HELM_VAR_RHS | |
| rhs(3,5,1) | ||
| ↓ | EG_BICGSTAB | |
| exner_prime_term(1,1,1) | ||
| ↓ | EG_ADD_PRESSURE | |
| u_np1(1,1,1) | U at next timestep | |
| ↓ | ATM_STEP_4A | Standard update of timestep |
| u(1,1,1) |
