& theta_Fe_max_hi,
_RL theta_Fe_max_hi
& (theta_Fe_max_hi-theta_Fe_max_lo)*Fe_lim(i,j,k)
theta_Fe_max_hi = 0.04 _d 0