SEPHIS Operating Manual

^.
^s^'f ORNL/CSD/TM-70
User's Guide to the SEPHIS Computer Code

for Calculating the Thorex Solvent
Extraction System
S. B. Watson
R. H. Rainey
WSTBl
Any further distributi r of the data
t l W ^ i ^ t D t h i r d parties repreS govern-
mentsTfoWl^pmpanies and 1 ons
of U.S. c o m p a m l ^ ^ ^ u l d be coordinated! irector, Divisior
Reactor Researcfi a n ^ l f c h n o l o g y , Departmenf
^-fVigours
Printed in the United States of America. Available from
the Department of Energy,
Technical Information Center
P.O. Box 62, Oak Ridge, Tennessee 37830
Price: Printed Copy $5.25: Microfiche $3.00
This report was prepared as an account of work sponsored by an agency of the United
States Government. Neither the United States Government nor any agency thereof, nor
any of their employees, contractors, subcontractors, or their employees, makes any
warranty, express or implied, nor assumes any legal liability or responsibility for any
third party's use or the results of such use of any information, apparatus, product or
process disclosed in this report, nor represents that its use by such third party would
not infringe pnvately owned rights.
i
DISCLAIMER
This report was prepared as an account of work sponsored by an

agency of the United States Government. Neither the United States
Government nor any agency Thereof, nor any of their employees,
makes any warranty, express or implied, or assumes any legal
liability or responsibility for the accuracy, completeness, or
usefulness of any information, apparatus, product, or process
disclosed, or represents that its use would not infringe privately
owned rights. Reference herein to any specific commercial product,
process, or service by trade name, trademark, manufacturer, or
otherwise does not necessarily constitute or imply its endorsement,
recommendation, or favoring by the United States Government or any
agency thereof. The views and opinions of authors expressed herein
do not necessarily state or reflect those of the United States
Government or any agency thereof.
DISCLAIMER
Portions of this document may be illegible in

