2007 Journal
of Volcanology & Geothermal Research (JVGR)
Dear
Colleague:
Thank you for your interest in my 2007 Journal of Volcanology and Geothermal Research (JVGR) paper:
Resmini, R.G., (2007). Modeling of crystal size distributions (CSDs) in sills.
Journal of Volcanology and Geothermal Research, v. 161, pp.
118-130, doi: 10.1016/j.jvolgeores.2006.06.023.
The code used in the modeling, as well as input and output files, is provided
here. A single case is highlighted; that for a position within the sill
located 3 meters from the sill/wallrock contact.
The code, though “tight”, is somewhat rudimentary in
that input/output file names and model parameters are hardcoded. Thus,
for each run, a recompile is required (except for the BASIC program). I
wouldn’t write code like this today but this method was convenient in that it
served as a conduit for the documentation for each model run.
I will periodically update the documentation on this
page. Thus, this site is a work-in-progress. I am happy, however,
to answer your questions. Please feel free to contact me at: rresmini@gmu.edu or at 703-470-3022 (please
leave voice-mail).
Please also see the American Geophysical Union (AGU)
briefings located elsewhere on this site. They are relevant to the topic
of the JVGR paper.
l1.dat – lithology type
t1.dat –
initial temperature
h1.dat – heat
production (i.e., radiogenic sources); always zero
state1.dat
– initial state (solid, liquid)
pctmlt1.dat
– initial percent melt present
Main
output file name is:
coolXXX.T1
where XXX is node number in the sill/wallrock system.
Note, within the code, that other files may be generated but the procedure
calls are commented out.
Provided here is file: cool412.T1 for the 412th node in the model
corresponding to a position within the sill located 3 meters from the
sill/wallrock contact.
csdcalc2.for
(a FORTRAN program) inputs the coolXXX.T1 files produced by melt1d3.for and
calculates crystal populations. Output file names are sizeXXX.TMY where
XXX is node number in the sill/wallrock system and Y is an arbitrary number
indicating output file version. The file sizeXXX.TMY has to be trimmed so
that the very first temperature is the first subliquidus temperature.
Similarly, the file must also be trimmed so that the last temperature entry is the
first subsolidus temperature. csdcalc2.for also produces an output file
named grthXXX.TMY recording the growth rates calculated by the code.
Here, as above, XXX is node number in the sill/wallrock system and Y is an
arbitrary number indicating output file version.
Provided here is file: size412.TM1 also for the 412th node in the model
corresponding to a position within the sill located 3 meters from the
sill/wallrock contact.
Also provided is file: GRTH412.TM1
csdsynq.bas
(a BASIC program) inputs file sizeXXX.TMY produced by csdcalc2.for and
calculates a CSD. Output file names are chosen by user at run-time but
are usually: sizeXXX.csd for the total CSD results files and sizeXXX.reg for a
file of two columns containing crystal size and the natural log of crystal
population density, ln(n). Here, too, XXX is node number in the
sill/wallrock system.
Provided here are files: SIZE412.CSD, SIZE412.REG, and SIZE412.REG.xls
(an MS Excel version of SIZE412.REG).
For additional information on the model, please see: Note_to_Colleagues.doc.
Thank you.
Ron Resmini
College of Science
Department of Geography and Geoinformation Science
George Mason University
4400 University Drive
MSN 6A2
Fairfax, VA 22030-4444
v: 703-470-3022
e: rresmini@gmu.edu