Ejemplo n.º 1
0
        static void Main(string[] args)
        {
            BusMovement testBM    = new BusMovement("BusMovement.txt");
            HydroZones  testHZ    = new HydroZones("CityZones.txt");
            River       testR     = new River("Riverinput.txt");
            Rainfall    testRf    = new Rainfall("TestRF.out");
            RiverFlow   testRFlow = new RiverFlow("TestRFlow.out");

            // Test output
            Console.Write(" Rainfall is: " + testRf.GetcumRF() + "\n");
            Console.Write(" Riverflow at Ouse at Skelton is: " + testRFlow.GetFlowOuseSkelton() + "\n");
            Console.Write(" Riverflow at Foss at Huntington is: " + testRFlow.GetFlowFossHuntington() + "\n");
            Console.Write(" Flow at the furthest downstream node is: " + testR.calcRiverFlow(testHZ, testBM, testRf, testRFlow) + "\n");
        }
Ejemplo n.º 2
0
        public double calcRiverFlow(HydroZones HZ, BusMovement BM, Rainfall Rf, RiverFlow RFlow)
        {
            double cumArea = 0, cumConc = 0;

            for (int i = 1; i <= numRiverReaches; i++)
            {
                // Get upstream node and check if node exists or not
                iUpstream = riverReach[i].GetUpstreamReach();
                if (iUpstream != 0)
                {
                    flowUpstream = riverReach[iUpstream].GetFlowAtNode();
                }
                else
                {
                    // Check to see if need to add other inflows
                    if (i == 1)
                    {
                        // Add in flow for Ouse at Skelton
                        flowUpstream = RFlow.GetFlowOuseSkelton() * 86400;
                    }
                    else if (i == 6)
                    {
                        // Add in flow for Foss at Huntington
                        flowUpstream = RFlow.GetFlowFossHuntington() * 86400;
                    }
                    else
                    {
                        flowUpstream = 0;
                    }
                }
                // Get tributary node and check if it exists or not
                iUpstream = riverReach[i].GetTributaryReach();
                if (iUpstream != 0)
                {
                    flowTributary = riverReach[iUpstream].GetFlowAtNode();
                }
                else
                {
                    flowTributary = 0;
                }
                // Calc inflow from CityZone associated with river reach (if any)
                numCityZones = riverReach[i].GetNumCZs();
                if (numCityZones > 0)
                {
                    inflowCityZone = 0;
                    for (int iCZs = 1; iCZs <= numCityZones; iCZs++)
                    {
                        int IDCZ = riverReach[i].GetIDCZ(iCZs);
                        inflowCityZone = inflowCityZone + (Rf.GetcumRF() / 1000) * HZ.GetOverallAreaHZ(IDCZ);
                        cumArea        = cumArea + HZ.GetUrbanAreaHZ(IDCZ);
                        cumConc        = cumConc + BM.GetROEmissions(IDCZ) / HZ.GetUrbanAreaHZ(IDCZ);
                    }
                }
                else
                {
                    inflowCityZone = 0;
                }
                // calc flow based on flow from upstream node, tributary, inflow from city zone and any WWTW inflow
                flowAtNode = flowUpstream + flowTributary + inflowCityZone + riverReach[i].GetFlowWWTW();
                riverReach[i].SetFlowAtNode(flowAtNode);
            }
            Console.Write(" Cum. Area is : " + cumArea + "\n");
            Console.Write(" Cum. Conc is : " + cumConc + "\n");
            return(riverReach[numRiverReaches].GetFlowAtNode());
        }