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.


melt1d3.for (a FORTRAN program) produces the temperature vs. time history for a position in a 40-meter thick sill.  Input files are:


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