public void TestCTL(TestData example1, TestData example2, TestData example3, TestData example4, TestData example5, TestData example6,
                            TestData example7, TestData example8, TestData example9)
        {
            Calcs calc = new Calcs();

            double CTL = calc.GetCTLLiqGas(example1.TempF, example1.Api60);

            Assert.True(EqualsToPrecision(example1.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example2.TempF, example2.Api60);
            Assert.True(EqualsToPrecision(example2.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example3.TempF, example3.Api60);
            Assert.True(EqualsToPrecision(example3.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example4.TempF, example4.Api60);
            Assert.True(EqualsToPrecision(example4.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example5.TempF, example5.Api60);
            Assert.True(EqualsToPrecision(example5.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example6.TempF, example6.Api60);
            Assert.True(EqualsToPrecision(example6.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example7.TempF, example7.Api60);
            Assert.True(EqualsToPrecision(example7.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example8.TempF, example8.Api60);
            Assert.True(EqualsToPrecision(example8.CTL, CTL, 1.0e-5));
            CTL = calc.GetCTLLiqGas(example9.TempF, example9.Api60);
            Assert.True(EqualsToPrecision(example9.CTL, CTL, 1.0e-5));
        }
        public void TestTankRoofVolumeCorrection(TestData example1, TestData example2)
        {
            Calcs calc = new Calcs();

            double BBL = calc.GetBarrelsDueToTankRoof(example1.Grp, example1.Api60, example1.TempF, bblPerApi: example1.FRABblPerApi, refApi: example1.FRAApi);

            Assert.True(EqualsToPrecision(example1.BBL, BBL, 1.0e-2));
            BBL = calc.GetBarrelsDueToTankRoof(example1.Grp, example2.Api60, example2.TempF, roofWgtLb: example2.RoofWtLb);
            Assert.True(EqualsToPrecision(example2.BBL, BBL, 1.0e-2));
        }
Example #3
0
        static void Main(string[] args)
        {
            // Section 12.2 Example 5 in MPSM 11.2.2
            double api60    = Conversions.SGtoAPI(0.5050);
            double temp     = 88.1;
            double vapPress = 159.5;  // At temp=88.1
            double press    = 425;

            // Get compressibility factor
            Calcs calc = new Calcs();

            ComputeCompress(calc, api60, temp, press, vapPress, cplCheck: 1.0133, fCheck: 4.9285e-5);

            // Section 11.2.2 - 11.2.2.8.3M Example
            api60 = Conversions.SGtoAPI(0.5300);
            temp  = Conversions.DegCtoDegF(5.0);
            ComputeCompress(calc, api60, temp, press, vapPress, aCheck: 281093, bCheck: 5.504);
        }
        public void TestVolumeCorrectionFactor(TestData example1, TestData example2, TestData example3, TestData example4, TestData example5)
        {
            Calcs calc = new Calcs();

            double CTPL = calc.GetCTPLFromApiDegFPsig(example1.Grp, example1.Api60, example1.TempF, example1.PressPsig);

            Assert.True(EqualsToPrecision(example1.CTPL, CTPL, 1.0e-3));

            CTPL = calc.GetCTPLFromApiDegFPsig(example2.Grp, example2.Api60, example2.TempF, example2.PressPsig);
            Assert.True(EqualsToPrecision(example2.CTPL, CTPL, 1.0e-3));


            CTPL = calc.GetCTPLFromApiDegFPsig(example3.Grp, example3.Api60, example3.TempF, example3.PressPsig);
            Assert.True(EqualsToPrecision(example3.CTPL, CTPL, 1.6e-3));             // Customer data seems suspect in this case.  Pressure seems rather low.

            CTPL = calc.GetCTPLFromApiDegFPsig(example4.Grp, example4.Api60, example4.TempF, example4.PressPsig);
            Assert.True(EqualsToPrecision(example4.CTPL, CTPL, 1.0e-3));

            CTPL = calc.GetCTPLFromApiDegFPsig(example5.Grp, example5.Api60, example5.TempF, example5.PressPsig);
            Assert.True(EqualsToPrecision(example5.CTPL, CTPL, 1.6e-3));  // Customer data seems suspect in this case.  Pressure seems rather low.
        }
Example #5
0
        static void ComputeCompress(Calcs calc, double api60, double tempF, double pressPsi, double vapPressPsi, double cplCheck = 0, double fCheck = 0, double aCheck = 0, double bCheck = 0)
        {
            double A;
            double B;
            double fp = calc.GetCompressFactorLiqGas(tempF, api60, pressPsi, vapPressPsi, out A, out B);

            Console.WriteLine("Compressibility Factor={0} for Commodity Group {1}", fp, COMMODITY_GROUP.LPG_NGL);
            // TODO NOT ABLE TO REPRODUCE THIS RESULT: fp = 4.9285e-5;
            Console.WriteLine("Print A={0} and B={1}", A, B);
            Console.WriteLine("Print A={0} for Dp in kPa", Conversions.PSItokPa(A));

            double CPL = calc.GetCPLLiquidGas(fp, pressPsi, vapPressPsi);

            CPL = calc.RoundUpAPI11_1(CPL, "cpl");
            Console.WriteLine("Pressure Correction CPL={0} for Commodity Group {1}", CPL, COMMODITY_GROUP.LPG_NGL);
            // TODO NOT ABLE TO REPRODUCE THIS RESULT: cpl = 1.0133


            // Checks
            if (cplCheck != 0)
            {
                Console.WriteLine("CPL Expected {0}, Actual {1} differ by {2}%", cplCheck, CPL, 100 * Math.Abs(CPL - cplCheck) / cplCheck);
            }
            if (fCheck != 0)
            {
                Console.WriteLine("Fp Expected {0}, Actual {1} differ by {2}%", fCheck, fp, 100 * Math.Abs(fp - fCheck) / fCheck);
            }
            if (aCheck != 0)
            {
                Console.WriteLine("A Expected {0}, Actual {1} differ by {2}%", aCheck, Conversions.PSItokPa(A), 100 * Math.Abs(Conversions.PSItokPa(A) - aCheck) / aCheck);
            }
            if (aCheck != 0)
            {
                Console.WriteLine("B Expected {0}, Actual {1} differ by {2}%", bCheck, B, 100 * Math.Abs(B - bCheck) / aCheck);
            }
        }
Example #6
0
        static void Main(string[] args)
        {
            // Check roundups
            double t    = 5.34;
            string uom  = "degC";
            Calcs  calc = new Calcs();

            Console.WriteLine("For temp of {0} {1} roundup is {2} {1}", t, uom, calc.RoundUpAPI11_1(t, uom));
            t = -5.34;
            Console.WriteLine("For temp of {0} {1} roundup is {2} {1}", t, uom, calc.RoundUpAPI11_1(t, uom));
            t   = 10.05;
            uom = "degF";
            Console.WriteLine("For temp of {0} {1} roundup is {2} {1}", t, uom, calc.RoundUpAPI11_1(t, uom));

            // Check other roundups
            t = -0.05;
            int pow10 = -1;

            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUpAPI11_1(t, Math.Pow(10, pow10)), pow10);
            t = -0.049;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUpAPI11_1(t, Math.Pow(10, pow10)), pow10);
            t = -0.051;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUpAPI11_1(t, Math.Pow(10, pow10)), pow10);
            t = 0.35555;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUpAPI11_1(t, Math.Pow(10, pow10)), pow10);
            t = 0.40289;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUpAPI11_1(t, Math.Pow(10, pow10)), pow10);
            t = 24.95;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUp(t, pow10), pow10);
            t     = 2015.0;
            pow10 = 1;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUp(t, pow10), pow10);
            t     = 2017.5;
            pow10 = 0;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUp(t, pow10), pow10);
            t     = 1.374174158511;
            pow10 = -5;
            Console.WriteLine("Value of {0} rounds up to {1} at precision 10e{2}", t, calc.RoundUp(t, pow10), pow10);


            // Check temperature conversions
            t = 32;
            Console.WriteLine("Temp {0} in degF converts to {1} in degC", t, Conversions.DegFtoDegC(t));
            Console.WriteLine("Temp {0} in degF converts to {1} in degK", t, Conversions.DegFtoDegK(t));
            t = -40;
            Console.WriteLine("Temp {0} in degF converts to {1} in degC", t, Conversions.DegFtoDegC(t));
            Console.WriteLine("Temp {0} in degF converts to {1} in degK", t, Conversions.DegFtoDegK(t));
            t = 60;
            Console.WriteLine("Temp {0} in degF converts to {1} in degC", t, Conversions.DegFtoDegC(t));
            Console.WriteLine("Temp {0} in degF converts to {1} in degK", t, Conversions.DegFtoDegK(t));
            t = 212;
            Console.WriteLine("Temp {0} in degF converts to {1} in degC", t, Conversions.DegFtoDegC(t));
            Console.WriteLine("Temp {0} in degF converts to {1} in degK", t, Conversions.DegFtoDegK(t));
            t = 0;
            Console.WriteLine("Temp {0} in degC converts to {1} in degF", t, Conversions.DegCtoDegF(t));
            t = -40;
            Console.WriteLine("Temp {0} in degC converts to {1} in degF", t, Conversions.DegCtoDegF(t));
            t = 15;
            Console.WriteLine("Temp {0} in degC converts to {1} in degF", t, Conversions.DegCtoDegF(t));
            t = 100;
            Console.WriteLine("Temp {0} in degC converts to {1} in degF", t, Conversions.DegCtoDegF(t));
            t   = 32;
            uom = "degF";
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68 - rounded up", t, uom, calc.RoundUpAPI11_1(Conversions.TempITS90toITPS68(t, uom), uom));
            t = -40;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68 - rounded up", t, uom, calc.RoundUpAPI11_1(Conversions.TempITS90toITPS68(t, uom), uom));
            t = 60;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68 - rounded up", t, uom, calc.RoundUpAPI11_1(Conversions.TempITS90toITPS68(t, uom), uom));
            t = 212;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68", t, uom, Conversions.TempITS90toITPS68(t, uom));
            t   = 0;
            uom = "degC";
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68", t, uom, Conversions.TempITS90toITPS68(t, uom));
            t = -40;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68", t, uom, Conversions.TempITS90toITPS68(t, uom));
            t = 15;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68", t, uom, Conversions.TempITS90toITPS68(t, uom));
            t = 100;
            Console.WriteLine("Temp {0} in {1} ITS90 converts to {2} {1} ITPS68", t, uom, Conversions.TempITS90toITPS68(t, uom));

            // Check Density Conversions
            double d = 10;

            uom = "API";
            Console.WriteLine("Density {0} in {1} converts to {2} in SG", d, uom, Conversions.APItoSG(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in kg/m3", d, uom, Conversions.APItoKgm3(d));
            d = 1.8400000333786;
            Console.WriteLine("Density {0} in {1} converts to {2} in SG", d, uom, Conversions.APItoSG(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in kg/m3", d, uom, Conversions.APItoKgm3(d));
            d   = 1;
            uom = "SG";
            Console.WriteLine("Density {0} in {1} converts to {2} in API", d, uom, Conversions.SGtoAPI(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in kg/m3", d, uom, Conversions.APItoKgm3(Conversions.SGtoAPI(d)));
            d   = 999.016;
            uom = "kg/m3";
            Console.WriteLine("Density {0} in {1} converts to {2} in API", d, uom, Conversions.Kgm3toAPI(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in SG", d, uom, Conversions.APItoSG(Conversions.Kgm3toAPI(d)));
            Console.WriteLine("Density {0} in {1} converts to {2} in lb/gal", d, uom, Conversions.KgM3ToLbGal(d));
            d = 610.6;
            Console.WriteLine("Density {0} in {1} converts to {2} in API", d, uom, Conversions.Kgm3toAPI(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in SG", d, uom, Conversions.APItoSG(Conversions.Kgm3toAPI(d)));
            d = 800;
            Console.WriteLine("Density {0} in {1} converts to {2} in API", d, uom, Conversions.Kgm3toAPI(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in SG", d, uom, Conversions.APItoSG(Conversions.Kgm3toAPI(d)));
            d   = 0.5250;
            uom = "SG";
            Console.WriteLine("Density {0} in {1} converts to {2} in API", d, uom, Conversions.SGtoAPI(d));
            Console.WriteLine("Density {0} in {1} converts to {2} in kg/m3", d, uom, Conversions.SGtoKgm3(d));

            // Get kCoeffs and print them
            var grp = COMMODITY_GROUP.CRUDE_OIL;
            var kc  = calc.GetKCoeffs(grp);

            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);
            grp = COMMODITY_GROUP.FUEL_OILS;
            kc  = calc.GetKCoeffs(grp);
            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);
            grp = COMMODITY_GROUP.GASOLINES;
            kc  = calc.GetKCoeffs(grp);
            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);
            grp = COMMODITY_GROUP.JET_FUELS;
            kc  = calc.GetKCoeffs(grp);
            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);
            grp = COMMODITY_GROUP.LUBRICATING_OIL;
            kc  = calc.GetKCoeffs(grp);
            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);
            grp = COMMODITY_GROUP.TRANSITION_ZONE;
            kc  = calc.GetKCoeffs(grp);
            Console.WriteLine("Commodity Group {0} has k0={1}, k1={2}, k2={3}", kc.CommodityGroup, kc.k0, kc.k1, kc.k2);

            // Get GasLiqProperties and print them
            var gasLiq = LIQ_GAS_FLUID.EE_68_32;
            var lgp    = calc.GetLiqGasProps(gasLiq);

            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.ETHANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.EP_65_35;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.EP_35_65;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.PROPANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.iBUTANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.nBUTANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.iPENTANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.nPENTANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.iHEXANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.nHEXANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);
            gasLiq = LIQ_GAS_FLUID.nHEPTANE;
            lgp    = calc.GetLiqGasProps(gasLiq);
            Console.WriteLine("Liquid Gas {0} has relDens60={1}, tempCritK={2}, compCrit={3}, densCrit={4}, k1={5}, k2={6}, k1={7}, k2={8}", lgp.Fluid, lgp.relDens60, lgp.tempCritK, lgp.comprFactCrit, lgp.densCrit, lgp.k1, lgp.k2, lgp.k3, lgp.k4);

            // Get Vapor Pressure Correlation Parameters
            var vpCorrParms = RulesLoader.GetVapPressCorrParams();

            foreach (var vpCorrParm in vpCorrParms)
            {
                Console.WriteLine("{0}{1}Relative Density{2}{3} - A0={4},A1={5},A2={6},B0={7},B1={8},B2={9}",
                                  vpCorrParm.RelDensityRange.Min, vpCorrParm.RelDensityRange.MinCompare == COMPARE.INCLUDE?"<=":"<",
                                  vpCorrParm.RelDensityRange.MaxCompare == COMPARE.INCLUDE ? "<=" : "<", vpCorrParm.RelDensityRange.Max,
                                  vpCorrParm.A0, vpCorrParm.A1, vpCorrParm.A2, vpCorrParm.B0, vpCorrParm.B1, vpCorrParm.B2);
            }

            // Validate range checks
            double pres = 0;
            double temp = 0;
            double dens = 800;
            double api  = Conversions.Kgm3toAPI(dens);

            grp = COMMODITY_GROUP.CRUDE_OIL;
            Console.WriteLine("Testing for commodity group {0}", grp);
            double vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);  // All good

            Console.WriteLine("Values in range test completed");
            // Temp to low
            temp = -58.1;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low temp for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Temp {0} is too low", temp);
            }
            // Temp to high
            temp = 302.1;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high temp for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e) {
                Console.WriteLine("Temp {0} is too high", temp);
            }
            temp = 0;
            // Press to low
            pres = -0.1;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low press for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Press {0} is too low", pres);
            }
            // Press to high
            pres = 1500.1;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high press for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Press {0} is too high", pres);
            }
            // Density too low
            dens = 610.5;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 1163.5;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            pres = 0;
            temp = 0;
            dens = 900;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.FUEL_OILS;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);              // All good
            Console.WriteLine("Values in range test completed");
            // Density too low
            dens = 818.3126;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 1163.6;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            dens = 800;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.JET_FUELS;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);              // All good
            Console.WriteLine("Values in range test completed");
            // Density too low
            dens = 787.5194;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 838.3127;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            dens = 780;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.TRANSITION_ZONE;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);              // All good
            Console.WriteLine("Values in range test completed");
            // Density too low
            dens = 770.351;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 787.5195;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            dens = 650;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.GASOLINES;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);              // All good
            Console.WriteLine("Values in range test completed");
            // Density too low
            dens = 610.5;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 770.352;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            dens = 900;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.LUBRICATING_OIL;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);              // All good
            Console.WriteLine("Values in range test completed");
            // Density too low
            dens = 800.8;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density too high
            dens = 1163.5;
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres);
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }
            // LPG and NGLs
            pres = 150;
            double vapPres = 40;

            temp = -40;
            dens = 500;
            api  = Conversions.Kgm3toAPI(dens);
            grp  = COMMODITY_GROUP.LPG_NGL;
            Console.WriteLine("Testing for commodity group {0}", grp);
            vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres, vapPres);  // All good
            Console.WriteLine("Values in range test completed");
            // Temp to low
            temp = -50.9;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres, vapPres);                 // All good
                Console.WriteLine("Test for low temp for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Temp {0} is too low", temp);
            }
            // Temp to high
            temp = 199.5;
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres, vapPres);                 // All good
                Console.WriteLine("Test for high temp for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Temp {0} is too high", temp);
            }
            // Density to low
            dens = Conversions.SGtoKgm3(0.349);
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres, vapPres);                 // All good
                Console.WriteLine("Test for low density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too low", dens);
            }
            // Density to high
            dens = Conversions.SGtoKgm3(0.689);
            api  = Conversions.Kgm3toAPI(dens);
            try
            {
                vcf = calc.GetCTPLFromApiDegFPsig(grp, api, temp, pres, vapPres);                 // All good
                Console.WriteLine("Test for high density for {0} failed", grp);
            }
            catch (ArgumentOutOfRangeException e)
            {
                Console.WriteLine("Density {0} is too high", dens);
            }

            Console.WriteLine("Press ENTER to terminate");
            Console.ReadLine();
        }