Returns the derivative of analytic saturation function.
Input: psi0 unsaturated flux pu Output: der_psi derivative of saturating flux with respect to psi0The analytic expression for the saturating flux function and derivative is
psi = (c[1]/c[2])*tanh(c[2]*psi0) + c[3]*psi0 d_psi/d_psi0 = c[1]/cosh(c[2]*psi0)^2 + c[3]For xratio < < 1 the coefficients c are related to
xratio ratio saturated/unsaturated coupling impedance psi_sat saturation value of flux puin the following way
c[1] = 1-xratio c[2] = (1-xratio)/(psi_sat-xratio) c[3] = xratio
function der_saturationAnalytic input Real[:] psi0 "unsaturated flux pu"; input Real[3] c; output Real[size(psi0, 1)] der_psi "derivative d_psi/d_psi0"; end der_saturationAnalytic;