Fortran code used in 'A fractal model for effective excess charge density in variably saturated fractured rocks'

This code is uploaded to support the research study 'A fractal model for effective excess charge density in variably saturated fractured rocks' by L. Guarracino and  D. Jougnot (submitted to JGR: Solid Earth, 2021).

Files:
a) Fortran source code (qvfrac.f) for estimating the effective excess charge density in fractured rocks. The calculation is based on model equations described in the research study.
b) Input data (network1.dat) to calculate the effective excess charge density for fracture network 1 described in Section 3 (Figure 5a).