5 REMARKABLE SMOG SIMULATION PROGRAM BY DAVE AHL
10 & "SMOG CITY SIMULATION PROGRAM":&:T=7
20 & "THIS PROGRAM REPRESENTS A MODEL OF A CITY WITH CERTAIN"
30 & "CONDITIONS WHICH YOU MAY SPECIFY.  THE PROGRAM WILL THEN"
40 & "PRINT OUT THE RESULTING SMOG PATTERNS BY HOUR OF THE DAY.":&
42 & "A PLOT FOR 5 DAYS IS PRINTED.":&
45 &:& "CARS PER MILE OF ROAD"
50 & "EXAMPLES:  LOS ANGELES   200 CARS PER MILE"
55 & "           DETROIT       100"
60 & "           TULSA          25"
65 INPUT "HOW MANY CARS PER MILE OF ROAD IN SMOG CITY";C
70 IF C>275 THEN &"TOKYO IS THE MOST DENSE CITY IN THE WORLD WITH 275/MILE.";
	:& "  AGAIN...":GOTO 65
80 &:& "SMOG GENERATION CONSTANT"
85 & "EXAMPLES:  1950 AUTO   .005"
90 & "           1972 AUTO   .001"
95 & "           BUS         .010"
100 INPUT "WHAT IS THE AVERAGE VALUE IN SMOG CITY";L1
105 IF L1>.01 THEN &"EVEN GARY, IND. HAS AN OVERALL RATE OF ONLY .008.";
	:& "  AGAIN...":GOTO 100
110 &:& "DAYTIME SMOG BREAKDOWN RATE (% PER HOUR)"
115 & "EXAMPLES:   CLEAR   0.02"
120 & "            CLOUDY  0.01"
125 INPUT "DAYTIME RATE IN SMOG CITY";R1
130 IF R1<.008 THEN &"EVEN IN HEAVY CLOUDS THE RATE IS 0.008 OR GREATER.";
	:& "  AGAIN...":GOTO 125
140 &:& "NIGHTTIME SMOG BREAKDOWN RATE (% PER HOUR)"
145 & "EXAMPLES:   HOT   0.05"
150 & "            COOL  0.10"
155 INPUT "NIGHTTIME RATE IN SMOG CITY";R2
160 IF R2<.04 THEN &"ARE YOU ON EARTH?  TRY AGAIN WITH 0.04 OR GREATER."
	:GOTO 155
170 &:& "DISPERSION FROM WEATHER (% PER HOUR)"
175 & "EXAMPLES:   HIGH WIND AND RAIN   1.00"
180 & "            NO WIND AND DRY      0.01"
185 INPUT "DISPERSION FROM WEATHER IN SMOG CITY";R3
190 IF R3<.008 THEN &"EVEN IN A THERMAL INVERSION THE RATE IS 0.008 OR MORE";
	:& ".  AGAIN...":GOTO 185
195 Z=Z+L1*C*10-R1*Z-R3*Z FOR I=1 TO 13:Z=30/Z
200 &:&:&"HOUR","SMOG LEVEL","PLOT OF SMOG LEVEL":&
220 R=R1:K1=L1:T1=INT((T-6)/12)
240 IF T1/2=INT(T1/2) THEN 270
250 K1=0:R=R2
270 S=S+K1*C*10-R*S-R3*S
275 IF S>0 THEN 280 ELSE S=0
280 T=T+1:T3=INT(T/12):T2=T-T3*12+1
300 IF T3/2=INT(T3/2) THEN A$=" AM" ELSE A$=" PM"
310 S1=28+S*Z:IF S1>72 THEN S1=72
320 & T2;A$,S,TAB(S1);"X":H=H+1:IF H<121 THEN 220
330 INPUT "ANOTHER RUN (Y OR N)";A$:IF A$="Y" THEN &:&:GOTO 45
999 END