GitHub - kanzakiy/HT-Oiso: 2D hydrothermal + oxygen isotope model

kanzakiy / HT-Oiso Public

Notifications You must be signed in to change notification settings
Fork 0
Star 0

2D hydrothermal + oxygen isotope model

Notifications

Name		Name	Last commit message	Last commit date
Latest commit History 9 Commits
HT1_v8uw.f90		HT1_v8uw.f90
HT1_v8uw_irr.f90		HT1_v8uw_irr.f90
HT1_v8uw_irr_off.f90		HT1_v8uw_irr_off.f90
HT1_v8uw_irr_off_btm.f90		HT1_v8uw_irr_off_btm.f90
LICENSE		LICENSE
MEMO.txt		MEMO.txt
Oiso+HT1_v6.f90		Oiso+HT1_v6.f90
Oiso+HT1_v6_irr.f90		Oiso+HT1_v6_irr.f90
README.txt		README.txt
steam.f90		steam.f90

Repository files navigation

Hydrothermal circulation + O isotope model

Oiso+HT1_v6_irr.f90 --> calculate steady state O isotopes of rock and porewater; including dispersion 
HT1_v8uw_irr.f90    --> Hydrothermal circulation without considering lateral free flow

Need UMFPACK & STEAM libraries
See https://github.com/PetterS/SuiteSparse/tree/master/UMFPACK for UMFPACK
See MEMO.txt for making STEAM library

First run hydrothermal model: 
    gfortran HT1_v8uw_irr_off.f90 umf4_f77wrapper.o -lumfpack -lamd -lcholmod -lcolamd -lsuitesparseconfig -lopenblas -lsteam -g -fcheck=all
    ./a
Then run isotope model:
    gfortran Oiso+HT1_v6_irr.f90 umf4_f77wrapper.o -lumfpack -lamd -lcholmod -lcolamd -lsuitesparseconfig -lopenblas -g -fcheck=all
    ./a