function correlation_open_expander_epsilon_s input Real N_rot; input Real rho; input Real log_rp; output Real epsilon_s; end correlation_open_expander_epsilon_s;