assistance with 'subscript out of range' error

[satan165]

   Hello.  I am historian/archeologist in Chicago studying an abandoned coke plant (coke is a processed version of coal used in iron making). 

I have been exploring the laboratory and came across the hard copy of what I recognized as a BASIC program used to evaluate coal. 

I am far from a computer programmer but I figured out enough to input the whole thing into Notepad, then I got it to run on  QB45 on Dosbox.  I used that to debug my many typos.  Then I installed QB64 which is obviously much easier to work with.

The first issue I had was that the program was written to output it's results to a printer.  I don't have a printer and even if I did, that isn't efficient.  So after some experimentation I was able to make the changes so that it would output to a file instead.  That seems to be working well now.  The program does start up and it is starting to output to the file before it crashes.

I am getting another error and this one is a bit over my head.  I am getting a 'subscript out of range' error on the following lines:

   

Code: QB64: [Select]

1. 210 E(1, 5) = E(1, 5) + A(N, 1) * B(N, 5) / 100
2.    220 E(1, 1) = E(1, 1) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 1) / 100
3.    230 E(1, 2) = E(1, 2) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 2) / 100
4.    240 E(1, 3) = E(1, 3) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 3) / 100
5.    250 E(1, 4) = E(1, 4) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 4) / 100
6.    260 E(1, 21) = E(1, 21) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 21) / 100
7.    270 E(1, 22) = E(1, 22) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 22) / 100

I suspect this is related to the way the 'E' variable is defined in the DIM statement on line 10:

Code: QB64: [Select]

1. 10 Dim A(40, 40), E(20, 50), B(20, 50), Z(20, 20), C(20, 20), A$(50), B$(50)

I saw another post where someone got this error and the 'OPTION _EXPLICIT" command was recommended.  I tried that, but then other variables start to have issues.

I have attached the entire program and I plan on scanning the hard copy.  I don't see any typo here so hopefully it isn't something that dumb.

I am way outside my depth and quietly proud that I made it this far.  If anyone thought they'd like to help that would be much appreciated.  Thank you very much.

[Dimster]

Are you able to provide the Looping statement ?

[satan165]

sorry i dont know what that means, but i did attach the entire program to my initial post

[bplus]

Can you provide test data with known correct solution?

[satan165]

I’m not sure what you are looking for, I must not understand because if I had the “correct solution” I would not need any help.

I guess what I am asking is for someone to review the code, if there is a glaring error and can point me in that direction it would be helpful.

[bplus]

I’m not sure what you are looking for, I must not understand because if I had the “correct solution” I would not need any help.

I guess what I am asking is for someone to review the code, if there is a glaring error and can point me in that direction it would be helpful.

I am looking to not have to answer a million questions about coal which I know nothing about, just to check out some code for you.

[bplus]

Code: QB64: [Select]

1. 109 Print "ENTER THE PREDICTED COKE STABILITY"
2. 110 Input S
3. 118 Print #1, "                         COAL YIELD EVALUATION                      "
4. 119 Print #1,
5. 123 Print #1,
6. 124 Print #1, "                           COAL ANALYSIS                            "
7. 127 Print #1,
8. 129 Print #1, " ***************************DRY BASIS****************************** "
9. 130 Print #1, "COAL  NAME   WET COAL   H20    F.C.    ASH     VOL. MAT    S     O2 "
10. 140 Print #1, "IN BLEND      IN MIX     %      %       %          %       %     %  "
11. 190 For N = 1 To S
12.    200 A(10, 1) = A(10, 1) + A(N, 1)
13.    210 E(1, 5) = E(1, 5) + A(N, 1) * B(N, 5) / 100
14.    220 E(1, 1) = E(1, 1) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 1) / 100
15.    230 E(1, 2) = E(1, 2) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 2) / 100
16.    240 E(1, 3) = E(1, 3) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 3) / 100
17.    250 E(1, 4) = E(1, 4) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 4) / 100
18.    260 E(1, 21) = E(1, 21) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 21) / 100
19.    270 E(1, 22) = E(1, 22) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 22) / 100
20. 275 Next N
21.  

If you input S too high or low, easy to get out of bounds error, what is safe S?

[Dimster]

The Loops in the program are in a For:Next. I was just wondering, where the coding is calling for an input (ie Coal Blends input = C then the Loop Statement is For n= 1 to c, the c being the number just inputted - or The Coal Stability = S and the loop statement for this input is For N= 1 to S) I was thinking the out of range error could be associated with the value being inputted for S.

