INPUT PARAMETERS:	PARAMETER EXPLANATIONS:
1			icalc+
74.0 			initial pressure (MPa) (used only if icalc=1)
28.5			conduit pressure gradient, Mpa/km (used only if icalc=2)
2	   		iteration number++
1.0   	 		initial velocity (m/s)
2800. 			unvesiculated magma density (kg/m3)
1200. 			initial temperature (C)
0.27,  0.05,   0.07	initial h2o, co2, S content (wt%)
2			vesiculation parameter*
-3000.			initial depth (m)
10.0		        conduit diameter (m)**
0.0025			fo (wall rock roughness term)

NOTES ON INPUT PARAMETERS:
+icalc		= 1 if specifying conduit of constant x-sectional area, and
			having the program calculate the pressure profile,
		  2 if specifying a constant pressure gradient in the conduit,
			and having the program calculate x-sectional area.
++iteration #=2 if the velocity is to be adjusted automatically to reach
                sonic velocities at the exit (valid only if icalc=1), or
             1 if no adjustment is desired.
*vesiculation p.= 2 if gas exsolution is to stop after fragmentation
		   1 if not
**diameter	= diameter throughout the conduit if icalc=1
		  diameter at base of conduit if icalc=2

OUTPUT PARAMETERS:
List of variables to be written out. Enter a number in the
first column indicating the column # where this variable will be written
in the output file.  You can write out up to seven variables.
 	1  x-sectional area (m2)
5	2  Mach number
7	3  pressure (MPa)
 	4  Reynolds number
 	5  mixture density
2	6  time (seconds) since entering the conduit
4	7  velocity (m/s)
3	8  volume fraction gas
 	9  log viscosity (Pa s)
1	10  z (depth, meters)
 	11  dadz (change in x-sectional area with depth, m2/m)
 	12  log(pressure)
 	13  dpdz (pressure gradient, Pa/m)
 	14  dz   (vertical step size, meters)
 	15  exco2 (wt % exsolved co2)
 	16  exh2o (wt % exsolved h2o)
 	17  exsulfur (wt % exsolved s)
 	18  f      (friction factor)
 	19  gamma  (Cp/Cv for gas phase)
 	20  mf    (mass fraction exsolved gas)
 	21  mm    (mass fraction exsolved magma)
 	22  r     (Universal Gas const. * n)
 	23  rhof  (gas density)
  	24  sv    (sonic velocity (m/s)
 	25  xco2  (mole fraction co2 in gas)
 	26  xh2o  (mole fraction h2o in gas)
 	27  xsulfur  (mole fraction sulfur species in gas)
 	28  cph2o (sp. heat at const. p. of H2O gas, kJ/kmol K)
 	29  cpco2 ("     "  "    "    "  "  CO2  "  )
 	30  cps   ("     "  "    "    "  "  S species)
	31  cvh2o ("     "  "    "   v. "  H2O gas,  kJ/kmol K)
	32  cvco2 ("     "  "    "   "  "  CO2  " )
	33  cvs   ("     "  "    "    "  "  S species)
 	34  temperature (C)
6	35  enthalpy of mixture (kJ/kg)
 	36  cp (sp. heat at const. p) of gas phase (kJ/kg C)
        37  conduit radius (m)