CCC 5

C PBOGBAH SEPHIS 10
C 15
C TBABSIIHT STAGE CALCOLATIOBS TOR IHOREI SOLVEBT EXTRACTIOB 20
C PBOCESS. CALCBLAIIOBS TALIS EOB 3 TBABSFEBBIB6 COBPOBEBTS, 25
C IBCIBDIBG SAITIBG EPPECT PEOB BOB-TBABSPEBBIBG BITRATE SAITS. 30
C COnPOBEBT 1 " DRABIOB COHPOBBBT 3 = BITBIC ACID 35
C COKPOHEBT 2 ' TRORIDB COBPOREBT » ' ALDHIBOn BITBATE DO
C STAGE TOLDBES EQOAL OBLESS BOTED OTHEBBISE 45
CCC 50
C
C PBOGRAB BBITTEB BI S. B. WATSOB AT OAK BIOGE BATIOBAl lABORATORI
C JOBB 1978 TBRSIOB
C <
BEAl*8 TITIE(10,10) 55
DIBEBSIOB AODT(100|, OODT(IOO). STTOL(IOO), DP(1I|, S0I.IB(3| , SOLOB 60
1T(3), AQOOT(IOO), ORGOOT(IOO) 65
DIBEBSIOB AA(IOO). TT(100), 00(100), RATIO (100) 70
COHBOB /PIBS / IRATIO,IO0T,I»Ol,RECT(100) ,ICHBGB,AtTRT(100),OLV»T(1 75
100) 80
COdBOB /STGS/ BSOLO,CODnfl(100),AR(100),OR(100),AT{100),OT(100) 85
COHBOB /PEED/ BTOST.APDBT(100),OPDRT(100) 90
COHBOB /PROP/ 1PBO,I(«,100,2) ,T(«,100,2),IPD(«,100),TPD(I1,100 ) ,TPH 95
10P(100,2),APDTIB(100) ,OPDTEB(100) 100
COBBOB /DEC/ A(100),0(100),RAIOdOO) ,RCTCA(100) ,RCICO(100) 105
CONBOB /OIST/ CTBP.ART(II) ,DTRT(I|) ,T0 110
COHBOB /BOBTE/ BSTR,ISTR,JSTR 115
C 120
BS0I(I<l| 125
C 130
C REIT CASE IS BE! 135
c mo
10 COBTIBBE 1»5
READ SaO, BTTL 150
DO 20 ITTL«1,BTT1 155
READ 550 , (TITLE(1TT1,I),I«1,10) 160
20 COBTIBOE 165
READ 560, CTBP,BTOST,IPRO,IRATIO,ICALC,IO0T,ITOL,IPR,IPBCH,ICALPL, 170
1BSTB 175
C 180
C TITLE IS TBE PROBLEB TITLE. 10 SPACES ARE ALIOBED. 185
C BTOST " TOTAL BOBBER OP STAGES, BDST BOT EXCEED 25 190
C BSOLB • BOBBER OF SOLOTES, BAXIBOB « 195
C CTBP = TOLOBB PBACIIOB OP DBT TBP 200
C IPRO - 0 POR A EERO IBITIAL COBCEBTBATIOB PROPILB 205
C 1 POR A BOB-IERO IBITIAL PEOPILE 210
C IRATIO - 0 IP PRASE RATIO - PLOB RATIO 215
C 1 IP PRASE RATIO CARDS POLLOV 220
C ICALC - 0 FOB SRDTDOfB CALCOLATIORS (BO D OR TH IB PEED STREABS) 225
C 1 POR START-DP OR IBTERBBPTIOB 230
C lOOT • 0 IP TRERE ARE BO EXTRA OBTGOIBG STREABS 235
C 1 IP TRERE ARE OOT GOIBG STREABS IB ADDITION TO EBD STREABS. 2a0
C ITOL - 0 IP STAGE TOLDBES ABE ALL EQOAL 245
C 1 IP STAGE TOUBES ARE DBIOIIAL BUT ACCORDIRG TO TOTAL FLOB 250
C 2 IF STAGE TOLDBES ARE IBDEPEBDEBTLT DBESOAL 255
C 3 IF STAGE TOLOHES ARE EODAL ABD SPECIFIED. (OPTIOBS 2 OR 3 260
C SHODLD BE BSED FOR HAB CALCOLATIORS TO CORRECT FOR STAGE 265
C RESIDEBCE TIBE.) 270
C IPR « 0 TO PBIBT ODT ALL TIBE IBTERTAIS 275
C 1 TO PRIBT OBLI LAST IBTERTAL 280
C -1 TO PRIBT ETERT 20TB TIBE IBTERTAL 285
C ICRHGE » 0 IP BEXT CASE IS BEB 290
C 1 IP BEIT CASE IBTOLTES OBLI A CHARGE IB FEED PIOBS, COHP, AB 295
C 2 STOP 300
C IPBCH = 0 IP BO POHCHED CARD OBTPOT DESIRED 305
C 1 FOR FIBAL PROFILE OOTPBT OB PBHCBED CARDS 310
C ICALPL « 0 BO CALCOBP PLOTS 315
C = 1 COBCEBTBATIOB TS TIBE IBTERTAL 320
C » 2 PERCEHT OF FEED TS TIBE IBTERTAL 325
C - 3 PLOT BOTH 1 AHD 2 330
C BSTR « 0 BO 0BOSOAL ROOTIBG PATTERB 335
C = 1 ROOTIRG PATTERB OTHER THAB BORHAL 340
C 345
C " EEBOIHG OF ARRATS 350
C 355
DO 30 J=1,BTOST 360
COD0H(J)=0.0 365
STTOL(a)'0.0 370
DO 30 K-1,2 375
TPROF(J,!t)>=25.0 380
DO 30 I-1,BS0L0 365
DP (I)-0.0 390
X(I,J,K)>=0.0 395
T(I,J,R)»0.0 400
30 COBTIBOE 405
ICHBGE-0 410
C 415
C ISTR ABO JSTR ASSIGBED TALDES TO ATOID COHPILER PBOBLEHS 420
C 425
ISTR=0 430
JSTR=0 435
C 440
C BEXT CASE IHTOLTES OHLT A CHARGE IK PEED FLOBS, COHP.TEHP 445
C 450
24
40 COBTIBOE 455
DO 50 J"1,BTOST 460
A(J)-0.0 465
O(J)«0.0 470
AFORT(J)-0.0 475
OFDBT(J)'0.0 480
ALTBT(J)-0.0 485
OLTBT(J)-0.0 490
AFDTEB(J)'25.0 495
OFDTER(J)'2S.O 500
50 COBTIBOE 505
DO 60 I<1,BS0LD 510
DO 60 J-1,BTOST 515
XFD(I,J)«0.0 520
TFD(I,J)-0.0 525
60 COBTIBOE 530
C 535
C PBOFIL SETS IBITIAL COHCEBTRATIOBS FOR TIBE ZERO, ALSO DEFIHES FEED S 540
C 545
CALL PROFIL 550
C 555
C BRITE FEED STREAB DATA 560
C 565
PRIBT 570 570
DO 70 ITTL»1,BTTL 575
PRIBT 580, (TITLE(ITTL,I),1-1,10) 580
70 COBTIBOE 585
BRITE (6,590) 590
DO 80 I«1,BTOST 595
IF (AFDRT(I).LE.O.) GO TO SO 600
PRIBT 600, I,XPD(1,I) ,XFD(2,I),XFD(3,I) ,XFD(a,I) ,AFDRT(I) ,AFDTEB(I 60S
1) 610
C 615
C COBTEBT TC BOLARITT 620
C 625
XFD(1,I)-XFD(1,I)/238 . 630
IPO (2,1)-IPO (2,1)/232 . 635
80 COBTIBOE 640
DO 90 1-1,BTOST 645
C 650
C SET AfTID IB ALL STAGES TO SCBOB ACID COBCEBTBATIOB 655
C 660
IF (IPRO.E(1.0.ABO.ICBHGE.Ea.O) X (3, I, 1) - (IPD (3, 1) •AFDRT( 1) tXPO (3,4 665
1) -AFORT (4) ) /(APDRT (1) •AFORT (4) ) 670
C 675
IF (OFDRT(I).LE.O.) GO TO 90 680
PRIBT 610, CTBP,I,TFD(1,I),IFD(2,I),IFD(3,I),TFD(4,I),0F0RT(I),0F0 685
ITEB(I) 690
C 695
C COBTEBT TC HOLARITI 700
C 705
TFD(1,I)-TFD(1,I)/238. 710
TFD(2,I)-TFD(2,I)/232. 715
C 720
C FLOBS SETS ALL IBTERSTAGE FLOB RATES, DEFIRES ODTGOIBG STREABS ABD PL 725
C RATIOS, ABD CALCOLATES BECESSART RECTCLE PLOHS 730
C 735
90 COBTIBOE 740
CALL FLOBS (STTOL) 745
C 750
C PRIBT ODTGOIBG STBEAHS (ALTBT ABD OLTRT) 755
C 760
C IP BO EXTRA ODTGOIBG STBEAHS, ATOID PRIBT OF ALTRT ABD OLTRT 765
C 770
IP (lOOT.EQ.O) 00 TO 120 775
C 780
PRIBT 620 785
DO 100 1-1,BTOST 790
IF (ALTRT(I).LE.O.) GO TO 100 795
PRIBT 630, I,ALTRT(I) 800
100 COBTIBOE 805
00 110 1-1,BTOST 810
IF (OLTBT(I).LE.O.) GO TO 110 815
PRIBT 640, I,OLTRT(I) 820
110 COBTIBOE 825
C 830
C PRIBT OBHSOAL ROOTIBG PATTERB 835
C 840
IF (BSTR.EQ.1) PRIBT 650, ISTR,JSTR 845
C 850
C PRIBT IBITIAL PROFILE BEFORE OBIT COBTERSIOH 855
C 860
120 COBTIBOE ' 865
PRIBT 660 870
DO 130 J-1,BTOST 875
PRIBT 670, J,I(1,J,1) ,I(2,J,1) ,X(3, J,1) ,T(1,J,1) ,T(2, J,1) ,I(3,J,1) 880
1,TPH0F(J,1) 885
130 COBTIBOE 890
C 895
C CBECK FOE BOB-ZERO IBITIAL PROFILE 900
C 905
IF (IPRO.EQ.O.ABD.ICHHGE.EQ.O) GO TO 150 910
C 915
DO 140 J-1,BTOST 920
C 925
25
C COBTERT TC BOLARITT 930
C 935
X(1,J,1)«X(1,J,1)/238. 940
T(1,J,1)-T(1,J,1)/23e. 945
X(2,J,1)-X(2,J,1)/232. 950
T(2,J,1)-T(2,J,1)/232. 955
140 COBTIBOE 960
150 COBTIBOE 965
C 970
C PRIBT APPROPRIATE TITLE FOR IRATIO 975
C 980
IP ([RATIO.EO-O) GO TO 160 985
PBIBT 680 990
GO TO 170 995
160 COBTIBOE 1000
PRIBT 690 1005
170 COBTIBOE 1010
C 1015
C PRIBT PHASE RATIOS 1020
C 1025
DO 200 1-1,BTOST 1030
IF ([. Eg. 1) GO TO 180 1035
IF (PRATIO.EQ.RAIO(I)) GO TO 190 1040
llt-I-1 1045
PBIBT 700, HST,IB,PHATI 0 1050
180 COBTIBOE 1055
PRATtO-HAIO(I) 1060
BST-I 1065
190 COBTIBOE 1070
IF (I.IT.BTOST) GO TO 200 1075
PRIBT 700, BST,I,PRATIO 1080
200 COBTIBOE 1085
PRIBT 710 1090
C 1095
C 1100
C A IS AOOEODS IBTERSTAGE PLOB 1105
C 0 IS ORGABIC IBTERSTAGE PLOB 1110
C AA IS TOTAL AQOEOOS FLOB LEATIBG STAGE 1115
C 00 IS TOTAL ORGABIC FLOB LEATIBG STAGE 1120
C ECICA IS AQDEODS RECTCLE TO SATISFY PHASE RATIO 1125
C RCICO IS ORGABIC RECYCLE TO SATISFY PHASE RATIO 1130
C AR IS AQOEOOS RECYCLE TO FILL STAGE 1135
C OR IS ORGABIC RECYCLE TO FILL STAGE 1140
C AT IS TOTAL AQOEODS FLOB HITRIB A STAGE 1145
C OT IS TOTAL ORGABIC FLOB BITRIH A STAGE 1150
C TT IS TOTAL FLOB OH TO A STAGE 1155
C 1160
DO 210 1=1,BTOST 1165
AA(I)=ALTBT(I)*A(I) 1170
00(I)=OLTBT(I)tO(I) 1175
RATI0(1)=AA(I)/00(I) 1180
AB (I) -RECY (I) •RAIO (I) / (1. 0»RAI0 (1) ) 1185
OR(l)=RECY(I)/(1.0tRAIO(I)) 1190
AT(I)-AA(I)»AR(I)*RCYCA(I) 1195
OT(I)-00(I)*OR(I)tRCTCO(I) 1200
TT(I)-AT(I)»OT(I) 1205
210 COBTIBOE 1210
C 1215
C PRIBT STAGE FLOBS 1220
C 1225
DO 220 J-1,BTOST 1230
PRIBT 720, J,RATIO(J),AA(J),AR(J),RCYCA(J) ,AT(J) 1235
220 COBTIBOE 1240
PRIBT 730 1245
DO 230 J=1,BTOST -» 1250
PBIBT 740, J,TT(J) ,STTOL(J) ,00(J) ,OR(J) ,RCICO(J) ,OT(J) 1255
230 COBTIBOE 1260
BRITE (6,750) 1265
C 1270
C H - IS 100 OH LAST TIHE IBTERTAL 1275
C ICEB - COOBTS BY 100"S HDHBER OF ITERATIOHS 1280
C ADF - nAX(DF(1),DF(2).DF(3)) FOR A GITEH TIBE IBTERTAL 1285
C BDF - ADF FOR PSETIOOS TIHE IBTERTAL 1290
C CDF - SET EQDAL TO ADF IP ADF-BDF (COBTEEGEHCE HAS OCCORRED) 1295
C DDF - DETERHIHES PDP,PPDF,PPPDF,PPPPDF (PERCEBT STEADY STATE) 1300
C 1305
H-100 1310
ICEB-0 1315
BOP-0. 1320
CDF-0. 1325
DDF-0.10 1330
C 1335
C IBITIAIIZATIOB OF TOTAL BOLES IH (SOLIH) AHD TOTAL HOLES OOT (SOLOOT) 1340
C 1345
DO 240 1=1,3 1350
SOL1B(I)=0.0 1355
SOLOOT(I)=0.0 1360
240 COBTIBOE 1365
DO 250 1-1, 3 1370
DO JSO J-1,BTOST 1375
SOLlH(I)='iOLIB(I)*APDRT(J)*XFD(I,J) •OFDRT (J) *IFD (I, J) 1380
C IF(JFD(I,J).LT.O) SOLIB(I) = SOLIB(I) - AFDRT(J) •XPD (I, J) 1385
250 COBTIBOE 1390
C 1395
SPH=0. 399*CTBP»0.321*(1. O-CTBP) 1400
26
C 1«05
C K - TIBE IBTESTAL CODBTEB (BOD 100) 111 10
C 1015
260 COBTIBOE 1«20
B-1 1*25
C 11130
C STAGES PERrOBBS STAGE CAICOIATIOBS FOB EACH TIHE IBTBBTAl 1«35
C IMO
270 COBTIBOE 11105
CALL STAGES (SPH) 1050
DO 280 J'l.BTOST 1055
TPB0r(J,1)-TPS0r(J,2) 1060
DO 280 I>1,0 1*65
I(I,J,1)<I(I,J,2) 1070
T(I,J,1)-T(I.J,2| 1075
280 COITIBOE 1080
C 1085
C ir BOBBEB or TIBE IBTEBTALS IS LESS THAB THE BDBBEB OE STAGES, 1090
C OBIT CALCOLATIOB Or SOLIB.SOLODT, ABD PEBCEBT STEADT STATE 1095
C 1500
IP (ICEB.,EQ.0.AB0.K.LT.BTOST| GO TO 390 1505
C 1510
C ir SBDTDOBI (ICALC'O), CRECB rOB BEAB ZEBO COBCEBTBATIQB 1515
C 1520
ir (ICALC.EC. 1) 60 TO 300 1525
C 1530
C rOE SHDTDOBB, ICALC-0, TEST EOB COBPOBEBTS 1 ABD 2 » 0, HOT COBPOBEBI 1535
C 1500
DO 290 J-1,BT0ST 1505
ir (I(1,J,1|.GT.0.000001) GO TO 390 1550
IF (T(1,J,1).GT.0.000001) GO TO 390 1555
ir (I(2,J,1).GT.0.000001) GO TO 390 1560
IF (T(2,J,1).6T.0.OO0001) GO TO 390 1565
290 COBTIBOE 1570
00 TO 380 1575
C 1580
C CALCDLATIOB OF PEBCEBT STEADT STATE 1585
C 1590
C 1595
C CALCULATE TOTAL HOLES IB (SOLIB) ABD TOTAL BOLES OOT (SOLODT) 1600
C 1605
300 COBTIBOE 1610
DO 310 1-1,3 1615
SOLODT(I)-A(BTOST)•!(I.BTOST,1)tO(1)•!(I,1,1) 1620
DO 310 J-1,BTOST 1625
S0L00T(I)-SOl00T(I)tALTBT(J)*I(I,J,1)t0L?RT(J)«I(I,J,1) 1630
310 COBTIBOE 1635
C 1600
DO 330 I«1,3 1605
IF (SOLIB(I).GT.O.) GO TO 320 1650
Dr(I)-0. 1655
GO TO 330 1660
320 COBTIBOE 1665
DF(I)<ABS((SOLIB(I)-SOLODT(I))/SOLIB(I)) 1670
330 COBTIBOE 1675
ADF-ABAX1(DF(1),DF(2),DF(3)) 1680
C 1685
C IF ADF-BDr (PEBCEBT STEADT STATE POE COBPEBT TIHE IBTEB7AL 1690
C SAHE AS PBETIOOS TIBE IBTEBTAL), TEBBIBATE CALCDLATIOBS 1695
C 1700
C IP(ADr.BE.BDr) 00 TO 175 1705
i r (ABS((ADr-BDr)/ADr).GT.5.E-0.OB.ADr.GT.DDr) GO TO 300 1710
CDr-AOP 1715
GO TO 380 1720
C 1725
300 COBTIBOE 1730
BDF-ADF 1735
IF (ICBB.OE. 3 . ABD.DDF.LT. 0.000) DDF-2.«DDF 1700
IF (ADF.GT.DDF) GO TO 390 1705-
C 1750
C 90% STEADI STATE 1755
C 1760
IF (DDF.LT.0.100) GO TO 350 1765
J«-K»100»ICEB 1770
PDF-100.«(1.0-DDF) 1775
DDF-0.05 1780
00 TO 390 1785
C 1790
C 95t STEADT STATE 1795
C 1800
350 COBTIBOE 1805
IF (DDF.LT.0.05) GO TO 360 1810
JJK-K*100»ICEB 1815
PPDr-100.«(1.0-DDr) 1820
DDr»0.010 1825
GO TO 390 1830
C 1835
C ir 99« STIADT STATE BAS BEEB BEACHED ABD ICEB>0, PBIBT TIBE IHTEBVAL 1800
C ABD PEBCEBT STEADT STATE 1805
C 1850
360 COBTIBOE 1855
IP (ICEH.GT.0) GO TO 380 1860
C 1865
C 1870
C 99» STEADI STATE ABD ICEB«0 1875
27
C 1880
i r (DDF.IT.0.010) GO TO 370 1885
JJJK"=Kt100*ICEB 1890
PPPDP»1O0.*(1.0-DDr) 1895
DDF-0.005 1900
GO TO 390 1905
C 1910
C IF 99.511 STAEDT STATE HAS BEEH REACHED ABD ICEH>0, PRIHT 1915
C TIBE IHTERVAL ABD PERCEBT STEADI STATE 1920
C 1925
370 COBTIBOE 1930
IF (][CEB.GT.O) GO TO 380 1935
C 1900
C 99.S» STEADT STATE ABD ICEB-0 1905
C 1950
IF (DDF.LT.0.005) GO TO 380 1955
JJJJI!"-B»100»ICEB 1960
PPPPDF-100.«(1.0-DDP) 1965
DDF-0.001 1970
GO TO 390 1975
C 1980
380 COBTIBOE 1985
H»1t 1990
390 COBTIBOE 1995
C 2000
C DETEEHIBE IF TIBE IBTERVAL IS TO BE PHIBTED 2005
C 2010
IF (^[PR.EO.O) GO TO 000 2015
IF (UPS. Eg..-1. ABD. (K.EQ.20.OR.B.Eg.00.aB.K.E<).60.0R.K.E9.80.OR.K.E 2020
10.HI) SO TO 400 2025
IF (XPR.EC. 1.ABD.K.EQ. B) 60 TO 000 2030
GO TO O70 2035
C 2000
000 COBTtBDE 2005
DO 020 J-1,BTOST 2050
AOBT(J)«A(J(*ALTRT(J) 2055
OO0T(J)-0(J)*OLVRT(J) 2060
IP (IOOI.EQ.0) 60 TO 010 2065
AQOOT(J)—ALTET(J| 2070
OROOOT (J) —OLVRT (J) 2075
010 COBTIBOE 2080
C 2085
C COBVEBT TO G/L 2090
C 2095
X(1,J,2)-T(1,J,2)*23e. 2100
T(1,J,2)'I(1,J,2)»238. 2105
I(2,J,2)=X(2,J,2)»232. 2110
T(2,.J,2)-T(2,J,2)*232. 2115
C 2120
020 COBTIBOE 2125
C 2130
C PRIST t ABD T OH DISK (OBIT 10) FOR CALCOBP PLOTS 2135
C 2100
L=K+100»ICEH 2105
DO 030 J=1,BTOST 2150
WRITE (11) L, (r(I,J,2 ) ,1=1,0) ,(I(I,J,2 ) ,I = 1,0),AODT(J),OO0T(J) 2155
030 COBTIBOE 2160
C 2165
C PBIHTOOT FOR A GIVEN TIHE IBTEBTAL 2170
C 2175
PRIHT 760 2180
DO 050 J=1,BTOST 2185
IF (J.GT.1) GO TO 000 2190
L=IC»100»ICEH 2195
PRIHT 770, L,J,(X(I,J,2) ,1=1,3) ,AaaT(J) 2200
IF (ALVBT(J).6T.0.) PRIBT 780, AQOOT(J) 2205
60 TO 050 2210
000 COBTIBOE 2215
PSIBT 790, J,(1[(I,J,2) ,1=1,3) ,A00T(J) 2220
IF (ALTHT(J).6T.O.) PRIHT 780, AQODT(J) 2225
050 COBTIBOE 2230
PRIHT 800 2235
DO 060 J=1,BTOST 2200
IF (OL»8T(J)-GT-0.) PRIHT 780, ORGODT(J) 2205
PRIHT 810, J,(T(I,J,2) ,1=1,3) ,OO0T(J) ,CODDB(J) ,TPR0P(J,2) 2250
C 2255
C CHECK FOR THIRD PHASE 2260
C 2265
IF (T(2,J,2)..GT.S0RT(. 0311727E-1-.3317607E-0*(TPROP(a,2)-55.03193) 2270
1**2t.36719*(I(3,J,2)-.3055787)«*2)*232.) PRIHT 820 2275
060 COBTIBOE 2280
C 2285
C 2290
C IBCREHEHT TIBE IBTEBTAL COOHTBR 2295
C 2300
070 COBTIBOE 2305
!t=Kt1 2310
C 2315
C 2320
C IT K LESS THAH OR EODAL TO H,PERFOEH STAGE CALCDLATIOBS FOR AH0THE8 2325
C TIHE IHTERTAL 2330
C 2335
IF (K-IE. S) 60 TO 270 2300
C 2305
C 2350
28
C IF SHDTDOBB (ICALC-0), BYPASS PRIBTOOT OF PERCEBT STEADI STATE 2355
C 2360
IF (ICALC.EO.O) GO TO 080 2365
C 2370
C 2375
C PRIBT PEBCEBT STEADI STATE AFTER BTERY 100 TIBE IBTERTALS 2380
C ABD OTBEB CASES WBEB K-B 2385
C 2390
APDF-100.«(1.0-ADF) 2395
PRIBT 830 2000
PRIBT 800, APDF 2005
IF (DDF.LT.0.1) PRIBT 850, PDF,JK 2010
IF (DDF.LT.0.05) PBIBT 850, PPDF,JJK 2015
IF (ODF.LT.0.01) PRIHT 850, PPPDF,JJJ« 2020
IF (DDF.LT.0.005) PRIBT 850, PPPPDF,JJJJB 2025
IF (lOOT.GT.O) PRIBT 860 2030
C 2035
C CHECK FOR C0BTER6EBCE 2000
C 2005
IF (AEP.BE.CDF) GO TO 080 2050
C 2055
PRIBT 870 2060
080 COBTIBDE 2065
IF (B.LT.100) 60 TO 090 2O70
ICEB-ICEB»1 2075
IF (ICEB.EQ.15) 60 TO 090 2080
60 TO 260 2085
C 2090
C COBTERT TOTAL BOLES IB AHD ODT TO G/L ABD PBIBT 2095
C 2500
090 COBTIBOE 2505
SOLIB(1)-SOLIB(1)•238. 2510
SOL0DT(1)-S0L0DT(1)*238. 2515
SOIIB (2)-SOLIB (21*232. 2520
S0L0DT(2)>S010IIT(2)*232. 2525
PRIHT 880, SOLIB(1),S0L0DT(1) 2530
PRIBT 890, SOLIB(2).SOLODT(2) 2535
PRIHT 900, S0LIB(3),S0L0DT(3) 2500
C 2505
C PRIBT BIBDTES PER TIHE IHTERTAL 2550
C 2555
TIBIBT-STT0l(1)/(AT(1)»0T(1)) 2560
IF (IT0L.6E.2) PRIHT 910, TIBIBT 2565
C 2570
C CALCOBP PLOTTER OOTPDT 2575
C 2580
IF (ICALPL.EO.O) GO TO 500 2585
CALL PLCTR (R«100*ICEH,ICALPL,BTOST,SOLIB) 2590
500 COBTIBDE 2595
C 2600
C PDBCBED CARD ODTPDT 2605
C 2610
IF (IPBCH.LE.O) 60 TO 520 2615
DO 510 J-1,BTOST 2620
POBCH 920, (X(I,J,2),I-1,3),(I(I,J,2),I-1,3) ,TPR0F(J,2) ,X(0,J,2) 2625
510 COBTIBDE 2630
C 2635
C 2600
C BEAD BEXT CASE 2605
C 2650
520 COBTIBDE 2655
READ 930, ICRBGE,IDIF 2660
IF (lOIF.EQ.I) READ 930, ICALC,IO0T,IPH,IPBCH,ICALPL 2665
IF (ICHBGE.E0.0) 60 TO 10 2670
IF (ICBB6E.Ea.2) STOP 2675
C 2680
C IF ICHB6E - 1, BEW FEED STREAH AHD ODT CARDS ARE REODIRED FOR EACH 2685
C FEED ABD EXIT STRIAE,ETEB IF BO CHAHGE IS BADE. HO PROFILE CARDS 2690
C ARE REODIRED. 2695
C 2700
C ICBBGE-1 2705
C 2710
C IBITIAL PROFILE IS SET AS TBAT AT EHD OF LAST CALCDLATIOB. 2715
C 2720
IPBO-O 2725
DO 530 1-1,0 2730
DO 530 J-1,BTOST 2735
I(I,J,1)-X(I,J,2 ) 2700
1(1,J,1)-Y(I,J,2 ) 2705
530 COBTIBDE 2750
60 TO 00 2755
C 2760
C 2765
C 2770
500 FOBHAT (12) 2775
550 FORBAT (10A8) 2780
560 FORBAT (F8.0,1012) 2785
570 FORBAT (•1CALCDLATIOBS FOR A SOLTEHT EXTRACTIOB PROCESS',' HATIHG 2790
13 IBTERACTIBG SOLDTES'/) 2795
580 FORBAT (• ',10A8) 2800
590 FORHAT ('OFEED STREAB DATA',16X,'STAGE*,3X,'DRABIDB',6X,'TBORIDB', 2805
17X,•Ba03',5X,'ALDBIBOB B03',2X,'FLOW RATE',6X,*TEBP',/,30X,*B0.', 5 2810
2X,'(6/1)',81,'(G/L)',8X,*(B)•,10X,*(B)',6X,'(TOL.OBITS)',5T,"DE 6 C 2815
3'/) 2820
600 FOBHAT (20X,'AQDEOOS ',10,0(2X , 1PE11.0),3X,E9.3,3X,0Pr7.1 ) 2825
29
610 rOIBAT (19X.2PF6.1,* - TBP *,I0.0(3I,1PB10.0),3X,E9.3,3X,0PF7.1) 2830
620 FORBAT ('O*) 2835
630 FORBAT (• AODEODS STIEAH BIROTBD AT STAGE*,13,*. FLOI •*,1PE9.3) 2800
600 FORHAT (• OBGAHIC STREAH IIHOTED AT STA6E*,I3,*. FLOB -•,1PE9.3) 2805
650 FORBAT (*00R6AHIC STREAH IXITIB6 AT STA6E*,I0,' FEEDS STA6I*,10) 2850
660 FORHAT (•OIBITIAl PROFILE*,171,*STAaE',13X,*AODEODS PBASI*,26X,*0R 2855
16AIIC PHASE*,101,/,30X,*BO.*,OX,*DRABIDH',6X.'TH0«IDH*,7X,'HBa3',8 2860
2I,*0RABIDB*,6I,'THORI0H*,7X,'HHO3',7X,'TEBP',/,31t,'(ODTTLOW)',2X, 2865
3*(a/L)*,8X,' (6/L)',SX,*(H|*,10X,' (G/L) ',8X,* (G/L) * ,8X,'(H) • ,8X,'DE 2870
OG C*,/) 2875
670 FOIBAT (321,10,6(3X,1PE10.0),3x,0PFe.1) 2880
680 FOIBAT ('OPRASE RATIOS DO HOT EQDAL FIOH RATIOS',3X,'STAGE B0.',2X 2885
1,'A/O PBASI RATIO*,/) 2890
690 FORHAT ('OPRASE RATIOS EODAL FLOI RATIOS*,11X,*STAGE B0.*,2X,*A/0 2895
1PHASE RATIO*,/) 2900
700 FORBAT (00X,I«,* -',I0,0X,1PE10.3) 2905
710 FORHAT (/,' STAGE FLOWS;',/,51X,'AODEODS PHASE',/,07X,'RBCICLE 2 2910
1RECTCLE 1 TOTAL',/," STAGE FLOW',22X,'FLOW ODT TO FILL 2915
2F01 PHASE*,* AODEODS',/,' HO. RATIO',21X,'OF STAGE STAG 2920
31',' RATIO FLOW') 2925
720 FORHAT (2X,I2,OX.1PE9.3.17X,0(2X,E9.3)) 2930
730 FOBHAT (/,9X,*0TIRALL*,36X,'0RGAHIC PHASE',/,10X,'TOTAL',32X,'RECI 2935
ICII 2 RECYCLE 1 TOTAL',/,* STAGE STAGE*,21X,* FLOW ODT TO 2900
2 FILL roB PRASE ORGAHIC',/,' BO. FLOW STAGE T0L.',7X, 2905
3*0F STAGE*,' STA6E RATIO FLOW') 2950
700 rOBHAT (2X,I2,2I,2(2X,1PE9.3),ex,0(2X,E9.3)) 2955
750 rORBAT (27B1TRABSXEBT BIHATIOB RESOLTS,10X,09HSD10TI COHCEBTBATIOB 2960
IS IB ErrLDIBT raoR EACR STAOE) 2965
760 rORBAT (///,311,'AQOEOOS PHASE',/,3I,'TIBE",OX,'STAGE*,31,*ORARI0H 2970
1*,0X,*TH0RIDH',5X,*HB03*,5X,'FL0« OOT',/,' IBTIRTAL',3X,'BO.',5X,' 2975
2(G/L)',6I,* (G/1)*,7I,* (H)*,5X,'0F STAGE') 2980
770 PORHAT (1X,IO,7X,I2,2X,0(1X,1PB10.3)) 2985
780 FORBAT (121,3(10X,'«'),5X,1PE10.3) 2990
790 FORBAT (12I,I2,2X,5(II,1PE10.3),12X,E10.3) 2995
800 FOBBAT (/.31X,'086ABIC PHASE',/,1IX,'STAGE',3X,'DRABIDB',OX,'THORI 3000
10B',5X,'HB03',5X,*FLOW ODT',6X,'ITIRA-',2X,*TEBP*,/,12X,'BO.*,5X,* 3005
2(G/l)',6X,* (6/l)',7X,'(H)',5X,'0F STA6E', 6X,'HOBS', 3X,'DEG C ) 3010
810 FORBAT (12X,I2,21,0(2X,1PE9.3),ex.0PF3.0,OX,FO.1) 3015
820 FORHAT (*•*,37X,'t*) 3020
830 FORBAT ('Ot IBDICATES THIRD PHASE WITBIH 10V TOLEBABCE*) 3025
800 FORBAT (*0*,F7.2,* PEBCEBT STEADT STATS BY OTEBALL BATERIAL',' BAL 3030
1ABCI') 3035
850 FORHAT (' ',F6.1,2X,'PFRCEBT OF STEADT STATE AFTER TIBE',' IBTIBTA 3000
1L',I6) 3005
860 FOBBAT (1H0,' • SIDI STREAB REBOTED FROB INTERSTAGE PLOW DETWEEB 3050
1 TBI IBSICATID STAGES.') 3055
870 rOHBAT (1H0,'CALCDLATIOB TERHIHATED BY C0HTER6EHCI AT IHDICATED 3060
1STEADT STATE.') 3065
880 rOBBAT (/,' TOTAL G/L DBABIDH IB:',2X,1PE11.0,', OOT:',Ell.0) 3070
890 rOBBAT (• TOTAL G/L THORIDR IB!',2X,1PI11.0,', OOT:',Ell.0) 3075
900 FORHAT (' TOTAL HOLES HB03 IB:*,5X,1PE11.0,*, OOT:*,Ell.0) 3080
910 FORBAT (/,* HIBDTIS/TIHI IHTERTAL:',F7. 3) 3085
920 FORBAT (1P8E10.0) 3090
930 FOBBAT (512) 3095
EHD 3100-
SDBROBTIBE PROFIL PRO 5
COHHOB /PBOF/ IPRO,X(0,100,2),T(0,100,2),XFD(0,100),YFD(0,100),TPRPRa 10
1OF(100,2),AFDT!B(100),OFDTEH(100) PRO 15
COHHOM /FEED/ HTOST,ArDRT(100),OFDRT(100) PRO 20
C PRO 25
C DEPIHI FEED STREAHS PRO 30
C 1 - STAGE BOBBER THAT FEED EHTERS JHAS = 1 FOR AQOEODS FEED PRO 35
C FDRT - FEED FLOW RATE (TOLDHE OBITS) 0 FOR ORGAHIC FEED PRO 00
C C0H1 - DBABIDH COHTEBT IB BOLABITY IBDEX = 1 IF ROBE FEED CARDS FOPEO 05
C C0B2 - THORIDR COHTEBT IB HOLABITT 0 FOR LAST CARD PRO 50
C C0B3 - ACID HOLARITY TEHP = TEHPERATDEE IB DEGREES PRO 55
C COBO - BITRATE BOLARITY FROH IHEITRACTABLE SALTS, IBCLDBIHG THE REOOCPRO 60
C PRO 65
10 COBTIIOE PBO 70
READ (5,90) I,JBAS,FDBT,COH1,COB2,COH3,TEBP,IHDEX,COHO PRO 75
IF (JHAS.EQ, 1) GO TO 30 PRO 80
C PRO 85
C ORGAHIC FIID PRO 90
C PRO 95
OFDHT(I)=FDHT PBO 100
IF (TEBP.EQ.O.) 60 TO 20 PRO 105
0FDTia(I)=TEBP PRO 110
20 COBTIBDE PRO 115
TFD(1„I)-C0B1 PRO 120
YFD(2,I)=C0H2 PRO 125
TPD(3,I)=COB3 PRO 130
IFD(0,I)-COBO PRO 135
60 TO 50 PRO 100
C PBO 105
C AQOEODS FEED PRO 150
30
C PRO 155
30 COBTIBDE PBO 160
AFBBT(I)-FDBT PBO 165
IF (TEBP.IO.O.) GO TO 00 PBO 170
AFDTER(I)-TIHP PRO 175
00 COBTIBDE PBO 160
XPD(1,I)-C0B1 PBO 185
XFD(2,I)-COB2 PBO 190
XFD(3,I)-CaB3 PBO 195
XFD (0,1)-COBO PRO 200
C PBO 205
C BORE CABDS7 PRO 210
C PBO 215
50 COBTIBDE PRO 220
IF (IBDIX.IO.1) 60 TO 10 PRO 225
C PRO 230
C PRO 235
C CHICK FOB BOB-IIBO IBITIAL PBOFILI PRO 200
C PRO 205
IF (IPRO.GT.0) 60 TO 70 PRO 250
C PBO 255
C FOR IPBO - 0, TBI IIBO PBOFILI IS ALBIADY II X ABD I ABBAIS PRO 2(0
C PRO 265
DO 60 J-1,IT0ST PRO 270
TPR0F(J,1)-APDTEB(1) PRO 275
60 COBTIBDE PBO 280
lETDBB PRO 285
C PRO 290
C FOR A BOB-XERO IBITIAL PBOFILE, TALDES ARE BOW READ PRO 29 5
C PRO 30 0
70 COBTIBDE PRO 305
DO 80 J-1,BTOST PBO 310
READ (5,100) (X(I,J,1) ,1-1,3) ,(Y(I,J,1),I-1,3),TPR0F(J,1) ,X(0,J,1)Pia 315
IF (TPtar(J,1).EQ.0.| TPIOF(J,1)-25.0 PRO 320
80 COITIBOE PRO 325
C PRO 330
C PRO 335
RITORB PRO 300
C PRO 305
90 FORBAT (2I2,5F8.0,I2,F8.0) PIG 350
100 FOBBAT (8F10.0) PBO 355
BID PRO 360-
SDBBOOTIIE FLOIS (STTOL) FLO 5
COBBOB /FLIS/ IRATIO,IODT,ITOl,RICY(100) ,ICHHGI,ALTRT(IOO),0LTRT(1FL0 10
100) FLO 15
COHHOB /PUD/ ITOSI,AFDBT(100),OFOBT(100) FLO 20
COBBOB /BIC/ A(100),0(100),RAIO(IOO) ,RCYCA(100) ,RCYCO(100) PLO 25
COHHOB /BODTI/ BSTB,ISTB,JSTB FLO 30
DIBEBSIOB STTOL(100) FLO 35
C FLO 00
C DEFIBE ODTGOIBG STBEARS OTHER THAB AODEODS RAFTIBATE ABD ORGAIIC PRODFLO 05
C I - STAGE BDBBIR THAT STREAH IIATIS JHAS - 1 FOR AODIODS STREAH PLO 50
C IBDEX - 1 IF ROBE CARDS FOLLOI 0 FOR ORGAIIC STREAH FLO 55
C 0 FOR LAST CARD OTRT • EXIT FLOI RATE (TOLDRE FLO 60
C FLO 65
C FLO 70
C CHECK FOR OOTGOIIG STREAHS (IODT»0-IO,IODT-1-IES) FLO 75
C FLO 80
IF (lODT.EQ.O) GO TO 60 FLO 85
C FLO 90
C READ ODTGCIIG STREAB DATA FLO 95
C FLO 100
C FLO 105
IBA-1 FLO 110
lEO-ITOST FLO 115
10 COBTIIOE FLO 120
READ 290, I,JHAS,OTRT,IBDEX FLO 125
IF (JHAS.EQ. 1) GO TO 30 FLO 130
C FLO 135
C ORGAHIC STBEAR FLO 100
C FLO 105
OLTRT(I)-0TRT FLO 150
IF (OLTRT(I).BE.O.) GO TO 50 FLO 155
DO 20 J-I,ItO FLO 160
OLTBT(I)-0LTBT(I)»OFDRT(J) FLO 165
20 COITIIDI FLO 170
110-1-1 FLO 175
GO TO 50 FLO 180
C FLO 185
C AQDEOOS STREAB FLO 190
C FLO 195
30 COBTIIOE FLO 200
ALTBT(I)-OTRT FLO 205
IF (ALTRT(I).lE.O.) GO TO 50 FLO 210
DO 00 J-IBA,I FLO 215
31
ALTRT(I)-ALTRT(I)•AFDRT(J)
00 COBTIBDE
IBA-I,»1
BORE CARDS?
50 COBTIBDE
IF (IBDEX.IQ.1) 60 TO 10
IE DO lOT PROTIDE FOR RITDRIIBG STREAHS
SET AODIODS ABD ORGABIC IITEBSTAGB FLOI (A ABD 0)
60 COBTIIOE
READ DIDSDAL E0DTIB6 PATTEBB
0B6AHIC STREAH FROH STAGE ISTR FEEDS IBTO STAGE JSTB
ir (BSTR.EO.1) READ 300, ISTR,JSTB
DO 80 1-1,BTOST
IP (I.GT.1) 60 TO 70
A(1)»APDRT(1)-ALTBT(1)
0(BTOST) -OTDRT (BTOST) -OLTRT (HOST)
GO TO 80
70 COBTIBOE
A (1) -A (1-1) •AFDRT (I) -ALTRT (I)
IST-HT0ST»1-I
0(1ST)«0(IST»1)•OFDRT(BST)-OLTRT(BST)
i r (nSTR.ia.l.AID.BST.IO.ISTB-1) 0 (BST)-OFDRT (BST)-OLTBT (BST)
IF (HSTR.ig.1.AB0.IST. IQ. JSTB) O(BST)-O(ISTB) •OPDRT(IST)-OLTRT
l)^O(HST^I)
80 COITIBOE
THE IBTBBSTAGE FLOWS ABE BOW SET
DEFIBE PHASE RATIOS
DO 90 1-1,BTOST
RCTC»(I)-0.0
BCICO(1)-0.0
90 COBTIBDE
FlOBAX-0.0
K-1
IF (IBATIO.E0.1) GOTO 110
DEFIBE PHASE BATIO WHII IRATIO-O,PHASE RATIO-FLOW RATIO
DO 100 1-1,BTOST
RAIO (1) • (ALTRT (I) •A (I) ) /(OLTRT(I) •O (I) )
100 COBTIIOE
60 TO 130
DO BOT RIDEFIHE PHASE BATIO IF ICHB6E-1
110 COBTIBOE
IF (ICHI6E.Ea. 1) 60 TO 150
DEFIBE PHASE RATIOS WREB IBATIO-1, PBASI BATIOS DO BOT EQDAL
FLOW RATIOS
I-BDHBER OF STAGES HATIHG PRASE RATIO SHOII
BEAD 310, I,RATIO
I-K^I-1
DO 120 J=K,I
RAIO(J)-RATIO
120 COBTIBDE
IF (I.EQ.BTOST) 60 TO 150
K-1^1
GO TO 110
CALCOLATE IICESSART RECYCLE FLOIS
TYPE 1 IS RECYCLE OF OBE PHASE OILY, TO SATISFY PHASE BATIO
DSED WHEI IRATIO = 1
TYPE 2 IS RECYCLE OF BOTH PHASES IITH PHASI RATIO RAIO(I)
TO ALLOW TOTAL STAGE FLOW TO EQDAL FLOBAX
DSED IBIB ITOL =0,2
CALCOLATE FLOTOT.BECY,ABD FLOBAX FOR IBATIO-O; ITOL-0,2,3
FOR ITOL-1, RECY SET TO 0.0
FLOTOT - TOTAL PLOW LIATIHG STA6E
FLOBAX - GREATEST TOTAL FLOW FROH A STAGE
130 COBTIBOE
IF (ITOL.EQ. 1) GO TO 210
DO 100 1=1,BTOST
FLOTOT=ALTRT (I) •A (I) •OLTRT (I) •O (I)
RICI(I)-FLOTOT
PL0HAX-A8AX1(PL0T0T,FLORAX)
100 COBTIBDE
IF (ITOL.EQ.O) 60 TO 190
IP (ITOL.E0.1) 60 TO 210
IF (ITOL.EO. 2.0S.IT0L.EQ.3) 60 TO 230
FLO 220
FLO 225
FLO 230
FLO 235
FLO 200
FLO 205
FLO 250
FLO 255
FLO 260
FLO 265
FLO 270
FLO 275
FLO 280
FLO 285
FLO 290
FLO 295
FLO 300
FLO 305
FLO 310
FLO 315
FLO 320
FLO 325
FLO 330
FLO 335
FLO 300
FLO 305
FLO 350
FLO 355
FLO 360
FLO 365
(BSTFLO 370
FLO 375
FLO 380
FLO 385
FLO 390
FLO 395
FLO 000
FLO 005
FLO 010
PLO 015
FLO 020
FLO 025
FLO 030
FLO 035
FLO 000
FLO 005
FLO 050
FLO 055
FLO 060
FLO 065
FLO 070
FLO 075
FLO 080
FLO 085
FLO 090
PLO 095
FLO 500
FLO 505
FLO 510
FLO 515
FLO 520
FLO 525
FLO 530
FLO 535
FLO 500
FLO 505
PLO 550
FLO 555
FLO 560
FLO 565
FLO 570
FLO 575
FLO 580
FLO 585
FLO 590
FLO 595
FLO 600
FLO 605
FLO 610
FLO 615
FLO 620
FLO 625
FLO 630
FLO 635
FLO 600
FLO 605
FLO 650
FLO 655
FLO 660
FLO 665
FLO 670
FLO 675
FLO 680
FLO 685
FLO 690
32
C CALCOLATE FLOTOT,RCYCA,RCYCO,RECI,FLOHAX FOR IRATIO-1
C ABO ITOL-0,1,2,3
C AFL - TOTAL AQDEOOS FOLI LEATIBG A STAGE
C OFL - TOTAL ORGABIC FLOI LEATIB6 A STAGE
C RFL - BATIO OF AQOEODS TO OBGABIC FLOI LEATIIG A STAGE (AFL/OFL)
C RCYCA - AQDEOOS BECTCII TO SATISFY PHASI BATIO
C HCYCO - ORGABIC BICYCLE TO SATISFY PHASE RATIO
C
150 COBTIBDE
DO 180 1-1, BTOST
AFL-ALTRT(I)^A(I)
OFL-OLVRT(I)»0(I)
FLOTOT-AFL^OFL
RFL-APL/OFL
C
C COBPABISOB OF FLOB BATIO (RFL) TO PRASE RATIO (RAIO)
C
IF (BFL.GT.RAIO(I)) GO TO 160
IF (RFL.EQ.RAIO(I)) GO TO 170
FLOTOT-FLOTOT- (RAIO(I) •I.O)/(RFL^I. 0)
RCYCA(I)-PLOTOT-APL-OFL
GO TO 170
160 COBTIBOE
FL0T0T-F10T0T»RFL» (RAIO (I) • 1. 0) / (RAIO (I) • (RFl^l. 0) )
RCYCO(I)-PLOTOT-APL-OFL
170 COBTIBDE
RECI(I)-FLOTOT
FL0BAX-ARAX1(FLOTOT,FLOHAX)
180 COBTIBDE
FLOBAX COBTAIBS TBE LABGEST TOTAL FLOI FROB A STAGE - STAGES BATIBG
TOTAL FLOI LESS TBAB FLOBAX HOST EXBIBIT RECYCLE
OF TYPE 2
IF (ITOL.EQ.0) GO TO 190
IF (ITOL.EQ.1) GO TO 210
IF (IT0L.Ea-2.0I.IT0L. EQ. 3) GO TO 230
ITOL-0; IRATIO-O,1
190 COBTIBDE
DO 200 1-1,BTOST
RECY(I)-FLORAX-RECY(I)
200 COITIBDE
EETDBB
ITOL-1; IRATIO-O,1
210 COITIBDE
DO 220 1-1,BTOST
RECI(I)-0.0
220 COBTIBOE
BITDBH
DITERRIIE STAGE TOLDRE
STAGE TOLDRES ARE ZERO FOR ITOL-0,1
IF ICHB6E-0, READ STAGE TOLDBES
IF ICHIGE-1, BYPASS BEADIBG OF STAGE TOLDBES
230 COITIBOE
IF (ICBBGE.EQ.1) 60 TO 260
IF (IT0L.EQ.2) GO TO 250
ITOL-3; STAGE TOLDBES ARE EQDAL ABD SPECIFIED
IRATIO-O,1
READ 320, TOL
DO 200 J-1,BTOST
STTOL(J)-TOL
200 COBTIBDE
60 TP 260
ITOL-2: STAGE TOLOHES ARE IBDIPIBDEITLY DIEQDAL
250 COITIBDE
READ 330, (STTOL(J),J-1,BTOST)
CALCOLATE LEIOTH OF TIHE IBTEaTAL(TIBIB) ABD RECY FOB ITOL-2,3
260 COBTIIOE
TIHI1=0.0
DO 270 1-1,BTOST
RECY (I) -BICY (I) /STTOL (I)
TIBIB-ABAI1(RICY(I).TIBIB)
270 COBTIBOE
FLOBAX-TIBII
DO 280 1-1,BTOST
BECYfI)=(FLOBAX-RECI(I))*STTOL(I)
280 COBTIBDE
FLO 695
FLO 700
FLO 705
FLO 710
FLO 715
FLO 720
FLO 725
FLO 730
FLO 735
FLO 700
FLO 705
FLO 750
FLO 755
FLO 760
FLO 765
FLO 770
FLO 775
FLO 780
FLO 785
FLO 790
FLO 795
FLO 800
FLO 805
FLO 810
FLO 815
FLO 820
FLO 825
FLO 830
FLO 835
FLO 800
AFLO 805
FLO 850
FLO 855
FLO 860
FLO 865
FLO 870
FLO 875
FLO 880
FLO 885
FLO 890
FLO 895
FLO 900
FLO 905
FLO 910
FLO 915
FLO 920
FLO 925
FLO 930
FLO 935
FLO 900
FLO 905
FLO 950
FLO 955
FLO 960
FLO 965
FLO 970
FLO 975
FLO 980
FLO 985
FLO 990
FLO 995
FLO1000
FL01005
FLO1010
PLO1015
FL01020
FLO1025
PLO1030
FL01035
FL01000
FL01005
FLO1050
FLO1055
FL01060
FL01065
FL01070
FL01075
FL01080
FLO1085
FLO1090
FL01095
FL01100
FL01105
FLO1110
FL01115
FL01120
FL01125
PLO1130
FL01135
FL01100
FL01105
FL01150
FL01155
FL01160
FL01165
33
290 FORHAT (2I2.F8.0,I2)
300 FORHAT (215)
310 PORRAT (12.F8.0)
320 FORBAT (18.0)
330 FORBAT (10F8.0)
IHD
FLO1170
FL01175
PLO1180
FLO1185
FLO1190
FLO1195
FL01200-
SDBBODTIBE STA6ES (SPH) STA
DIHEBSIOI AII(O), 011(0) STA
CORROI /STGS/ ISOLD,CODDB(IOO),AB(IOO) ,OR (100) ,AT(100),OT(100) STA
CORROI /FEED/ ITOST,AFDRT(100),OFDBT(100) STA
COHHOB /PROF/ IPBO,X (1,100, 2) ,Y (0,100,2) , XFD(0,100) ,YFD(0,100) ,TPRSIA
10F(100,2) ,AFDTIH(10a) ,0FDTBH(100)
CORHOH /lie/ A(100),0(100),BAIO(IOO) ,RCYCA(100),RCYCO(100)
COHROB /DIST/ CTBP,ARI(0),DTRY(0),T0
COBBOB /BODTE/ BSTR,ISTR,JSTB
C
C PEBFORR STAGE CALCDLATIOBS FOB CDRRIBT TIHE IBTEBTAL - COBSIDIR
C IBPDTS FIOB AOJACEBT STAGES, FEED STBEARS, ABD HOST IBPORTABTIT
C BECYCII WITHIB THE STAGS - RECYCLE A BESDLT OF STAGS FLOW BSIW6
C LESS THAB FLORA! (TYPE 2) AID PBASI RATIO lOT EQDAL TO FLOI
C RATIO (TYPE 1)
C
C AR - AQOEODS RECYCLE TO FILL STAOE
C OR - ORGABIC BICYCLE TO FILL STAGE
C ECICA - AQDEOOS RECYCLE TO SATISFY PHASE RATIO
C BCYCO - ORGABIC BECYCLE TO SATISFY PHASE RATIO
C AT - TOTAL AQOEODS FLOI IITHIH A STAGE =
C TOTAL AQOEODS FLOI LEATIBG A STAGE • AQOEODS BECYCLE
C TO FILL STAGE • AQOIODS BICYCLE TO SATISFY PRASI BATIO
C OT - TOTAL OBGABIC FLOW WITHIB A STA6I -
C TOTAL 0R6ABIC FLOI LEATIBG A STAGE • OBGAHIC BECYCLE
C TO FILL STAGE • ORGABIC RECYCLE TO SATISFY PHASE RATIO
C XFD - AQDEOOS PEED COBTEBT OF STAOE
C YFD - 0R6ABIC FEED COBTEBT OF STAOE
C AFDRT - AODIODS STA6I FLOW EATF
C OFDRT - OReABIC STA6I FLOW RATE
C X - SOLDTI COICIBTBATIOB IB AQDEOOS PHASE FOR A 6ITEI STAGE
C Y - SOLDTE COBCIBTEATIOH IB ORGAHIC PHASE FOR A GITEI STAGE
C AFDTIH - TIRPIIATDRE OF AODEODS FEED
C OFDTEH - TEHPIBATDRI OF ORGAIIC FEED
C TPROF - TEHPIBATDRI OF SOLDTE IB A GITEB STAGE
C AEY - TRIAL AQDEOOS PHASE COBPOSITIOH
C DTRY - DISTBIBOTIOH RATIOS
STA
STA
STA
STA
STA
STA
STA
STA
STA
STA
STA
STA
STA
STA
5
10
15
20
25
30
35
00
05
50
55
60
65
70
75
80
85
90
95
STA 100
STA 105
STA 110
STA 115
STA 120
STA 125
STA 130
STA 135
STA 100
STA 105
STA 150
STA 155
STA 160
STA 165
STA 170
STA 175
STA 180
STA 185
STA 190
STA 195
STA 200
STA 205
STA 210
STA 215
STA 220
STA 225
STA 230
DO 100 J-1,BTOST
JJ«J-1
JJJ=J^1
DO 50 1-1,ISOLD
IF (J.GT.1) GO TO 10
AIH (I) -AFDRT (J) •XFD (I, J) • (AR (J) •RCYCA (J) ) •X (I, J, 1)
60 TO 20
10 COBTIBOE
AIH(I)=AFDRT(J)»XF0(1,J)^(AR(J)+BCYCA(J))^X(I,3,1)^A(JJ)^X(I,JJ,1) STA 235
20 COHTIIDE STA 200
IF (J.EQ.ITOST) GO TO 30 STA 205
OIH (I) -OFDHT (J) • YFD (I, J) • (OR (J) •RCYCO (J) ) •Y (I, J, 1) •O (JJJ) •! (I,
JJJ , STA 250
11) STA 255
IF (ISTR.EQ.1.AID.J.EQ.ISTR-1) OII(I) =OFDRT(J) *IFD (I, J) • (OB (J) •BCISTA 260
1C0(J))»Y(I,J,1) STA 265
IF (BSTB.E0.1.AHD.J.EQ.JSTR) OIH (I) =0TDRT (J) *IFD (I. J) • (OR (J) •RCICOSTA 270
1 (J))^Y(I,J,1)+0(ISTa)^Y(I,ISTE,1)^0(JJJ)^T(I,JJJ,1) STA 275
60 TO 00 STA 280
30 COBTIBOE STA 285
OII(I)=OFDBT(J)»YPD(I,J)^(OE(J) •BCYC0(J))^Y(1,J,1) STA 290
00 COITIBDE STA 295
A1H(I)=AIH(I)/AT(J) STA 300
0IB(1)=0II(I)/0T(J) STA 305
ARY(I)=X(I,J,1) STA 310
IF (ABY(2).LT.0) ARY (2) = (All (2) •RAIO (J) •OIH (2))/RAIO (J) STA 315
ARy(0)=AIB(0) STA 320
50 COBTIBOE STA 325
IF (J. 5T. 1) 60 TO 60 STA 330
AHTIB=AFDET(J)^AFDTEH(J)^(AS(J) •RCICA(J))*TPR0F(J,1) STA 335
GO TO 70 STA 300
60 COBTIBDE STA 305
AHTIB-AFDRT (J) •AFDTEH (J) • (AR (J) •RCYCA (J)) *TPROF (J, 1) •A (JJ) *TPROF(JSTA
350
1J, 1) STA 355
70 COBTIBOE STA 360
IF (J. EQ. ITOST) GO TO 80 STA 365
OHTIH-SPH*(OFDRT(J)•OFDTEH(J)•(OB(J)•RCYCO(J))*TPROF(J.1)•O(JJJ)'TSTA 370
1PBOF(JJJ,1)) STA 375
IP (BSTB.ig.1..ABD.J.EQ.ISTB-1) OHTIS=SPH^ (OFDBT(J) •OFDTEH (J) • (OB (JSTA 380
l)^RCICO(J)) •TPROF (J.I)) STA 385
IF (HSTR.EQ. 1.ABD.J.Eg.JSTR) OHTIN=SPH» (OFDRT (J) •OFDTEH (J) • (OB (J) •STA 390
34
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
1 RCYCO (J)) •TPROF (J, 1) •O (ISTR) •TPROF (ISTR, 1) •O (JJ J) •TPROF (JJJ, 1))
60 TO 90
80 COBTIBDE
OHTIB-SPB*(OFDRT(J)•OFDTIB(J)•(OB(J)•BCYCO(J))•TPBOF (J,1))
90 COBTIBDI
TPROF (J, 2) • (ARTIB^OBTIB) / (AT (J) •SPB^OT (J) )
TIRPC-TPR0F(J,2)
CCD • 0.001
CCD-0.005
BCK IS TBE ITERATIOH CODBTEB
BDEX IS TBE FIT CODBTEB
RCK-0
100 COBTIBDE
BDIX-0
DCOR CALCDLATES DISTBIBOTIOH RATIOS (DTBY) FOR EACH TBIAL AQOEODS
CORPOSITIOB (ABY) AT TBE STAOE EODILIBBIDR TERPIBATDRE (TEBPC)
CALL OCOB (ABY,DTBY,TEBPC,CTBP)
DO 110 I<1,BS0LD
TSIB-AIB(I)-BAIO(J)•Oil(I)
IF (TSII.LT.O.) DTRY(I)-0.
DTIY(0)-0.
X(I,J,2)-TSIB/ (RAIO(J)•DTBY(I))
IF (ABS(X(I,J,2)).LT.1.E-6 ) GO TO 110
01-(X (I, J, 2)-ARY (I) )/X (I, J, 2)
DIF-ABS(DI)
ARI(I)-X(I,J,2 )
IF (DIF.LE.CCD) GO TO 110
BDEX CODBTS lORBER OF SOLDTES THAT HOT C0ITER6ED
HDEX-BDEX^I
110 COBTIBOE
BCK-BCK^I
ALLOI OILY 20 ITERATIOIS
IT (HCK. 10.20) 60 TO 120
IF ALL SOLDTES RATI lOT C0ITER6ID, RIPIAT CALL TO DCOR
IITH III TALDIS OF ARY
IF (RDIX.GT.O) 60 TO 100
120 COBTIBDI
DO 130 1-1.ISOLD
Y(I.J.2)-I(I,J,2)^DTEY(I)
130 COBTIBDI
CODaB(J)-HCK
100 COITIIDI
RITORI
IBD
STA 395
STA 000
STA 005
STA 010
STA 015
STA 020
STA 025
STA 030
STA 035
STA 000
STA 005
STA 050
STA 055
STA 060
STA 065
STA 070
STA 075
PBASTA 080
STA 085
STA 090
STA 095
STA 500
STA 505
STA 510
STA 515
STA 520
STA 525
STA 530
STA 535
STA 500
STA 505
STA 550
STA 555
STA 560
STA 565
STA 570
STA 575
STA 580
STA 585
STA 590
STA 595
STA 600
STA 605
STA 610
STA 615
STA 620
STA 625
STA 630
STA 635
STA 600
STA 605
STA 650
STA 655
STA 660
STA 665
STA 670
SDBRODTIIE DCOR (ARY,DTRY,TEBPC,CTBP)
BBAL*0 COIP(12,3),KD(3),KR(3),KTB(3),ABY(0),DTRY(0),A(0) ,TBP(3)
REAL^O B03
C DCOR RITORIS A TALDI OF TBI DISTBIBOTIOH COBFPICIIIT -
C DTRY
C 0 - TOTAL lOIIC STREI6TB
C TEBPC - TEHPEHATDEE IB DE6BEES CEBTI6BADE
C ABY - AODEODS PBASE COflPOSITIOB
C T - HOLABITY OP TBP IB DRY, SOLOTE-FREE SOLTEHT
C KD,KTa,KH - BASS ACTIOI EQDILIBBIOB COISTAITS.
C THESE EQDILIBBIOB COBSTAITS BAY BE EXPBESSID AS
C A^BX^CY^DE
C IRIBI X,Y, ABD Z ABE 1ST,2ID ABD 3RD POWERS OF
C lOBIC STBEieTH
C FOR KD(I) - A,B,C,D ARE C0EF(1,I), C0EF(2,I), C0EF(3,I).
C AID COIF(O.I) - D COIPFICIEITS OPTIBIZED
C FOB 1-1 D ACCDRACY
C 1-2 TH ACCDRACY
C 1-3 H ACCDRACI
C FOB KTB(I) - A.B.C.D ARE C0EF(5,I), C0EF(6,I), C0EF(7,I),
C ABD C0EF(8,I) - TB COEFFICIIBTS OPTIBIZID
C FOR 1-1 0 ACCDRACY
C 1-2 TH ACCDRACY
C 1-3 HACCORACY
C FOR KH(I) - A,B,C,D ARE C0EF(9,I| , COEF(10,I), C0EF(11,I),
C ABD C0EF(12,I) - H COEFFICIEBTS OPTIHIZED
C FOR 1-1 0 ACCDRACY
C 1=2 TH ACCDRACY
C 1=3 H ACCDRACY
OCO
DCO
OCO
DCO
OCO
OCO
DCO
DCO
DCO
OCO
DCO
DCO
OCO
OCO
OCO
OCO
OCO
OCO
DCO
OCO
OCO
DCO
DCO
DCO
OCO
OCO
DCO
OCO
OCO
5
10
15
20
25
30
35
00
05
50
55
60
65
70
75
80
85
90
95
100
105
110
115
120
125
130
135
100
105
35
C DCO 150
C DCO 155
DATA COIF/5.16760,12.7385,-6.10939,0.825665 , DCO 160
1 0.895310,-0.258202,0.0628035,0.00105090 , DCO 165
2 0.080975.-0.222063,0.0320680,-0.000857809 , DCO 170
3 0..98513,13.7875,-7.02178, .883583, OCO 175
0 1.52821,-.527579,.0320897,.00016738, DCO 180
5 .538601.-.207956,.0300380,-.00118109, DCO 185
6 2.55963,-7.66961,2.03295,1.58978. DCO 190
7 1,31060,-0.075190,0.0370571,0.00577316, DCO 195
8 0.398105,-0.158715.0.0312360,-0.00160366/ DCO 200
D-3.0»ARY(1)^10.0*ABI(2)^ARY(3) DCO 205
IO3«2.0*ARY(1)^O.0*ARY(2) •ARY(3) DCO 210
T-(CTBP*0.973»1000.0)/26«.3 DCO 215
C DCO 220
C CALCDLATI BASS ACTIOB EODILIBBIDR COISTAITS DCO 225
C OCO 230
DO 50 1-1, 3 OCO 235
KD(I)-COIF (1,1)•COIF(2,1)•D+COIF(3,1)•O*^ 2^C0EF(0,1)•0^« 3 OCO 200
KTR(I) -COEF (5,1) •COIF (6,1) •D^COEF (7,1) •D^»2^C0IF(8,1) *<S**3 OCO 205
KH(1)-COST(9.1)•COIF(10,1)•D^COEF(11,1)•0*^2^C0IF(12,1)•D^«3 DCO 250
C OCO 255
C SOLTE CDBIC EQDATIOI FOR TBP OCO 260
C DCO 265
A(1)--T OCO 270
A(2)-1.0^KR(I)»ARY(3)^BO3 OCO 275
A(3)-2.a^KD(I)^ABY(1) •B03^^2 DCO 280
A(0)-3.0*KTR(I)^AIY(2)*IO3*^0 OCO 285
T-ABII1(-A(1),-A(1)/A(21) DCO 290
IF (Y.LT.O.) Y—A(1) DCO 295
DO 20 J-1,20 DCO 300
J J-J OCO 305
D«A(0) DCO 310
C-A(O) DCO 315
DO 10 11-1,2 DCO 320
K-O-II OCO 325
0»A(K)^T*D OCO 330
C-D^Y«C DCO 335
10 COITIIDE OCO 300
D-A(1)^Y^D DCO 305
DELTY-D/C DCO 350
IF (ABS(DELTY) .LT. 1.E-5.AID.ABS(D).LT.1.E-5) 60 TO 30 DCO 355
Y=Y-DELTY OCO 360
20 COITIIDE DCO 365
30 COITIIDI OCO 370
TBP(I)-Y OCO 375
IF (Y.LT.O.) GO TO 00 OCO 380
60 TO 50 DCO 385
00 COITIIDE DCO 390
PEIIT 70 , TBP(I) ,I,K0(I).KTH(1) ,18(1 ) ,ABI, (C0EP(J,I) ,J-1,12 ) DCO 395
Y=0. DCO 000
50 COITIIDE DCO 005
0TBY(1)-K0(1)^BO3^*2^TBP(1)^^2 DCO 010
DTBY(2)=KTH(2) •BO3^^0*TBP(2)^^3 DCO 015
DTRI(3)-KB(3)^B03^TBP(3) OCO 020
DTRY(0)=0. DCO 025
C DCO 030
C TEBPEBATDRI ADJDSTHEBT OCO 035
C OCO 000
IF (TERPC.E0.25.) 60 TO 60 DCO 005
IT (TEHPC.LT. 15..0E.TERPC.6T. 70.) GO TO 60 DCO 050
XTIHP-1./(273.^TEHPC)-1./298. OCO 055
DTRY(1)-DTRY(1)^EXP(2500.^XTIHP) DCO 060
DTRY (2)-DTBY (2)^IXP( (891. 660^879. SBS^ABY (2) -1070 . 38^A»Y (2) ••2 ) •XTEDCO
065
1HP) OCO 070
0TtlT(3)-DTBY(3)*IXP(-81B.802^XTIBP) DCO 075
60 COBTIBDI DCO O80
C OCO 085
C DCO 090
REIORB DCO 095
C DCO 500
70 FORHAT (//// ' TALDE OF TBP IS '.1PE11.3. ' FOE COBPOBEBI '.12./. * OCO 505
1PSE0D0-EQDIL1BRIDB COBSTAITS <,3E11.3,/, ' AQOIODS PHASE CORPOSITIODCO 510
2B '.0E11.3,/, * PARABETERS *,0E11.3,/,( ' •,1IX.OSII.3)) OCO 515
IBB DCO 520-
SOBBODTIBI PLOTR (KOOIT,ICALPL,HT0ST,S01IB)
C
C THIS HODTIBE ISEHERATES CALCOHP PLOTS OF TWO KIBDS:
C ICALPL = 1 COBCEITRATIOB TS TIHE IBTEBTAL
C ICALPL - 2 PEBCEBT OF PIED TS TIHE IBTEBTAL
C ICIILPL = 3 PLOT BOTH
C
C KODIT - TOTAL HDHBE8 OF TlHI IBTERTALS
C HIOST - TOTAL HDHBES OF STAGES
C SOLII - TOTAL ARODIT IB PEED
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
5
10
15
20
25
30
35
00
05
50
36
A - AQOIODS COICIITRATIOIS
0 - OBGABIC COICIITIATIOIS
AODT - AODEODS FLOI OOT OP STAGE
OODT - ORGAIIC FLOI ODT OF STAOE
DIBEBSIOB A(0,100), 0(0,100) , AOnT(IOO), OoaT(IOO), 1(1500)
10,3) , SOLIB (100), 1C0R(3), 8DFFER(0000)
REAL«8 YLAB(8)/'C0ICEIIB',*ATI0I ( *,*0,TB IB ','0/L-, HBO',
1 )','«','PIBCII T ','OF FIIDS*/
BIAL^B XIAB(2)/'TIBE IBT', ' EBTALt'/
REAL»0 C0B(3)/'D * , * TR ',' H •/,LL(2)/'LIB ','LOG ' /
L0GICAL»1 AQ0/*A*/,OIG/*O'/,IAO
LOGICAL^I DESCB(3a)
LOGICAL*1 DISC 1(30)/'D','I','A','B','I', * 0*.*H*,*»*,*T*,'H'
1','I','0','R*,*t*.*B*,*I*,'T','B*,*I*,'C', ' ','A',*C*,*I', '
IF ICALPL - 3 , PLOT BOTB COICIITRATIOI TS TIBE ABD
PERCEBT OF FEED TS THE
IC-0
IF (ICALPL. EQ. 3) IC-3
IF (ICALPL.EO.3) ICALPL-1
CALL QQBIIS (0)
IRITI (9,190)
IBITI (9,200)
IPL-0
BIAD SIT OP PLOT SIIICBES
BSTB - STIIAB BDBBEB TO BE PLOTTED
BCOR - BDRBII OF CORPOBEBTS TO BE PLOTTED
lAO - A PLOT AODEODS STREAH
0 PLOT ORGAIIC STREAH
ILL - 0 LIIIAI X-AXIS, IIBEAR Y-AXIS
1 LIBIAR I-AIIS , LOG Y-AXIS
BAIT - 0 PLOT BTIBT POUT
1 PLOT ITER! OTHER POUT
ILAST - 0 LAST PLOT
1 OTHER PLOTS FOLLOI
BEAD IIDICES OF COBPOIEITS TO BE PLOTTED
1 - DBABIDB
2 - TBOBIDB
3 - BITBIC ACID
10 COITIIDE
READ (5,210) ISTR,ICOB,IAO,ILL,IALT,ILAST
READ (5,220) (ICOB (I) ,1-1,ICOB)
20 COBTIBDE
IPL-IPL^I
SET OP POIBTS TO BE PLOTTED
YRAX-0.
YBII-1.E10
EEIIID 11
I-O
DO 50 lK-1,KODIT
BEAD COICEITBATIOBS (A ABD 0) AID FLOIS (AODT ABD OOOT)
DO 30 J-1,BTOST
BEAD (11) ITB, (A(K,J),K-1,0), (0(K,J).K<1,0),AODT(J),OODT(J)
30 COBTIBDE
IF (lALT.EQ.1.AID.IK/2^2.EQ.IK) GO TO 50
1-1^1
DO 00 K-1,IC0H
IF (ICALPL. ig . LAID. lAO.IQ.AQD) Y (I ,K)-A (ICOB (K) , BSTB)
IF (ICALPL. ig.l . ABD. IA0.ta.ORG) Y (I,K)- 0 (ICOB (K) ,BSTR)
IF (ICALPL. IQ. 2. AID. IAO.EQ.AQD) T (I , K)-A (ICOB (K) , ISTR) •AODT
1 SOLII (ICOB (K))^100.
IF (ICALPL.EO. 2 . ABD.BAO.E0.ORG) Y (I,K)- 0 (ICOB (K) , BSTB) •OODT
1SOLIB(ICOR(K))•100.
YRAX-AHAX1(Y(I,K) ,YRAX)
IF (BLL.I0.1.ABD.Y(I,K).LE.O.) GO TO 00
YBII-ARII1 (Y(1,K),YHII)
00 COITIBDE
X(I)-I
50 COBTIBDE
YBAX=IHAX^YHAX^.05
XRII=0.
XBAX-KCDBT
CALL IBTL (•BECB',BDFFER,0000,'S.B. lATSOBS')
CALL QOCCPA (3,XRIB)
CALL QQQCPA (0,IBAX)
CALL QQQCPA (13,YBIB)
CALL QOCCPA (10,YBAX)
IF (KO0IT.LE.100) CALL SCALE (1 , 10. , 0. ,'LIES')
IF (KODIT.61.100) CALL SCALE (1,KOOBT/10.,0.,'LIKE')
IF (BLI.ECO) CALL SCALE (-1, 10. ,0. ,'LIIE')
IF (BLI..EC.1) CALL SCALE (-1. 10. . 0. ,'L062')
CALL AXIS (1,-15,.15)
PLO
PLO
PLO
PLO
PLO
, Y(150PLa
PLO
'3 II RPLO
PLO
PLO
PLO
PLO
PLO
,*0*,*BPLO
D*,*$*/PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
(BSTB)/PLO
PLO
(BSTR)/PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
PLO
55
60
65
70
75
80
85
90
95
100
105
110
115
120
125
130
135
100
105
150
155
160
165
170
175
180
185
190
195
200
205
210
215
220
225
230
235
200
205
250
255
260
265
270
275
280
285
290
295
300
305
310
315
320
325
330
335
300
305
350
355
360
365
370
375
380
385
390
395
000
005
010
015
0 20
025
030
035
000
005
050
055
060
065
070
075
080
085
090
095
500
505
510
515
520
525
37
IF (ILL.EQ.O) CALL AXIS (-1,.15,-.15) PLO 530
IF (ILL.SQ.1) CALL AXIS (-1.. 1.-. 15) PLO 535
CALL BOX (3) PLO 500
IF (ILL.EQ.O) CALL BOX (0,1.5,-1.5) PLO 505
IF (ILL.SQ.1) CALL BOX (0.-1.5.1.5) PLO 550
IF (ILL.EQ.O) CALL GBID (.35,1.5,-1.5) PLO 555
IF (BLL.IQ.O) CALL GRID (-.05,1.5,-1.5) PLO 560
IF (ILL.Ea.1) CALL GRID (.35,-1.5,1.5) PLO 565
IF (ILL.EQ.1) CALL GRID (-.05,-1.5,1.5) PLO 570
C PLO 575
C LABEL AXES PLO 580
C PLO 585
IF IICALPUEQ. 1) CALL TITLE (0. 0,*EIIT STREAB COBPOSITIOH TS TIHESPLO 590
1*,XI.AB(1),YLAB(1)) PLO 595
ir (ICALPL.EO.2) CALL TITLE (0.0,'EXIT STREAB COBPOSITIOB TS TIBESPLO 600
1*,XLAB(1),YLAB(7)) PLO 605
C PLO 610
C PLOT GRAPHS PLO 615
C PLO 620
DO 80 K-1,BC0H PLO 625
IS-ICOB(K) PLO 630
lB-1 PLO 635
BPTS-KODIT PLO 600
IP (ILL.EQ.O) 60 TO 70 PLO 60S
DO 60 J-1,KODIT PLO 650
IB-J PLO 655
IF (Y(J,K).6T.O.) 60 TO 70 PLO 660
BPTB-BPTS-1 PLO 665
60 COBTIBDE PLO 670
70 COITIIDE PLO 675
IF (IB.EQ.KODIT) 60 TO 80 PLO 680
CALL LIIPLT (X(IB),Y(IB,K),IPTS,IS,-1) PLO 685
80 COITIIDE PLO 690
C PLO 695
C SET DP LEGEID PLO 700
C PLO 705
lP-1 PLO 710
DO 1150 K-1,BC0a PLO 715
IP (ICCH(K).EO. 1) GO TO 90 PLO 720
IF (IC0H(K).EQ.2) 60 TO 110 PLO 725
IF (ICCR(K).Eg.3) GO TO 130 FLO 730
GO TO 150 PLO 735
90 COITIBDE PLO 700
DO 1100 1-1,8 PLO 705
DESCR (IP)-DESCI(l) PLO 750
IP-IP^1 PLO 755
100 COBTIIOE PLO 760
GO TO 150 PLO 765
110 COBTIBDE PLO 770
DO 120 1=9,16 PLO 775
DISCR(IP)=DISC1(I) PLO 780
IP-1P^1 PLO 785
120 COBTIBDI PLO 790
GO TO 150 PLO 795
130 COBTIBDI PLO 800
DO 100 1-17,28 PLO 805
DESCR(IP)-DISCI(1) PLO 810
IP-1P^1 PLO 815
100 COBTIBDI PLO 820
150 COBTIBDI PLO 825
IF (BLL.EO.O) CALL LEGEID (*COHPOBEITS«',ICOB,ICOB,DESCR.2,1.5,-1.PLO 830
15) PLO 835
IF (ILl.tg.1) CALL LEGEID ('COHPOIEBTSS*,ICOH,HCOn,DESCR,2,-1.5,1.PLO 800
15) PLO 805
C PLO 850
WBITE (9,230) IPL,HSTR,BAO,LL(BLL^I),(COH(ICOB(K)),K=1,BCOR) PLO 855
C PLO 860
IF (IC.EQ. 3. ABD.ICALPL. Sg. 1) GO TO 160 PLO 865
IF (BLAST.Sg.O) GO TO 180 PLO 870
IF (IC.Ig.3.ABD.ICALPL.I0.2) GO TO 170 PLO 875
IT (BLAST.BI.O) 00 TO 10 PLO 880
160 COBTIBDI PLO 885
ICALPL=2 PLO 890
GO TO 20 PLO 895
170 COBTIBDI PLO 900
ICALPL-1 PLO 905
60 TO 10 PLO 910
180 COBTIBDI PLO 915
CALL ADTAIS PLO 920
C PLO 925
C PLO 930
RITDRI PLO 935
C PLO 900
190 PORHAT (' ',///, ' ',21X,'PL0T SDBHARY'/) PLO 905
200 FORHAT (' ','PLOT*,0X,*STA6E',0X,*AQDEODS OR*,0X,*LIIEAR OH'.OI.'CPLO 950
10HP0IEITS*/* *,1X,'I0.'.SX.'IO . ',6X,'0RGAIIC?',7X,'L06?' ,BX,'PLOTTPLO 955
2ID'/) PLO 960
210 FOBHAT (212,IX,AI,312) PLO 965
220 FORHAT (1012) PLO 970
230 FORRAT (' ',I2,7X,I2,10X,AI,12X,A0,6X,3A0) PLO 975
BID PLO 980