[satan165]

I am looking to not have to answer a million questions about coal which I know nothing about, just to check out some code for you.

I apologize - I thought you were referring to the code itself. 

I was remiss not to have offered this data initially....in my mind it didn't matter, it seemed that the program would crash regardless of what I entered.

For the first question the program asks for the number of coals in the blend.  Answer '2'.

The following question asks for nine values regarding the coal quality.  Here are some suitable values (these are based on real world values, so I assume the code will accept these):
5,50,10,50,15,16,3,15,12000

the unit of measurement for 'S' is not given by the program, but typically this is expressed as a % (coal stability factor) so 50 should work

If you input S too high or low, easy to get out of bounds error, what is safe S?

what do you see as a safe value based on the code?

I was thinking the out of range error could be associated with the value being inputted for S

i tried a number of different values but the program crashed regardless.  i did not consider that my input is what was crashing it, i assumed it was a problem with the code but you may be right.  What range do you see is acceptable based on the code?

I am not at home right now to continue to test but I will later today for sure.  if you are able to get the program to run beyond line 44 let me know what values you used.

Thank you both so much!

[Dimster]

Satan, I am, like bplus, out of my depth when comes to my knowledge in coal. Hope you can help with one observation on the algorythm. If you look at lines 45 and 50.
Line 45 :The question is asking for 10 data items (1-Cost per Ton, 2-%wet Petcoke, 3-%coal H2O ... 10-BTU/LB)
Line 50 : The arrays being filled - 1 Cost/ton at C(n1), 2 - %Wet Petcoke at A(N,1), 3- %coal H2O at (B(n,5) but as you go along there are only 9 array items. I'm not sure where the BTU's per pound is being used. Likely there but it's not jumping out for me.

The looping control for this data input is generated by the answer to ENTER THE NUMBER of COALS IN THE BLEND? The the qualifying remarks indicate the program can deal with 7 coals plus petcoke. That would seem to imply a total of 8? But 8 can't be the upper limit of your B array which is taking the bulk of the inputted data - the B array is dimensioned (20,50). If answer to "the number of coals in the blend" can exceed 20 you should generate the range error limit.

The Coal Stability is similar. The A array is (40,40). an S greater than 40, (which might be the % Wet Petcoke in Blend) could generate the error
                                        - the E array is (20,50) but seems this array is not affected by the value of S,
                                         - then the B array, which we just addressed is (20,50) so S greater than 20 should generate the range error.

I have no idea if this is helpful or not as I'm a Pepsi drinker.

[TempodiBasic]

Hi evil165
ops satan165

I think that the error you got depends from your input.
Why? Because this old fashion BASIC program doesn't check the input given from user.
Changing the input I have got different errors at runtime at different lines of code.
Cause it is a simpler huge math task, you can port those calculations into a spreedsheet like excel or similar.
But if you like play with OLD BASIC the first suggestion is
give a readable name at those orrible variable with one character as name!  :-)
2 write a check input data entry routine with manage the 10-20 variables taken by a single INPUT command.
3 kill all those GOTOs
Here a screenshot of that I have got typing value for all  variables asked
  [ You are not allowed to view this attachment ]  
good luck for your archeologist study

[satan165]

