FORTRAN 77 programs for conductive cooling of dikes with temperature-dependent thermal properties and heat of crystallization

Temperature histories obtained from transient heat-conduction theory are applicable to most dikes despite potential complicating effects related to magma flow during emplacement, groundwater circulation, and metamorphic reaction during cooling. Here. machine-independent FORTRAN 77 programs are presented to calculate temperatures in and around dikes as they cool conductively. Analytical solutions can treat thermal-property contrasts between the dike and host rocks, but cannot address the release of magmatic heat of crystallization after the early stages of cooling or the appreciable temperature dependence of thermal conductivity and diffusivity displayed by most rock types. Numerical solutions can incorporate these additional factors. The heat of crystallization can raise the initial temperature at the dike contact, ??c1, about 100??C above that which would be estimated if it were neglected, and can decrease the rate at which the front of solidified magma moves to the dike center by a factor of as much as three. Thermal conductivity and diffusivity of rocks increase with decreasing temperature and, at low temperatures, these properties increase more if the rocks are saturated with water. Models that treat these temperature dependencies yield estimates of ??c1 that are as much as 75??C beneath those which would be predicted if they were neglected. ?? 1988.