# '#' 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