Line 45 :The question is asking for 10 data items (1-Cost per Ton, 2-%wet Petcoke, 3-%coal H2O ... 10-BTU/LB)
Line 50 : The arrays being filled - 1 Cost/ton at C(n1), 2 - %Wet Petcoke at A(N,1), 3- %coal H2O at (B(n,5) but as you go along there are only 9 array items. I'm not sure where the BTU's per pound is being used. Likely there but it's not jumping out for me.

There are nine.  I added numbers to each to illustrate

(1) COST/TON, (2) %WET COAL IN BLEND, (3) %COAL H20, AND (DRY BASIS) (4) %FIXED CARBON, (5) %ASH, (6) %VOLATILE MATTER, (7) %SULFUR, (8) %OXYGEN, (9) BTU/LB

The "And (Dry Basis)" is distinguishing the remaining 5 variables as being measured while dry.  Superfluous comment for sure.

If answer to "the number of coals in the blend" can exceed 20 you should generate the range error limit.

In reality, I think it is doubtful that anyone would have a coal blend of more than 4-5 types.

You may be right that there is a small window of acceptable answers for 'S' I perhaps I just need to continue to experiment with that

Because this old fashion BASIC program doesn't check the input given from user.
Changing the input I have got different errors at runtime at different lines of code.
Cause it is a simpler huge math task, you can port those calculations into a spreedsheet like excel or similar.

I agree completely.  As far as recreating the program in Excel, I have done this once before already.  I found some documentation that described a 'program' but I did not have the code, but it was simpler so I just built it in Google Sheets https://docs.google.com/spreadsheets/d/12MwtYzWKHrRA_lldHyX30KwqV4_MCTDckFx7eApC8wg/edit?usp=sharing

I think I may be out of my depth in trying to institute the changes you described but perhaps I will try.  If I can get this program to output finished data one time, I will be satisfied.

I actually found a second version of this program (another hard copy), but that one is a bit better because rather than asking you to input the nine variables, it simply lets you type the name of the coal, the program already has the nine variables for a variety of coals so no need to input manually.  However, that hard copy was not printed correctly, so some of the longer lines are cut off.  I may try to input it anyways and figure out what the rest of the missing line is.

Thank you all for your patience and help, if anyone can output a file and get the program to run all the way through, or can just help in specifying with the lower/upper limits for the input on these variables is please share.

Thank you all so much.

[bplus]

Thanks for S.

Code: QB64: [Select]

1. 190 For N = 1 To S
2.     '' dimmed   A(40, 40), E(20, 50), B(20, 50), Z(20, 20), C(20, 20), A$(50), B$(50)
3.  
4.     ' where you see N that is where S is going to go
5.     ' A(N, )  A has top limit of 40 so S at 50 is going out of bounds
6.     ' B(N, )  B has top limit of 20 so B is going to beat A for going out of bounds
7.  
8.  
9.     200 A(10, 1) = A(10, 1) + A(N, 1)
10.    210 E(1, 5) = E(1, 5) + A(N, 1) * B(N, 5) / 100
11.    220 E(1, 1) = E(1, 1) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 1) / 100
12.    230 E(1, 2) = E(1, 2) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 2) / 100
13.    240 E(1, 3) = E(1, 3) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 3) / 100
14.    250 E(1, 4) = E(1, 4) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 4) / 100
15.    260 E(1, 21) = E(1, 21) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 21) / 100
16.    270 E(1, 22) = E(1, 22) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 22) / 100
17. 275 Next N
18.  

The use of S as percentage looks fishy to me, unless S never is expected to exceed 20 (20%).

EDIT Oh we cross posted sorry.

[satan165]

Just got home and immediately went back to just inputting data.  Honestly, sorry for wasting everyone's time on checking code when I just didn't spend enough time trying different input values.

I got it to output on the first try

number of coals = 2
coal quality data = 5,50,10,50,15,16,3,15,12000
petcoke = no
predicted coke stability = 4

so as I just said an hour ago, if I can get the program to calculate and output the data, I am satisfied.  now I can have some fun fixing the formatting and messing around and trying to understand the upper/lower limits for these values

I think I will work on the second set of code I found, I may return with some questions on the missing bits there

you guys are awesome - thanks for humoring an idiot like me.  If you are interested in the iron making industry, or urban exploration, or chicago history you can check out my site www.acmecoke.com . ill be putting up a post about this program in the next week or so and ill be sure to give the forum a shout out.

thanks again, so much! made my day. 

[OldMoses]

The program is almost entirely linear in execution, with just a couple of GOTO jumps. which I circumvented with IF...END IF blocks and then striped out the old line numbers and REMs. I did that to get rid of the clutter of old style programming which makes it hard for me to see issues sometimes. It's hard to believe I used to bash my head against such things back in my Applesoft daze.

I can't see much that jumps out at me as a problem, and like you I got it to give output entering your numbers, but one line looks like a possible typo. I left the old line number in to help locate it:

Code: QB64: [Select]

1. 107 A$(I) = A$(N)  
2.  

I say it looks like a typo because I don't see any declaration or assignment for 'I' prior to that line, and all the rest of the string array is referenced by literals. Maybe I just don't understand what the original programmer was trying to accomplish. In fact, I'm sure I don't. ;)

I include my changes to the original code below, it shows the 'I' in lower case, which was what caught my eye. I also added some comments that highlight those lines that are to be changed if altering the 1986 standards. Happy bug hunting...

Code: QB64: [Select]