electronic image products. Images are produced
from the best available original document.
[ This report contains information of a preliminary nature and
i was prepared primarily for internal use at the ongmatmg
installation It a subject to revision or correction and therefore
does not represent afinalreport It is passed to the reapient in
I confidence and should not be abstracted or further disclosed
i without the approval of the ongmatmg mstallation or USDOE ORNL/CSD/TM-70
I Technical Information Center, Oak Ridge, TN 37830
D i s t . Category UC-79c
Contract No. W-7405-.eng-26

Hiis report was prepared as an account of work
sponsored by the United States Government Neither the
United States nor the United States Department of
Energy, nor any of their employees, nor any of their
COMPUTER SCIENCES DIVISION contractors, subcontractors, or their employees, makes
any warranty express or implied, or assumes any legal
liability or responsibility for the accuracy, completeness
or usefulness of any information, apparatus, product or
process disclosed, or represents that its use would not
infrmge pnvately owned nghts
USER'S GUIDE TO THE SEPHIS COMPUTER CODE

FOR CALCULATING THE THOREX SOLVENT
EXTRACTION SYSTEM
S. B. Watson and R. H. Rainey*
(Sponsor: R. H. Rainey; O r i g i n a t o r : S. B. Watson)

This document is
PUBLICLY RELEASABLE
Date Published - May 1979
Âuthorizing Official
Date_„J^^JtQf
NOTICE
This document contains information of a preliminary nature.
I t is subject to r e v i s i o n or correction and therefore does
not represent a f i n a l r e p o r t .
This work was funded by

HTGR Fuel Recycle Program
(189a 0H045) Reprocessing Development.
(now part of the Consolidated Fuel Reprocessing Program)
UNION CARBIDE CORPORATION, NUCLEAR DIVISION

Operating the
Oak Ridge Gaseous D i f f u s i o n Plant Oak Ridge National Laboratory
Oak Ridge Y-12 Plant Paducah Gaseous D i f f u s i o n Plant
f o r the
DEPARTMENT OF ENERGY
*Oak Ridge Associated U n i v e r s i t i e s

CONTENTS
ABSTRACT 1
I. INTRODUCTION . 1
II. PROGRAM DESCRIPTION 3
III. USER INFORMATION 8
Input 8
Output 16
REFERENCES 19
APPENDICES
A. A LISTING OF SEPHIS-THOREX ROUTINES 21
B. SAMPLE INPUT AND CORRESPONDING COMPUTER OUTPUT 39
iii
USER'S GUIDE TO THE SEPHIS COMPUTER CODE
FOR CALCULATING THE THOREX SOLVENT
EXTRACTION SYSTEM
S. B. Watson and R. H. Rainey*
ABSTRACT
The SEPHIS computer program was developed to

simulate the countercurrent solvent extraction process.
The code has now been adapted to model the Acid Thorex
flow sheet. This report represents a practical user's
guide to SEPHIS-Thorex containing a program description,
user information, program listing, and sample input and
output.
I. INTRODUCTION
SEPHIS, a computer program for Solvent Extraction P^rocesses Having
interacting Solutes, was developed by W. S. Groenier^ for the purpose of
studying the application of a ^5% TBP solvent extraction system for
recovered irradiated LMFBR fuel. G. L. Richardson of HEDL has adapted
the SEPHIS program to the conventional Purex process which uses 30% TBP
and has made valuable changes in program structure.^'^ The authors
have made additional modifications to the Richardson version of SEPHIS-
Purex to make the code more efficient and adaptable for Oak Ridge National
Laboratory (ORNL) uses.'* Research to improve SEPHIS-Purex is being
continued by A. D. Mitchell who has written a code entitled SEPHIS-M0D4.^»6
This improved model of the Purex solvent extraction system is better
able to simulate run conditions during transient periods than previous
versions of SEPHIS and also better represents the operation of mixer-
settlers.
*Chemical Technology Division
1
2
A modification of the SEPHIS program has now been made for the
simulation of the Acid Thorex flowsheet,'' and this report represents a
documentation of the SEPHIS-Thorex computer code for the user. It
should be noted that a version of SEPHIS-M0D4 for Thorex solvent extrac-
tion is in preparation.^
The SEPHIS program simulates the countercurrent multi-stage ex-
traction of three interacting solutes. In the Purex version of the
program, these solutes are uranium, plutonium IV and nitric acid. In
the Thorex version thorium is substituted for the plutonium.
The SEPHIS computer program is useful for obtaining estimates of the

concentration of components throughout a solvent extraction contactor
system under either transient or equilibrium conditions or in estimating
the time required to reach steady-state operation. The program is based
on an empirical fit to data which had considerable scatter. Comparison
of calculated data with data from a limited number of countercurrent
batch experiments indicate good agreement in the middle range of concen-
trations. In the low concentration range, the calculated data do not fit
experimental data. At low concentrations of metal ions, trace quantities
of impurities strongly influence extraction characteristics. These effects
are not included in the program. In spite of a considerable effort, the
calculations for nitric acid are less accurate than for the metal ions.
Experience has shown that analyses of free nitric acid have been erratic
in solvent extraction systems of the type simulated by this program.
Use of the computer program requires an input of the number of theoretical
stages in the contacted system. In practical operating systems, the
3
number of theoretical stages varies with operating conditions and, there-

fore, is usually not accurately known. The program does not simulate the
extraction of fission products.
Many of the technical aspects of the SEPHIS-Thorex computations

are discussed elsewhere.'' These aspects include calculation of the
distribution coefficients and partial molar volumes, explanation of uses
and limitations of the code, and comparisons between SEPHIS-Thorex calcu-
lations and experimental data. In addition, the study of temperature
effect on the extraction of thorium and nitric acid and the formation of
third phase are discussed in another publication.^ Thus, this report
includes only a general program description and user information nec-
essary to execute the program. Also included are a program listing and
sample input and output.
The Thorex version of the code would not have been possible without
the valuable guidance of W. S. Groenier. This work was funded by the
High-Temperature Gas-Cooled Reactor (HTGR) program.
II. PROGRAM DESCRIPTION
The description of the flow sheet and the contactor model for the
Thorex version parallel those of Groenier.^ in addition the calculational
procedure used for the stage computations and used to determine the distri-
bution coefficient are patterned after the technique of Groenier and
discussed in Ref. 7.
In general the calculation simulates the extraction of the three
interacting solutes by a cascade of mixer-settler units with 100%
stage efficiency. Flow between stages is batchwise. This flow can either
4
be proportional to the feedstream flow rates or can be varied to simulate

internal recycle or to increase stage holdup. The number of stages, the
number of incoming and outgoing streams, the flow rate or volume of each
stream, the flow patterns through the system, and the composition of
each feedstream can be designated. Also, the composition of each phase
in the contactor at the beginning of the calculation can be stipulated,
that is, the initial profile of the system can be specified.
The simulation starts by computing the concentration of the solutes

for each phase in each stage based on the distribution coefficients which
are calculated as part of the program. These conditions are then used
as the basis of a second calculation, and this iteration procedure is
continued until the amount of material entering and leaving the system
reaches a pre-determined material balance - normally 99.9%. This system
gives an estimate of the composition of each phase in each stage of the
system.
The program SEPHIS-Thorex consists of six routines that are
programmed in a well-structured manner and are heavily augumented
with program comment cards. The basic methodology and mathematical
formulations have been described previously;^•'^ thus, we include in this
section only the diagram of program structure and a description of each
routine.
Figure 1 shows the overall structure of the SEPHIS-Thorex program.
Arrows indicate the calling sequence of routines. All routines are
written in the IBM System/360 and System/370 FORTRAN IV language. A
listing of SEPHIS-Thorex routines is given in Appendix A.

5
MAIN
'
r 'r f f
PROFIL FLOWS STAGES PLOTR
r
1 UCOR 1
Fig. 1. Overall Structure of SEPHIS-Thorex.
The MAIN routine reads switches which determine program flow, directs
the calling of subroutines according to these switches, determines steady
state, and prints requested output. The routine functions as follows;
1. Read title information.

2. Read program switches (see Section III, Table 1 ) .
3. Zero arrays and constants.
4. Call PROFIL to set initial concentrations for time zero and
define feed streams.
5. Write feed stream data.
6. Convert feed streams concentrations to molarity.
7. Call FLOWS to set interstage flow rates, define outgoing streams

and flow ratios and calculate necessary recycle flows.
8. Print outgoing streams.
9. Print initial profile in the stages of the contactor system.
10. Print phase ratios and stage flows.
11. Initialize total molarity coming in (SOLIN) and total molarity

coming out (SOLOUT) of system.
12. Call STAGES to perform stage calculation.

6
13. Calculate percent of material balance which has been reached.

14. Convert concentration back to g/1.
15. Print concentrations for the given time interval.
16. Check for convergence. If not converged, return to 12.
17. Print total moles of each component in and out.
18. Call PLOTR for Calcomp plots, if requested.

19. Punched card output, if requested.
20. Read next case information. Go to 1., 3. or stop depending on
switches read.
The routine PROFIL defines the feed streams and the initial profile.
PROFIL proceeds as follows:

1. Read aqueous and organic feed concentrations and stream flow rates.
2. Read nonzero initial profile, if requested. Otherwise zero

profile is assumed.
The manner in which 1 and 2 are executed depends on switches read by the
MAIN routine.
FLOWS defines outgoing streams, unusual routing patterns, interstage
flows, phase ratios, flow ratios, recycle volumes and ratios, and stage
volumes. FLOWS proceeds thusly:
1. Define outgoing streams, if there are any.
2. Set aqueous and organic interstage flows. Determine unusual

routing pattern for the flows, if there is one.
3. Define phase ratios.
4. Calculate necessary recycle flows.
5. Determine stage volumes.
The above listing is merely a general diagram of the flow of this routine.
The manner in which 1-5 are executed (or possibly omitted) depends on
certain switches read by MAIN.
7
STAGES performs stage calculations for the current time interval.
It considers inputs from adjacent stages, input streams, and recycle
within the stage. It executes in the following manner.

1. Compute aqueous and organic concentrations coming into the stage
for the current time interval.
2. Zero iteration counter.
3. Calculate distribution coefficients for the current aqueous phase
composition.
4. Using information calculated in 1 and 3, determine new aqueous
phase concentration of each solute.
5. Compare new aqueous phase concentration with previous aqueous phase

concentrations for each solute.
6. If new and previous values are "close" by a fixed criterion for
each solute, proceed to 7. If not, increment the iteration counter
and return to 4 with current aqueous phase composition for each
solute redefined as the calculated value for this iteration.
7. Calculate the organic phase stage concentrations for the current

time interval.
Steps 1-7 are performed in sequence for each stage in the contactor system.
UCOR calculates distribution coefficients for a given aqueous phase
composition at the stage equilibrium temperature. UCOR embodies the
equations for the mass action equilibrium constants which have been
deterTTiined by the least squares fits described in Ref. 7. In addition,
UCOR contains the temperature correction equations published in Ref. 9.
PLOTR may be used to plot concentration versus time and/or percent

of feed versus time. A set of plot switches read by this routine allows
the user to determine the stage number, the component (either aqueous or
organic), the type of graph, and the number of points to be plotted. A
list of these switches is given in Section III.
8
III. USER INFORMATION
SEPHIS-Thorex is written in the IBM System/360 and System/370

FORTRAN IV language. The program has been compiled using the H-level
compiler and has been executed on both the IBM 360/91 and 360/75 at
ORNL. The program requires about 21 OK (K = 1024 bytes; byte = 8 bits)
of core storage for execution. Approximately 1/20 second of computing
time on the IBM 360/91 is required for each time interval when calculating
10 stages.
The remainder of this section describes the content of the card
input which must be prepared in order to execute the program and the
output produced by the program. Sample input and the corresponding
output are presented in Appendix B. The example is an Acid Thorex Flow
sheet for extraction of uranium and thorium from solutions of Consolidated
Edison fuel.^°
Input
The input data are to be prepared in the order shown in Table 1.
Referring to Table 1, "Subroutine" indicates the subroutine in which
the card is read. "Card Number" denotes the order within the group
(Program Cards, Feed Stream Cards, etc.) or quantity of cards for each
group to be input. "Columns" refers to the actual columns on the card
in which the data must be punched. "Remarks, Variables, etc." gives
an explanation as to the variable name, definition, and form of the
data to be punched on the card. In the FORMAT descriptor the letter "A"
denotes alphabetic data; "I" refers to integer data (no decimal point)
which must always be right-justified in the field; "F" indicates a real

9
ORNL-DWG. 78-124I9R
Table 1. Input to SEPHIS Code
Subroutine Card Number Columns Remarks, Variables, etc.
Program Cards
MAIN 1 FORMAT (12)
1-2. NTTL - number of title cards
2 to NTTL+1 F0RMAT (10A8)
1-80 TITLE - problem title

There are NTTL cards
of this form.
NTTL+2 F0RMAT (F8.2, 1012)
1-8 CTBP - volume fraction of dry TBP
9-10 NTf5ST - total number of stages -

must not exceed 25
11-12 IPR0 - 0 for a zero initial concen-

tration profile
1 for a non-zero initial
profile
13-14 IRATI0 - 0 if phase ratio = flow

ratio
1 if phase ratio values are
to be input
15-16 ICALC - 0 for shutdown calculations

(no U or Pu in feed
streams)
1 for start-up or interrup-
tion
17-18 I0UT -0 if there are no extra out

going streams
1 if there are outgoing
streams in addition to
end streams
10
ORNL-DWG. 78-12420
Table 1. Input to SEPHIS Code (cont'd)
MAIN 19-20 IV0L - 0 stage volumes are equal

(cont'd) 1 stage volumes are unequal
but according to total
flow
2 stage volumes are inde-
pendently unequal
3 stage volumes are equal
and specified
(Options 2 or 3 should be used for
HAN* calculations to correct for stage
residence time)(*Hydroxylamine
Nitrate - Ref. 3)
21-22 IPR - 0 print out all time inter-
vals
1 print only last interval
-1 print every 20th time
interval
23-24 IPNCH - 0 no punched card output
desired
1 final profile output
punched on cards
25-26 ICALPL - 0 no Calcomp plots
1 plot concentration vs
time
2 plot percent of feed vs
time
3 plot both 1 and 2
If plots are desired, i.e., ICALPL > 0
IPR must be set equal to 0.
27-28 NSTR - 0 no unusual routine pattern
1 routing pattern other than
normal
Feed Stream Cards
PR0FIL One card for F55RMAT ( 2 1 2 , 5 F 8 . 0 , 1 2 . F 8 . 0 )
each input
stream
1-2 stage number that feed enters
11
ORNL-DWG.78-12421
PR0FIL 3-4 JHAS - 1 for aqueous feed

0 for organic feed
5-12 FORI - feed flow rate (volume units)
13-20 C0N1 - uranium content (g/1)
21-28 C0N2 - plutonium content (g/1)
29-36 C0N3 - acid molarity
37-44 TEMP - temperature in degrees C
45-46 INDEX - 1 more cards to follow
0 last card
47-54 C0N4 - nitrate normality fron inex-
tractable salts, including
the reductant
(Enter Pu reductant normality as
negative Pu concentration^
Initial Profile Cards
(required only for non-zero initial profile)
PR0FIL One card for F0RMAT (8F8.0)
each stage
1-8 X(1,J,1) - aqueous-phase uranium con-
centration in stage J (g/1)
9-16 X(2,J,1) - aqueous-phase plutonium
concentration in stage
J (g/1)
17-24 X(3,J,1) - aqueous-phase nitric acid
J (M)
25-32 Y(1,J,1) - organic-phase uranium con-
centration in stage
J (g/1)
12
ORNL-DWG. 78-12422
PR0FIL 33-40 Y(2,J,1) - organic-phase plutonium

(cont'd) concentration in stage
J (g/1)
41-48 Y(3,J,1) - organic-phase nitric acid
J (M)
49-56 TPR0F(J,1) - temperature in stage J
57-64 X(4,J,1) - aqueous-phase nitrate con-

centration in stage J (M)
Outgoing Stream Data

(other than aqueous raffinate and organic product)
FL0WS One card for F0RMAT (212, F8.0, 12)

each out-
going stream 1-2 I - stage number that stream
leaves
3-4 JHAS - 1 for aqueous stream

0 for organic stream
5-12 0TRT - exit flow rate (volume

units) - total volume
must exit
13-14 INDEX - 1 more cards follow

0 last card
Unusua I Routing Data

[NSTR = T) •
FL0WS 1 F0RMAT (215)
1-5 ISTR - organic stream exits this

stage and feeds stage JSTR
6-10 JSTR - organic stream from stage

ISTR feeds this stage
13
ORNL-DWG.78-12423
• - -
Phase Ratio Cards

(used only if phase ratio f flow ratio)
FL0WS One card for F0RMAT (12, F8.0)
each region
of constant
phase ratio
1-2 I - number of stages having phase
ratio shown
3-12 RATI0 - aqueous/organic phase ratio

stage Volumes
(for IV0L = 2 and IV0L = 3)
FLOWS For IV0L = 2 stage volumes are unequal
and not flow related
1 F0RMAT (10A8)
1-8 STV0L (1)
72-80 STV0L (10)

2 F0RMAT (10A8)
1-8 STV0L (11)
72-80 STV0L (20)

3 F0RMAT (10A8)
1-8 STV0L (21)
32-40 STV0L (25)
Only NT0ST ;jtage volumes lead to be entered; that is, 1 to 3 cards are needed
for input.
14
ORNL-DWG 78-12424R
FL0WS For IV0L = 3 stage volumes are equal and

specified
F0RMAT (F8.0)
1-8 V0L - stage volume for all stages
Graph Cards
(ICALPL = 1, 2, 3)
For each plot include cards of the following form:

PL0TR 1 F0RMAT (212, IX, Al, 312)
1-2 NSTR - number of outgoing stream to be

plotted
3-4 NC0M - number of components to be

plotted (£ 3)
6 NA0 - A aqueous stream

0 organic stream
7-8 NLL - 0 linear plot on both axes

1 linear abscissa axis, log
ordinate axis
9-10 NALT - 0 plot every point

1 plot every other point
n-12 NLAST - 0 last plot
1 more plot"; to come
F0RMAT (1012)
st
1-2 IC0M(i; 1 component to be plotted
(1-uranium, 2-thorium, or
3-acid)
3-4 IC0M(2) - Z"'' component to be plotted
(1, 2, or 3)
5-6 IC0M(3) - 3'"'^ component to be plotted
(1, 2, or 3)
15
ORNL-DWG. 78-12425
For ICALPL- 3 include only one set of cards for each group of two plots that is
generated.
Next Case
MAIN 1 F0RMAT (212)
1-2 ICHNGE - 0 next case is new
1 new Feed Stream Cards and
Outgoing Stream Data required
for each feed and exit stream-
no Program Cards, Initial Pro-
file Cards, Phase Ratio Cards,
or Stage Volume Cards are
required
2 stop
3-4 IDIF - 0 no changes in Program Card
parameters
1 changes in at least one of
these Program Card parameters -
ICALC, I0UT, IPR, IPNCH,
ICALPL (these parameters are
explained in the section
Program Cards)
If IDIF = 1, input a card of the foMowing form; if IDIF = 0, omit this card.
2 F0RMAT (512)
1-2 ICALC
3-4 I0UT
5-6 IPR
7-8 IPNCH
9-10 ICALPL
See Program Cards for explanation of
these parameters .
16
number with or without signed exponent (decimal point should be punched).

The variable names are those used in the FORTRAN program. An example
of sample input is included in Appendix B.
Output
An example of program printout is given in Appendix B. A general
explanation of this example and hence all output are herein described.
The first page of printed output provides a synopsis of input
and various initial conditions. Included are the feed stream data,
initial profile, phase ratios, and stage flows. The remainder of the
printout consists of transient behavior results from each time interval.
Printout for each time interval includes aqueous- and organic-phase
concentrations, aqueous and organic flows, temperature, and number of
iterations required for convergence in the stage calculations; all of
which are printed for each stage. It should be noted that a "#" printed
to the right of an organic thorium concentration indicates third phase
within accuracy of about 10%. Experimental evidence indicate that third-
phase formation did not occur during this study,^° but the SEPHIS program
indicated that the third phase occurred starting at time interval 15.
This discrepancy can be due either to the fact that the third phase is
extremely difficult to detect experimentally or the mathematical represen-
tation of the third-phase surface in SEPHIS is only accurate to within
10%.
After every time interval which is a multiple of 100 and upon
convergence of the transient calculations, the percent steady state and
the time intervals at which 90, 95, 99, and 99.5% steady state have
been reached and printed.

17
If desired, output may also be obtained in the form of Calcomp
plots. Graphs for concentration versus time interval and/or percent of
feed versus time may be produced by (1) requesting plots on the "Program
Cards" and (2) entering the correct input as described in "Graph Cards"
(see Input).
19
REFERENCES
1. W. S. Groenier, CaZoulation of Transient Behavior of a Dilute-
Purex Solvent Extraction Process Having Application to the
Reprocessing of LMFBR FuelSy ORNL-4746, April 1972.
2. G. L. Richardson, Effect of High Solvent Irradiation Exposures
on TBP Processing of Spent LMFBR Fuels^ Hanford Engineering
Development Laboratory, HEDL-TME 73-51, June 1973.
3. G. L. Richardson and J. L. Swanson, Plutonium Pax'titioning in the
Purex Process with Hydrazine-Stabilized Hydroccylamine Nitrate,
Hanford Engineering Development Laboratory, HEDL-TME 75-31,
June 1975.
4. S, B. Watson, and R. H. Rainey, Modifications of the SEPHIS
Computer Code for Calculating the Purex Solvent Extraction
System, ORNL-TM-5123, December 1975.
5. A. D. Mitchellji A Comparison Between SEPHIS-M0D4 and Previous
Models of the Purex Solvent Extraction Model, ORNL/TM-6565,
February, 1979.
6. A. D. Mitchell, SEPHIS-M0D4: A User's Manual to a Revised Model
of the Purex Solvent Extraction System, ORNL-5471, (in press).
7. R. H. Rainey and S. B. Watson, interim Report of the Modifications
of the SEPHIS Computer Code for Simulating the Acid Thorex Solvent
Extraction System, ORNL/TM (in press).
8. A. D. Mitchell, Modification of the SEPHIS-M0D4 Computer Program
to Simulate the Thorex Solvent Extraction Process, ORNL/TM-6825,
(in preparation).
20
REFERENCES (Cont'd)
9. S. B. Watson and R. H. Rainey, Modeling the Effect of Temperature

on Thorium and Nitric Acid Extraction and the Formation of Third
Phase for Modification of the SEPHIS — Thorex Computer Program,
Oak Ridge National Laboratory, ORNL/CSD/TM-69, (in press).
10, R. H. Rainey and J. G. Moore, "Laboratory Development of the Acid
Thorex Process for Recovery of Thorium Reactor Fuel," Nuclear
Science and Engineering, JO^, 367—71 (1961).
^
APPENDIX A
A LISTING OF SEPHIS-THOREX ROUTINES
23
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|, S 0 I . I B ( 3 | , SOLOB 60
1 T ( 3 ) , AQOOT(IOO), ORGOOT(IOO) 65
DIBEBSIOB AA(IOO). T T ( 1 0 0 ) , 0 0 ( 1 0 0 ) , RATIO (100) 70
COHBOB / P I B S / 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/ 1 P B O , I ( « , 1 0 0 , 2 ) , T ( « , 1 0 0 , 2 ) , I P D ( « , 1 0 0 ) , T P D ( I 1 , 1 0 0 ) ,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 5 5 0 , ( T I T L E ( 1 T T 1 , I ) , I « 1 , 1 0 ) 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 5 8 0 , ( T I T L E ( I T T L , I ) , 1 - 1 , 1 0 ) 580
70 COBTIBOE 585
BRITE ( 6 , 5 9 0 ) 590
DO 80 I«1,BTOST 595
IF (AFDRT(I).LE.O.) GO TO SO 600
PRIBT 6 0 0 , I , X P D ( 1 , I ) , X F D ( 2 , I ) , X F D ( 3 , I ) , X F D ( 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 ) - I P O ( 2 , 1 ) / 2 3 2 . 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,PHATI0 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 , 7 5 0 ) 1265
C 1270
C H - I S 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 I P 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 I F ( J F D ( I , J ) . L T . 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 C O B T I B O E 11105
CALL S T A G E S (SPH) 1050
DO 2 8 0 J'l.BTOST 1055
TPB0r(J,1)-TPS0r(J,2) 1060
DO 2 8 0 I>1,0 1*65
I(I,J,1)<I(I,J,2) 1070
T(I,J,1)-T(I.J,2| 1075
280 C O I T I B O E 1080
C 1085
C ir BOBBEB o r TIBE I B T E B T A L S IS LESS THAB T H E BDBBEB O E S T A G E S , 1090
C O B I T C A L C O L A T I O B O r S O L I B . S O L O D T , A B D P E B C E B T STEADT S T A T E 1095
C 1500
IP (ICEB.,EQ.0.AB0.K.LT.BTOST| G O TO 390 1505
C 1510
C ir SBDTDOBI ( I C A L C ' O ) , CRECB rOB BEAB ZEBO C O B C E B T B A T I Q B 1515
C 1520
ir (ICALC.EC. 1) 6 0 TO 300 1525
C 1530
C rOE S H D T D O B B , I C A L C - 0 , T E S T EOB C O B P O B E B T S 1 ABD 2 » 0, HOT C O B P O B E B I 1535
C 1500
DO 2 9 0 J - 1 , B T 0 S T 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 C O B T I B O E 1570
00 TO 380 1575
C 1580
C C A L C D L A T I O B O F P E B C E B T STEADT S T A T E 1585
C 1590
C 1595
C C A L C U L A T E TOTAL H O L E S IB (SOLIB) A B D TOTAL B O L E S O O T (SOLODT) 1600
C 1605
300 C O B T I B O E 1610
DO 310 1-1,3 1615
SOLODT(I)-A(BTOST)•!(I.BTOST,1)tO(1)•!(I,1,1) 1620
DO 3 1 0 J-1,BTOST 1625
S0L00T(I)-SOl00T(I)tALTBT(J)*I(I,J,1)t0L?RT(J)«I(I,J,1) 1630
310 C O B T I B O E 1635
C 1600
DO 330 I«1,3 1605
IF (SOLIB(I).GT.O.) GO TO 3 2 0 1650
Dr(I)-0. 1655
GO TO 330 1660
320 C O B T I B O E 1665
DF(I)<ABS((SOLIB(I)-SOLODT(I))/SOLIB(I)) 1670
330 C O B T I B O E 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 3 8 0 1720
C 1725
300 COBTIBOE 1730
BDF-ADF 1735
IF (ICBB.OE. 3 . A B D . D D F . L T . 0 . 0 0 0 ) 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 C O B T I B O E 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 B E A C H E D A B D I C E B > 0 , P B I B T T I B E IHTEBVAL 1800
C ABD P E B C E B T S T E A D T S T A T E 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 I S TO BE PHIBTED 2005
C 2010
IF (^[PR.EO.O) GO TO 000 2015
IF (UPS. E g . . - 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 I S 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 ( F 8 . 0 , 1 0 1 2 ) 2785
570 FORBAT (•1CALCDLATIOBS FOR A SOLTEHT EXTRACTIOB PROCESS',' HATIHG 2790
13 IBTERACTIBG SOLDTES'/) 2795
580 FORBAT (• ' , 1 0 A 8 ) 2800
590 FORHAT ('OFEED STREAB DATA',16X,'STAGE*,3X,'DRABIDB',6X,'TBORIDB', 2805
17X,•Ba03',5X,'ALDBIBOB B03',2X,'FLOW R A T E ' , 6 X , * T E B P ' , / , 3 0 X , * B 0 . ' , 5 2810
2X,'(6/1)',81,'(G/L)',8X,*(B)•,10X,*(B)',6X,'(TOL.OBITS)',5T,"DE6 C 2815
3'/) 2820
600 FOBHAT (20X,'AQDEOOS ' , 1 0 , 0 ( 2 X , 1 P E 1 1 . 0 ) , 3 X , E 9 . 3 , 3 X , 0 P r 7 . 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 E A C R 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 ) ' , 6 I , * (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 295
C PRO 300
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 / F L I S / IRATIO,IODT,ITOl,RICY(100) ,ICHHGI,ALTRT(IOO),0LTRT(1FL0 10
100) FLO 15
COHHOB / P U D / 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) FLO 220

00 COBTIBDE FLO 225
IBA-I,»1 FLO 230
FLO 235
BORE CARDS? FLO 200
FLO 205
50 COBTIBDE FLO 250
IF (IBDEX.IQ.1) 60 TO 10 FLO 255
FLO 260
IE DO lOT PROTIDE FOR RITDRIIBG STREAHS FLO 265
FLO 270
FLO 275
SET AODIODS ABD ORGABIC IITEBSTAGB FLOI (A ABD 0) FLO 280
FLO 285
60 COBTIIOE FLO 290
FLO 295
READ DIDSDAL E0DTIB6 PATTEBB FLO 300
0B6AHIC STREAH FROH STAGE ISTR FEEDS IBTO STAGE JSTB FLO 305
FLO 310
ir (BSTR.EO.1) READ 300, ISTR,JSTB FLO 315
DO 80 1-1,BTOST FLO 320
IP (I.GT.1) 60 TO 70 FLO 325
A(1)»APDRT(1)-ALTBT(1) FLO 330
0(BTOST) -OTDRT (BTOST) -OLTRT (HOST) FLO 335
GO TO 80 FLO 300
70 COBTIBOE FLO 305
A (1) -A (1-1) •AFDRT (I) -ALTRT (I) FLO 350
IST-HT0ST»1-I FLO 355
0(1ST)«0(IST»1)•OFDRT(BST)-OLTRT(BST) FLO 360
i r ( n S T R . i a . l . A I D . B S T . I O . I S T B - 1 ) 0 (BST)-OFDRT (BST)-OLTBT (BST) FLO 365
IF (HSTR.ig.1.AB0.IST. IQ. JSTB) O(BST)-O(ISTB) •OPDRT(IST)-OLTRT (BSTFLO 370
l)Ô(HSTÎ) FLO 375
80 COITIBOE FLO 380
FLO 385
THE IBTBBSTAGE FLOWS ABE BOW SET FLO 390
FLO 395
DEFIBE PHASE RATIOS FLO 000
FLO 005
RCTC»(I)-0.0 PLO 015
BCICO(1)-0.0 FLO 020
90 COBTIBDE FLO 025
FlOBAX-0.0 FLO 030
K-1 FLO 035
IF (IBATIO.E0.1) GOTO 110 FLO 000
FLO 005
DEFIBE PHASE BATIO WHII IRATIO-O,PHASE RATIO-FLOW RATIO FLO 050
FLO 055
RAIO (1) • (ALTRT (I) •A (I) ) /(OLTRT(I) •O (I) ) FLO 065
100 COBTIIOE FLO 070
60 TO 130 FLO 075
FLO 080
DO BOT RIDEFIHE PHASE BATIO IF ICHB6E-1 FLO 085
FLO 090
110 COBTIBOE PLO 095
IF (ICHI6E.Ea. 1) 60 TO 150 FLO 500
FLO 505
DEFIBE PHASE RATIOS WREB IBATIO-1, PBASI BATIOS DO BOT EQDAL FLO 510
FLOW RATIOS FLO 515
I-BDHBER OF STAGES HATIHG PRASE RATIO SHOII FLO 520
FLO 525
BEAD 310, I,RATIO FLO 530
I-KÎ-1 FLO 535
DO 120 J=K,I FLO 500
RAIO(J)-RATIO FLO 505
120 COBTIBDE PLO 550
IF (I.EQ.BTOST) 60 TO 150 FLO 555
K-1^1 FLO 560
GO TO 110 FLO 565
FLO 570
CALCOLATE IICESSART RECYCLE FLOIS FLO 575
TYPE 1 IS RECYCLE OF OBE PHASE OILY, TO SATISFY PHASE BATIO FLO 580
DSED WHEI IRATIO = 1 FLO 585
TYPE 2 IS RECYCLE OF BOTH PHASES IITH PHASI RATIO RAIO(I) FLO 590
TO ALLOW TOTAL STAGE FLOW TO EQDAL FLOBAX FLO 595
DSED IBIB ITOL = 0 , 2 FLO 600
FLO 605
FLO 610
CALCOLATE FLOTOT.BECY,ABD FLOBAX FOR IBATIO-O; ITOL-0,2,3 FLO 615
FOR ITOL-1, RECY SET TO 0.0 FLO 620
FLOTOT - TOTAL PLOW LIATIHG STA6E FLO 625
FLOBAX - GREATEST TOTAL FLOW FROH A STAGE FLO 630
FLO 635
130 COBTIBOE FLO 600
IF (ITOL.EQ. 1) GO TO 210 FLO 605
DO 100 1=1,BTOST FLO 650
FLOTOT=ALTRT (I) •A (I) •OLTRT (I) •O (I) FLO 655
RICI(I)-FLOTOT FLO 660
PL0HAX-A8AX1(PL0T0T,FLORAX) FLO 665
100 COBTIBDE FLO 670
IF (ITOL.EQ.O) 60 TO 190 FLO 675
IP (ITOL.E0.1) 60 TO 210 FLO 680
IF (ITOL.EO. 2.0S.IT0L.EQ.3) 60 TO 230 FLO 685
FLO 690
32
C CALCOLATE FLOTOT,RCYCA,RCYCO,RECI,FLOHAX FOR IRATIO-1 FLO 695

C ABO I T O L - 0 , 1 , 2 , 3 FLO 700
C AFL - TOTAL AQDEOOS FOLI LEATIBG A STAGE FLO 705
C OFL - TOTAL ORGABIC FLOI LEATIB6 A STAGE FLO 710
C RFL - BATIO OF AQOEODS TO OBGABIC FLOI LEATIIG A STAGE (AFL/OFL) FLO 715
C RCYCA - AQDEOOS BECTCII TO SATISFY PHASI BATIO FLO 720
C HCYCO - ORGABIC BICYCLE TO SATISFY PHASE RATIO FLO 725
C FLO 730
DO 180 1-1, BTOST FLO 700
AFL-ALTRT(I)Â(I) FLO 705
OFL-OLVRT(I)»0(I) FLO 750
FLOTOT-AFLÔFL FLO 755
RFL-APL/OFL FLO 760
C FLO 765
C COBPABISOB OF FLOB BATIO (RFL) TO PRASE RATIO (RAIO) FLO 770
C FLO 775
IF (BFL.GT.RAIO(I)) GO TO 160 FLO 780
IF (RFL.EQ.RAIO(I)) GO TO 170 FLO 785
FLO 790
FLOTOT-FLOTOT- (RAIO(I) • I . O ) / ( R F L ^ I . 0) FLO 795
RCYCA(I)-PLOTOT-APL-OFL FLO 800
GO TO 170 FLO 805
FL0T0T-F10T0T»RFL» (RAIO (I) • 1. 0) / (RAIO (I) • (RFl^l. 0) ) FLO 815
RCYCO(I)-PLOTOT-APL-OFL FLO 820
RECI(I)-FLOTOT FLO 830
FL0BAX-ARAX1(FLOTOT,FLOHAX) FLO 835
AFLO 805
FLOBAX COBTAIBS TBE LABGEST TOTAL FLOI FROB A STAGE - STAGES BATIBG FLO 850
TOTAL FLOI LESS TBAB FLOBAX HOST EXBIBIT RECYCLE FLO 855
OF TYPE 2 FLO 860
FLO 865
IF (ITOL.EQ.0) GO TO 190 FLO 870
IF (ITOL.EQ.1) GO TO 210 FLO 875
IF ( I T 0 L . E a - 2 . 0 I . I T 0 L . EQ. 3) GO TO 230 FLO 880
FLO 885
ITOL-0; IRATIO-O,1 FLO 890
FLO 895
FLO 910
RECY(I)-FLORAX-RECY(I) FLO 915
200 COITIBDE FLO 920
EETDBB FLO 925
FLO 930
ITOL-1; IRATIO-O,1 FLO 935
FLO 900
210 COITIBDE FLO 905
RECI(I)-0.0 FLO 955
BITDBH FLO 965
FLO 970
DITERRIIE STAGE TOLDRE FLO 975
STAGE TOLDRES ARE ZERO FOR ITOL-0,1 FLO 980
FLO 985
FLO 990
IF ICHB6E-0, READ STAGE TOLDBES FLO 995
IF ICHIGE-1, BYPASS BEADIBG OF STAGE TOLDBES FLO1000
FL01005
230 COITIBOE FLO1010
PLO1015
IF (ICBBGE.EQ.1) 60 TO 260 FL01020
IF (IT0L.EQ.2) GO TO 250 FLO1025
PLO1030
FL01035
ITOL-3; STAGE TOLDBES ARE EQDAL ABD SPECIFIED FL01000
IRATIO-O,1 FL01005
FLO1050
READ 320, TOL FLO1055
DO 200 J-1,BTOST FL01060
STTOL(J)-TOL FL01065
200 COBTIBDE FL01070
60 TP 260 FL01075
FL01080
ITOL-2: STAGE TOLOHES ARE IBDIPIBDEITLY DIEQDAL FLO1085
FLO1090
250 COITIBDE FL01095
READ 330, (STTOL(J),J-1,BTOST) FL01100
FL01105
CALCOLATE LEIOTH OF TIHE IBTEaTAL(TIBIB) ABD RECY FOB ITOL-2,3 FLO1110
FL01115
260 COBTIIOE FL01120
TIHI1=0.0 FL01125
DO 270 1-1,BTOST PLO1130
RECY (I) -BICY (I) /STTOL (I) FL01135
TIBIB-ABAI1(RICY(I).TIBIB) FL01100
FL01105
270 COBTIBOE FL01150
FLOBAX-TIBII FL01155
DO 280 1-1,BTOST FL01160
BECYfI)=(FLOBAX-RECI(I))*STTOL(I) FL01165
280 COBTIBDE
33
FLO1170
290 FORHAT (2I2.F8.0,I2) FL01175
300 FORHAT (215) PLO1180
310 PORRAT (12.F8.0) FLO1185
320 FORBAT (18.0) FLO1190
330 FORBAT (10F8.0) FLO1195
IHD FL01200-
SDBBODTIBE STA6ES (SPH) STA 5

DIHEBSIOI AII(O), 011(0) STA 10
CORROI /STGS/ ISOLD,CODDB(IOO),AB(IOO) ,OR (100) ,AT(100),OT(100) STA 15
CORROI /FEED/ ITOST,AFDRT(100),OFDBT(100) STA 20
COHHOB /PROF/ IPBO,X (1,100, 2) ,Y (0,100,2) , XFD(0,100) ,YFD(0,100) ,TPRSIA 25
10F(100,2) ,AFDTIH(10a) ,0FDTBH(100) STA 30
CORHOH /lie/ A(100),0(100),BAIO(IOO) ,RCYCA(100),RCYCO(100) STA 35
COHROB /DIST/ CTBP,ARI(0),DTRY(0),T0 STA 00
COBBOB /BODTE/ BSTR,ISTR,JSTB STA 05
C STA 50
C PEBFORR STAGE CALCDLATIOBS FOB CDRRIBT TIHE IBTEBTAL - COBSIDIR STA 55
C IBPDTS FIOB AOJACEBT STAGES, FEED STBEARS, ABD HOST IBPORTABTIT STA 60
C BECYCII WITHIB THE STAGS - RECYCLE A BESDLT OF STAGS FLOW BSIW6 STA 65
C LESS THAB FLORA! (TYPE 2) AID PBASI RATIO lOT EQDAL TO FLOI STA 70
C RATIO (TYPE 1) STA 75
C STA 80
C AR - AQOEODS RECYCLE TO FILL STAOE STA 85
C OR - ORGABIC BICYCLE TO FILL STAGE STA 90
C ECICA - AQDEOOS RECYCLE TO SATISFY PHASE RATIO STA 95
C BCYCO - ORGABIC BECYCLE TO SATISFY PHASE RATIO STA 100
C AT - TOTAL AQOEODS FLOI IITHIH A STAGE = STA 105
C TOTAL AQOEODS FLOI LEATIBG A STAGE • AQOEODS BECYCLE STA 110
C TO FILL STAGE • AQOIODS BICYCLE TO SATISFY PRASI BATIO STA 115
C OT - TOTAL OBGABIC FLOW WITHIB A STA6I - STA 120
C TOTAL 0R6ABIC FLOI LEATIBG A STAGE • OBGAHIC BECYCLE STA 125
C TO FILL STAGE • ORGABIC RECYCLE TO SATISFY PHASE RATIO STA 130
C XFD - AQDEOOS PEED COBTEBT OF STAOE STA 135
C YFD - 0R6ABIC FEED COBTEBT OF STAOE STA 100
C AFDRT - AODIODS STA6I FLOW EATF STA 105
C OFDRT - OReABIC STA6I FLOW RATE STA 150
C X - SOLDTI COICIBTBATIOB IB AQDEOOS PHASE FOR A 6ITEI STAGE STA 155
C Y - SOLDTE COBCIBTEATIOH IB ORGAHIC PHASE FOR A GITEI STAGE STA 160
C AFDTIH - TIRPIIATDRE OF AODEODS FEED STA 165
C OFDTEH - TEHPIBATDRI OF ORGAIIC FEED STA 170
C TPROF - TEHPIBATDRI OF SOLDTE IB A GITEB STAGE STA 175
STA 180
C AEY - TRIAL AQDEOOS PHASE COBPOSITIOH STA 185
C DTRY - DISTBIBOTIOH RATIOS STA 190
STA 195
DO 100 J-1,BTOST STA 200
JJ«J-1 STA 205
JJJ=J^1 STA 210
DO 50 1-1,ISOLD STA 215
IF (J.GT.1) GO TO 10 STA 220
AIH (I) -AFDRT (J) •XFD (I, J) • (AR (J) •RCYCA (J) ) •X (I, J, 1) STA 225
60 TO 20 STA 230
10 COBTIBOE
AIH(I)=AFDRT(J)»XF0(1,J)^(AR(J)+BCYCA(J))^X(I,3,1)Â(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)ÂFDTEH(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
1 RCYCO ( J ) ) •TPROF ( J , 1) •O (ISTR) •TPROF (ISTR, 1) •O (JJ J) •TPROF ( J J J , 1))

STA 395
60 TO 90 STA 000
80 COBTIBDE STA 005
OHTIB-SPB*(OFDRT(J)•OFDTIB(J)•(OB(J)•BCYCO(J))•TPBOF ( J , 1 ) ) STA 010
90 COBTIBDI STA 015
TPROF ( J , 2) • (ARTIBÔBTIB) / (AT (J) •SPBÔT (J) ) STA 020
TIRPC-TPR0F(J,2) STA 025
c CCD • 0.001 STA 030
CCD-0.005 STA 035
c STA 000
c BCK IS TBE ITERATIOH CODBTEB STA 005
c BDEX IS TBE FIT CODBTEB STA 050
c STA 055
RCK-0 STA 060
100 COBTIBDE STA 065
BDIX-0 STA 070
c STA 075
c DCOR CALCDLATES DISTBIBOTIOH RATIOS (DTBY) FOR EACH TBIAL AQOEODS PBASTA 080
c CORPOSITIOB (ABY) AT TBE STAOE EODILIBBIDR TERPIBATDRE (TEBPC) STA 085
c STA 090
CALL OCOB (ABY,DTBY,TEBPC,CTBP) STA 095
DO 110 I<1,BS0LD STA 500
TSIB-AIB(I)-BAIO(J)•Oil(I) STA 505
IF (TSII.LT.O.) DTRY(I)-0. STA 510
DTIY(0)-0. STA 515
X ( I , J , 2 ) - T S I B / (RAIO(J)•DTBY(I)) STA 520
IF ( A B S ( X ( I , J , 2 ) ) . L T . 1 . E - 6 ) GO TO 110 STA 525
01-(X ( I , J , 2)-ARY (I) ) / X ( I , J , 2) STA 530
DIF-ABS(DI) STA 535
ARI(I)-X(I,J,2) STA 500
IF (DIF.LE.CCD) GO TO 110 STA 505
c STA 550
c BDEX CODBTS lORBER OF SOLDTES THAT HOT C0ITER6ED STA 555
c STA 560
HDEX-BDEXÎ STA 565
110 COBTIBOE STA 570
BCK-BCKÎ STA 575
c STA 580
c ALLOI OILY 20 ITERATIOIS STA 585
c STA 590
IT (HCK. 1 0 . 2 0 ) 60 TO 120 STA 595
c STA 600
c IF ALL SOLDTES RATI lOT C0ITER6ID, RIPIAT CALL TO DCOR STA 605
c IITH I I I TALDIS OF ARY STA 610
c STA 615
IF (RDIX.GT.O) 60 TO 100 STA 620
120 COBTIBDI STA 625
DO 130 1-1.ISOLD STA 630
Y(I.J.2)-I(I,J,2)^DTEY(I) STA 635
130 COBTIBDI STA 600
CODaB(J)-HCK STA 605
100 COITIIDI STA 650
RITORI STA 655
c STA 660
c STA 665
IBD STA 670
SDBRODTIIE DCOR (ARY,DTRY,TEBPC,CTBP) OCO 5

BBAL*0 COIP(12,3),KD(3),KR(3),KTB(3),ABY(0),DTRY(0),A(0) ,TBP(3) DCO 10
REALÔ B03 OCO 15
C DCOR RITORIS A TALDI OF TBI DISTBIBOTIOH COBFPICIIIT - DCO 20
C DTRY OCO 25
C 0 - TOTAL lOIIC STREI6TB OCO 30
C TEBPC - TEHPEHATDEE IB DE6BEES CEBTI6BADE DCO 35
C ABY - AODEODS PBASE COflPOSITIOB DCO 00
C T - HOLABITY OP TBP IB DRY, SOLOTE-FREE SOLTEHT DCO 05
C KD,KTa,KH - BASS ACTIOI EQDILIBBIOB COISTAITS. OCO 50
C THESE EQDILIBBIOB COBSTAITS BAY BE EXPBESSID AS DCO 55
C A^BX^CY^DE DCO 60
C IRIBI X,Y, ABD Z ABE 1ST,2ID ABD 3RD POWERS OF OCO 65
C lOBIC STBEieTH OCO 70
C FOR KD(I) - A,B,C,D ARE C 0 E F ( 1 , I ) , C 0 E F ( 2 , I ) , C 0 E F ( 3 , I ) . OCO 75
C AID COIF(O.I) - D COIPFICIEITS OPTIBIZED OCO 80
C FOB 1-1 D ACCDRACY OCO 85
C 1-2 TH ACCDRACY OCO 90
C 1-3 H ACCDRACI DCO 95
C FOB KTB(I) - A.B.C.D ARE C 0 E F ( 5 , I ) , C 0 E F ( 6 , I ) , C 0 E F ( 7 , I ) , OCO 100
C ABD C0EF(8,I) - TB COEFFICIIBTS OPTIBIZID OCO 105
C FOR 1-1 0 ACCDRACY DCO 110
C 1-2 TH ACCDRACY DCO 115
C 1-3 HACCORACY DCO 120
C FOR KH(I) - A,B,C,D ARE C 0 E F ( 9 , I | , C O E F ( 1 0 , I ) , C 0 E F ( 1 1 , I ) , OCO 125
C ABD C0EF(12,I) - H COEFFICIEBTS OPTIHIZED OCO 130
C FOR 1-1 0 ACCDRACY DCO 135
C 1=2 TH ACCDRACY OCO 100
C 1=3 H ACCDRACY OCO 105
35
C DCO 150
C DCO 155
DATA C O I F / 5 . 1 6 7 6 0 , 1 2 . 7 3 8 5 , - 6 . 1 0 9 3 9 , 0 . 8 2 5 6 6 5 , 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)ÂRY(3) DCO 205
IO3«2.0*ARY(1)Ô.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 ) • C O I F ( 2 , 1 ) • D + C O I F ( 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)ÂBY(1) •B03^^2 DCO 280
A(0)-3.0*KTR(I)ÂIY(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 . A I D . A B S ( D ) . L T . 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 7 0 , TBP(I) , I , K 0 ( I ) . K T H ( 1 ) , 1 8 ( 1 ) ,ABI, (C0EP(J,I) , J - 1 , 1 2 ) DCO 395
Y=0. DCO 0 0 0
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 0 3 0
C TEBPEBATDRI ADJDSTHEBT OCO 035
C OCO 000
IF (TERPC.E0.25.) 60 TO 60 DCO 005
IT (TEHPC.LT. 15..0E.TERPC.6T. 7 0 . ) GO TO 60 DCO 050
XTIHP-1./(273.^TEHPC)-1./298. OCO 055
DTRY(1)-DTRY(1)ÊXP(2500.^XTIHP) DCO 060
DTRY (2)-DTBY (2)ÎXP( ( 8 9 1 . 660^879. SBSÂBY (2) - 1 0 7 0 . 38Â»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 ' . 1 P E 1 1 . 3 . ' FOE COBPOBEBI ' . 1 2 . / . * OCO 505
1PSE0D0-EQDIL1BRIDB COBSTAITS < , 3 E 1 1 . 3 , / , ' AQOIODS PHASE CORPOSITIODCO 510
2B ' . 0 E 1 1 . 3 , / , * PARABETERS * , 0 E 1 1 . 3 , / , ( ' • , 1 I X . O S I I . 3 ) ) OCO 515
IBB DCO 520-
SOBBODTIBI PLOTR (KOOIT,ICALPL,HT0ST,S01IB) PLO 5

C PLO 10
C THIS HODTIBE ISEHERATES CALCOHP PLOTS OF TWO KIBDS: PLO 15
C ICALPL = 1 COBCEITRATIOB TS TIHE IBTEBTAL PLO 20
C ICALPL - 2 PEBCEBT OF PIED TS TIHE IBTEBTAL PLO 25
C ICIILPL = 3 PLOT BOTH PLO 30
C PLO 35
C KODIT - TOTAL HDHBE8 OF TlHI IBTERTALS PLO 00
C HIOST - TOTAL HDHBES OF STAGES PLO 05
C SOLII - TOTAL ARODIT IB PEED PLO 50
36
A - AQOIODS COICIITRATIOIS PLO 55

0 - OBGABIC COICIITIATIOIS PLO 60
AODT - AODEODS FLOI OOT OP STAGE PLO 65
OODT - ORGAIIC FLOI ODT OF STAOE PLO 70
PLO 75
DIBEBSIOB A ( 0 , 1 0 0 ) , 0 ( 0 , 1 0 0 ) , AOnT(IOO), OoaT(IOO), 1 ( 1 5 0 0 ) , Y ( 1 5 0 P L a 80
1 0 , 3 ) , SOLIB ( 1 0 0 ) , 1C0R(3), 8DFFER(0000) PLO 85
REAL«8 YLAB(8)/'C0ICEIIB',*ATI0I ( *,*0,TB IB ' , ' 0 / L - , HBO', '3 II RPLO 90
1 ) ' , ' « ' , ' P I B C I I T ' , ' O F FIIDS*/ PLO 95
BIAL^B XIAB(2)/'TIBE I B T ' , ' EBTALt'/ PLO 100
REAL»0 C 0 B ( 3 ) / ' D * , * TR ' , ' H • / , L L ( 2 ) / ' L I B ','LOG ' / PLO 105
L0GICAL»1 AQ0/*A*/,OIG/*O'/,IAO PLO 110
LOGICALÎ DESCB(3a) PLO 115
LOGICAL*1 DISC 1 ( 3 0 ) / ' D ' , ' I ' , ' A ' , ' B ' , ' I ' , * 0 * . * H * , * » * , * T * , ' H ' ,*0*,*BPLO 120
1','I','0','R*,*t*.*B*,*I*,'T','B*,*I*,'C',' ' , ' A ' , * C * , * I ' , ' D*,*$*/PLO 125
PLO 130
IF ICALPL - 3 , PLOT BOTB COICIITRATIOI TS TIBE ABD PLO 135
PERCEBT OF FEED TS THE PLO 100
PLO 105
IC-0 PLO 150
IF (ICALPL. EQ. 3) IC-3 PLO 155
IF (ICALPL.EO.3) ICALPL-1 PLO 160
PLO 165
CALL QQBIIS (0) PLO 170
IRITI ( 9 , 1 9 0 ) PLO 175
IBITI ( 9 , 2 0 0 ) PLO 180
IPL-0 PLO 185
PLO 190
BIAD SIT OP PLOT SIIICBES PLO 195
BSTB - STIIAB BDBBEB TO BE PLOTTED PLO 200
BCOR - BDRBII OF CORPOBEBTS TO BE PLOTTED PLO 205
lAO - A PLOT AODEODS STREAH PLO 210
0 PLOT ORGAIIC STREAH PLO 215
ILL - 0 LIIIAI X-AXIS, IIBEAR Y-AXIS PLO 220
1 LIBIAR I - A I I S , LOG Y-AXIS PLO 225
BAIT - 0 PLOT BTIBT POUT PLO 230
1 PLOT ITER! OTHER POUT PLO 235
ILAST - 0 LAST PLOT PLO 200
1 OTHER PLOTS FOLLOI PLO 205
PLO 250
PLO 255
BEAD IIDICES OF COBPOIEITS TO BE PLOTTED PLO 260
1 - DBABIDB PLO 265
2 - TBOBIDB PLO 270
3 - BITBIC ACID PLO 275
PLO 280
10 COITIIDE PLO 285
READ ( 5 , 2 1 0 ) ISTR,ICOB,IAO,ILL,IALT,ILAST PLO 290
READ ( 5 , 2 2 0 ) (ICOB (I) ,1-1,ICOB) PLO 295
20 COBTIBDE PLO 300
IPL-IPLÎ PLO 305
PLO 310
SET OP POIBTS TO BE PLOTTED PLO 315
PLO 320
YRAX-0. PLO 325
YBII-1.E10 PLO 330
EEIIID 11 PLO 335
I-O PLO 300
DO 50 lK-1,KODIT PLO 305
PLO 350
BEAD COICEITBATIOBS (A ABD 0) AID FLOIS (AODT ABD OOOT) PLO 355
PLO 360
DO 30 J-1,BTOST PLO 365
BEAD (11) ITB, (A(K,J),K-1,0), (0(K,J).K<1,0),AODT(J),OODT(J) PLO 370
30 COBTIBDE PLO 375
IF (lALT.EQ.1.AID.IK/2^2.EQ.IK) GO TO 50 PLO 380
1-1^1 PLO 385
DO 00 K-1,IC0H PLO 390
IF (ICALPL. i g . L A I D . lAO.IQ.AQD) Y (I ,K)-A (ICOB (K) , BSTB) PLO 395
IF (ICALPL. i g . l . ABD. IA0.ta.ORG) Y ( I , K ) - 0 (ICOB (K) ,BSTR) PLO 000
IF (ICALPL. IQ. 2 . AID. IAO.EQ.AQD) T (I , K)-A (ICOB (K) , ISTR) •AODT(BSTB)/PLO 005
1 SOLII (ICOB ( K ) ) ^ 1 0 0 . PLO 010
IF (ICALPL.EO. 2 . ABD.BAO.E0.ORG) Y ( I , K ) - 0 (ICOB (K) , BSTB) •OODT(BSTR)/PLO 015
1SOLIB(ICOR(K))•100. PLO 0 20
YRAX-AHAX1(Y(I,K) ,YRAX) PLO 025
IF ( B L L . I 0 . 1 . A B D . Y ( I , K ) . L E . O . ) GO TO 00 PLO 030
YBII-ARII1 (Y(1,K),YHII) PLO 035
00 COITIBDE PLO 000
X(I)-I PLO 005
50 COBTIBDE PLO 050
PLO 055
PLO 060
YBAX=IHAX^YHAX^.05 PLO 065
XRII=0. PLO 070
XBAX-KCDBT PLO 075
CALL IBTL (•BECB',BDFFER,0000,'S.B. lATSOBS') PLO 080
CALL QOCCPA (3,XRIB) PLO 085
CALL QQQCPA (0,IBAX) PLO 090
CALL QQQCPA (13,YBIB) PLO 095
CALL QOCCPA (10,YBAX) PLO 500
IF (KO0IT.LE.100) CALL SCALE (1 , 10. , 0. ,'LIES') PLO 505
IF (KODIT.61.100) CALL SCALE (1,KOOBT/10.,0.,'LIKE') PLO 510
IF (BLI.ECO) CALL SCALE (-1, 10. ,0. ,'LIIE') PLO 515
IF (BLI..EC.1) CALL SCALE (-1. 10. . 0. ,'L062') PLO 520
CALL AXIS (1,-15,.15) PLO 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
DO 100 1-17,28 PLO 805
DESCR(IP)-DISCI(1) PLO 810
IP-1P^1 PLO 815
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 , 2 3 0 ) IPL,HSTR,BAO,LL(BLLÎ),(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
ICALPL=2 PLO 890
GO TO 20 PLO 895
ICALPL-1 PLO 905
60 TO 10 PLO 910
CALL ADTAIS PLO 920
C PLO 925
C PLO 930
RITDRI PLO 935
C PLO 900
190 PORHAT (' ' , / / / , ' ' , 2 1 X , ' P L 0 T SDBHARY'/) PLO 905
200 FORHAT ( ' ','PLOT*,0X,*STA6E',0X,*AQDEODS OR*,0X,*LIIEAR OH'.OI.'CPLO 950
10HP0IEITS*/* * , 1 X , ' I 0 . ' . S X . ' I O . ' , 6 X , ' 0 R G A I I C ? ' , 7 X , ' L 0 6 ? ' ,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-
APPENDIX B
SAMPLE INPUT AND CORRESPONDING COMPUTER OUTPUT
SAMPLE INPUT TO SEPHIS-TWREX
—1
—7
Ma [71 4 s 8 9 10 n la Fl;? IS 16 17 1» 19|2D|2I|22 23l2«|2s|2E 27y29l3oUl]u{3)lM 3S 36 3,|3.|39 40|4l|42|43|44|4»|4«|47{M|49Js0]sl]32|s3|s«|ssl96|57|H|s9|co|tt|c2|«3|«4|6^|«|H|s9|;0J7l|72|73|74|7«|7f|77|;a|»|nj
l 2 .• ^ • : : • ; ! : : - !! M 1
AGIO T H9 REX FLOWSHEET wm EX T R A C,T I » H| 9 Fj URAN I UM AN0 T H9 R I UM j
FRIiN S0LUTI0NS 0F C9N sifUi 0 A T iE 0 E O I S J H
' ! i
FUEL ( jS E E :R E F . 7 , F I S . 2) ; i j
j . 3 1 1 0 0 1 0 1 0 0 lOi I• 1; . ]
1 1 . . 8 0 . . 1 .0
!:':,:; 1
1 4. . 1 , ;. . \-'.2[ : . : : : o : . ! .| \ , ; : . ! 0 ; . i \ , i \ \ 's[.\ ; :.; i ; 4 ; o : . ; . : i , M .1 : i : 0 , - : ; . ; ::; i : ; : . : ' . • ; I ; ! ; . . . : 1
1 T 1 ; . . ,1.-, ^ i . 2 i 0 i . :0i 1, ; , ;2;6:5Li i ' r l - i^ls, ; : , i .3:5:. '^T 1 ;. . ; ; ; :0.:.. i . ; : , . ' . : . i';: • / : ; ; ; : 1, . . |
10 1 -.3 ,0 . , ,0 . 1 ' 3 ; . , 1 : i i ;2^5:. : ^' ^ ] : • \ -0 . • ': ] . : \ j
1 1 0 7 . 0 . Oi. ;8j-| h i ; :2:51.i ,0: i ; i , , :<>:. : ; _ ; ; , . ! . ; : . , / ! ; _ . : , . |
1 2 ' i ! i 1 • T i t• i ! i ' • i ! • : t - ( • • • s ^ • ^ • i • : " ! ' 1 • i -j i • : • ' . - 1
• : I ' •
1
• i •
•
... ; .
"1 ^
•
. ':
•
•
! • ' i .
1 , i : 1 . 1 f .
' :
. 1 . j :
\ : •
:
^ i
1
\
1 r
' '. '
: • • f
:
!
;1 ; ;
l!
; •
:
h ^ . 1 1 t • i :' 1 i ' : . i r •
• ; '
- r •
•
i i i- 1 • • ;
i i ' f • ; ;
• i • ' ! ' '
1 ;i i i
;
i i 1 :
! ; t
; i : i i
— '•M
! ! • j !" !1
I i
'i : I ' M
i • ' ;
; ! • 1
1 " • : i ;
i1
., • t - 1! ;
i
1
•
I
1
• ! ' ; ! 1 1 . • ; : , i
1 j' -J.i -... ; i' i[

1 ' ' ' i !
• 1 - -i i J
[ ; ; ; 1
1 i. 1 • • - 1 ;
ii.H
. 1 . i i I i i
! i1
il L
DATA SHEET (SOcaiiNM)
_ 1 i!
t i
11 ln\ •
• i M ! ! h
; i 1M . i ;
M i ,i
42
C»1CD1»II0»S FOR A SOLTERt EXTRACTIOR PROCESS HATIRS 3 IRtERACTIRG SOIOTES

ACID THOREX FlOtSREET FOR EXTRACTIOR OF ORARIDR ARD TgORIDR
FROa SOLOTICRS OF CORSOLIDATED EDISOR FUEL (SEE REF. 7 , FIG. 21
FIED STREAK DATA STAGS DRARIOn AIDRIHDR R03 FLOa RATE

RO. (S/1) (G/l) («) (") (rOL.ORITS) DEG C
AOOEOOS 0 0.0 0.0 0.0 8.000E-01 • 0.0
AOOEODS 0 0.0 S.OOOOE 00 0.0 2.000E-01 • 0.0
AOOEOnS 7 OOOOE 01 2.t500E 02 -1.5000E-01 0.0 1.000E 00 35.0
AQOEOnS 10 0 0.0 1.3000E 01 0.0 3.000E-01 25.0
3 0 . 0 - TBP 11 0 0.0 0.0 0.0 7.000E 00 25.0
IRITIAL PROFILE STAGE AQOEOOS PHASE OBGARIC PHASE
RO. ORARIOn THORIOa RR03 ORARIOR TRORIOB HH03 TEB?
(OOTFIOR) (6/1.) (») (G/l) (G/1) DEG :
1 0.0 0.0
(G/l)
1.OOOOE 00 0.0 0.0 0.0
m • 0.0
2 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
3 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
« 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
5 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
6 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 •0.0
7 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
8 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
9 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
10 0.0 0.0 1.. OOOOE 00 0.0 0.0 0.0 • 0.0
11 0.0 0.0 1.OOOOE 00 0.0 0.0 0.0 • 0.0
PRASE RATIOS tgOAI. FLOR RATIOS STAGE HO. A/O PBASE RATIO
1 - 3 1.1«3E-01
• - 6 1.»29E-01
7 - 9 2. 857E-01
10 - 11 3.286E-01
STAGE FtORS:
AOOEOOS PHASE
RECICIE 2 RECTCIE 1 TOTAL
STAGE FIOR FLOR OOT FOR PHASE AQOEOns
RO. RATIO OF STAGE STAGE RATIO FLOR
1 1.143E-01 8.000E-01 0.0 0.0 8.OOOE-01
2 1.U3E-01 8.000E-01 0.0 0.0 8.OOOE-01
3 1.U3E-01 8.OOOE-01 0.0 0.0 S.OOOE-OI
H29E-01 1.000E 00 0.0 0.0 1.000E 00
•
5 1.»29B-01 1.000E 00 0.0 0.0 1.000E 00
6 1.*29E-01 1. OOOE 00 0.0 0.0 1.000E 00
7 2.857E-01 2.000E 00 0.0 0.0 2.000E 00
8 2.857E-01 2.OOOE 00 0.0 0.0 2.OOOE 00
9 2.857E-01 2. OOOE 00 0.0 0.0 2.OOOE 00
10 3.28SI-01 2.300E 00 0.0 0.0 2 . 3 0 0 e 00
11 3.286E-ai 2.300E 00 0.0 0.0 2.300E 00
OTEBAU ORGARIC PRASE

TOTAL RECICLE 2 RECYCLE 1 TOTAL
STAGE STAGE FLOW OOT TO FILL FOR PHASE ORGARIC
RO. FLOI STAGE TOL. OF STAGE STAGE RATIO FLOR
1 7.800E 00 0.0 7.OOOE 00 0.0 0.0 7.000E 00
2 7.800E 00 0.0 7. OOOE 00 0.0 0.0 7.000E 00
3 7.800E 00 0.0 7.OOOE 00 0.0 0.0 7.OOOE 00
8.000E 00 0.0 7.OOOE 00 0.0 0.0 7.000E 00
«
5 S.OOOE 00 0.0 7.OOOE 00 0.0 0.0 7.000E 00
6 8.000E 00 0,0 7.OOOE 00 0.0 0.0 7.OOOE 00
7 9.000E 00 0,0 7.OOOE 00 0.0 0.0 7.OOOE 00
8 9.000E 00 0,0 7.OOOE 00 0.0 0.0 7.000E 00
9 9 . 0 0 0 E 00 0.0 7.OOOE 00 0.0 0.0 7.OOOE 00
10 9.300E 00 0.0 7.OOOE 00 0.0 0.0 7.OOOE 00
11 9.300R 00 0.0 7.OOOE 00 0.0 0.0 7.OOOE 00
43
TRARSIERT BEHATIOB RESOLTS SOLOTI CONCERTRATIORS I R EFFLDERT FROB EACH STAGE
AQOEOOS PHASE
TIBE STAGE ORARIOB THORIO B HR03 FLOR OOT
IHTERTAL RO. (GA) (G/L) ("1 OF STAGE
1 1 0.0 0.0 0.0 8.000S-01
2 0.0 0.0 «.108E-01 8.OOOE-01
3 0.0 0.0 •.108E-01 8.OOOE-01
t 0.0 0.0 6.913E-01 1.000E 00
5 0.0 0.0 «.«87E-01 1.000E 00
6 0.0 0.0 •.•87E-01 1.000E 00
7 5.208E-01 5..976E 01 2.722E-01 2.OOOE 00
8 0.0 0.0 5.79SE-01 2. OOOE 00
9 0.0 0.0 5.79«E-01 2.OOOE 00
10 0.0 0.0 1. 391E 00 2.300E 00
11 0.0 0.0 6.201E-01 2 . 3 00E 00
ORGARIC PBASE
STAGE ORARIOB THOSIOH HR03 FLOR OOT ITEHA- TEHP
RO. (GA) (G/1) OF STAGE TIORS DEG C
1 0.0 0.0 0.0 (") 7.000E 00 1. • 0.0
2 0.0 0.0 6.73aE-02 7.OOOE 00 8. • 0.0
3 0.0 0.0 6.73«E-02 7.000E 00 8. •0.0
0.0 0.0 1. 580E-01 7.000E 00 5. • 0.0
«5 0.0 0.0 7.875E-02 7.000B 00 7. •0.0
6 0.0 0.0 7.875E-02 7. OOOE 00 7. • 0.0
7 2.708E 00 2.07aE 01 • . 3 6 5 E - 0 2 7.000E 00 5. 38.9
8 0.0 0.0 1. 202E-01 7.000E 00 5. • 0.0
9 0.0 0.0 1.202E-01 7.OOOE 00 5. •0.0
10 0.0 0.0 3.a58E-01 7. OOOE 00 3. 39.0
11 0.0 0.0 1.2lt8E-01 7.OOOE 00 4. 32.3
AQOEOOS PHASE
TIBE STAGE ORARIOB TBORIOH HR03 FLOa OOT
IRTERFAL RO. (G/L) (G/L) (B) OF STAGE
2 1 0.0 0.0 2.803E-01 8.OOOE-01
2 0.0 0.0 2.838E-01 8.000E-01
3 0.0 0.0 6. 211E-01 8.000E-01
0.0 0.0 7.160E-01 1.000E 00
•
5 0.0 0.0 5. 261E-01 1.000E 00
6 0 . 3 3 5 1 - 01 «.368E 01 3.340E-01 1.000E 00
7 5 . 1 1 8 E - 01 5 . 7 6 6 E 01 3.6^0E-01 2. OOOE 00
8 2.332R- 02 2.91«E 01 •.088E-01 2.OOOE 00
9 0.0 0.0 1.003E 00 2.000E 00
10 0.0 0.0 1.»20E 00 2.300E 00
11 0.0 0.0 8. 173E-01 2.300E 00
ORGARIC PHASE
STAGE ORARIOB THORIOB HR03 FLOR OOT ITERA- TEHP
HO. (G/l) (G/L) OF STAGE TIORS DEG C
1 0.0 0.0 m
3.»8»E-02 7 . OOOE 0 0 9. • 0.0
2 0.0 0.0 3.«90E-02 7.000E 00 7. •0.0
3 0.0 0.0 1. 3O3E-01 7.OOOE 00 6. • 0.0
0.0 0.0 1.663E-01 7.000E 00 3. •0.0
•5 0.0 0.0 1.023E-01 7.000E 00 5. • 0.0
6 2.696E 00 1.«5«E 01 6.003B-02 7.OOOE 00 7. 39.2
7 2.711E 00 2.138E 01 5.886E-02 7.OOOE 00 0. 38.9
8 1.«21R- 01 8.7«9E 00 8.117E-02 7.OOOE 00 7. 39.5
9 0.0 0.0 2.6S9E-01 7. OOOE 00 39.5
10 0.0 0.0 3.811E-01 7.000E 00 2•.. 35.1
11 0.0 0.0 1. 8B5E-01 7.OOOE 00 31.9
•.
AQOEOOS PHASE
TIHE STAGE ORARIOB THORIOB HR03 FLOH OOT
IRTERTAL RO. (G/L) (GA) (H) OF STAGE
3 1 0.0 0.0 1.781E-01 8.OOOE-01
2 0.0 0.0 5. 361E-01 8.000E-01
3 0.0 0.0 6.064E-01 8.000B-01
a 0.0 0.0 8. 1178-01 1.000E 00
5 lt.089E- 01 3.0»OE 01 •.758E-01 1.000E 00
6 » . 2 2 5 B - 01 «.237E 01 «. 130B-01 1.000E 00
7 7 . 3 2 9 E - 01 8.38eE 01 2.567E-01 2.OOOE 00
8 1.e9«E- 02 2.202E 01 7. 359E-01 2. OOOE 00
9 7.963Ii- 04 9.95«E 00 9.512E-01 2.OOOE 00
10 0.0 0.0 1.738E 00 2.300E 00
11 0.0 0.0 8.36"tE-01 2.300E 00
ORGARIC PHASE
STAGE ORARIOB THORIOB HH03 FLOH OOT ITERA- TEHP
RO. (G/L) (GA) OF STAGE TIORS DEG C
1 0.0 0.0 1.456E-02 (") 7.000E 00 7. •0.0
2 0.0 0.0 1.056E-01 7.000E 00 7. • 0.0
3 0.0 0.0 1.29aE-01 7.000E 00 3. • 0.0
0.0 0.0 2.002E-01 7.OOOE 00 u. • 0.0
•
5 2 . 5 8 8 1 00 1.019E 01 9 . » 3 5 E - 0 2 7 . 0 0 0 E 00 7. 39.4
6 2.651E 00 1.533E 01 7 . 5 0 3 E - 0 2 7.OOOE 00 u. 39.2
7 2.852E 00 2.889E 01 3. « 1 1 ! - 0 2 7.OOOE 00 t. 38.4
8 1.4 08E- 01 1.018E 01 1.586E-01 7.OOOE 00 5. 39.2
9 6.U35E- 03 5.U8 3E 00 2. 261E-01 7.000E 00 5. 37.1
10 0.0 0.0 a.611E-01 7. OOOE 00 2. 34.7
11 0.0 0.0 1.916E-01 7.000E 00 2. 29.9
44
AQOEOOS PHASE
TIHE STAGE ORARIOB THORIOB HR03 FLOH OOT
IHTERTAL HO. (G/L) (G/L) (B) OF STAGE
• 1 0.0 0.0 3.880E-01 8 . 0 0 0 E-01
2 0.0 0.0 •.955E-01 8.OOOE-01
3 0.0 0.0 7.4^8t-01 8.000E-01
4 3.531E-01 1.672E 01 8.301E-01 1.000E 00
5 3.909ET01 2.963E 01 5.558E-01 1.000E 00
6 5.997E-01 e.528E 01 3.^67E-01 1.000E 00
7 7.305E-01 8.190E 01 •.71«E-01 2.000E 00
8 2.903E-02 3.899E 01 6.329E-01 2.OOOE 00
9 5.2S8E-0^ S.783E 00 1.2S3E 00 2.000E 00
10 1.432E-05 i.ee«E 00 1.761E 00 2.300E 00
11 0.0 0.0 1.000E 00 2.300E 00
ORGARIC PHASE
STAGE ORARIOR THORIOB HR03 FLOH OOT ITERA- TEBP
RO. (G/L) (G/L) (B) OF STAGE TIORS DEG C
1 0.0 0.0 6. 123E-0 2 7 . OOOE 0 0 8. • 0.0
2 0.0 0.0 9.312E-02 7.OOOE 0 0 • 0.0
3 0.0 0.0 1.76Ê-01 7.OOOE 0 0 • 0.0
4 2.S38E 00 7.806E 00 1.879E-01 7 . OOOE 0 0 39.6
5 2.595E 00 1.110E 01 1. 116E-01 7 . OOOE 0 0 39. «
6 2.82;E 00 2.391E 01 5 . 2 5 5 E - 0 2 7.000E 00 38.7
7 2.850E 00 3.069E 01 6 . 1 5 3 E - 0 2 7.OOOE 00 38.2
8 2.075E-01 1.830E 01 1.186B-01 7.000E 00 37.7
9 5.262E-03 4 . 6 3 9 E 00 3 . 1 3 ^ E - 0 1 7. OOOE 00 36.7
10 2. 228E-0^ 2.297E 00 • . • 1 8 E - 0 1 7.OOOE 00 32.7
11 0.0 0.0 2.^23E-01 7.000E 00 3'. 29.7
AQOEOOS PHASE
TIHE STAGE ORARIOR THORIOB HR03 FLOR DOT
IHTERTAL RO. (GA) (GA) (B) OF STAGE
5 1 0.0 0.0 3.55^B-01 8.OOOE-01
2 0.0 0.0 6.553E-01 8.OOOE-01
3 3.922E-01 1.581E 01 7.186B-01 8.OOOE-01
•5 3.337E-01 i . e e 2 E 01 9. 1 8 « ( - 0 1 1.000E 00
4.64Ê-01 • . 7 0 4 E 01 5.4UE-01 1.000E 00
6 6.150B-01 S.521E 01 • . 831E-01 1.000E 00
7 6.974E-01 1.005E 02 3.607B-01 2.OOOE 00
8 2.,e^2E-02 3.217E 01 9.355E-01 2.000E 00
9 7.9«6E-0^ 1.257E 01 1.217B 00 2.OOOE 00
10 8.084E-0e 8.320E-01 2.003E 00 2.300E 00
11 5.254E-07 6.203E-01 1.017E 00 2.300E 00
ORGARIC PHASE
STAGE ORARIOB TRORIOB HR03 FLOH OOT ITERA- TEBP
RO. (GA) (GA) (H) OF STAGE TIORS DEG C
1 0.0 0.0 5.250E-02 7.000E 00 5. • 0.0
2 0.0 0.0 1.^58E-01 7.000E 00 • 0.0
3 2.493E 00 5.999E 00 1.62Ê-01 7.000E 00 39.7
4 2.547E 00 8.723E 00 2 . 0 8 « E - 0 1 7. OOOE 00 39.6
5 2.809E 00 1.958E 01 9 . 3 8 1 E - 0 2 7.000E 00 39.0
6 2.818E 00 2.561E 01 7 . 1 9 0 E - 0 2 7. OOOE 00 38.6
7 2.951E 00 3.67SE 01 • . 3 6 5 S - 0 2 7.000E 00 37.3
8 2.084E-01 1.885E 01 1. 808E-01 7.000E 00 37.*
9 8. 289E-03 9.8^7E 00 2 . 7 5 0 B - 0 1 7. OOOE 00 3^.9
10 1.«76E-04 1.379E 00 • . 9 9 3 E - 0 1 7.OOOE 00 32. «
11 •.531E-06 3.436E-01 2.«^6E-01 7. OOOE 00 28.7
AQOEOOS PHASE
TIBE STAGE ORARIOB TBORIOB HR03 FLOR OOT
IHTERTAL RO. (G/L) (G/L) (H) OF STAGE
e 1 0.0 0.0 • . 872E-01 8.OOOE-01
2 •.572E-01 1.501E 01 6.05»E-01 8.0008-01
3 3.577E-01 1.538E 01 8. 2«3E-01 8.000E-01
4 3.812E-01 3.060E 01 9 . « 6 6 E - 0 1 1.000E 00
5 •.719E-01 • . 7 ^ 8 E 01 6. • 8 2 E - 0 1 1.000E 00
e 7.^22E-01 8.^53E 01 • . « 6 « E - 0 1 1.000E 00
7 6.800E-01 9.811E 01 5. 631E-01 2.000E 00
e 3. U S E - 0 2 • . 7 6 5 E 01 8 . 2 ^ 0 8 - 0 1 2.000E 00
9 6.336E-0^ 8.297E 00 1.«76E 00 2.000E 00
10 1.260E-05 2.01tE 00 2.025E 00 2.300E 00
11 2.523E-07 2.6^0E-01 1. 1«3E 00 2.300E 00
ORGARIC PHASE
STAGE ORARIOB THORIOB HR03 FLOR OOT ITERA- TEBP
RO. (G/l) (G/L)
(") OF STAGE TIORS DEG C
1 0.0 0.0 9.0ieE-02 7.000E 00 7. • 0.0
2 2.440E 00 •.283E 00 1.338E-01 7.000E 00 6. 39.8
3 2 . 5 0 e E 00 6 . 9 e 5 E 00 1. 890E-01 7.000E 00 6. 39.7
« 2.799E 00 1.701E 01 1.836B-01 7.000E 00 5. 39.2
5 2.798E 00 2.120E 01 1. 105E-01 7.000E 00 5. 38.9
e7 2.911E 00 3.140E 01 5.7218-02 7.000E 00 37.8
2-957E 00 3.799E 01 6.7^9E-0 2 7.000E 00 37.1
8 1.985E-01 2.^95E 01 1.«26E-01 7.OOOE 00 36.0
9 7.516E-03 8.201E 00 3.^«9E-01 7. OOOE 00 3«.7
10 2.280E-04 3.272E 00 «. 8«0E-01 7. OOOE 00 31.1
11 2.573E-06 1.866E-01 2. 825E-01 7.OOOE 00 28.6
45
AQOEOOS PHASI
TIBE STAGE ORARIOB THORIO n HR03 FLOR OOT
IHTERTAL RO. (GA) (GA) (H) OF STAGE
7 1 6.146E-01 1 . 5 4 « E 01 4.225E-01 8.OOOE-01
2 4.009E-01 1 . 4 5 6 1 01 7. 1 5 3 1 - 0 1 8.000E-01
3 3.869E-01 3 . 0 4 1 E 01 8.1071-01 8.000E-01
3.8S5E-01 3 . 1 4 7 E 01 1.046E 00 1 . 0 0 0 1 00
•5 6.275E-01 e . 2 8 0 E 01 6.6»4E-01 1.000E 00
6 7.367E-01 8 . 4 0 1 E 01 5.952E-01 I.OOOE 0 0
7 5.906E-01 1 . 1 1 2 E 02 4.5898-01 2.OOOE 0 0
8 2.854E-02 4 . 0 5 4 E 01 1.094E 00 2 . 0 0 0 1 00
9 7.908E-04 1 . 4 0 0 E 01 1.439E 00 2 . 0 0 0 1 00
10 &,798E-06 1. 172E 00 2 . 2 2 5 1 00 2.300E 00
11 3.760E-07 6.2161-01 1 . 1 5 9 E 00 2 . 3 0 0 1 00
ORGARIC PHASE
STAGE ORARIOB THORIOB RH03 FLOH OOT ITERA- TEHP
RO. OF STAGE TIORS DEG C
1
(GA)
2 . 3 7 0 E 00
(GA)
2 . 5 1 9 1 00 («)
8.556E-02 7 . 0 0 0 E 00 12. 39.8
2 2 . 4 6 0 E 00 5 . 3 0 1 1 00 1.630E-01 7.OOOE 0 0 6. 39.8
3 2 . 8 0 7 E 00 1 . 5 2 5 1 01 1.601E-01 7.000E 00 6. 39.3
4 2 . 7 8 4 E 00 1.8^6E 01 1. 9 8 ' W - O I 7.OOOE 0 0 5, 39.1
5 2 . 8 7 6 E 00 2 . 6 8 0 E 01 9.7531-02 7.000E 00 5. 38.2
2 . 9 1 9 E 00 3 . 2 7 7 E 01 7.507E-02 7 . 0 0 0 1 00 37.6
6
7 2 . 9 9 3 1 00 • . 3 1 2 1 01 S.386E-02 7 . OOOE 0 0 •.
5. 36.2
8 1.938E-01 2 . 4 e 5 E 01 \. 9311-01 7 . 0 0 0 1 00 5. 35.8
9 8.991E-03 1 . 2 8 9 E 01, 3 . 0 8 2 1 - 0 1 7.OOOE 0 0 5. 33.3
10 1.807E-04 2., 1 7 2 E 0 0 5.3011-01 7.000E 00 30.9
11 4.015E-06 4.575E-01 2.8»»1-01 7.000E 00 «.
3. 28.0
AQOEOOS PHASI
TIBE STAGE ORARIOB THORIOB HR03 FLOR OOT
IHTERTAL HO. (G/L) OF STAGE
8 1
(G/l)
5.203E-01 1 . 5 8 8 E 01 («)
•. 999B-01 8.000E-01
2 4.061E-01 3 . 0 9 » E 01 6.6591-01 8.0001-01
3 3.811E-01 3 . 0 7 3 E 01 9.051E-01 8.000E-01
4 4.978E-01 « . 3 7 8 E 01 1.085E 00 I.OOOE 0 0
5 6.550E-01 6 . 3 7 1 1 01 7.8801-01 I.OOOE 0 0
6 6.805E-01 9 . 7 9 ^ E 01 5.672E-01 I.OOOE 0 0
7 5.804E-01 1.082E 02 6. 3 8 0 E - 0 1 2 . OOOE 0 0
8 3.079E-02 5 . ^ 5 0 E 01 9.775E-01 2 . 0 0 0 1 00
9 6.443E-04 9 . 9 1 6 E 00 1.666E 00 2 . 0 0 0 1 00
10 1.138E-05 2.092E 00 2.2»1E 00 2.300E 00
11 2.339E-07 3.207E-01 1. 2 e e E 0 0 2.300E 00
ORGARIC PHASE
STAGE ORARIOB THORIOB HR03 FLOR OOT ITERA- TEBP
RO. (G/L) (G/L) (H) OF STAGE TIORS DEG C
1 2 . S 0 1 E 00 3 . « 8 6 E 00 1. 0 5 8 E - 0 1 7.OOOE 0 0 8. 39.8
2 2.831E 00 1 . 3 4 8 E 01 1.323E-01 7.000E 00 5. 39.5
3 2 . 7 8 6 E 00 1 . 6 6 2 1 01 1.76«E-01 7.OOOE 0 0 4. 39.3
4 2 . 8 4 9 E 00 2 . 4 0 3 1 01 1.781E-01 7.000E 00 5. 38.6
5 2 . 8 8 1 E 00 2 . 8 1 6 1 01 1.1201-01 7.OOOE 0 0 4. 38.1
2.985E 00 3 . 8 1 0 1 01 6. 7 7 ^ 1 - 0 2 7 . 0 0 0 1 00 5. 36.8
«
7 2.990E 00 • . 3 6 0 1 01 7.••01-02 7.000E 00 4. 36.0
8 1.690E-01 2 . 9 0 8 E 01 1,6 00B-01 7.OOOE 0 0 5. 34.6
9 8.151E-03 1 . 0 9 2 E 01 3.667B-01 7,OOOE 0 0 4. 33.1
10 2. 2 6 2 I - 0 * 3 . 7 7 0 E 00 5.167E-01 7 , 0 0 0 1 00 3. 30.1
11 2„81Ê-06 2.796E-01 3. 1 5 1 E - 0 1 7 , 0 0 0 1 00 4. 27.9
AQOEOOS PHASE
TIBE STAGE ORARIOB THORIOB HR03 FLOR OOT
IHTERTAL RO. (GA) (GA) (H) OF STAGE
9 1 4.5»»B-01 3 - 1 5 5 E 01 4.265E-01 8.OOOE-01
2 3.930E-01 3 . 1 2 2 E 01 7. 5 0 0 E - 0 1 8.000E-01
3 4.765E-01 « . « 0 6 E 01 8.961E-01 8.OOOE-01
4 5.287E-01 • . 5 2 4 1 01 1 . 1 9 0 1 00 I.OOOE 0 0
5 7.304E-01 7 . 8 7 2 E 01 8.300E-01 I.OOOE 0 0
6 6.6^5E-01 9 . 6 3 4 E 01 7.155E-01 I.OOOE 0 0
7 5. 100E-01 1 . 1 7 6 E 02 5.474E-01 2.OOOE 0 0
8 2.816E-02 4 . 7 2 2 E 01 1. 215E 0 0 2.000E 00
9 7.5631-04 1 . 5 0 0 E 01 1.625E 00 2.000E 00
10 8.291E-06 1 . 3 3 8 1 00 2.414E 00 2.3001 00
11 2.942E-07 5.630E-01 1.2781 00 2 . 3 0 0 1 00
ORGARIC PHASE
STAGE ORARIOB TRORIOB HR03 FLOH OOT ITEHA- TIBP
HO. (GA) (G/L) ("I OF STAGE TIOHS DEG C
1 2 . 7 7 9 1 00 9 . 8 7 5 1 00 8.355E-02 7.OOOE 0 0 5. 39.6
2 2 . OOOE 0 0 1„«86E 01 1.078E-01 7.000E 00 5. 39. •
3 2 . 8 4 1 E 00 2 . 2 5 3 E 01 1.517E-01 7,OOOE 0 0 5. •38.8
4 2 . 8 4 9 E 00 2 . 5 2 1 E 01 1.8831-01 7 , 0 0 0 1 00 3. 38.5
5 2 . 9 5 2 E 00 3 . 3 1 1 E 01 1.0«21-01 7.000E 00 4. 37.3
6 2..989E 00 3 . 8 9 4 E 01 8. • 7 6 E - 0 2 7 , 0 0 0 1 00 5. 36.6
7 2 . 9 7 8 E 00 4 . 7 3 2 E 01 6.325E-02 7.000E 00 4. 35.2
8 1.659E-01 2 . 8 3 4 1 01 2.018E-01 7 . 0 0 0 1 00 5. 34.5
9 8.806E-03 1 . 5 0 5 E 01 3.318E-01 7.OOOE 0 0 4. 32.1
10 1.842E-04 2 . 6 7 3 E 00 5.553E-01 7.OOOE 0 0 29.9
11 3.642E-06 5,0251-01 3. 162E-01 7,OOOE 0 0 «.
3. 27.5
46
AQOEOOS PRASE
TIBl STAGl ORARIOR TBORIOB HH03 FLOR OOT
IHTERTAL RO. (GA) (GA) (B) OF STAGE
10 1 •.313E-01 3 . 2 3 4 E 01 •.766B-01 8.0008-01
2 4.6911-01 4 . 5 4 1 E 01 6. 9 9 1 E - 0 1 8.OOOE-01
3 •.9501-01 • . 5 0 0 1 01 9.804E-01 8.0008-01
4 6.7011-01 5 . 8 7 1 1 01 1 , 251E 0 0 I.OOOE 0 0
5 7,2151-01 7 . 8 7 5 8 01 9.562E-01 I.OOOE 00
6 5.8261-01 1 . 0 6 2 1 02 7,017E-01 I.OOOE 0 0
7 5.0731-01 1 . 1 • • E 02 7.055E-01 2.OOOE 0 0
8 2,8761-02 5.928E 01 1 . 102E 0 0 2,OOOE 0 0
9 6.287E-0^ 1.116E 01 1.822E 00 2.OOOE 0 0
10 1.021E-05 2.115E 00 2.420E 00 2 , 3 0 0 1 00
11 1.957E-07 3.2701-01 1.373E 00 2 . 3 0 0 1 00
ORGARIC PHASE
STAGE DRARIDH THORIOB HR03 FLOH OOT ITERA- TEHP
HO. (GA) (GA) (8) OF STAGE TIORS DEG C
1 2 . 7 5 1 E 00 1 . 1 1 7 E 01 9.334E-02 7 . 0 0 0 E 00 39.6
2 2 . 8 « 0 E 00 2.09Ê 01 1. 2061-01 7 . 0 0 0 E 00 •
5 39.0
3 2 , 8 3 7 E 00 2 . 3 6 ^ E 01 1.6201-01 7 . 0 0 0 8 00 38.7
4 2 . 9 1 1 E 00 2.976E 01 1.7071-01 7 , 0 0 0 1 00 •
5 37.8
5 2 . 9 6 1 E 00 3 . » 1 5 E 01 1.181E-01 7.OOOE 00
6 2 . 9 9 9 E 00 « . 3 ^ 0 E 01 8. 158E-02 7 . 0 0 0 E 00 •5 37.2
35.8
7 2 . 9 7 3 1 00 « . 7 2 8 E 01 8.1048-02 7 . OOOE 0 0 35.1
8 1. •631-01 3 . 1 7 3 1 01 1.735E-01 7 . 0 0 0 1 00 •
5 33.5
9 8,0511-03 1.298E 01 3.819E-01 7.OOOE 0 0 5 32.0
10 2.16^1-0^ • . 0 9 3 1 00 5.424E-01 7 . 0 0 0 E 00 3 29.3
11 2.660E-06 3.3201-01 3.420E-01 7 . 0 0 0 E 00 3 27.4
AQOEOOS PHASE
TIHE STAGE ORARIOB TBORIOB HR03 FLOR OOT
IHTERTAL HO. OF STAGE
11 1
(GA)
•.4961-01
(GA)
(")
• . 3 8 0 E 01 4 . 1 0 0 1 - 0 1 8.000E-01
2 4.806E-01 • . 6 3 1 E 01 7.631E-01 8.0008-01
3 6.4151-01 5 . 9 1 7 E 01 9.7251-01 8.000E-01
4 6.8661-01 6 . 0 2 0 F 01 1 . 357E 0 0 I.OOOE 0 0
5 6.5651-01 8 . 9 9 8 E 01 9.987E-01 1 . 0 0 0 F 00
6 5.7221-01 1 . 0 3 9 E 02 8.4211-01 1 . 0 0 0 1 00
7 •.5821-01 1 . 2 1 6 E 02 6.3131-01 2 . 0 0 0 1 00
8 2.6871-02 5.212E 01 1 . 310E 0 0 2 , 0 0 0 1 00
9 7.0631-04 1 . 5 6 4 E 01 1.775E 00 2 . 0 0 0 1 00
10 7.6911-06 1 . 4 4 9 E 00 2 . 5 7 0 1 00 2.3001 00
11 2.3751-07 5.138E-01 1.379E 00 2 . 3 0 0 1 00
ORGARIC PHASE
STAGE ORARIOB THORIOB HR03 FLOR OOT ITERA- TEHP
RO. (GA) (GA) (H) OF STAGl TIORS DIG C
1 2 . 7 8 8 1 00 1 . 5 9 4 E 01 7.3721-02 7 . 0 0 0 1 00 39.2
2 2 . 8 3 1 1 00 2.20Ê 01 1.2931-01 7 . 0 0 0 1 00 •
4 38.9
3 2 . 8 9 1 1 00 2 . 8 1 9 E 01 1.3951-01 7 . 0 0 0 1 00 5 38.1
4 2 . 9 2 0 E 00 3 . 0 7 0 E 01 1.791E-01 7.OOOE 0 0 3 37.7
5 3 . 0 0 1 E 00 3 . 8 9 3 E 01 1. 176B-01 7,00OE 00 4 36.4
6 2 . 9 9 4 E 00 • . 3 6 9 E 01 9.734E-02 7 . 0 0 0 E 00 5 35.7
7 2 . 9 5 6 1 00 5 . 0 0 1 E 01 7. 195E-02 7.OOOE 0 0 4. 34.4
8 1.453E-01 3 . 0 7 6 E 01 2.090E-01 7.000E 00 5 33.4
9 8. 2 3 2 1 - 0 3 1.6S6E 01 3. 5 0 1 E - 0 1 7,OOOE 0 0 4 31.2
10 1.7971-04 3 . 0 ^ 6 E 00 5,754E-01 7 . OOOE 00 4. 29.2
11 3.2771-06 5.263E-01 3. 4 1 8 8 - 0 1 7.OOOE 0 0 2 27.1
AQOIOOS P H A S I
TIBE STAGE ORARIOB TBORIOH HH03 FLOR OOT
IHTERTAL RO. (G/l) (G/l) (B) OF STAGE
12 1 4.4941-01 • . • 7 1 1 01 •.•19E-01 8 . OOOE-01
2 6.098E-01 6 . 0 U I 01 7.06Ê-01 8,0001-01
3 6.647E-01 6 . 0 5 3 1 01 1,050E 00 8.0001-01
4 6.8791-01 7 . 3 1 2 1 01 1.411B 00 I.OOOE 00
5 6.^02E-01 8 . 9 1 2 1 01 1.121E 00 I.OOOE 0 0
6 5.087E-O1 1 . 1 1 2 1 02 8.3468-01 1 . 0 0 0 1 00
7 •.596E-01 1 . 1 8 2 1 02 7. 7 1 4 8 - 0 1 2 . 0 0 0 1 00
8 2.647E-02 6 . 2 1 8 1 01 1 . 2 0 9 1 00 2 . 0 0 0 1 00
9 6.040E-04 1 . 2 0 7 1 01 1.9481 00 2.000E 00
10 9.130E-06 2 . 1 1 6 1 00 2.5668 00 2 . 3 0 0 1 00
11 1.660E-07 3.2521-01 1 . 4 6 3 1 00 2 . 3 0 0 1 00
ORGARIC PBASE
STAGE ORARIOR THORIOB HH03 FLOH OOT ITIRA- TEHP
RO. (G/1) (G/L) (H) OF STAGE TIORS DEG C
1 2 . 7 8 0 E 00 1 . 6 9 3 1 01 7. 8 7 6 1 - 0 2 7.OOOE 0 0 4 39.2
2 2 . 8 7 3 E 00 2 . 6 3 2 1 01 1.056B-01 7 . 0 0 0 1 00 5 38.4
3 2 . 8 9 9 1 00 2 . 9 0 7 1 01 1 . 46 3 8 - 0 1 7.OOOE 0 0 3 38.0
4 2 . 9 7 6 1 00 3 . 5 2 5 1 01 1.701E-01 7.OOOE 00 5 37,0
5 3 . 0 0 1 E 00 3 . 9 5 6 E 01 1. 311E-01 7.OOOE 0 0 5 36.3
6 2 . 9 7 7 E 00 4 . 6 9 8 E 01 9.5^0B-02 7 . 0 0 0 E 00 5 35.0
7 2 . 9 5 3 E 00 4 . 9 6 9 F 01 8.751E-02 7 . 0 0 0 E 00 4 3^.3
8 1.316E-01 3 . 3 5 » E 01 1.851E-01 7 . 0 0 0 F 00 5 32.7
9 7.683E-03 1.449E 01 3.933E-01 7 . 0 0 0 1 00 5 31.1
10 2.021E-0^ « . 3 0 0 E 00 5.6 288-01 7 . 0 0 0 1 00 3 28-7
11 2.472E-06 3.693E-01 3.638E-01 7 . 0 0 0 E 00 3 27.1
47
AOOEODS PHASE
Tim ST»GE OEA1II0B TBOKIBll BUDS flOB OUT
IIITE!IV»1 SO. (G/1) (G/L)
(") OF STAGE
13 1 S.0II5E-01 S.II10E 01 3.77IIE-01 8.000E-01
2 S.312E-01 6 . 1 I I 6 E 01 7.558E-01 8.000E-01
3 7.2II0E-01 7.518E 01 1.051E 00 8.0001-01
« 6.706E-01 7.356E 01 1.50«E 00 1.000E 00
5 5.630E-01 9 . 7 5 5 E 01 1. 153E 00 1.000E 00
«7 5.0»OE-01 1.086E 02 9.651E-01 1.000E 00
«.239E-01 1.2II1E 02 7. 103B-01 2.000E 00
8 2.527E-02 5.537E 01 1.393E 00 2.000E 00
9 «„533E-0a 1.S92E 01 1.900E 00 2 . 0 0 0 1 00
10 7.1«8E-0e 1.526E 00 2.697E 00 2 . 3 0 0 1 00
11 1.95eE-07 «.7«8E-01 1. I16«E 00 2 . 3 0 0 1 00
OIGAHIC BRASE
STkCE nUAEIOH TROEIDH HII03 PLOB OOT ITERA- TEHP
II0. «G/l) (G/M ("1 OP STAGE TIOBS DBG C
1 2.815E 00 2„01«B 01 6. 2I18E-02 7.000E 00 a. 38.8
2 2.878E 00 2.716E 01 1.105E-01 7.000E 00 3. 38.3
3 2.963E 00 3.353E 01 1. 308E-01 7.000E 00 II. 37. »
4 2 . 9 8 1 1 00 3.597E 01 1 . 7 9 2 E - 0 1 7 . 0 0 0 E 00 3. 36.9
5 2.995E 00 ».3119E 01 1. 322E-01 7.000E 00 5. 35.6
e 2.972E 00 II.691E 01 1 . 0 9 8 8 - 0 1 7.000B 00 5. 34.9
7 2.9aOE 00 5 . 1 8 1 E 01 7 . 9 9 6 E - 0 2 7.000E 00 a. 33.7
8 1.318E-01 3.2«»E 01 2 . 1 5 9 B - 0 1 7.000E 00 5. 32.5
9 7.579E-03 1.752E 01 3.65IIE-01 7.000E 00 5. 30.5
10 1.7271-01 3.318F 00 5 . 9 1 1 B - 0 1 7.000E 00 3. 28.6
11 2.93SE-06 5„39«E-01 3 . 6 2 3 E - 0 1 7.000E 00 2. 26.8
AQOEOOS PHASE
TIHE STAGE DRAIIIIH THOBIUB HII03 PLOa OUT
mTERTAl HO. (GA) (G/L) ("t OP STAGE
H 1 5.130E-01 5.509E 01 3 . 9 7 1 E - 0 1 8.0001-01
2 7.310E-01 7 . 5 6 1 E 01 7. 106E-01 8.0001-01
3 7.160E-01 7.596E 01 1.120B 00 8.0001-01
II 6.005E-01 8.35IIE 01 1.S30E 00 1.000E 00
5 S.517E-01 9 . 6 08E 01 1.263E 00 1.000E 00
6 II.571E-01 I.IHllE 02 9 . 5 7 9 E - 0 1 1.000E 00
7 l|_275E-01 1.207E 02 8 . 3 5 1 B - 0 1 2.000E 00
8 2.II38E-02 6.362E 01 1. 30I1E 00 2 . 0 0 0 1 00
9 5.7««E-0« 1.265E 01 2.0S1E 00 2 . 0 0 0 1 00
10 8.172E-06 2.090E 00 2 . e 8 e E 00 2 . 3001 00
11 1. llllOE-07 3..207E-01 1.537E 00 2.300E 00
ORGAHIC PHASE
STtSE ORAHIOH THOHIOH BROS PLOB OOT ITIBA- TEBP
HO. (G/L) (G/l) (B) OP STAGE TIOBS DEG C
1 2 . 8 1 9 E 00 2.08eE 01 6.510E-02 7.000E 00 4. 38.7
2 2 . 9 3 7 E 00 3.107E 01 9. 268E-02 7 . 0 0 0 1 00 5. 37.7
3 2.971E 00 3.11311 01 1.375E-01 7.000E 00 3. 37.3
U 2.992E 00 • . 0 1 5 E 01 1.7651-01 7 . 0 0 0 1 00 II. 36.3
5 2.989E 00 II.369E 01 1.11II1E-01 7 . 0 0 0 1 00 5. 35.5
6 2.955E 00 «.9ll0E 01 1.079E-01 7 . 0 0 0 1 00 4. 34.3
7 2.939E 00 5.133E 01 9.3701-02 7 . 0 0 0 1 00 4. 33.6
8 1-217E-01 3.II81E 01 1. 958E-01 7 . 0 0 0 1 00 a. 32.0
9 7. 227E-03 1.552E 01 •1.0311-01 7 . 0 0 0 1 00 5. 30.4
10 1.869E-0II e.llOOE 00 5.7951-01 7 . 0 0 0 1 00 3. 28.3
11 2.301E-Oe 3.960E-01 3.8111-01 7 . 0 0 0 1 00 3. 26.8
AQOEOnS PHASE
TIHE STAGE OBAIIIDH THOltOa HB03 PlOB OOT
IHTEBTAl HO. (G/l( (G/l) (") OP STAGE
15 1 5.937E-01 6.1I38E 01 3. »78B-01 8.aO0E-O1
2 7.356E-01 7 . 6 8 5 E 01 7. 5 3 5 1 - 0 1 8.000E-01
3 6.555E-01 8.76«E 01 1.112B 00 8.000E-01
» 5.852E-01 8.313E 01 1. eiOE 00 1.000E 00
5 11.925E-01 1.029E 02 1.278E 00 1.000E 00
e I1.570E-01 1. 117E 02 1 . 0 7 5 1 00 1 . 0 0 0 1 00
7 3.998E-01 1.258E 02 7 . 8 3 7 1 - 0 1 2 . 0 0 0 1 00
8 2.373E-02 5.733E 01 1. I167E 00 2 . 0 0 0 1 00
9 e.028E-0« 1.588E 01 2 . 0 0 5 1 00 2 . 0 0 0 1 00
10 6.635E-06 1.570E 00 2. 800E 00 2.300E 00
11 1.611 HE-07 4.405E-01 1.5311B 00 2.300E 00
OSGAHIC PHASE
STAGE nBAmon TBOKID n HB03 PLOB OOT ITEBA- TIHP
HO. (G/ll (G/l) (B) OP STAGE TIOBS DEG C
1 2.869E 00 2.371E 01 5 . 2 9 3 E - 0 2 7.000E 00 4. 38.3
2 2.9«6E 00 3 . 182E 01 9 . 6 7 2 1 - 0 2 7.000E 00 3. 37.6
3 3.000E 00 3.877E 01 1.306E-01 7 . 0 0 0 1 00 4. 36.6
a 2 . 9 8 8 E 00 I1.050E 01 1 . 8 5 0 E - 0 1 7.000E 00 3. 36.1
5 2.971E 00 a.eetE oi 1. U39E-01 7 . 0 0 0 1 00 a. 34.8
6 2.952E 00 •I.910E 01 1 . 2 0 6 1 - 0 1 7 . 0 0 0 1 00 5. 34.2
7 2-930E 00 5 . 3 0 7 E OK 8 . 7 2 7 E - 0 2 7 . 0 0 0 1 00 a. 33.2
8 1.226E-01 3.362E 01 2 . 2 2 6 1 - 0 1 7 . 0 0 0 1 00 a. 31.9
9 6.981E-03 1.80IIE 01 3 . 7 9 1 1 - 0 1 7 . 0 0 0 1 00 5. 29.9
10 1.6«2E-0'I 3.ll9eE 00 6.OUOl-O1 7 . 0 0 0 1 00 3. 28.2
11 2.631E-06 5.U19E-01 3. 7 8 7 1 - 0 1 7.000E 00 3. 26.6
48
AQDEODS PBASE
TIBE STAGE DRAIIOa TBOBIDB BB03 FLOI ODT

IBTERT AL 10. (6/U (G/L) (B) OF STAGE
16 1 6.0701-01 6.506E 01 3.655B-01 8.0001-01
7.1361-01 8 . 9 0 5 1 01 7. 205E-01 e.OOOE-01
6.019E-01 8.760E 01 1.172B 00 8.000E-01
5.208E-01 9.082E 01 1.618E 00 I.OOOE 00
0.882E-01 1.010E 02 1.370E 00 I.OOOE 00
0.219E-01 1. 169E.02 1.063B 00 I.OOOE 00
0.051E-01 1.223S 02 8.938E-01 2.000E 00
2.268E-02 6 . 0 1 1 E 01 1. 3911 00 2,000E 00
5.039E-00 1.290E 01 2 . 1 3 8 1 00 2 , 0 0 0 1 00
7.3351-06 2.037E 00 2. 786E 00 2.300E 00
1. 2681-07 3.138E-01 1.599E 00 2.3aOE 00
ORGABIC PBASI
STAGE DBABIDB THORIOH RI03 FLOI ODT ITEBA- TEBP
10. (GA) (GA) (HI OF STAGE TIOBS DEG C
2 . 8 7 7 1 00 2.030E 01 5 . 0 9 5 1 - 0 2 7,O0OE 00 3, 38.2
2 . 9 8 7 1 00 3.S96I 01 8 . 7 5 9 1 - 0 2 7.000E 00 37.0
2.998E 00 3.9271 01 1 . 3 7 2 1 - 0 1 7 . 0 0 0 1 00 36.5
2.971E 00 0.3681 01 1 . 8 2 7 E - 0 1 7 . 0 0 0 1 00 35.6
2.966E 00 0.6501 01 1.503E-01 7 . 0 0 0 1 00 30.7
2.900E 00 5.1071 OH 1 . 1 8 0 E - 0 1 7.000E 00 33.7
2.929E 00 5.2081 OH 9 . 9 3 2 E - 0 2 7.000E 00 33.1
1.107E-01 3.567E 01 2 . 0 5 7 E - 0 1 7 . 0 0 0 1 00 31.0
6.788E-03 1.618E 01 0 . 1 2 0 E - 0 1 7.000E 00 29.9
1.720E-00 0.0091 00 5 , 9 3 5 B - 0 1 7.000E 00 27.9
2.1381-06 0.1271-01 3.90SB-01 7 . 0 0 0 1 00 26.6
AQOIODS PRASE
TIBE STAGE DRABIDB TBOBIDB HB03 FLOI ODT
IBTEBTAL BO. (G/L) (6/1) (R) OF STAGE
17 1 6.9181-01 7 . 5 1 8 1 01 3.060E-01 8.000E-01
2 7.0391-01 8.909E 01 7,626E-01 8.000B-01
3 5.736E-01 9 . 6 2 5 E 01 1.158E 00 B.OOOE-OI
0 5.1601-01 8.987E 01 1.687S 00 I.OOOE 00
5 0.058E-01 1.068E 02 1 . 3 7 6 1 00 I.OOOE 00
6 0.202E-01 1.1OOE 02 1.165E 00 I.OOOE 00
7 3.8261-01 1.269E 02 8.0871-01 2.000E 00
8 2.2351-02 5.833E 01 1. 530E 00 2.000E 00
9 5.591E-00 1.561E 01 2.096E 00 2.000E 00
10 6.1581-06 1.583E 00 2.S85E 00 2.300E 00
11 1.0001-07 0,093E-01 1. 592E 00 2.300E 00
0B6AIIC PHASE
STA6I DBABIDH TBOBIDB HB03 FLOW OOT ITERA- TEHP
BO. (G/L) (G/L) (B) OF STAGE TIOIS DIG C
1 2 . 9 0 8 1 00 2.736E 01 0. 8 0 0 1 - 0 2 7.0001 00 0. 37.8
2 2.987E 00 3.652E 01 9 . 1 7 8 B - 0 2 7.0001 00 36.9
3 2.987E 00 0.286E 01 1 . 3 3 2 E - 0 1 7.0001 00 35.9
0 2 . 9 6 6 1 00 0.372E 01 1.902E-01 7.0001 00 35.5
5 2.951E 00 0.878E Olf 1, 525E-01 7.0001 00 30.2
6 2.93SE 00 5.063E 0 1 * 1. 291E-01 7.000E 00 33.6
7 2.923E 00 5.396E O l t 9 . 3 6 3 B - 0 2 7.000E 00 32.7
8 1. 161E-01 3.006E 01 2 . 2 9 3 E - 0 1 7.000E 00 31.3
9 6.093E-03 1.8271 01 3 . 9 1 9 E - 0 1 7.000E 00 29.5
10 1.555E-00 3 . 5 9 1 1 00 6. 108E-01 7.000E 00 3 27.9
11 2.360E-06 5.309E-01 3.921B-01 7.000E 00 3 26.0
AQDEODS PBASE
TIRI STA6I DBABIDB TRORIDB BB03 FLOI OOT
IBTIRTAL 10. (6/L) (6/L) (R) OF STAGE
18 1 6. 9 82 E-01 7.593E 01 1.605E-01 8.OOOE-01
2 6.012E-01 9.802E 01 7.008E-01 8.0001-01
3 5.657E-01 9.560E 01 1. 210E 00 8.000E-01
0 0.730E-01 9.600E 01 1.683B 00 I.OOOE 00
5 0.051E-01 1.008E 02 1.060E 00 I.OOOE 00
6 3.969E-01 1 . i e 7 E 02 1.109E 00 l.OOOB 00
7 3.887E-01 1.235E 02 9.056E-01 2.000E 00
8 2.135E-02 6.001E 01 1.068E 00 2.000E 00
9 5.100E-00 1.300E 01 2.212E 00 2-OOOE 00
10 6.630E-06 1.966E 00 2.869E 00 2.300E 00
11 1.129E-07 3.003E-01 1.650E 00 2.300E 00
0B6AIIC PRASE
STAGE DRAHIDR TBOBIDB HB03 FLOW ODT ITERA- TEHP
10. (G/L) (G/L) (B) OF STAGE TIOIS DEG C
1 2.907E 00 2 . 7 8 0 1 01 5 . 0 1 2 B - 0 2 7.0001 00 3. 37.7
2 2.993E 00 0 , 0 2 0 1 01 8 . 7 6 5 1 - 0 2 7.0001 00 5. 36.0
3 2.982E 00 0.302E 01 1. 390B-01 7.0001 00 3. 35.8
0 2.909E 00 0.607E OH 1.873E-01 7.000T 00 35.0
5 2.909E 00 0.850E OH 1.616E-01 7.000E 00 30.1
6 2.930E 00 5 , 2 2 6 1 OH 1 . 2 6 0 1 - 0 1 7.000E 00 33.2
7 2.923E 00 5 . 3 3 2 1 OH 1 . 0 0 2 1 - 0 1 7.000E 00 32.7
8 1.097E-01 3 . 6 2 0 1 01 2 . 1 0 9 1 - 0 1 7.000E 00 31.0
9 6.395E-03 1.6501 01 0 . 2 1 2 E - 0 1 7,000E 00 29.0
10 1.599E-00 0.350E 00 6 . 0 5 5 E - 0 1 7.000E 00 3I 27.6
11 1.986E-06 O.203T-01 0 . 0 5 9 E - 0 1 7.000E 00 3 26.0
49
TIMF INTERVALS 19 THROUGH 112 CHITTED TO SAVE SPACE
AQDEODS PHASE
TIBE STAGE DBABIDB THOBIDH RB03 FLOW ODT
IBTEBTAL BO. (G/l) (6/1) (R) OF STAGE
113 1 6.3111-01 9 . 8 5 3 1 01 0. 250E-01 8.000S-01
2 0.5701-01 1 . 1 2 6 1 02 8.807E-01 8.000E-01
3 0.0091-01 1.089E 02 1. 392E 0 0 8.000E-01
0 3.725E-01 1 . 0 2 1 1 02 1.979E 00 I.OOOE 0 0
5 3.517E-01 1.0991 02 1. 817E 0 0 I.OOOE 0 0
6 3.091E-01 1.108E 02 1.590B 00 I.OOOE 0 0
7 3.569E-01 1.195E 02 1. 260S 0 0 2.000E 00
8 1.755E-02 0 . 9 6 1 1 01 1.981B 00 2 . 0 0 0 1 00
9 3.506E-00 9 . 1 3 3 1 00 2.e22E 00 2.000E 00
10 1.2311-06 1.017E 00 3.308S 00 2.300E 00
11 0.905E-08 1.602E-01 1.90BI 00 2.300E 00
0 R 6 A B I C PHASE
STAGE DBABIDB TBOBIDB BI03 FLOW ODT ITIRA- TEHP
BO. (6/L) (G/L) (R) OF STAGI TIOBS DIG C
1 2 . 8 5 7 E 00 3.7B0E 01 5.121B-02 7.D00E 00 1 30.6
2 2 . 9 2 9 E 00 0 . 9 1 0 E 01 9,9808-02 7.000E 00 1 32.9
3 2 . 9 0 9 E 00 5 . 0 7 0 1 0 1 * 1. 5 2 3 1 - 0 1 7 . 0 0 0 E 00 1 32.3
0 2 . 9 0 3 E 00 5.0281 OH 2.100S-01 7.000E 00 1 32.1
5 2 . 9 1 0 E 00 5.2021 O H 1,911E-01 7.000E 00 1 31.0
6 2 . 9 0 7 E 00 5.3531 OH 1.680E-01 7.000E 00 1 31.1
7 2 . 9 0 7 E 00 5 . 0 2 2 1 O H 1. 3 5 6 E - 0 1 7.000E 00 1 30.9
8 1.020E-01 3 . 0 1 0 1 01 2.911E-01 7.000E 00 1 29.2
9 5.0101-03 1 . 0 1 2 1 01 0.961E-01 7.000E 00 1 27.8
10 1.0021-00 2 . 5 5 7 1 00 6 . 7 9 3 1 - 0 1 7.000E 00 1 26.6
11 1.0061-06 2.8031-01 0.5991-01 7 . 0 0 0 1 00 1 25.8
AODEODS PHASE
TIRE STAGE DRABIOR TBORIDB BB03 FLOW ODT
IBTEBTAL 10. (6A) (GA) (B) OF STAGE
110 1 6.3111-01 9 . 8 5 3 E 01 0.251E-01 8.0001-01
0.5701-01 1.125E 02 8. 807E-01 8.0001-01
0.050E-01 1 . 0 8 9 E 02 1.393E 00 8.000E-01
3.7251-01 1,021E 02 1.979E 00 I.OOOE 0 0
3.5171-01 1 . 0 9 9 1 02 1.8171 00 I.OOOE 0 0
3.091E-01 1.108E 02 1. 590E 0 0 I.OOOE 0 0
3.570E-01 1.195E 02 1.260E 00 2.000E 00
1.755E-02 0 . 9 5 9 E 01 1. 981E 0 0 2.000E 00
3.5061-00 9 . 1 2 9 1 00 2.622E 00 2 . 0 0 0 E 00
10 3.2301-06 1.016E 00 3 . 308E 0 0 2.300E 00
11 0.9001-08 1.6021-01 1.908E 00 2.300E 00
OBGABIC PRASE
STAGE DBABIDH TRORIDB HBa3 FLOW OOT ITIBA- TIRP
BO. (6A) (GA) (H) OF STAGI TIOBS DIG
2.857E 00 3 . 7 8 0 1 01 5. 122E-02 7.000E 00 1 30.6
2 . 9 2 9 1 00 0 . 9 1 0 B 01 9.9791-02 7.000E 00 1 32.9
2 . 9 0 9 1 00 5.0701 OH 1,520E-01 7.000E 00 1 32.3
2 . 9 0 3 1 00 5.028E OH 2.100E-01 7.000E 00 1 32.1
2.910E 00 5.202E OH 1.9111-01 7.000E 00 1 31.0
2.907E 00 5.352E OH 1.680E-01 7.000E 00 1 31.1
2,907E 00 5.022E 011 1.356E-01 7.000E 00 1 30.9
1.0 201-01 3 . 0 1 0 E 01 2.911B-01 7.000E 00 1 29.2
5.010E-03 1 . 0 1 2 E 01 0. 962E-0 1 7.000E 00 1 27.8
10 1.002E-00 2 . 5 5 6 E 00 6,793S-01 7 . 0 0 0 1 00 1 26.6
11 1.005E-06 2.802E-01 0. 599E-01 7.000E 00 1 25.8
# IBDICATES THIRD PRASE WITHIB 10« TOLIRABCE

99.90 PIRCEBT STEADY STATI BY OTIBALL RATBRIAL BALAICI
90.0 PIRCEIT OF STEADT STATI AFTER TIRE IITEBTAL 23
95.0 PERCEIT OF SHADY STATI AFTIR TlHI IBTIRTAL 30
99.0 PERCEIT OF STEADY STATE AFTER TIRE IITISTAL 50
99.5 PERCEIT OF STEADY STATE AFTER TIHE IHTERTAL 76
TOTAL 6 A DRAIIDH IB: 2.0000E 0 1 , ODT: 1.9999E 01

TOTAL G/L THOBIDH I B : 2 . 6 5 0 0 E 0 2, OOT: 2 . 6 5 26E 02
TOTAL HOLES HI03 I I : 0.7500E 0 0 , OOT! 0 . 7 0 731 00
^
51
ORNL/CSD/TM-70
Dist. Category UC-79c
INTERNAL DISTRIBUTION
1. W. D. Burch 10-19. R. H. Rainey

2. H. P. Carter/A. A. Brooks/ 20-29. S. B. Watson
CSD Library 30-31. Central Research Library
3. W. L. Carter 32. Document Reference Section
4. D. A. Costanzo Y-12
5. W. S. Groenier 33-35. Laboratory Records
6. M. V. Keigan 36. Laboratory Records - Record
7. R. E. Leuze Copy
8. G. W. Morrison 37. ORNL Patent Office
9. K. J . Notz
EXTERNAL DISTRIBUTION
38. J . N. Rogers, Division 8324, Sandia Laboratories, Livermore,
California 94550.
39. Chief, Mathematics and Geoscience Branch, Department of Energy,
Mathematical and Computer Sciences Program, Molecular Sciences
and Energy Research, Division of Physical Research,
Washington, D.C. 20545.
40. Office of Assistant Manager for Energy Research and Development,
Department of Energy, Oak Ridge Operations, Oak Ridge,
Tennessee 37830.
41-281. Given distribution as shown in TID-4500 under UC-79c, Fuel
Recycle Category (Applied).

SEPHIS Operating Manual

Uploaded by

Copyright:

Available Formats

SEPHIS Operating Manual

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

SEPHIS Operating Manual

Uploaded by

Copyright:

Available Formats

^.

User's Guide to the SEPHIS Computer Code

This report was prepared as an account of work sponsored by an

Portions of this document may be illegible in

Contract No. W-7405-.eng-26

USER'S GUIDE TO THE SEPHIS COMPUTER CODE

S. B. Watson and R. H. Rainey*

(Sponsor: R. H. Rainey; O r i g i n a t o r : S. B. Watson)

This work was funded by

UNION CARBIDE CORPORATION, NUCLEAR DIVISION

*Oak Ridge Associated U n i v e r s i t i e s

S. B. Watson and R. H. Rainey*

The SEPHIS computer program was developed to

SEPHIS, a computer program for Solvent Extraction P^rocesses Having

interacting Solutes, was developed by W. S. Groenier^ for the purpose of

studying the application of a ^5% TBP solvent extraction system for

recovered irradiated LMFBR fuel. G. L. Richardson of HEDL has adapted

and has made valuable changes in program structure.^'^ The authors

have made additional modifications to the Richardson version of SEPHIS-

Laboratory (ORNL) uses.'* Research to improve SEPHIS-Purex is being

continued by A. D. Mitchell who has written a code entitled SEPHIS-M0D4.^»6

This improved model of the Purex solvent extraction system is better

able to simulate run conditions during transient periods than previous

versions of SEPHIS and also better represents the operation of mixer-

*Chemical Technology Division

The SEPHIS computer program is useful for obtaining estimates of the

number of theoretical stages varies with operating conditions and, there-

Many of the technical aspects of the SEPHIS-Thorex computations

II. PROGRAM DESCRIPTION

Thorex version parallel those of Groenier.^ in addition the calculational

bution coefficient are patterned after the technique of Groenier and

In general the calculation simulates the extraction of the three

interacting solutes by a cascade of mixer-settler units with 100%

be proportional to the feedstream flow rates or can be varied to simulate

each feedstream can be designated. Also, the composition of each phase

in the contactor at the beginning of the calculation can be stipulated,

that is, the initial profile of the system can be specified.

The simulation starts by computing the concentration of the solutes

as the basis of a second calculation, and this iteration procedure is

reaches a pre-determined material balance - normally 99.9%. This system

gives an estimate of the composition of each phase in each stage of the

The program SEPHIS-Thorex consists of six routines that are

programmed in a well-structured manner and are heavily augumented

with program comment cards. The basic methodology and mathematical

formulations have been described previously;^•'^ thus, we include in this

section only the diagram of program structure and a description of each

Figure 1 shows the overall structure of the SEPHIS-Thorex program.

Arrows indicate the calling sequence of routines. All routines are

written in the IBM System/360 and System/370 FORTRAN IV language. A

listing of SEPHIS-Thorex routines is given in Appendix A.

PROFIL FLOWS STAGES PLOTR

1. Read title information.

5. Write feed stream data.

6. Convert feed streams concentrations to molarity.

7. Call FLOWS to set interstage flow rates, define outgoing streams

9. Print initial profile in the stages of the contactor system.

10. Print phase ratios and stage flows.

11. Initialize total molarity coming in (SOLIN) and total molarity

12. Call STAGES to perform stage calculation.

13. Calculate percent of material balance which has been reached.