double arshup(double g, double p, int *ierr);
Compute pressure factor based upon pressure factor .
ENTRY: Specific heat constant .
ENTRY: Pressure fraction .
EXIT: (optional) Return status code.
Specific heat ratio .
Given out of range.
Return Value
Density fraction or NAN on error.