# # scenario with true sources for generating observations # # lines with a hash-sign are ignored # # input for 1d pollution model # # grid x = [0.0, 60.0, 3600.0] # stationary flow u = [1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0] # cross sectional area a = [1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0] # initial concentrations c = [0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0] # simulation timespan refdate = '01 dec 1999' #unit is always seconds unit = 'seconds' # step is 1min total 5hours, while a particle travels through domain in 1hour time = [ 0,60,18000] # sources mass/m^3/s source_locations = [5, 30] source_labels = ['factory1','factory2'] # generated as source_values['factory1'] = [23.05, 26.02, 22.55, 15.11, 27.15, 39.25, 32.86, 49.47, 30.62, 72.28, 71.86, 97.23,120.07, 88.43,137.25,174.95,124.66,135.25, 80.32, 69.87, 36.49, 28.24, 25.99, 41.12, 38.30, 27.36, 27.49, 17.04, 13.21, 25.93, 27.59, 19.87, 19.27, 14.78, 11.20, 10.44, 8.86, 9.63, 19.08, 17.08, 14.97, 21.91, 31.19, 40.20, 40.90, 25.58, 34.27, 30.70, 33.04, 27.23, 33.10, 29.68, 23.52, 13.16, 8.63, 18.50, 45.16, 29.44, 28.26, 29.20, 24.46, 53.03, 66.71, 51.12, 30.29, 31.69, 32.73, 18.54, 24.90, 25.72, 45.01, 40.80, 50.59, 38.43, 38.70, 64.26,114.89,117.14,332.78,310.26,513.25,545.96,455.61,944.82,1011.43,857.65,756.17,826.17,784.33,767.53,681.89,436.05,414.60,209.87,191.29,193.76,217.37,310.28,184.18,108.38,247.26,741.34,1035.45,774.87,944.15,971.81,389.59,386.42,353.96,412.86,413.06, 92.33, 56.86, 30.68, 22.16, 10.32, 9.46, 6.59, 8.69, 6.66, 8.77, 5.04, 2.35, 1.63, 2.22, 3.76, 3.03, 6.07, 5.95, 6.47, 5.91, 4.84, 5.60, 7.35, 6.94, 2.33, 6.00, 4.15, 9.62, 15.06, 5.02, 5.19, 5.96, 8.56, 11.25, 14.63, 16.80, 20.23, 41.15, 51.01, 47.23, 47.32, 61.45, 60.12, 42.94, 27.80, 31.70, 47.86, 82.32,118.54, 85.89,119.45, 77.25, 72.59, 40.91, 71.27, 85.25,124.02,115.53,128.94, 80.71, 46.33, 30.08, 64.15, 74.21, 64.37, 51.81, 33.18, 57.29, 85.76,102.35,121.44,204.87,320.11,569.32,359.30,265.71,425.75,646.24,974.35,978.64,1033.42,855.04,800.68,1352.03,1451.31,1514.85,1691.59,795.66,1775.87,1335.32,3070.43,2829.30,3445.95,2612.57,4783.44,4634.51,4321.40,5068.21,7447.13,4452.18,2345.82,2959.58,3116.83,1400.84,1970.76,1439.79,1662.80,2546.66,3517.33,2482.50,2286.72,3146.44,4902.50,5800.36,10148.97,13676.49,15490.18,10417.54,10755.03,13450.50,18993.69,17098.02,14845.72,21480.42,33390.04,30626.23,20567.50,18762.49,18495.96,21765.84,4988.93,8563.92,3747.14,2417.25,1172.48,796.50,435.65,495.70,867.07,568.66,410.43,236.53,307.04,180.56,169.25,115.01, 96.22,107.92,141.14, 78.89, 49.65, 35.80, 26.13, 33.33, 26.01, 18.22, 6.81, 5.96, 7.39, 11.00, 15.12, 11.71, 8.55, 15.27, 5.20, 8.89, 9.65, 6.58, 8.47, 6.21, 10.48, 22.36, 45.84, 42.98, 29.86, 20.97, 13.68, 10.64, 4.33, 5.33, 7.08, 6.70, 4.84, 5.06, 4.08, 3.71, 4.91, 6.57, 3.95, 6.11] source_values['factory2'] = [5.97, 4.83, 6.15, 11.59, 9.59, 7.82, 13.18, 13.39, 18.25, 30.27, 33.84, 52.12, 69.97, 39.19, 31.07, 38.54, 16.71, 12.10, 4.37, 3.37, 3.28, 3.80, 3.03, 0.86, 1.72, 1.50, 2.81, 1.64, 2.46, 3.17, 2.67, 7.53, 17.61, 24.40, 24.42, 30.23, 39.58, 30.34, 26.79, 49.74, 73.92, 81.38, 56.07, 44.52, 62.54, 70.71, 66.39, 88.18, 40.41, 65.66, 62.23, 51.41, 64.69, 82.25,121.84,173.07,276.31,236.37,167.48,125.14,107.88,174.62, 87.56, 72.09, 78.47,107.33, 68.79, 63.97, 65.62, 65.47, 74.31, 69.02, 50.14, 41.66, 41.08, 71.21, 82.25, 81.97, 89.74,105.68,214.32,227.73,407.12,438.22,577.07,689.41,1353.06,846.11,538.00,590.81,809.25,804.75,689.81,1199.72,857.73,1119.85,2321.20,2509.06,1603.51,4075.19,2606.20,2679.43,3448.63,5238.50,4612.56,5469.32,2361.28,2217.71,1481.42,675.82,1001.40,1712.13,1740.11,1013.87,2578.24,1994.57,2682.50,2317.24,3673.02,3821.06,2847.19,2684.28,2348.73,3885.12,4742.11,5006.64,3523.41,2247.86,2708.35,4569.63,2953.73,1583.67,2994.57,2851.28,1863.97,1017.70,1583.91,2762.20,5294.55,6226.99,5812.42,12803.77,19316.79,8238.24,7348.79,18488.97,12386.38,12304.80,9392.35,14573.67,9867.78,9061.01,5973.77,5001.28,7113.00,6774.72,6059.87,3696.04,1814.92,1941.97,2178.32,2285.18,957.54,1182.55,1207.84,1875.91,1817.94,1802.86,1289.17,774.25,638.83,1097.45,2158.83,1135.95,1237.92,1135.44,1102.57,1435.52,870.63,533.42,563.57,596.56,504.34,432.65,592.93,425.21,371.71,473.62,505.24,905.68,827.57,761.95,606.49,451.08,626.49,338.96,282.78, 90.77,121.50,150.31,213.52,132.82,124.05,123.46, 56.85, 26.36, 19.75, 34.40, 82.36, 85.50, 40.64, 44.69, 51.53, 48.27, 50.75, 23.03, 12.17, 17.12, 15.43, 7.78, 6.86, 4.97, 6.72, 6.67, 7.70, 9.58, 6.26, 6.60, 6.34, 8.54, 18.83, 21.40, 20.37, 20.45, 15.69, 7.59, 9.33, 7.77, 9.20, 9.44, 7.74, 8.63, 11.70, 12.47, 13.18, 8.65, 14.99, 24.61, 8.30, 3.50, 4.19, 2.98, 1.82, 5.66, 3.93, 2.40, 2.80, 2.11, 1.35, 0.78, 0.63, 0.83, 0.46, 0.47, 0.49, 0.45, 0.72, 0.98, 3.47, 2.10, 2.31, 3.30, 2.65, 3.38, 6.93, 9.75, 14.14, 12.84, 14.52, 11.64, 10.91, 10.49, 15.93, 18.99, 18.65, 15.79, 14.21, 8.71, 12.01, 19.50, 12.62, 25.63, 20.59, 21.93, 25.16, 22.42, 20.87, 12.70, 12.57, 16.25, 35.19] #output (index based and 0 based) output_file = 'pollution_model_generate_observations.output' output_locations = [10, 20, 40] output_labels = ['locA','locB','locC'] # boundaries # only left and right at locations 0 and -1 allowed at the moment bound_labels=['left', 'right'] bound_locations=[0, -1] bound_values['left']=[0.01] bound_values['right']=[0.0]