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