#
# 
#
# lines with a hash-sign are ignored
#
# input for 1d reactive 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]
# reation time
reaction_time = [600.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
c1 = [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]
c2 = [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_substance = [1, 1]
source_labels = ['c1_factory1','c1_factory2']
# generated as 
source_values['c1_factory1'] = [120.02, 119.00, 110.58, 119.29, 119.17, 123.05, 117.80, 113.07, 115.45, 116.28, 129.71, 127.54, 143.43, 153.95, 139.29, 126.78, 122.50, 121.10, 121.13, 127.15, 121.79, 116.52, 115.10, 115.02, 117.11, 124.80, 129.21, 116.36, 121.46, 127.30, 131.77, 141.05, 143.05, 137.78, 136.39, 118.37, 121.80, 122.61, 118.27, 113.54, 105.80, 105.67, 108.64, 106.42, 103.01, 103.01, 102.01, 95.36, 105.12, 108.54, 117.99, 114.40, 106.38, 112.44, 112.84, 107.47, 107.01, 92.66, 100.95, 94.14, 89.57, 99.76, 93.49, 90.92, 87.72, 99.97, 100.87, 93.34, 85.62, 98.75, 113.11, 125.05, 115.86, 114.36, 113.00, 115.33, 111.24, 104.29, 101.30, 109.43, 126.77, 128.31, 126.22, 131.19, 127.45, 140.80, 148.46, 139.08, 133.90, 143.37, 145.78, 148.18, 135.33, 136.71, 141.39, 136.41, 128.86, 127.35, 119.58, 119.05, 121.08, 130.99, 132.09, 127.95, 139.64, 145.28, 146.51, 150.93, 146.72, 152.46, 153.86, 157.81, 149.58, 160.20, 162.66, 148.06, 150.98, 159.68, 163.64, 172.90, 181.57, 173.95, 156.34, 168.68, 170.69, 169.99, 166.20, 154.31, 155.33, 144.77, 138.27, 136.44, 141.63, 143.98, 133.82, 144.31, 144.13, 157.26, 147.87, 141.73, 135.79, 134.74, 127.41, 120.29, 103.03, 97.23, 99.45, 98.45, 97.57, 97.08, 98.66, 95.37, 90.65, 94.80, 92.16, 94.26, 87.72, 87.18, 93.10, 87.08, 88.61, 87.41, 87.38, 85.31, 70.82, 70.58, 81.47, 86.81, 84.61, 80.39, 82.83, 80.65, 77.46, 78.55, 88.30, 70.56, 75.59, 69.01, 76.68, 59.82, 59.18, 60.59, 58.48, 58.01, 48.06, 41.08, 36.46, 40.48, 37.54, 50.02, 52.18, 39.79, 37.81, 34.43, 36.58, 36.34, 35.69, 31.75, 32.73, 23.70, 29.71, 21.60, 18.42, 19.02, 22.47, 23.51, 37.68, 50.77, 56.44, 58.14, 53.30, 51.04, 41.98, 40.54, 38.43, 52.45, 59.00, 60.47, 59.31, 60.90, 67.41, 84.47, 87.40, 94.18, 99.51, 103.85, 111.36, 120.61, 120.90, 118.59, 120.24, 127.44, 136.11, 131.90, 138.29, 136.73, 134.00, 137.68, 144.12, 139.57, 128.38, 127.86, 131.72, 132.12, 143.33, 140.53, 145.93, 139.80, 130.31, 134.96, 124.61, 115.09, 110.78, 100.13, 104.49, 102.72, 93.62, 87.61, 80.98, 81.09, 64.23, 60.11, 51.01, 54.52, 59.92, 61.96, 62.17, 44.03, 45.47, 40.34, 39.72, 38.66, 31.60, 37.25, 40.33, 39.09, 46.88, 59.50, 52.46, 53.37, 46.41, 53.05, 57.76, 55.03, 52.04, 47.69, 55.57, 54.82, 54.80, 46.75, 50.21, 70.83, 73.72, 64.17, 66.86, 57.33, 59.42, 48.15, 54.53, 63.86, 76.01]
source_values['c1_factory2'] = [105.08, 99.42, 102.56, 109.03, 101.02, 102.38, 107.79, 98.25, 91.29, 96.99, 96.91, 100.68, 110.47, 113.71, 107.70, 112.94, 111.43, 105.25, 120.25, 120.41, 113.20, 117.51, 129.32, 133.05, 138.15, 138.20, 136.34, 141.96, 139.78, 139.58, 135.03, 129.41, 126.90, 117.76, 119.19, 119.53, 118.83, 102.16, 98.34, 100.09, 102.68, 103.86, 103.44, 91.69, 94.20, 93.27, 96.81, 101.12, 100.40, 102.66, 100.11, 97.69, 94.88, 87.98, 89.79, 98.82, 94.57, 90.94, 96.25, 96.51, 107.65, 95.12, 87.63, 82.23, 98.20, 101.28, 97.96, 98.13, 95.79, 102.86, 102.86, 108.66, 113.52, 111.54, 110.45, 111.12, 105.06, 103.17, 102.39, 94.85, 101.65, 93.34, 88.22, 79.87, 87.83, 83.03, 70.73, 77.02, 87.76, 92.75, 93.79, 87.59, 82.26, 91.38, 99.77, 95.29, 101.13, 102.60, 99.77, 93.67, 95.77, 82.58, 71.40, 86.38, 86.01, 105.20, 103.27, 104.39, 110.92, 109.36, 105.35, 115.55, 108.91, 114.53, 110.51, 115.64, 121.44, 114.99, 118.52, 125.62, 117.49, 119.74, 119.07, 109.78, 107.54, 106.18, 98.14, 83.94, 77.54, 77.79, 65.59, 75.22, 71.20, 56.65, 58.09, 69.47, 70.28, 71.86, 72.05, 73.69, 80.43, 81.99, 87.68, 91.48, 92.49, 94.57, 87.46, 80.74, 80.89, 83.69, 77.39, 83.57, 84.55, 104.51, 107.33, 97.66, 100.40, 103.93, 100.08, 101.70, 97.26, 97.38, 100.37, 108.91, 114.27, 94.80, 96.87, 102.83, 102.64, 108.73, 114.12, 123.26, 131.71, 138.05, 125.36, 117.68, 122.27, 114.76, 110.04, 114.77, 114.59, 121.97, 111.81, 115.01, 102.84, 92.27, 94.44, 93.43, 85.27, 83.32, 79.29, 73.06, 75.22, 64.76, 67.54, 53.32, 57.84, 60.95, 67.00, 62.45, 67.92, 60.60, 67.69, 52.76, 55.52, 50.87, 45.99, 47.91, 46.25, 54.08, 53.98, 62.00, 59.11, 60.20, 63.66, 72.27, 77.14, 75.35, 65.74, 61.26, 64.22, 57.50, 60.23, 68.83, 71.04, 72.63, 75.94, 85.06, 95.72, 88.28, 89.88, 94.23, 92.33, 110.45, 116.40, 109.88, 115.52, 120.25, 125.34, 112.43, 123.28, 111.15, 118.92, 110.49, 113.02, 119.72, 125.24, 125.01, 117.81, 116.62, 110.45, 98.13, 91.55, 95.87, 99.87, 96.79, 97.21, 100.48, 101.83, 103.59, 118.52, 101.64, 103.98, 105.93, 110.55, 106.09, 112.34, 113.30, 119.14, 125.72, 134.79, 133.67, 142.57, 143.43, 132.28, 130.58, 126.34, 113.30, 109.97, 105.27, 110.32, 112.53, 118.63, 120.30, 118.83, 117.79, 119.51, 127.54, 132.84, 122.40, 122.82, 120.75, 124.47, 121.81, 129.58, 127.71, 116.55, 126.30, 128.97, 132.51, 126.86]
#output (index based and 0 based)
output_file = 'reactive_pollution_model.output'
matlab_output_file = 'reactive_pollution_model_output.m'
output_map_times = [0, 60, 180]
output_locations = [10, 20, 40, 10, 20, 40]
output_substance = [1, 1, 1, 2, 2, 2]
output_labels = ['c1_locA','c1_locB','c1_locC','c2_locA','c2_locB','c2_locC']
# boundaries
# only left and right at locations 0 and -1 allowed at the moment
bound_labels=['c1_left', 'c1_right','c2_left','c2_right']
bound_locations=[0, -1, 0, -1]
bound_substance=[1, 1, 2, 2]
bound_values['c1_left']=[0.01]
bound_values['c1_right']=[0.0]
bound_values['c2_left']=[0.01]
bound_values['c2_right']=[0.0]