ROck Physics Toolbox
1.0
A microgeodynamicsbased toolkit for rock physics.

This module contains data types and functions applicable to the region of interest. Each location in the region is assigned a unit cell, which contains the physical properties. The parameters for the equation of state of the melt phase are stored in the derived type MAGMA. It is assumed that only one kind of melt exists in the region. Three dimensional location of each point in the region are stored in the array LOC. Typically this data is stored as lat,lon, depth (km). The variable composition is also regional, it contains the basalt fraction of the bulk ocmposition and the average potential temperature of the region. Created by Saswata HierMajumder August, 2012.!<. More...
Data Types  
type  region 
Public Member Functions  
type(region) function  load_seismo (c, dim, fname, depth_in) 
type(region) function  melt_calculate (reg) 
Reads in input data from seismological models. The format of the input data should be the following: Lat Lon Depth Vs Vp If the data set is two dimensional, then an optional argument, depth_in (km) needs to be supplied. This feature is not yet fully implemented. !<.  
subroutine  vtk_write (reg, mod) 
This function returns the bulk modulus and pressure from the properties of melt set in function melt_set in module microgeodynamics. Either the Vinet or a third order Birch Murnaghan EOS can be used for the calculation. Just replace vinet with bm3 in the following line. Since both EOS are implicit equations, the density at depth must be known. If not known ahead of time, run a few trials and read in the pressure. Once the pressure is satisfactory for the region of interest, then use the prefactor. This process of iteration can also be automated either within the function or by calling an external function. Once the bulk modulus of the melt is known, it uses a nonlinear root finder to calculate the melt volume fraction based on the discrepancy between the predicted and observed Shear wave velocities. Currently, two solvers are implemented, a bisection algorithm and a combined bisection/ NewtonRaphson algorithm. The combined algorithm is called nonlinear_solver. The bisection algortihm is called bisection. They both take the same arguments. See more on these routines in the module microgeodynamics.!<. 
This module contains data types and functions applicable to the region of interest. Each location in the region is assigned a unit cell, which contains the physical properties. The parameters for the equation of state of the melt phase are stored in the derived type MAGMA. It is assumed that only one kind of melt exists in the region. Three dimensional location of each point in the region are stored in the array LOC. Typically this data is stored as lat,lon, depth (km). The variable composition is also regional, it contains the basalt fraction of the bulk ocmposition and the average potential temperature of the region. Created by Saswata HierMajumder August, 2012.!<.
Definition at line 13 of file regional.f90.
type(region) function regional::load_seismo  (  type(composition), intent(in)  c, 
integer, intent(in)  dim,  
character(len=*)  fname,  
real(sp), intent(in), optional  depth_in  
) 
Definition at line 29 of file regional.f90.
Reads in input data from seismological models. The format of the input data should be the following: Lat Lon Depth Vs Vp If the data set is two dimensional, then an optional argument, depth_in (km) needs to be supplied. This feature is not yet fully implemented. !<.
Here we load the observed variables into the region !< calculate potential temperature from depression of 410 Using function depth2dT in module microgeodynamics!< Convert the temperature to Shear wave velocity using the function temp2vs in module microgeodynamics!< Convert the temperature to P wave velocity using the function temp2vp in module microgeodynamics!< Convert the temperature to density using the function temp2rho in module microgeodynamics!< Calculate observed Shear wave velocity from impedance contrast observation!<
Definition at line 113 of file regional.f90.
subroutine regional::vtk_write  (  type(region), intent(in)  reg, 
type(mantle), intent(in)  mod  
) 
This function returns the bulk modulus and pressure from the properties of melt set in function melt_set in module microgeodynamics. Either the Vinet or a third order Birch Murnaghan EOS can be used for the calculation. Just replace vinet with bm3 in the following line. Since both EOS are implicit equations, the density at depth must be known. If not known ahead of time, run a few trials and read in the pressure. Once the pressure is satisfactory for the region of interest, then use the prefactor. This process of iteration can also be automated either within the function or by calling an external function. Once the bulk modulus of the melt is known, it uses a nonlinear root finder to calculate the melt volume fraction based on the discrepancy between the predicted and observed Shear wave velocities. Currently, two solvers are implemented, a bisection algorithm and a combined bisection/ NewtonRaphson algorithm. The combined algorithm is called nonlinear_solver. The bisection algortihm is called bisection. They both take the same arguments. See more on these routines in the module microgeodynamics.!<.
Definition at line 157 of file regional.f90.