Modify RHS parameter functions

[kmdeck]

Description

for robustness, replace theta with
theta = max(eps(FT), theta)
in all functions of theta on the RHS, but do not change the state value.

Figure out if we want to issue a warning if theta <0, though.