ortho-para H2 phase calculation

[GZALAT]

Description

Phase and quality are determined incorrectly for ortho-para Hydrogen mixtures. For some supercritical regions, "Phase" returns "Two-phase" and "QMass" returns a numerical value. This issue is not present for pure paraHydrogen or normal Hydrogen. It also occurrs only when using the Excel wrapper, not in the Refprop GUI.

Steps to Reproduce

1. Calculate "=REFPROP("Qmass";"OrthoHydrogen;ParaHydrogen|.25;.75";"TP";"SI";70;1.5)" in any MS Excel Cell...

Expected behavior: [What you expect to happen]
Should be "#State is single phase, quality not calculated", not a numerical value

Actual behavior: [What actually happens]
Numerical value is returned.

Versions

REFPROP Version: 10.0 ; file ParaHydrogen.fld adjusted in line 14 to work on common reference point with orthoHydrogen
Operating System and Version: Windows with MS Excel 2013
Access Method: Excel Wrapper

Additional Information

opH2 Phase issue.xlsx

[ianhbell]

What does your excel sheet give for Tc and pc? I can't run it on my machine.

[GZALAT]

Below the screenshot - critical properties for the mixture seem ok:

[GZALAT]

For normal H2, calculated critical properties are a bit off:

[ianhbell]

This seems to be a bug, I agree that this doesn't seem right. @EricLemmon can you please look into this?