Скрипт MasterScada

Добрый день. Очень нужна ваша помощь. У меня не выполняется скрипт, хотя компиляция пройдена успешно. Не могу понять в чем проблема. Заранее большое спасибо.

Код:

using System; using MasterSCADA.Script.FB; using MasterSCADA.Hlp; using FB; using System.Linq; public partial class ФБ : ScriptBase { public override void Execute() { //температура на тарелках double[] Temp = new double[57]; double t1 = 350;double t23 = 320;double t35 = 240; double t47 = 140;double t56 = 100; double p1 = 95;double p56 = 10;double pTR1 = 100;double pTR2 = 100;double pTLT = 1;double p360 = 100; double ro1 = 0.71;double ro2 = 0.76;double ro3 = 0.82;double ro4 = 0.9;double ro5 = 0.99;double roF = 1.1; double TF = 200;double TR1 = 200;double TR2 = 200;double TLI23 = 320;double TLI35 = 220;double TLI47 = 120; double F = (double)RawTB; double eTB = 0.02; double Part1TB = 18; double Part2TB = 14; double Part3TB = 14; double Part4TB = 7; double Part5TB = 45; double ZakTB = 20; double Zok1TB = 30; double Zok2TB = 40; double LTTB = 130; double LC1TB = 200; double LC2TB = 400; double LC3TB = 150; double TLTTB = 3; double TLC1TB = 2; double TLC2TB = 4; double TLC3TB = 3; double p1TB = 95; double p56TB = 10; double pTR1TB = 100; double pTR2TB = 100; double pTLTTB = 1; double t1TB = 350; double p360TB = 100; double t23TB = 320; double t35TB = 240; double t47TB = 140; double t56TB = 100; double e1 = eTB; //if (e1 <= 0 || F <= 0 || e1 > 1 ) //{ throw new Exception("Параметры должны быть положительными, а доля отгона не может быть больше 1"); } double part1 = Part1TB / 100; double part2 = Part2TB / 100; double part3 = Part3TB / 100; double part4 = Part4TB / 100; double part5 = Part5TB / 100; //if (part1 <= 0 || part2 <= 0 || part3 <= 0 || part4 <= 0 || part5 <= 0 || part1 + part2 + part3 + part4+part5 > 1) //{ throw new Exception("Доли фракций должны быть положительными, а их сумма не должна превышать 100%"); } double ro12 = (part1 + part2) / ((part1 / ro1) + (part2 / ro2)); //средняя плотность бензин+дизель double ro123 = (part1 + part2 + part3) / ((part1 / ro1) + (part2 / ro2) + (part3 / ro3)); //бензин+нафта+дизель double ro1234 = (part1 + part2 + part3 + part4) / ((part1 / ro1) + (part2 / ro2) + (part3 / ro3) + (part4 / ro4)); //бензин+нафта+дизель+газойль double Zak = ZakTB; double Zok1 = Zok1TB; double Zok2 = Zok2TB; // if (Zak <= 0 || Zok1 <= 0 || Zok2 <= 0) // { throw new Exception("Количество пара должно быть положительным"); } double LT = LTTB; double LC1 = LC1TB; double LC2 = LC2TB; double LC3 = LC3TB; //if (LT <= 0 || LC1 <= 0 || LC2 <= 0 || LC3 <= 0) //{ throw new Exception("Количество орошения должно быть положительным"); } double TLT = TLTTB; double TLC1 = TLC1TB; double TLC2 = TLC2TB; double TLC3 = TLC3TB; p1 = p1TB; p56 = p56TB; pTR1 = pTR1TB; pTR2 = pTR2TB; pTLT = pTLTTB; p360 = p360TB; // if (p1 <= 0 || p56 <= 0 || pTR1 <= 0 || pTR2 <= 0 || pTLT <= 0 || p360 <= 0) // { throw new Exception("Давление должно быть положительным"); } t1 = t1TB; t23 = t23TB; t35 = t35TB; t47 = t47TB; t56 = t56TB; double hz360 = Enth.Enthalpy_OHSteam(360, p360); //расчет энтальпии ПВП double p22 = Enth.Pressure(p1, p56, 22); double hzT22 = Enth.Enthalpy_OHSteam(Temp[22], p22); double p34 = Enth.Pressure(p1, p56, 34); double hzT34 = Enth.Enthalpy_OHSteam(Temp[34], p34); double p46 = Enth.Pressure(p1, p56, 46); double hzT46 = Enth.Enthalpy_OHSteam(Temp[46], p46); double hzT56 = Enth.Enthalpy_OHSteam(Temp[56], p56); double hzTLT = Enth.Enthalpy_OHSteam(TLT, pTLT); double hzTR1 = Enth.Enthalpy_OHSteam(TR1, pTR1); double hzTR2 = Enth.Enthalpy_OHSteam(TR2, pTR2); double LSfull = e1 * F * (Enth.Enthalpy_Vapors(TF, roF)) + (1 - e1) * F * (Enth.Enthalpy_Liquids(TF, roF)) + (Zak + Zok2 + Zok1) * hz360 - (Zak + Zok2 + Zok1) * hzT56; // hvTf пары нефти на входе hlTf нефть на входе double LS1 = part1 * LSfull / (Enth.Enthalpy_Liquids(Temp[56], ro1)); // hlT56 бензин на выходе LS1TB = LS1; double LS2 = part2 * LSfull / (Enth.Enthalpy_Liquids(Temp[47], ro2)); // hlT47 нафта на выходе LS2TB = LS2; double LS3 = part3 * LSfull / (Enth.Enthalpy_Liquids(Temp[35], ro3)); // hlT35 дизель на выходе LS3TB = LS3; double LS4 = part4 * LSfull / (Enth.Enthalpy_Liquids(Temp[23], ro4)); // hlT23 газойль на выходе LS4TB = LS4; double LS5 = part5 * LSfull / (Enth.Enthalpy_Liquids(Temp[1], ro5)); // hlT1 мазут на выходе LS5TB = LS5; double Z22 = Zak / 38, //пар по тарелкам Z34 = Z22, Z46 = (Zak + Zok2) / 12, Z56 = (Zok1 + Zak + Zok2) / 6; double Qlt = LT * ((Enth.Enthalpy_Vapors(Temp[56], ro1)) - (Enth.Enthalpy_Liquids(TLT, ro1))); //тепло верхнего орошения // hvT23 hlTLT double Qlc3 = LC3 * (Enth.Enthalpy_Liquids(Temp[47], ro12)-Enth.Enthalpy_Liquids(TLC3,ro12)); //тепло промежуточных орошений // hlT47 hlTLC3 double Qlc2 = LC2 * (Enth.Enthalpy_Liquids(Temp[35], ro123) - Enth.Enthalpy_Liquids(TLC2, ro123)); // hlT35 hlTLC2 double Qlc1 = LC1 * (Enth.Enthalpy_Liquids(Temp[23], ro1234) - Enth.Enthalpy_Liquids(TLC1, ro1234)); // hlT23 hlTLC1 double LI47 = (LS1 * (Enth.Enthalpy_Vapors(Temp[56], ro1) - Enth.Enthalpy_Vapors(Temp[46], ro1)) + LS2 * (Enth.Enthalpy_Liquids(Temp[47], ro2) - Enth.Enthalpy_Vapors(Temp[46], ro2)) + Z56 * hzT56 - Z46 * hzT46 - Zok2 * hz360 + Qlt+Qlc3) / (Enth.Enthalpy_Vapors(Temp[46], ro12) - Enth.Enthalpy_Liquids(Temp[47], ro12)); // hvT56 hvT46 hlT47 hvT46 hvT46 hlT47 LI47TB = LI47; double LI35 = (LS1 * (Enth.Enthalpy_Vapors(Temp[56], ro1) - Enth.Enthalpy_Vapors(Temp[34], ro1)) + LS2 * (Enth.Enthalpy_Liquids(Temp[47], ro2) - Enth.Enthalpy_Vapors(Temp[34], ro2)) + LS3 * (Enth.Enthalpy_Liquids(Temp[35], ro3) - Enth.Enthalpy_Vapors(Temp[34], ro3)) + Z56 * hzT56 - Z34 * hzT34 - Zok2 * hz360 - Zok1 * hz360 + Qlt + Qlc2 + Qlc3) / (Enth.Enthalpy_Vapors(Temp[34], ro123) - Enth.Enthalpy_Liquids(Temp[35], ro123)); // hvT56 hvT34 hlT47 hvT34 hlT35 hvT34 hvT34 hlT47 LI35TB = LI35; double LI23 = (LS1 * (Enth.Enthalpy_Vapors(Temp[56], ro1) - Enth.Enthalpy_Vapors(Temp[22], ro1)) + LS2 * (Enth.Enthalpy_Liquids(Temp[47], ro2) - Enth.Enthalpy_Vapors(Temp[22], ro2)) + LS3 * (Enth.Enthalpy_Liquids(Temp[35], ro3) - Enth.Enthalpy_Vapors(Temp[22], ro3)) + LS4 * (Enth.Enthalpy_Liquids(Temp[23], ro4) - Enth.Enthalpy_Vapors(Temp[22], ro4)) + Z56 * hzT56 - Z22 * hzT22 + Qlt + Qlc1 + Qlc2 + Qlc3) / (Enth.Enthalpy_Vapors(Temp[22], ro1234) - Enth.Enthalpy_Liquids(Temp[23], ro1234)); LI23TB = LI23; // hvT56 hvT22 hlT47 hvT22 hlT35 hvT22 hlT23 hvT22 hvT22 hlT23 //количество углеводородных паров из отпарных колонн double LR2 = (LS2 * ((Enth.Enthalpy_Liquids(Temp[47], ro2)) - (Enth.Enthalpy_Liquids(TLI47, ro12))) - Zok2 * (hz360 - hzTR1)) / ((Enth.Enthalpy_Liquids(TLI47, ro12)) - (Enth.Enthalpy_Vapors(TR1, ro1))); // hlT47 hlTLI47 Т флегмы T на выходе отп hlTLI13 hvTR1 double LR1 = (LS2 * ((Enth.Enthalpy_Liquids(Temp[35], ro3)) - (Enth.Enthalpy_Liquids(TLI35, ro123))) - Zok1 * (hz360 - hzTR2)) / ((Enth.Enthalpy_Liquids(TLI35, ro123)) - (Enth.Enthalpy_Vapors(TR2, ro12))); // hlT35 hlTLI35 Т флегмы T на выходе отп hlTLI13 hvTR2 LR1TB = LR1; LR2TB = LR2; } } public class Enth { public static double Enthalpy_Vapors(double t, double p420) { double h = 0; double p1515 = p420 + (0.035 / p420); //Перевод относительной плотности нефтепродукта из Советской системы 4-20 в европейскую 15-15 h = ((50.2 + 0.109 * t + 0.00014 * t * t) * (4 - p1515) - 73.8) * 4.184; // Расчет энтальпии по формуле Войнова с переводом из ккал в кдж return h; } public static double Enthalpy_Liquids(double t, double p420) { double h = 0; double p1515 = p420 + (0.035 / p420); //Перевод относительной плотности нефтепродукта из Советской системы 4-20 в европейскую 15-15 h = (1 / (Math.Sqrt(p1515)) * (0.403 * t + 0.000405 * t * t)) * 4.184; // Расчет энтальпии по формуле c переводом из ккал в кдж return h; } public static double Pressure(double p1, double p56, double i) //профиль давления { double pi; pi = p1 - i * (p1 - p56) / 56; return pi; } public static double Enthalpy_OHSteam(double t, double p) { double h = 0, TK=0; TK = (1/(0.00397-0.0004455*(Math.Log10(p))))-273; //расчет температуры конденсации пара h = (606.5+0.305*TK+0.5*(t-TK)) * 4.184; // Расчет энтальпии по формуле c переводом из ккал в кдж return h; } }