1. OPEN "OUTPUT", 1, "#1"
2. PRINT "                                         COAL EVALUATION PROGRAM"
3. PRINT
4. PRINT "NOTE: PROGRAM USES NOVEMBER 1986 STANDARDS FOR BY-PRODUCT VALUES, FIXED AND VARIABLE CONVERSION COSTS AND BULK DENSITY OIL USAGE."
5. PRINT "UPDATING REQUIRES CHANGES IN LINES 1410, 1420, 1430, 1440, 4150, 4160, 4170, 4410, 4690 AND 4720."
6. DIM A(40, 40), E(20, 50), B(20, 50), Z(20, 20), C(20, 20), A$(50), B$(50)
7. PRINT
8. PRINT "THIS PROGRAM IS DESIGNED FOR A BLEND CONTAINING UP TO SEVEN COALS PLUS PETCOKE.  ENTER THE NUMBER OF COALS IN THE BLEND."
9. INPUT C
10. FOR N = 1 TO C
11.     PRINT
12.     PRINT "ENTER NAME OF COAL"
13.     INPUT A$(N)
14.     PRINT
15.     PRINT "ENTER DATA IN THIS ORDER: COST/TON, %WET COAL IN BLEND, %COAL H20, AND (DRY BASIS) %FIXED CARBON, %ASH, %VOLATILE MATTER, %SULFUR, %OXYGEN, BTU/LB                               "
16.     INPUT C(N, 1), A(N, 1), B(N, 5), B(N, 1), B(N, 2), B(N, 3), B(N, 4), B(N, 21), B(N, 22)
17. NEXT N
18. PRINT
19. N = C
20. PRINT "DOES THE BLEND CONTAIN PETCOKE? ENTER YES OR NO "
21. INPUT D$: D$ = UCASE$(MID$(D$, 1, 1))
22. A$(8) = "PETCOKE"
23. IF D$ <> "N" THEN
24.     PRINT
25.     PRINT "ENTER DATA IN THIS ORDER: COST/TON, %WET PETCOKE IN BLEND, %PETCOKE H20, AND (DRY BASIS) %FIXED CARBON, %ASH, %VOLATILE MATTER, %SULFUR, %OXYGEN, BTU/LB"
26.     INPUT C(8, 1), A(8, 1), B(8, 5), B(8, 1), B(8, 2), B(8, 3), B(8, 4), B(8, 21), B(8, 22)
27. END IF
28. A$(9) = "COAL PORTION"
29. A$(10) = "BLEND"
30. A$(12) = "FCE COKE"
31. A$(i) = A$(N)
32. PRINT
33. PRINT "ENTER THE PREDICTED COKE STABILITY"
34. INPUT S
35. PRINT #1, "                         COAL YIELD EVALUATION                      "
36. PRINT #1,
37. PRINT #1,
38. PRINT #1, "                           COAL ANALYSIS                            "
39. PRINT #1,
40. PRINT #1, " ***************************DRY BASIS****************************** "
41. PRINT #1, "COAL  NAME   WET COAL   H20    F.C.    ASH     VOL. MAT    S     O2 "
42. PRINT #1, "IN BLEND      IN MIX     %      %       %          %       %     %  "
43. FOR N = 1 TO S
44.     A(10, 1) = A(10, 1) + A(N, 1)
45.     E(1, 5) = E(1, 5) + A(N, 1) * B(N, 5) / 100
46.     E(1, 1) = E(1, 1) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 1) / 100
47.     E(1, 2) = E(1, 2) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 2) / 100
48.     E(1, 3) = E(1, 3) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 3) / 100
49.     E(1, 4) = E(1, 4) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 4) / 100
50.     E(1, 21) = E(1, 21) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 21) / 100
51.     E(1, 22) = E(1, 22) + A(N, 1) * ((100 - B(N, 5)) / 100) * B(N, 22) / 100
52. NEXT N
53. FOR N = 1 TO C
54.     P$ = "\                             \###,##   ##,###   ##,###   ##,###        ##,###   ##,###   ##,###   "
55.     PRINT #1, USING P$; A$(N); A(N, 1); B(N, 5); B(N, 1); B(N, 2); B(N, 3); B(N, 4); B(N, 21)
56. NEXT N
57. IF D$ <> "N" THEN
58.     PRINT #1, USING P$; A$(8); A(8, 1); B(8, 5); B(8, 1); B(8, 2); B(8, 3); B(8, 4); B(8, 21)
59. END IF
60. N = C
61. PRINT #1,
62. PRINT #1,
63. PRINT #1, "                                 BLEND ANALYSIS"
64. PRINT #1,
65. ' COAL BLEND LESS PETCOKE ANALYSIS
66. A(9, 1) = A(10, 1) - A(8, 1)
67. E(11, 5) = (E(1, 5) - A(8, 1) * B(8, 5) / 100) * 100 / A(9, 1)
68. E(11, 1) = (E(1, 1) - A(8, 1) * B(8, 1) / 100) * 100 / A(9, 1)
69. E(11, 2) = (E(1, 2) - A(8, 1) * B(8, 2) / 100) * 100 / A(9, 1)
70. E(11, 3) = (E(1, 3) - A(8, 1) * B(8, 3) / 100) * 100 / A(9, 1)
71. E(11, 4) = (E(1, 4) - A(8, 1) * B(8, 4) / 100) * 100 / A(9, 1)
72. E(11, 21) = (E(1, 21) - A(8, 1) * B(8, 21) / 100) * 100 / A(9, 1)
73. E(11, 22) = (E(1, 22) - A(8, 1) * B(8, 22) / 100) * 100 / A(9, 1)
74. E(8, 1) = B(8, 1) * (100 - B(8, 5)) / 100
75. E(8, 2) = B(8, 2) * (100 - B(8, 5)) / 100
76. E(8, 3) = B(8, 3) * (100 - B(8, 5)) / 100
77. E(8, 4) = B(8, 4) * (100 - B(8, 5)) / 100
78. E(8, 21) = B(8, 21) * (100 - B(8, 5)) / 100
79. PRINT #1, "                                 **************************WET BASIS****************               "
80. PRINT #1, "                 WET COAL       H20      F.C.      ASH          VOL.MAT.              S        O2  "
81. PRINT #1, "                 IN MIX%         %        %         %               %                 %        %   "
82. Y = E(11, 4)
83. X = E(11, 21)
84. PRINT #1, USING P$; A$(9); A(9, 1); E(11, 5); E(11, 1); E(11, 2); E(11, 3); Y; X
85. Y = E(8, 4)
86. X = E(11, 21)
87. PRINT #1, USING P$; A$(8); A(8, 1); E(8, 5); E(8, 1); E(8, 2); E(8, 3); Y; X
88. Y = E(1, 4)
89. X = E(1, 21)
90. PRINT #1, USING P$; A$(10); A(10, 1); E(10, 5); E(10, 1); E(10, 2); E(10, 3); Y; X
91. PRINT #1,
92. ' PRODUCED COKE (COKE PUSHED FROM THE OVENS)=COAL+PETCOKE
93. Z(4, 1) = E(11, 1) + E(11, 2) + 9.000001E-02 * E(11, 3) + 1.18 * (100 - E(11, 5)) / 100
94. Z(4, 1) = (Z(4, 1) / 100) * 2000 * (A(9, 1) / 100)
95. Z(4, 2) = (B(8, 1) + B(8, 2)) * ((100 - B(8, 5)) / 100)
96. Z(4, 2) = (Z(4, 2) / 100) * 2000 * (A(8, 1) / 100)
97. Z(4, 3) = Z(4, 1) + Z(4, 2)
98. Z(5, 1) = E(11, 1) + E(11, 2)
99. Z(6, 1) = E(11, 3) - E(11, 21)
100. ' T IS THE COKING TEMPERATURE IN KELVIN
101. T = 1323
102. Z(6, 1) = (Z(6, 1) ^ (1.26 * (T / 1000) ^ 1.4)) / (9.5 * (T / 1000) ^ 7.3)
103. Z(5, 1) = Z(5, 1) + Z(6, 1)
104. Z(5, 1) = (Z(5, 1) / 100) * 2000 * (A(9, 1) / 100)
105. Z(5, 2) = (B(8, 1) + B(8, 2)) * ((100 - B(8, 5)) / 100)
106. Z(5, 2) = (Z(5, 2) / 100) * 2000 * (A(8, 1) / 100)
107. Z(5, 3) = Z(5, 1) + Z(5, 2)
108. Z(6, 1) = (1 / 2) * (Z(4, 1) + Z(5, 1))
109. Z(6, 2) = (1 / 2) * (Z(4, 2) + Z(5, 2))
110. Z(6, 3) = (1 / 2) * (Z(4, 3) + Z(5, 3))
111. ' TAR YIELD
112. Z(1, 1) = 8.399999 * (E(11, 3) - E(11, 21)) ^ 2 / T
113. Z(1, 1) = (Z(1, 1) / 100) * (A(9, 1) / 100) * 2000
114. Z(1, 2) = .95 * B(8, 3) * ((100 - B(8, 5)) / 100)
115. Z(1, 2) = (Z(1, 2) / 100) * (A(8, 1) / 100) * 2000
116. Z(1, 3) = Z(1, 1) + Z(1, 2)
117. ' LIGHT OIL YIELD
118. Z(2, 1) = -.284 + .0384 * E(11, 3)
119. Z(2, 1) = (Z(2, 1) / 100) * (A(9, 1) / 100) * 2000
120. Z(2, 2) = 0
121. Z(2, 3) = Z(2, 1) + Z(2, 2)
122. ' COKE OVEN GAS YIELD
123. Z(3, 1) = 14500 * E(11, 1) / 100
124. Z(3, 1) = .105 * (E(11, 22) - Z(3, 1))
125. Z(3, 1) = Z(3, 1) * (T / E(11, 3)) ^ (1 / 2)
126. Z(3, 1) = Z(3, 1) * 2000 / (500 * 35.6)
127. Z(3, 1) = E(3, 1) * (A(9, 1) / 100)
128. Z(3, 2) = .05 * B(8, 3) * (100 - B(8, 5)) / 100
129. Z(3, 2) = (Z(3, 2) / 100) * 2000
130. Z(3, 2) = Z(3, 2) * (A(8, 1) / 100)
131. Z(3, 3) = Z(3, 1) + Z(3, 2)
132. ' DETERMINE COKE FC, ASH, VM AND S ANALYSIS ON A DRY WT BASIS
133. A(12, 1) = (Z(6, 3) / 2000) * 100
134. E(12, 2) = (E(1, 2) / A(12, 1)) * 100
135. E(12, 3) = .9
136. E(12, 4) = .6 * (A(9, 1) / 100) * E(11, 4)
137. E(12, 4) = E(12, 4) + .8 * (A(8, 1) / 100) * E(8, 4)
138. E(12, 4) = (E(12, 4) / A(12, 1)) * 100
139. E(12, 1) = 100 - E(12, 2) - E(12, 3)
140. X = E(12, 4)
141. Q$ = "THE OXYGEN CONTENT OF THE WET COAL OPRTION OF THE BLEND IS BY WT ##,###%"
142. PRINT #1, USING Q$; E(11, 21)
143. R$ = "THE HEATING VALUE OF THE WET COAL PORTION OF THE BLEND IN BTU/LB IS #####"
144. PRINT #1, USING R$; E(11, 22)
145. ' COKE OVEN GAS FOR UNDERFIRING REQUIREMENT
146. ' 1,022 BTU/LB OF DRY COAL; 1,710 BTU/LB OF WATER
147. ' C.O.G. IS AT 35.6 CU FT/LB AND 500 BTU/LB OF WATER
148. Z(3, 4) = 1022 * 2000 * ((100 - E(1, 5)) / 100)
149. Z(3, 4) = E(3, 4) + 1710 * 2000 * (E(1, 5) / 100)
150. Z(3, 4) = Z(3, 4) / (500 * 35.6)
151. 'PRODUCED COKE TO FURNACE COKE, BREEZE AND BUCKWHEAT
152. 'L=L(2,1); M=E(2,2); S=E(1,17); N=L(1,1)
153. 'BREEZE; PRODUCED COKE ASH; STABILITY; FURNACE COKE
154. M = E(12, 2)
155. L = .79 + (1.824 * M - .152 * S)
156. Z(7, 3) = Z(6, 3) * (L / 100)
157. Z(7, 5) = Z(7, 3)
158. Z(8, 3) = Z(6, 3) - Z(7, 3)
159. FOR i = 1 TO 8
160.     Z(i, 5) = Z(i, 3) - Z(i, 4)
161. NEXT i
162. 'PRODUCTS PER NT OF DRY FURNACE COKE
163. FOR i = 1 TO 8
164.     Z(i, 6) = Z(i, 5) * (2000 / Z(8, 5))
165. NEXT i
166. 'UNIT VALUES
167. Z(1, 7) = .38 / 9.68 '                                              1410
168. Z(2, 7) = .5 / 7.41 '                                               1420
169. Z(3, 7) = (2.94 / 1000000!) * 500 * 35.6 '                          1430
170. Z(7, 7) = 39.7 / (2000 * (100 - 12) / 100) '                        1440
171. 'DETERMINE VALUE PER NT OF FURNACE COKE
172. FOR i = 1 TO 8
173.     Z(i, 8) = Z(i, 6) * Z(i, 7)
174. NEXT i
175. PRINT #1,
176. PRINT #1,
177. PRINT #1,
178. PRINT #1, "                                               PRODUCTS OF COKING "
179. PRINT #1,
180. B$(1) = "COAL"
181. B$(2) = "PETCOKE"
182. B$(3) = "BLEND"
183. B$(4) = "USED"
184. B$(5) = "PRODUCT"
185. B$(6) = "/NT COKE"
186. B$(7) = "$/LB"
187. B$(8) = "$/NT COKE"
188. PRINT #1, "ITEM                  TAR       L.OIL     COKE OVEN      COKE     COKE     FURNACE"
189. PRINT #1, "LB/NT                                        GAS         TOTAL    BREEZE     COKE "
190. FOR i = 1 TO 6
191.     S$ = "\                     \    ###.####  ###.####  ###.####  ###.####  ###.####  ###.####"
192.     PRINT #1, USING S$; B$(i); Z(1, i); Z(2, i); Z(3, i); Z(6, i); Z(7, i); Z(8, 1)
193. NEXT i
194. PRINT #1,
195. PRINT #1,
196. PRINT #1, "   COKE PRODUCED FROM COAL ONLY    TOTAL COKE     BREEZE &     FURNACE"
197. PRINT #1, "   FORMULA 1   FORMULA 2 AVERAGE    PRODUCED      BUCKWHEAT      COKE "
198. T$ = "    ####.##   ####.##   ####.##   ####.##   ####.##   ####.##"
199. PRINT #1, USING T$; Z(4, 1); Z(5, 1); Z(6, 1); Z(6, 3); Z(7, 3); Z(8, 3)
200. 'COAL COST
201. C(9, 2) = 2000 / Z(8, 5)
202. FOR i = 1 TO 8
203.     C(i, 2) = C(9, 2) * A(i, 1) / 100
204.     C(i, 3) = C(i, 2) * (100 - B(i, 5)) / 100
205.     C(i, 4) = C(i, 1) * C(i, 2)
206.     C(9, 3) = C(9, 3) + C(i, 3)
207.     C(9, 4) = C(9, 4) + C(i, 4)
208. NEXT i
209. C(9, 1) = C(9, 4) / C(9, 2)
210. PRINT #1,
211. PRINT #1,
212. PRINT #1,
213. PRINT #1, "  COAL NAME       UNIT WET     WET COAL NT/        DRY COAL NT/       COAL COST   "
214. PRINT #1, "   IN BLEND       COAL COST      NT DRY              NT DRY          $/NT OF DRY  "
215. PRINT #1, "                    $/NT         FURNACE COKE      FURNACE COKE     FURNACE COKE  "
216. FOR i = 1 TO N
217.     U$ = "\                \  ##.##            ##.##             ##.##               ###.###"
218.     PRINT #1, USING U$; A$(i); C(i, 1); C(i, 2); C(i, 3); C(i, 4)
219. NEXT i
220. IF D$ <> "N" THEN
221.     PRINT #1, USING U$; A$(8); C(8, 1); C(8, 2); C(8, 3); C(8, 4)
222. END IF
223. PRINT #1,
224. PRINT #1, USING U$; A$(10); C(9, 1); C(9, 2); C(9, 3); C(9, 4)
225. PRINT #1,
226. PRINT #1,
227. 'WET COAL BLEND, DRY COAL BLEND AND COKE ANALYSIS
228. 'WET COAL BLEND ANALYSIS
229. A = A(10, 1)
230. B = E(1, 5)
231. C = E(1, 1)
232. D = E(1, 2)
233. E = E(1, 3)
234. F = E(1, 4)
235. A$ = "WET BLEND"
236. PRINT #1, USING P$; A$; A; B; C; D; E; F; E(1, 21)
237. 'DRY COAL BLEND ANALYSIS
238. A = A(10, 1) - E(1, 5)
239. C = E(1, 1) * 100 / A
240. D = E(1, 2) * 100 / A
241. E = E(1, 3) * 100 / A
242. F = E(1, 4) * 100 / A
243. G = E(1, 21) * 100 / A
244. H = E(1, 22) * 100 / A
245. A$ = "DRY BLEND"
246. E(1, 5) = 0!
247. PRINT #1, USING P$; A$; A(10, 1); E(1, 5); C; D; E; F; G
248. PRINT #1,
249. PRINT #1,
250. V$ = "THE HEATING VALUE OF THE WET BLEND IN BTU/LB IS #####"
251. PRINT #1, USING V$; E(1, 22)
252. W$ = "THE HEATING VALUE OF THE DRY BLEND IN BTU/LB IS #####"
253. PRINT #1, USING W$; H
254. PRINT #1,
255. PRINT #1,
256. A$(12) = "FURNACE COKE"
257. A$(13) = "TOTAL COKE"
258. PRINT #1, "           DRY COKE    H20     F.C.    ASH     VOL. MAT.    S "
259. PRINT #1, '           YIELD%       %       %       %           %       % "
260. PRINT #1, USING P$; A$(13); A(12, 1); E(12, 5); E(12, 1); E(12, 2); E(12, 3); X
261. Y = (Z(8, 3) / 2000) * 100
262. PRINT #1, USING P$; A$(12); Y; E(12, 5); E(12, 1); E(12, 2); E(12, 3); X
263. PRINT #1,
264. PRINT #1,
265. PRINT #1, " DESCRIPTION                                   $NT DRY          $NT DRY   "
266. PRINT #1, "                                             FURNACE COKE    FURNACE COKE "
267. PRINT #1, " RAW MATERIAL COST"
268. PRINT #1, "            COAL  "
269. X$ = "                                                    ###.###"
270. PRINT #1, USING X$; C(9, 4)
271. PRINT #1, "         BULK DENSITY OIL"
272. PRINT #1, "         1.53 GAL/NT   FCE COKE" '                           4150
273. PRINT #1, "         AT $0.50/GALLON       " '                           4160
274. PRINT #1, "                                               0.770 " '     4170
275. PRINT #1, "             INVENTORY COST AT 1%"
276. C(11, 4) = Z(3, 8) * Z(3, 4) / Z(3, 5)
277. PRINT #1, USING X$; C(10, 4)
278. PRINT #1, "       COG FOR UNDERFIRING"
279. C(11, 4) = Z(3, 8) * Z(3, 4) / Z(3, 5)
280. PRINT #1, USING X$; C(11, 4)
281. PRINT #1, "           SUBTOTAL RAW MATERIALS"
282. C(12, 4) = C(9, 4) + .77 + C(10, 4) + C(11, 4) '                        4410
283. Y$ = "                                                       ###.###            ###.###"
284. PRINT #1, USING Y$; C(12, 4); C(12, 4)
285. PRINT #1,
286. PRINT #1, " BY-PRODUCT CREDITS"
287. PRINT
288. PRINT #1, "     EXCESS COKE OVEN GAS"
289. PRINT #1, USING X$; Z(3, 8) * (-1)
290. PRINT #1, "       COG FOR UNDERFIRING"
291. A = Z(3, 8) * Z(3, 4) / Z(3, 5) * (-1)
292. PRINT #1, USING X$; A
293. PRINT #1, "      TOTAL COKE OVEN GAS"
294. T = Z(3, 8) * Z(3, 3) / Z(3, 5) * (-1)
295. PRINT #1, USING X$; T
296. PRINT #1,
297. PRINT #1, "       TAR"
298. PRINT #1, USING X$; Z(1, 8) * (-1)
299. PRINT #1, "       LIGHT OIL"
300. PRINT #1, USING X$; Z(2, 8) * (-1)
301. PRINT #1, "              COKE BREEZE AND BUCKWHEAT"
302. PRINT #1, USING X$; Z(7, 8) * (-1)
303. B = T - Z(1, 8) - Z(2, 8) - Z(7, 8)
304. PRINT #1, "               TOTAL BY-PRODUCTS CREDITS"
305. PRINT #1, USING Y$; B; B
306. PRINT #1,
307. PRINT #1, " CONVERSION COST  "
308. PRINT #1,
309. PRINT #1, "              CONVERSION - VARIABLE"
310. C = 16.23 '                                                         4690
311. PRINT #1, USING X$; C
312. PRINT #1, "             CONVERSION - FIXED"
313. D = 193073! / (7 * 100 * 18 * (Z(8, 5) / 2000)) '                   4720
314. PRINT #1, USING X$; D
315. PRINT #1, "             TOTAL CONVERSION"
316. E = C + D
317. PRINT #1, USING Y$; E; E
318. PRINT #1,
319. PRINT #1, "  TOTAL PRODUCTION COST"
320. PRINT #1, " $/NT DRY FURNACE COKE"
321. G = C(12, 4) + B + E
322. PRINT #1, USING Y$; G; G
323. PRINT #1,
324. PRINT #1,
325. PRINT #1,
326. PRINT #1,
327. PRINT #1,
328. PRINT #1, "END OF RUN"
329. CLOSE COM1
330. END
331.