# '#' means comment line
# This is file for data input to run 'nbubble'
#
# Name of output file (must be first non-comment line)
#
rhy-st
#
#
# Information lines to be in your output file (up to 5 lines)
#
***start*of*info*input************************************************
RHYOLITE
rhy-st
Standard conditions
See table 2 and 4 of paper
***end*of*info*input**************************************************
#
# Input of precision parameters
#
# Dimensionless parameter "time" govern a computational time step for
# the integration. You can vary it to eliminate CPU time. Recommended
# range for "time" is between 1.0d-03 (slower) and 1.0d+00 (faster)
#
time=5.0d-02
#
# If you want to specify the nine dimentionless parameters, then change
# them all from 0.0d+00. Otherwise, leave them alone and the code will
# do it for you.
#
ypar(1)=0.0d+00
ypar(2)=0.0d+00
ypar(3)=0.0d+00
ypar(4)=0.0d+00
ypar(5)=0.0d+00
ypar(6)=0.0d+00
ypar(7)=0.0d+00
ypar(8)=0.0d+00
ypar(9)=0.0d+00
ypar(10)=0.0d+00
#
#
# If any of above parameters is 0.0d+00, then the following values
# will be taken. If you specified dimensionless parameters, only
# c0,temp0,wm,denmelt,enact,vcoef,dpower,tempvc,dTvc will be taken.
#
# It is possible to run a model of decompressive gas exsolution and
# bubble growth with a constant decompression (or magma rise) rate.
# To start it, make rise rate (vh) in m/s different from 0.0d+00.
# Initial depth of magma rise (depth0) in m must be specified, and
# final depth of magma rise is surface at atmospheric pressure.
# Initial concentration (c0) of the volatile is saturation
# concentration, if c0 is given below as 0.0d+00.
#
vh=0.0d+00
depth0=0.0d+00
#
# Initial volatile concentration (c0) in the melt
#
c0=1.0d-02
#
# (1) Initial bubble radius (r0) and sphere radius (s0) in m
#
r0=1.0d-05
s0=1.0d-03
#
# (2) Initial temperature (temp0) in Centigrade, pressure (pa) in MPa,
# molecular weight of gas (wm) in kg/mole, melt density (denmelt)
# in kg/m**3; and (3) surface tension (sigma) in N/m
#
temp0=1.0d+03
pa=0.1d+00
wm=18.0d-03
denmelt=2.2d+03
sigma=0.32d+00
#
# If vh=0.0d+00 and pa=0.0d+00, then "pa" will be calculated from
# the value of "depth0"
#
# (4) Viscosity formula for rhyolite is (Hess, Dingwell, 1996):
# log(eta)=[-3.545+0.833*ln(H2O)]+[9601-2368*ln(H2O)]/
# {T-[195.7+32.25*ln(H2O)]}
# - it is used if enact is set to 0.0d+00. Otherwise:
# eta=eta00*exp(Eact/BT)
# Eact=enact*(1-vcoef*con) and eta00=10**(-4.5)
# "enact" is about 304500 (rhyolite) or 150000 (basalt) in J/mole,
# "vcoef" (dimensionless) is about 11 (H2O in rhyolite) or 6 (basalt)
#
# Diffusion coefficient formula is
# lnD=ln(c)-dcoef-denact/RT
# "dcoef" is 12.574 and "denact" is 87300 in J/mole for rhyolite
# "dcoef" is 12.49 and "denact" is 15200 in J/mole for basalt
#
enact=0.0d+00
vcoef=11.0d+00
dcoef=12.574d+00
denact=87.3d+03
#
# (5) Henry's constant (hen) in 1/Pa; (6) temperature diffusivity (xi)
# in m**2/s; (7) heat capacity of melt (cm) in J/(kg*K); (8) heat
# capacity of gas (cp) in J/(kg*K)
#
hen=1.6d-11
xi=1.42d-07
cm=1.35d+03
cp=2.5d+03
#
# (9) If volatile is water (wm=18.0d-03), initial heat of gas
# vaporization (heat0) in J/mole is evaluated for initial pressure
# and then taken as pressure dependent.
# In this case "hea0" must be 0.0d+00. But if gas is not water or
# "heat0" is different from zero, then the "heat0" is assumed as
# pressure independent.
#
heat0=0.0d+00
#
# (10) Heat of vitrification or crystallization (VC) (heatvc) in J/kg,
# liquidus VC temperature (tempvc) in C, VC temperature interval (dTvc)
# in Centigrade. To avoid error dTvc must not be 0.0d+00
#
heatvc=7.15d+03
tempvc=8.02d+02
dTvc=5.0d+01
#
# If you want interactive messages on your monitor during computations,
# which is sometimes necessary to control precision of the computations,
# make "monitor=1.0d+00", otherwase make "monitor=0.0d+00"
#
monitor=1.0d+00
#
end