Gas Mixing Ratio Feedbacks

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).

gas_mix_ratio(1,85) (aka atm%gas_mix_ratio(1,85,4))
SCALE_ABSORB (SOCRATES)
gas_frac_rescaled(1,85,4) rescaled gas mass fraction
SOLVE_BAND_K_EQV (SOCRATES)
k_eqv(1,85) Equivalent extinction
SOLVE_BAND_K_EQV (SOCRATES)
k_grey_tot(1,85) (aka ss_prop%k_grey_tot(1,85,0)) Grey extinction
SINGLE_SCATTERING (SOCRATES)
ss_prop%tau(1,85,0) (aka tau(1,85)) Optical depth
TRANS_SOURCE_COEFF (SOCRATES)
source_coeff(1,85,2) Source coefficients in two-stream equations
IR_SOURCE (SOCRATES)
s_down(1,85) Downward source function
SOLVER_HOMOGEN_DIRECT (SOCRATES)
h(1,85) Working array
SOLVER_HOMOGEN_DIRECT (SOCRATES)
flux_total_incr(1,171) (aka flux_total_part, flux_total) Total flux
AUGMENT_RADIANCE (SOCRATES)
radout%flux_up(1,85,1) Upward flux
LW_RAD
lw_down(1,1) Downwards surface flux
SF_EXP_L (JULES)
radnet(1), radnet_tile(1,9) Surface net radiation on tiles
SF_FLUX (JULES)
dtstar(1) Change in TSTAR over timestep
SF_FLUX (JULES)
fqw(1,1,1) FQW(,1) is total water flux from surface, 'E'
BDY_IMPL4
q(1,1,1) (aka qw)
BDY_IMPL4 (later call from above)
q_latest(1,1,1)
ATMOS_PHYSICS2
q_star(1,1,1)
SET_STAR_ZERO_LEVEL Simple q_star(:,:,0) = q_star(:,:,1)
q(1,1,0) (aka qstar) I think qstar = -(w'q')/u_*
EG_Q_TO_MIX
m_star(1,1,0)
EG_R_S
R_m_v_a(1,1,1), R_m_cf_a(1,1,1) (aka S_m_v, S_m_cf)
UPDATE_MOISTURE
m_v(1,1,1), m_cf(1,1,1)
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)

Summary

Simplified feedback

Change GAS_MIX_RATIO
Alters absorption, diffusion of air (SOCRATES)
Alters radiation budget (SOCRATES)
Alters evapouration of water at surface (JULES) (Probably turbulence in general)
Alters moisture
Alters the other dynamical fields like U