This repository contains the codes for running ray tracing in 1D layered cake model, implemented in Fortran. There are multiple alternatives with better usage / more user-friendly interface to my codes; though they are also "heavier" to leverage MCMC seismic tomography.
The main purpose of Ray Tracing here is to be super quick for the particular case of 1D, so I could easily switch different velocity models, and still be able to quickly trace the rays with arrival times.
The repository contains both the scripts that invoke the "workhorse" functions, as well as the functions themselves.
- Compile the function codes with
R CMD SHLIB subroutineR-quiet.f90 - In R, load the shared object:
dyn.load("./subroutineR-quiet.so") - Call the function in this fashion:
timeP <- .Fortran("dff" ,vels=as.numeric(v), depths=as.numeric(d),NLayers=as.integer(NLayers), src_offset=as.numeric(srco), src_depth=as.numeric(srcd), NSrc=as.integer(NSrc), timeP = as.numeric(timeP))
I will populate later this section, after I "clean" this repository. For now, look at the details in rayTracerR.R.
- in the filebase.txt, first line is string length, second line is the string to be prepended to all filenames.
- mapfile : (transpose these rows into a long one row file)
k layer1-depth-to-interface layer1-vp layer2-depth-to-.. layer2-vp... .. .. vp -- for halfspace 0 0
iaccept into TEMPSWP for neighbours ( check the difference in beta equal to dTlog) running covariance matrix