~o|16 ~?!i&32768 ~o=295 ~$>end_of_copyright ~/!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ~/Copyright (C) Heinz Spiess, CH-2558 Aegerten, 1994. All rights reserved. ~/ ~/The right to use this macro is granted to all EMME/2 users, provided the ~/following conditions are met: ~/ 1) The macro cannot be sold for a fee (but it can be used and distributed ~/ without charge within consulting projects). ~/ 2) The user is aware that this macro is not a part of the EMME/2 software ~/ licence and there is no explicit or implied warranty or support ~/ provided with this macro. ~/ 3) The comments in this macros must not be removed and any additions or ~/ modification must be appropriately identified as such and give at least ~/ date, name and the reason of the modification. ~/!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ~:end_of_copyright ~/============================================================================= ~/MLUS92 1.3 - (C) Heinz Spiess, EMME/2 Support Center, CH-2558 Aegerten, 1995 ~/ ~/Berechnung der Luftverunreinigungen an Strassen ohne oder mit lockerer ~/Randbebauung gemaess MLuS-92 ******** ALPHA TEST VERSION !!! ******** ~/ ~/Aufruf: ~/ ~ [ ...] ~/ ~/wobei: ~/ zu berechnender Luftschadstoff, moegliche Werte: ~/ co Kohlenmonoxid Resultate in @co ~/ hc Kohlenwasserstoffe @hc ~/ nox Stickoxide @nox ~/ pb Blei @pb ~/ so2 Schwefeldioxid @so2 ~/ pm Russpartikel @pm ~/ c6h6 Benzol @c6h6 ~/ no2 Stickstoffdioxid @no2 ~/ alle alle oben beschriebenen Schadstoffe ~x=%0% ~?x>0 / check if enough parameters are provided ~$>start_macro ~t9=~*Mehr Informationen? ~t9=%t9.1% ~?t9=n ~$>end_of_macro ~/Die Berechnungen entsprechen der Publikation: ~/MLuS-92: Merkblatt ueber Luftverunreinigungen an Strassen, ~/Strassen ohne oder mit lockerer Randbebauung Ausgabe 1992, herausgegeben ~/von der Forschungsgesellschaft fuer Strassen- und Verkehrswesen ~/ ~/Achtung: Folgende Einschraenkungen sind zu beachten (siehe 3.1.1) ~/ - vPkw > 50 km/h ~/ - vLkw <= 80 km/h ~/ - pLkw <= 0.5 ~/ - Laengsneigung der Strasse n <= 6% ~/ ~/Die Extra-Streckenattribute (@co, @hc, @nox, ...) muesssen vor Aufruf ~/des Makros mittels Modul 2.42 erstellt werden. Skalare ms97, ms98 ~/und ms99 werden zur Speicherung von Temporaerdaten verwendet. ~$>end_of_macro >>> Inhalt der Register fuer Berechnungen: >>> t1: DTV (Netzrechnerausdruck, in Kfz/24h) >>> t2: pLkw (Netzrechnerausdruck, nicht in %) >>> t3: vPkw (Netzrechnerausdruck, in km/h) >>> t4: vLkw (Netzrechnerausdruck, in km/h) >>> t5: Abstand zu Fahrbahnrand (Netzrechnerausdruck, in m) >>> t6: mittlere Windgeschwindigkeit in m/s >>> t8: Schadstoff: co, hc, nox, pb, so2, pm >>> x Bits 0-3: Bezugsjahr (0=1985, 1=1990, 2=1995, 3=2000) >>> 4: Windverhaeltnisse (0=Fall I, 1=Fall II) >>> 5: 0=Jahresmittelwert, 1=98%Perzentile ~:start_macro ~/ ~+|~t1=volau|~t9=~*Eingabe: DTV in Kfz/24h (%%%t1%%%)=|~?!t9=|~t1=%%%t9%%% ~+|~t2=volad/(volau+volad)|~t9=~*Eingabe: LkW-Anteil (%%%t2%%%)=|~?!t9=|~t2=%%%t9%%% ~+|~t3=60*timau/len|~t9=~*Eingabe: PkW-Geschwindigkeit in km/h (%%%t3%%%)=|~?!t9=|~t3=%%%t9%%% ~+|~t4=60*timau/len.min.80|~t9=~*Eingabe: LkW-Geschwindigkeit in km/h (%%%t4%%%)=|~?!t9=|~t4=%%%t9%%% ~+|~t6=3|~t9=~*Eingabe: Mittlere Windgeschwindigkeit in m/s (%%%t6%%%)=|~?!t9=|~t6=%%%t9%%% ~+|~t5=50|~t9=~*Eingabe: Fahrbahnabstand in m (%%%t5%%%)=|~?!t9=|~t5=%%%t9%%% ~:Eingabe Bezugsjahr ~+|~z=1995|~t9=~*Eingabe: Bezugsjahr 1985, 1990, 1995 oder 2000 (%%%z%%%)=|~?!t9=|~z=%%%t9%%% ~+|~x=-1|~?z=1985|~x=0|~?z=1990|~x=1|~?z=1995|~x=2|~?z=2000|~x=3|~?x<0|~$Eingabe Bezugsjahr ~+|~z=1|~t9=~*Eingabe: 1=Jahresmittelwert, 2=98%-Perzentile (%%%z%%%)=|~?t9=2|~x+16 ~+|~z=1|~t9=~*Eingabe: 1=Fall I, 2=Fall II (%%%z%%%)=|~?t9=2|~x+32 ~t8=%1% ~?t8=alle ~t0=co hc nox pb so2 pm c6h6 no2 2.41 1 / Dummyausdruck fuer Streckenauswahl ~/ ~/Streckenauswahl fuer Schadstoffberechnungen ~+|y|tmpj3|len||~z=%o%|~o=6|4 ~+|~t8=~*|~?t8=|~o=%z%|~?t8=*|~o=%z%|~?t8=all|~o=%z%|%%%t8%%%|~?e|~$|~?q=0|~$|r|~/|~o=%z% ~/DTV: %t1% ~/pLkw: %t2% ~/vPkw: %t3% ~/vLkw: %t4% ~/W'geschw. u: %t6% ~/Abstand s: %t5% ~+|~z=%x%|~z&7|~z*5|~z+1985 ~/Bezugsjahr: %z% ~?!x&16 ~/Konzentration: Jahresmittelwert ~?x&16 ~/Konzentration: 98% Prozentile ~?!x&32 ~/Windverh.: Fall I ~?x&32 ~/Windverh.: Fall II ~/ ~:next_pollutant ~?t0= ~$>all_done ~t8=%1% ~/... Berechnung der Konzentration des Schadstoffs %t8% in @%t8% 1 y @%t8% ~?e ~+|~/Fehler: Extra-Attribut @%t8% nicht gefunden!|~%||~$next_pollutant y %t8% nach MLuS92(f%x%/s%t5%/u%t6%) ~?t8=no2 ~$>Spezialfall NO2 0*( ~:Konstanten al, Tabelle 2 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ~?t8=co put(150.500)+put(2.047)+put(252.876)+put(1.303)+put(-.0266)+put(14.287) ~?t8=hc put(38.046)+put(1.100)+put(210.508)+put(-.982)+put(0.0086)+put(1.360) ~?t8=nox put(24.940)+put(2.409)+put(67.824)+put(15.961)+put(0.0286)+put(6.248) ~?t8=pb put(0.107)+put(2.000)+put(0)+put(0)+put(0)+put(0.00935) ~?t8=so2 put(2.124)+put(1.117)+put(10.582)+put(1.775)+put(-.0041)+put(0.270) ~?t8=pm put(0)+put(0)+put(6.194)+put(0.173)+put(0)+put(0.0338) ~?t8=c6h6 put(38.046)+put(1.100)+put(210.508)+put(-.982)+put(0.0086)+put(1.360) ~:Bezugskonzentrationen Tabelle 1 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ~?t8=co +put(0.58)+put(1.80)+put(2.49)+put(9.66) ~?t8=hc +put(0.06)+put(0.18)+put(0.23)+put(0.77) ~?t8=nox +put(0.28)+put(0.87)+put(0.71)+put(2.15) ~?t8=pb +put(.00063)+put(.0019)+put(.0018)+put(.0048) ~?t8=so2 +put(0.011)+put(0.034)+put(0.049)+put(0.184) ~?t8=pm +put(0.0014)+put(.0043)+put(.0059)+put(.0224) ~?t8=c6h6 +put(0.06)+put(0.18)+put(0.23)+put(0.77) ~:Reduktionsfaktoren, Tabelle 3 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ~?t8=co +put(1.00)+put(0.56)+put(0.34)+put(0.24) ~?t8=hc +put(1.00)+put(0.71)+put(0.53)+put(0.43) ~?t8=nox +put(1.00)+put(0.82)+put(0.70)+put(0.62) ~?t8=pb +put(1.00)+put(0.56)+put(0.21)+put(0.00) ~?t8=so2 +put(1.00)+put(1.00)+put(1.00)+put(1.00) ~?t8=pm +put(1.00)+put(1.00)+put(1.00)+put(1.00) ~?t8=c6h6 +put(1.00)+put(0.71)+put(0.53)+put(0.43) )+ ~: Ki(s)*g(s)*( (1) get(%x%/16+7)*(1-0.166*ln(1+(%t5%)))*( ~: (ePkw (5,7,7a) (get(1)*exp(get(2)*put(%t3%)/100)/get(15) ~: *(1-p)*rij (7,9) *(1-put(%t2%))*get(11+(%x%&7)) ~: +eLkW*p)/e' (7b,7) +(get(3)/put(%t4%)+get(4)+get(5)*get(17))*get(16))/get(6) ~: )*fd (5) )*(%t1%)/46000 ~: *fu (1) *(3.6-1.9*((%x%&32)!=0))/(%t6%) ~?!t8=c6h6 ~$>Ende der Formel ~: Zusatz nur fuer Benzolkonzentration *KC6H6(s) (3,4) *((.0616-.000194*get(15)+(get(15)/780.63)^2) *(1-get(16))+0.046*get(16)) ~$>Ende der Formel ~:Spezialfall NO2 (0.052+0.058*((%x%&16)!=0))*(1-0.088*ln(1+(%t5%))) *.00447*((%t1%)^0.514)*exp((%t1%)/-241545) ~:Ende der Formel 5 3 99 mx%t8% Maximalkonzentration %t8% (MLuS92) 98 i-%t8% I-Knoten Max %t8% (MLuS92) 97 j-%t8% J-Knoten Max %t8% (MLuS92) r ~/ Maximalwert %t8%: %ms99% mg/m3 auf Strecke %ms98%-%ms97% ~% ~$next_pollutant ~:all_done q ~/ ~/Makro MLUS92 normal beendet. ~/============================================================================= ~:end_of_macro ~o=6