csqrtY computation and meaning

[elizamanelli]

Hi,

I'm currently trying to fully understand the pCqSDHC model and your implementation thereof. So while looking through the associated papers (http://dx.doi.org/10.1016/j.jqsrt.2013.06.015, http://dx.doi.org/10.1016/j.jqsrt.2013.05.034, http://dx.doi.org/10.1016/j.jqsrt.2013.10.015) I've been able to retrace almost every step of the computation. I'm only grappeling with your variable csqrtY which is defined in line 33861:

csqrtY = (Gam2 - iz * Shift2) / (2.0e0 * cte * (1.0e0-eta) * (Gam2 ** 2 + Shift2 ** 2))

I expect this variable csqrtY to be the square root of the variable Y defined just above based on its name and the further usage. Y is defined just above in line 33860:

Y = __ComplexType__(1.0e0 / ((2.0e0 * cte * c2t)) ** 2)

which would simply be

Y = 1/ ( 2 * cte * c2t ) = 1 / (2 * cte * (1 - eta) * c2) = 1 / (2 * cte * (1-eta) * Gam2 + i Shift2)

could you please explain the step I missed in your computation of csqrtY? Your help would be greatly appreciated, thanks!