FEHM Version 3.6.2 - October 2024
This page describes version updates for FEHM V3.6 release with major updates from George Zyvoloski. Previous release FEHM V3.4 September 2019
This release was compiled and tested primarily on Windows using VisualStudio, Linux with GCC 9.4 and 13.2, and mac OSX compiled on GCC 13.2
Summary
This release includes improved phase change algorithms, saved zone usage, and boundary conditions primarily related to new isothermal two-phase flow with dissolved gas. Also updated zone_saved capability.
Code enhancements are generally related to improvements in numerical performance for the AMANZI isothermal, 2-phase, noncondensible, gas-water physics with gas solubility. Included are subroutine changes to accommodate restarts with the new AMANZI physics.
-
Added variables for global gas and dissolved gas, added cnlof() for last times step dissolved gas
-
Added tolerances for phase change, improved phase change (dissolved gas) algorithm
-
Changed VG upper/lower cutoff parameters
-
Added code to read and output isothermal multiphase solubility variables
-
Added new subroutine geneq2_sol to generate equations for iso 2phase with sol
-
Added use module com_prop_data (for future use)
-
Added new subroutine inpress_solubility to read variables associated with isothermal two phase w solubility; added old mass fraction cnlof() set to cnlf()
-
Modified code so initial liquid mass fraction is not read in ‘air’ keyword but in ‘pres’ keyword
-
Added coding to set phase state(ieos) based on initial conditions (including iso 2phase w sol)
-
Modified boundary condition (kq = -202), corrected but in solubility calculation
-
Modified code so upper cutoff is passed through subroutine header
-
Some code corrections for simulations related to VV Test cases and how output files are written.
-
Documentation has been changed from sphinx to github pages for easier updates. Missing pages and sections have been fixed. fehmpytests descriptions and instructions have been added.
FEHM Testing Improvements
fehmpytests is included in this repo and is recommended for developers as a fast, easy to use set of test cases to verify correctness of the code. This does not replace the VV Test Suites used since the YMP QA program in 2003, but is easier to use and faster to run. Most the test cases were constructed from the VV Test Suite and have the same input and output files.
The VV Test Suite is still used for checking analytical solutions and comprehensive multi-run models. It uses 2003 perl scripts and fortran programs to run and then compare results to analytical solutions. The perl scripts and programs may have problems running, but the test problems can still be run and results evaluated manually. The Windows version V3.4 of the VV Test Suite has been improved and is successful with this release. The Linux version 3.3 has not been updated but was uses for testing this release. Both VV Test Suites are available as assets with this release, see VERIFICATION_V3.4win_v3.6_clean.tar.gz and VERIFICATION_V3.3.0lnx_V3.6_clean.tar.gz
These 30 test cases were successful using fehmpytests: avdonin baro_vel bodyforce cden cflxz colloid_filtration dryout dispersion evaporation fracture_aperture head heat2d heat2d_quad heat3d heat_pipe henrys_law mptr multi_solute perm_test potential_energy ramey richards salt_perm_poro sorption sptr_btc theis toronyi vapor_extraction wvtest rad_decay
These 84 Windows VV Test cases were successful for this release: nts_thermal avdonin baro_trans baro_vel boun_test cellbased.init cellbased chain colloid_filtration colloid_ptrk colloid_mptr colloid_stream cden_test cflxz_test convection dissolution div_coarse div_fine doe dryout dual dual-k erosion_test evaporation forward reverse fracture_aperture fracture_transport gdpm head heat2d heat3d heat3d.finv 1D_convection heat2d_quad avdonin84 heat_gdkm_2D heat_pipe henrys_law infiltration infiltration_rlpm infiltration_gdkm multi_solute particle_capture perm_anisotropy perm_test potential_energy pruess_test ramey richards recycle_test_init recycle_test rh_test salt_test sorption binmode_flux binmode streamline sptr_btc tensor tensor_perm stress_2D stress_3D stress_3Dbeam stress_bodyforce stress_perm stress_restart theis toronyi transferf1 transferf2 transferf3 mptr_test transport3D.init transport3D vapor_extraction water_table lost_part wtrise_fehm wtrise_ini wtsi_test uz_test wvtest
For the Linux version V3.3 of the VV Test Suite, all tests were successful but perl scripts failed for some. These test had to be checked manually and compared to the Windows results: barometric, infiltration_gdkm, divergence, forward, particle_capture, uz_test
Code Bugs Fixed
Issue #79 Fixed trac option to write specified time steps for trac. The option was being ignored and now works. See fehmpytests case dispersion for use of this option.
Issue #51 Windows to Linux issues fixed. When compiling code written on Windows, these are some common problems. Fortran line length too long, replace eq with eqv for logicals, and replace windows call close(lu,dispose = ‘delete’) with close(lu). The Makefile.depends was updated to include all the new code and modules.
Known Issues
Issue 69. See test transport3D. Linux and windows write different formats for the *.trc files. Linux v3.6 and old fehm v3.4 write 6 values in a single line. The newer Windows v3.6 writes 6 values in double lines. fehmpytests allows either format and will check 6 values even if they are on double lines. This will be addressed in later versions of FEHM.
Modern compilers have less tolerance for old fortran codes and may have issues compiling. Some of these issues may cause FEHM to interrupt a run. The compiler issues will be addressed in the future with help from collaborators.
There are problems compiling FEHM with Intel Visual Fortran for debugging. Checking out the flag “Check Routine Interfaces” in the “Diagnostics” option of the Fortran workspace properties can let compilation of FEHM go through.
The code may produce the following, which do not affect correctness of the run: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
Code Changes
Software changes Feb 14, 2024
Changes compared to Master V3.4 release on repo.
From George
1. air_cp.f –
2. airctr_2j.f – Added and modified phase change when a soluble gas is present
3. allocmem_2b.f – corrected minor bug with heat conduction arrays
4. co2ctr_2w.f –added some coding to allow for mass fraction CO2 to
5. com_prop_data_2b.f – added character variable zones_char(:) for zones used by a node
6. comai_2h.f – added parameter strd_mass for new phase change with solubility
7. combi_2c.f –
8. comfi_2de.f – added global variables associated with new isothermal gas flow with solubility
9. connectivity_symmmetry_test_2a.f– minor bug corrected
10. convctr_2h.f – corrected minor bug to enable saved zones to properly work when inputting head in two-phase or WH or AWH fluid physics
11. datchk_2e.f-
12. dated_a5.f –changed variable FEHM verno to FEHM3.50_VER1ab_021324
13. diskwrite_new_2a.f –
14. fehmn_pcx_4_2h.f – added/modified comment
15. flowrate_vectors_2c.f – disabled for now, will fix soon
16. geneq2_2m.f – corrected errors with gas diffusion
17. gensl2_2e.f – added minor coding for isothermal gas diffusion
18. indiff_2a.f– new subroutine to read in diffusivities for isothermal gas diffusion
19. infiles_2b.f – changed coding for saved zones (added ‘infi’ to saved zones call)
20. inflo3_2e.f –test for dissolved species
21. inhist_2b.f –small correction for printing rlp data
22. initdata2_2g.f – added subroutine zone_char_fill that fills character variable zone_char, added subroutine zone_saved that manages all the saved zones.
23. innode_2c.f –
24. input_2g.f – added coding to allocate diffusivities for isothermal gas diffusion
25. namefile2.f– minor modification, changed output format
26. nr_stop_ctr_2f.f – fixed bug when counting phase changes
27. plot_new_2c.f – added coding to output air pressure as well as water pressure
28. rlp_cap_2a.f90 – no changes
29. scanin_b6.f – added modifications to support saved zones
30. thermw_prop_2g.f – more corrected bugs for a heat conduction(only) problem
31. thrair_3m.f- fixed some bugs for two-phase flow with dissolved gas
32. vgrlp_2a.f- corrected some errors in gas relative permeabilities
33. vgrlps_2e.f- modified code to have variable cutoff tolerances
34. wrtout_2d.f – corrected bugs for a heat conduction(only) problem
Note: write_avs_node_mat_s.f replaces write_avs_node_mat.f
subroutines replaces write_avs_node_mat are commented out
new call is now write_avs_node_mat_s which is very different
old subroutine write_avs_node_mat(lu,ifdual,mout)
new subroutine write_avs_node_mat_s(icall,
2 neq,
3 nscalar,
4 lu,
5 ifdual,
6 iriver2, mout)
These VV Window test cases fail:
************ BOUN_TEST *************
************ EVAPORATION ************* (gaz will fix)
************ RICHARDS ************* (gaz will fix)
************ REVERSE_TRACKING *************
************ COLLOID_STREAM *************
************ TENSOR_PERM *************
************ WATER_TABLE ************* (gaz will fix)
Set of 17 files with fixes Sent Apr 30:
SRC_TO_LANL_043024
-rw-rw-r-- 1 tamiller sft 93280 Mar 8 20:42 scanin_b7.f
-rw-rw-r-- 1 tamiller sft 58289 Mar 9 00:05 input_2g.f
-rw-rw-r-- 1 tamiller sft 19050 Mar 18 21:51 diskwrite_2a.f
-rw-rw-r-- 1 tamiller sft 37459 Apr 7 14:21 gensl2_2f.f
-rw-rw-r-- 1 tamiller sft 7304 Apr 8 07:31 dated_a6.f
-rw-rw-r-- 1 tamiller sft 11942 Apr 12 11:57 vgrlps_2f.f
-rw-rw-r-- 1 tamiller sft 60655 Apr 15 10:14 rlperm_2h.f
-rw-rw-r-- 1 tamiller sft 2154 Apr 15 10:49 comrlp_param.f
-rw-rw-r-- 1 tamiller sft 72070 Apr 18 20:20 diskread_2d.f
-rw-rw-r-- 1 tamiller sft 25268 Apr 19 07:44 diskwrite_new_2b.f
-rw-rw-r-- 1 tamiller sft 13382 Apr 20 15:16 comfi_2h.f
-rw-rw-r-- 1 tamiller sft 49053 Apr 21 12:23 geneq2_sol_3a.f
-rw-rw-r-- 1 tamiller sft 9604 Apr 25 12:16 air_cp_2e.f
-rw-rw-r-- 1 tamiller sft 54191 Apr 25 16:54 thrair_3q.f
-rw-rw-r-- 1 tamiller sft 58322 Apr 25 16:54 startup_2i.f
-rw-rw-r-- 1 tamiller sft 18209 Apr 25 16:54 inpres_2d.f
-rw-rw-r-- 1 tamiller sft 49203 Apr 29 17:09 airctr_3o.f
These files are new:
comrlp_param.f
geneq2_sol_3a.f
SENT TO LANL --------------------------------------------------------
May 1 Set of current build files that pass VV Tests on Windows
/project/eesdev/FEHM/from_george/SRC_W_BUILD_050124
FEHM3.6_VER2.vfproj SRC_MASTER SRC_MASTER_UPDATE_1 SRC_MASTER_UPDATE_2
AT LANL May 9 -------------------------------------------------------
Added missing file geneq2.f
/project/eesdev/FEHM/from_george/src_050924
src_050924.zip (2.48 MB)
VERIFICATION_V3.4win_64bit_clean_CP.zip (398 MB)
All tests pass except BAROMETRIC
Sep 2024 The correct 'mainrip.f' file for Windows was in the ./src/PC folder, GFORTRAN has the UNIX versions of the code.
These were fixed so Windows uses all src files with correct mainrip.f file, linux uses files in GFORTRAN.
FEHM now passes all VV problems, including the Barometric.
Oct 2024 Bug fix for trac files using NPRTTRC to write specified times.
Group 3 - IACCMX, DAYCM, DAYCMM, DAYCMX, NPRTTRC
Fixed with new files: csolve.f tracout.f write_avs_node_con.f
This is now tested in fehmpytests/dispersion
The following is a copy of George’s compile file with comments regarding code changes.
Note: File updates are embedded in the relative source locations and names may differ from repo file names.
<File RelativePath=".\SRC_MASTER\accum_derivative _2d.f"/>
<File RelativePath=".\SRC_MASTER\accum_term_2d.f"/>
<File RelativePath=".\SRC_MASTER\active_nodes_ctr.f"/>
<File RelativePath=".\SRC_MASTER\add_gdpm_2a.f"/>
<File RelativePath=".\SRC_MASTER\air_combine_2d.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\air_cp_2e.f"/>
1. Added variables for global gas and dissolved gas, added cnlof()for last times step dissolved gas
<File RelativePath=".\SRC_MASTER\air_properties_new.f90"/>
<File RelativePath=".\SRC_MASTER\air_rdof.f"/>
<File RelativePath=".\SRC_MASTER\air_rdof_part.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\airctr_3o.f"/>
2. Added tolerances for phase change, improved phase change (dissolved gas) algorithm
<File RelativePath=".\SRC_MASTER\airctr_part_2d.f"/>
<File RelativePath=".\SRC_MASTER\allocmem_2e.f"/>
<File RelativePath=".\SRC_MASTER\angle3.f"/>
<File RelativePath=".\SRC_MASTER\anonp.f"/>
<File RelativePath=".\SRC_MASTER\anonp_ani.f"/>
<File RelativePath=".\SRC_MASTER\area2d_tri.f"/>
<File RelativePath=".\SRC_MASTER\area_flow_bc.f"/>
<File RelativePath=".\SRC_MASTER\area_length_calc.f"/>
<File RelativePath=".\SRC_MASTER\area_vol_tet.f"/>
<File RelativePath=".\SRC_MASTER\area_vol_weightctr.f"/>
<File RelativePath=".\SRC_MASTER\avs_io_2d.f"/>
<File RelativePath=".\SRC_MASTER\avs_write_cord.f"/>
<File RelativePath=".\SRC_MASTER\avs_write_struc.f"/>
<File RelativePath=".\SRC_MASTER\avsio_2e.f"/>
<File RelativePath=".\SRC_MASTER\bc_far_ctr.f"/>
<File RelativePath=".\SRC_MASTER\bcon_2a.f"/>
<File RelativePath=".\SRC_MASTER\bcrlp.f"/>
<File RelativePath=".\SRC_MASTER\bit.f"/>
<File RelativePath=".\SRC_MASTER\bnswer_2d.f"/>
<File RelativePath=".\SRC_MASTER\bnswer_part_2d.f"/>
<File RelativePath=".\SRC_MASTER\brooks_corey_cap.f"/>
<File RelativePath=".\SRC_MASTER\brooks_corey_cap3.f"/>
<File RelativePath=".\SRC_MASTER\brooks_corey_rlp.f"/>
<File RelativePath=".\SRC_MASTER\brooks_corey_rlp3.f"/>
<File RelativePath=".\SRC_MASTER\CalcPoros.f90"/>
<File RelativePath=".\SRC_MASTER\calculate_cdf.f"/>
<File RelativePath=".\SRC_MASTER\calculate_ret.f"/>
<File RelativePath=".\SRC_MASTER\cappr.f"/>
<File RelativePath=".\SRC_MASTER\cappr2.f90"/>
<File RelativePath=".\SRC_MASTER\cappr_hyd.f"/>
<File RelativePath=".\SRC_MASTER\cascade_sat.f"/>
<File RelativePath=".\SRC_MASTER\cden_cor.f90"/>
<File RelativePath=".\SRC_MASTER\cell_time.f"/>
<File RelativePath=".\SRC_MASTER\cfluxz.f"/>
<File RelativePath=".\SRC_MASTER\check_rlp_2c.f"/>
<File RelativePath=".\SRC_MASTER\check_rlp_carb.f"/>
<File RelativePath=".\SRC_MASTER\check_sx.f"/>
<File RelativePath=".\SRC_MASTER\chemod.f"/>
<File RelativePath=".\SRC_MASTER\cliffnodes.f"/>
<File RelativePath=".\SRC_MASTER\close_files.f"/>
<File RelativePath=".\SRC_MASTER\cnswer.f"/>
<File RelativePath=".\SRC_MASTER\cntlin.f"/>
<File RelativePath=".\SRC_MASTER\cntlio.f"/>
<File RelativePath=".\SRC_MASTER\co2_array.f"/>
<File RelativePath=".\SRC_MASTER\co2_fugacity.f"/>
<File RelativePath=".\SRC_MASTER\co2_properties.f90"/>
<File RelativePath=".\SRC_MASTER\co2_properties_new.f90"/>
<File RelativePath=".\SRC_MASTER\co2ctr_2w.f"/>
<File RelativePath=".\SRC_MASTER\co2h2o_combine.f"/>
<File RelativePath=".\SRC_MASTER\co2wh_properties_new.f90"/>
<File RelativePath=".\SRC_MASTER\coef_replace_ctr.f"/>
<File RelativePath=".\SRC_MASTER\coeff_management.f"/>
<File RelativePath=".\SRC_MASTER\coeffc.f"/>
<File RelativePath=".\SRC_MASTER\com_exphase.f"/>
<File RelativePath=".\SRC_MASTER\com_part.f"/>
<File RelativePath=".\SRC_MASTER\com_prop_data_2b.f"/>
<File RelativePath=".\SRC_MASTER\comai_2h.f"/>
<File RelativePath=".\SRC_MASTER\combi_2c.f"/>
<File RelativePath=".\SRC_MASTER\combine_a.f"/>
<File RelativePath=".\SRC_MASTER\comchem.f"/>
<File RelativePath=".\SRC_MASTER\comci.f"/>
<File RelativePath=".\SRC_MASTER\comco2.f"/>
<File RelativePath=".\SRC_MASTER\comcomp.f"/>
<File RelativePath=".\SRC_MASTER\comcouple.f"/>
<File RelativePath=".\SRC_MASTER\comdi_4b.f"/>
<File RelativePath=".\SRC_MASTER\comdti.f"/>
<File RelativePath=".\SRC_MASTER\comei.f"/>
<File RelativePath=".\SRC_MASTER\comevap.f"/>
<File RelativePath=".\SRC_MASTER\comfem.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\comfi_2h.f"/>
3. Changed Comments only
<File RelativePath=".\SRC_MASTER\comflow.f"/>
<File RelativePath=".\SRC_MASTER\comgi.f"/>
<File RelativePath=".\SRC_MASTER\comhi.f"/>
<File RelativePath=".\SRC_MASTER\comii_2.f"/>
<File RelativePath=".\SRC_MASTER\comji.f"/>
<File RelativePath=".\SRC_MASTER\comki.f"/>
<File RelativePath=".\SRC_MASTER\commass_AWH.f"/>
<File RelativePath=".\SRC_MASTER\commeth.f"/>
<File RelativePath=".\SRC_MASTER\comnire.f"/>
<File RelativePath=".\SRC_MASTER\compart.f"/>
<File RelativePath=".\SRC_MASTER\compfrac.f"/>
<File RelativePath=".\SRC_MASTER\compute_average_stress.f"/>
<File RelativePath=".\SRC_MASTER\compute_exit_new.f"/>
<File RelativePath=".\SRC_MASTER\compute_flow_residual.f"/>
<File RelativePath=".\SRC_MASTER\compute_ipsi.f"/>
<File RelativePath=".\SRC_MASTER\compute_permfactor.f"/>
<File RelativePath=".\SRC_MASTER\compute_permfactor_effstrs.f"/>
<File RelativePath=".\SRC_MASTER\compute_permfactor_effstrs_pp.f"/>
<File RelativePath=".\SRC_MASTER\compute_permfactor_vonMises.f"/>
<File RelativePath=".\SRC_MASTER\compute_strains.f"/>
<File RelativePath=".\SRC_MASTER\comriv.f"/>
<File RelativePath=".\SRC_MASTER\comrlp_2a.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\comrlp_param.f"/>
4. Changed VG upper/lower cutoff parameters
5.
<File RelativePath=".\SRC_MASTER\comrtd.f"/>
<File RelativePath=".\SRC_MASTER\comrxnb.f"/>
<File RelativePath=".\SRC_MASTER\comrxni.f"/>
<File RelativePath=".\SRC_MASTER\comsi.f"/>
<File RelativePath=".\SRC_MASTER\comsk.f"/>
<File RelativePath=".\SRC_MASTER\comsplitts.f"/>
<File RelativePath=".\SRC_MASTER\comsptr.f"/>
<File RelativePath=".\SRC_MASTER\comsteady_2.f"/>
<File RelativePath=".\SRC_MASTER\comtable.f"/>
<File RelativePath=".\SRC_MASTER\comuserc.f"/>
<File RelativePath=".\SRC_MASTER\comwellphys.f"/>
<File RelativePath=".\SRC_MASTER\comwt.f"/>
<File RelativePath=".\SRC_MASTER\comxi.f"/>
<File RelativePath=".\SRC_MASTER\comzeoli.f"/>
<File RelativePath=".\SRC_MASTER\comzone.f"/>
<File RelativePath=".\SRC_MASTER\concadiff.f"/>
<File RelativePath=".\SRC_MASTER\concen.f"/>
<File RelativePath=".\SRC_MASTER\coneq1.f"/>
<File RelativePath=".\SRC_MASTER\coneq1_ani.f"/>
<File RelativePath=".\SRC_MASTER\coneq1_gdpm.f"/>
<File RelativePath=".\SRC_MASTER\coneq1mdnode.f"/>
<File RelativePath=".\SRC_MASTER\connections_list.f"/>
<File RelativePath=".\SRC_MASTER\connectivity_symmmetry_test_2a.f"/>
<File RelativePath=".\SRC_MASTER\constant_value.f"/>
<File RelativePath=".\SRC_MASTER\contr_2b.f"/>
<File RelativePath=".\SRC_MASTER\contrc.f"/>
<File RelativePath=".\SRC_MASTER\contrj.f"/>
<File RelativePath=".\SRC_MASTER\convctr_2h.f"/>
<File RelativePath=".\SRC_MASTER\corey.f"/>
<File RelativePath=".\SRC_MASTER\cplxcalc.f"/>
<File RelativePath=".\SRC_MASTER\crdpdp.f"/>
<File RelativePath=".\SRC_MASTER\csolve.f"/>
<File RelativePath=".\SRC_MASTER\ctdpdp.f"/>
<File RelativePath=".\SRC_MASTER\darcy_derivatives_BLI.f"/>
<File RelativePath=".\SRC_MASTER\data_4c.f"/>
<File RelativePath=".\SRC_MASTER\datchk_2e.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\dated_a6.f"/>
6. Changed the FEHM banner to FEHM3.6VER2
<File RelativePath=".\SRC_MASTER\davidi.f"/>
<File RelativePath=".\SRC_MASTER\daycrl.f"/>
<File RelativePath=".\SRC_MASTER\degree.f"/>
<File RelativePath=".\SRC_MASTER\DeltaHeight.f90"/>
<File RelativePath=".\SRC_MASTER\DeltaPoros.f90"/>
<File RelativePath=".\SRC_MASTER\den_vis_spatial.f"/>
<File RelativePath=".\SRC_MASTER\determ.f"/>
<File RelativePath=".\SRC_MASTER\diagnostics_2b.f"/>
<File RelativePath=".\SRC_MASTER\directional_search.f"/>
<File RelativePath=".\SRC_MASTER\diskc.f"/>
<File RelativePath=".\SRC_MASTER\diskp.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\diskread_2d.f"/>
7. Added coding to output isothermal multiphase solubility variables
<File RelativePath=".\SRC_MASTER_UPDATE_1\diskwrite_2a.f"/>
8. No changes
<File RelativePath=".\SRC_MASTER_UPDATE_1\diskwrite_new_2b.f"/>
9. Added coding to read isothermal multiphase solubility variables
<File RelativePath=".\SRC_MASTER\dispersion_divergence.f"/>
<File RelativePath=".\SRC_MASTER\dispersion_divergence_omr.f"/>
<File RelativePath=".\SRC_MASTER\dispersion_gradient.f"/>
<File RelativePath=".\SRC_MASTER\dispersion_node.f"/>
<File RelativePath=".\SRC_MASTER\done_macro.f"/>
<File RelativePath=".\SRC_MASTER\DoubleDot.f"/>
<File RelativePath=".\SRC_MASTER\dpdp.f"/>
<File RelativePath=".\SRC_MASTER\dpdp3.f"/>
<File RelativePath=".\SRC_MASTER\dpdpfa.f"/>
<File RelativePath=".\SRC_MASTER\dpdpfh.f"/>
<File RelativePath=".\SRC_MASTER\dpdpta.f"/>
<File RelativePath=".\SRC_MASTER\dual.f"/>
<File RelativePath=".\SRC_MASTER\dualex.f"/>
<File RelativePath=".\SRC_MASTER\dualfa.f"/>
<File RelativePath=".\SRC_MASTER\dualfh.f"/>
<File RelativePath=".\SRC_MASTER\dualta.f"/>
<File RelativePath=".\SRC_MASTER\dualtx.f"/>
<File RelativePath=".\SRC_MASTER\dvacalc_a3.f"/>
<File RelativePath=".\SRC_MASTER\eflxz.f"/>
<File RelativePath=".\SRC_MASTER\eigenvectors.f90"/>
<File RelativePath=".\SRC_MASTER\elastic_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\elastic_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\elem_type_2a.f"/>
<File RelativePath=".\SRC_MASTER\end_macro.f"/>
<File RelativePath=".\SRC_MASTER\enrich_ctr.f"/>
<File RelativePath=".\SRC_MASTER\enthp.f"/>
<File RelativePath=".\SRC_MASTER\equal_array.f"/>
<File RelativePath=".\SRC_MASTER\eullag3.f"/>
<File RelativePath=".\SRC_MASTER\evaporation.f"/>
<File RelativePath=".\SRC_MASTER\explicit.f"/>
<File RelativePath=".\SRC_MASTER\explicit_phase_update_a1.f"/>
<File RelativePath=".\SRC_MASTER\exponential.f"/>
<File RelativePath=".\SRC_MASTER\fd_calc_heat.f"/>
<File RelativePath=".\SRC_MASTER\fehmn_pcx_4_2h.f"/>
<File RelativePath=".\SRC_MASTER\fem_DruckerPrager_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\fem_DruckerPrager_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\fem_elastic_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\fem_elastic_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\fem_material_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\fem_material_stiffness_2D.f"/>
<File RelativePath=".\SRC_MASTER\fem_material_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\fem_shapefunctions_2r.f"/>
<File RelativePath=".\SRC_MASTER\fem_shapefunctions_3r.f"/>
<File RelativePath=".\SRC_MASTER\fem_transverse_isotropy_elastic_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\fem_update_stress.f"/>
<File RelativePath=".\SRC_MASTER\fem_update_stress_2D.f"/>
<File RelativePath=".\SRC_MASTER\fem_vonMises_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\fem_vonMises_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\file_prefix.f"/>
<File RelativePath=".\SRC_MASTER\fimpf.f"/>
<File RelativePath=".\SRC_MASTER\find_closest_curve.f"/>
<File RelativePath=".\SRC_MASTER\find_v_3lsq.f"/>
<File RelativePath=".\SRC_MASTER\find_v_lsq.f"/>
<File RelativePath=".\SRC_MASTER\flow_boun_b0.f"/>
<File RelativePath=".\SRC_MASTER\flow_boundary_conditions_b0.f"/>
<File RelativePath=".\SRC_MASTER\flow_flux_perm_residual.f"/>
<File RelativePath=".\SRC_MASTER\flow_humidity_bc.f"/>
<File RelativePath=".\SRC_MASTER\flowrate_vectors_2c.f"/>
<File RelativePath=".\SRC_MASTER\fluid_prop_control_7h.f"/>
<File RelativePath=".\SRC_MASTER\flux_net.f"/>
<File RelativePath=".\SRC_MASTER\flxn.f"/>
<File RelativePath=".\SRC_MASTER\flxo.f"/>
<File RelativePath=".\SRC_MASTER\flxz_2.f"/>
<File RelativePath=".\SRC_MASTER\fnroot.f"/>
<File RelativePath=".\SRC_MASTER\ftime.f"/>
<File RelativePath=".\SRC_MASTER\gdkm_calc.f"/>
<File RelativePath=".\SRC_MASTER\gdkm_connect.f"/>
<File RelativePath=".\SRC_MASTER\gdkm_volume_fraction_interface_2c.f"/>
<File RelativePath=".\SRC_MASTER\gdpm_corr.f"/>
<File RelativePath=".\SRC_MASTER\gdpm_geneqh.f"/>
<File RelativePath=".\SRC_MASTER\gen_mixmodel.f"/>
<File RelativePath=".\SRC_MASTER\gencof.f"/>
<File RelativePath=".\SRC_MASTER\gencon_2a.f"/>
<File RelativePath=".\SRC_MASTER\gendat.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_2a.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_ani.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_stress.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_stress_coupl.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_stress_coupl_porosi.f"/>
<File RelativePath=".\SRC_MASTER\geneq1_well.f"/>
<File RelativePath=".\FEHM_SRC_BEFORE_JULY2020\geneq2.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_ani.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_ani_2p.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_part.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_rich_2d.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\geneq2_sol_3a.f"/>
10. Created new subroutine geneq2_sol .f to generate equations for iso 2phase with sol
<File RelativePath=".\SRC_MASTER\geneq2_switch.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_uz_wt.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_wellphysics.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_wtsi_2d.f"/>
<File RelativePath=".\SRC_MASTER\geneq2_wtsi_well_2d.f"/>
<File RelativePath=".\SRC_MASTER\geneq3.f"/>
<File RelativePath=".\SRC_MASTER\geneq_air_co2.f"/>
<File RelativePath=".\SRC_MASTER\geneq_co2.f"/>
<File RelativePath=".\SRC_MASTER\geneq_co2_wellphysics.f"/>
<File RelativePath=".\SRC_MASTER\geneq_flow_coupled.f"/>
<File RelativePath=".\SRC_MASTER\geneq_flowflux_perm.f"/>
<File RelativePath=".\SRC_MASTER\geneq_h2o.f"/>
<File RelativePath=".\SRC_MASTER\geneq_h2o_co2.f"/>
<File RelativePath=".\SRC_MASTER\geneq_hyd.f"/>
<File RelativePath=".\SRC_MASTER\geneq_hyd_equil.f"/>
<File RelativePath=".\SRC_MASTER\geneq_meth.f"/>
<File RelativePath=".\SRC_MASTER\geneq_plastic_3D.f"/>
<File RelativePath=".\SRC_MASTER\geneq_stress_fem_2D_XY.f"/>
<File RelativePath=".\SRC_MASTER\geneq_stress_fem_3D.f"/>
<File RelativePath=".\SRC_MASTER\geneq_stress_uncoupled_2D.f"/>
<File RelativePath=".\SRC_MASTER\geneq_stress_uncoupled_2D_FEM_rad.f"/>
<File RelativePath=".\SRC_MASTER\geneq_stress_uncoupled_3D.f"/>
<File RelativePath=".\SRC_MASTER\geneqc.f"/>
<File RelativePath=".\SRC_MASTER\geneqmdnode.f"/>
<File RelativePath=".\SRC_MASTER\generate_flow.f"/>
<File RelativePath=".\SRC_MASTER\genrcm.f"/>
<File RelativePath=".\SRC_MASTER\gensco2h2o.f"/>
<File RelativePath=".\SRC_MASTER\gensdp.f"/>
<File RelativePath=".\SRC_MASTER\gensdp3.f"/>
<File RelativePath=".\SRC_MASTER\gensdp_switch.f"/>
<File RelativePath=".\SRC_MASTER\gensl1.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\gensl2_2f.f"/>
11. Added use module com_prop_data not used yet
<File RelativePath=".\SRC_MASTER\gensl2_part.f"/>
<File RelativePath=".\SRC_MASTER\gensl2_switch_2d.f"/>
<File RelativePath=".\SRC_MASTER\gensl3_2a.f"/>
<File RelativePath=".\SRC_MASTER\gensl4_2d.f"/>
<File RelativePath=".\SRC_MASTER\gensl_stress_2D.f"/>
<File RelativePath=".\SRC_MASTER\gensl_stress_3D.f"/>
<File RelativePath=".\SRC_MASTER\gensl_stress_coupled_3D.f"/>
<File RelativePath=".\SRC_MASTER\gensmethh2o.f"/>
<File RelativePath=".\SRC_MASTER\gensmethh2o_equil.f"/>
<File RelativePath=".\SRC_MASTER\gentdp.f"/>
<File RelativePath=".\SRC_MASTER\geoin.f"/>
<File RelativePath=".\SRC_MASTER\getarea.f"/>
<File RelativePath=".\SRC_MASTER\getconc14.f"/>
<File RelativePath=".\SRC_MASTER\getconc36.f"/>
<File RelativePath=".\SRC_MASTER\GetSubsTime.f90"/>
<File RelativePath=".\SRC_MASTER\GetSubsWaveLength.f90"/>
<File RelativePath=".\SRC_MASTER\gncf2.f"/>
<File RelativePath=".\SRC_MASTER\gncf3.f"/>
<File RelativePath=".\SRC_MASTER\gp_global_coord.f"/>
<File RelativePath=".\SRC_MASTER\grad_array.f"/>
<File RelativePath=".\SRC_MASTER\gradctr_1.f"/>
<File RelativePath=".\SRC_MASTER\h2o_properties.f"/>
<File RelativePath=".\SRC_MASTER\h2o_properties_new_1.f90"/>
<File RelativePath=".\SRC_MASTER\harmonic_mean.f"/>
<File RelativePath=".\SRC_MASTER\head_2phase_2d.f"/>
<File RelativePath=".\SRC_MASTER\headctr_2d.f"/>
<File RelativePath=".\SRC_MASTER\heatloc.f"/>
<File RelativePath=".\SRC_MASTER\hstz.f"/>
<File RelativePath=".\SRC_MASTER\humidity.f"/>
<File RelativePath=".\SRC_MASTER\hyddiss.f"/>
<File RelativePath=".\SRC_MASTER\hydrate_array.f"/>
<File RelativePath=".\SRC_MASTER\hydrate_equil.f"/>
<File RelativePath=".\SRC_MASTER\hydrate_properties.f"/>
<File RelativePath=".\SRC_MASTER\icectr.f"/>
<File RelativePath=".\SRC_MASTER\icectrco2.f"/>
<File RelativePath=".\SRC_MASTER\implicit_well2.f"/>
<File RelativePath=".\SRC_MASTER\impsample.f"/>
<File RelativePath=".\SRC_MASTER\impsample_ptrk.f"/>
<File RelativePath=".\SRC_MASTER\inanpe.f"/>
<File RelativePath=".\SRC_MASTER\incden.f90"/>
<File RelativePath=".\SRC_MASTER\incond.f"/>
<File RelativePath=".\SRC_MASTER\incoord.f"/>
<File RelativePath=".\SRC_MASTER\inctrl.f"/>
<File RelativePath=".\SRC_MASTER\indices.f"/>
<File RelativePath=".\SRC_MASTER\indiff_2a.f"/>
<File RelativePath=".\SRC_MASTER\indpdp.f"/>
<File RelativePath=".\SRC_MASTER\infiles_2b.f"/>
<File RelativePath=".\SRC_MASTER\inflo2.f"/>
<File RelativePath=".\SRC_MASTER\inflo3_2e.f"/>
<File RelativePath=".\SRC_MASTER\infloa.f"/>
<File RelativePath=".\SRC_MASTER\inflogh.f"/>
<File RelativePath=".\SRC_MASTER\inflow_2d.f"/>
<File RelativePath=".\SRC_MASTER\inflow_wt.f"/>
<File RelativePath=".\SRC_MASTER\infrlp.f"/>
<File RelativePath=".\SRC_MASTER\ingdpm_2b.f"/>
<File RelativePath=".\SRC_MASTER\ingrowth.f"/>
<File RelativePath=".\SRC_MASTER\inhflx.f"/>
<File RelativePath=".\SRC_MASTER\inhist_2b.f"/>
<File RelativePath=".\SRC_MASTER\inhstz.f"/>
<File RelativePath=".\SRC_MASTER\inhyco.f"/>
<File RelativePath=".\SRC_MASTER\initchem.f"/>
<File RelativePath=".\SRC_MASTER\initdata2_2g.f"/>
<File RelativePath=".\SRC_MASTER\initfc.f"/>
<File RelativePath=".\SRC_MASTER\inmentat.f"/>
<File RelativePath=".\SRC_MASTER\inmptr.f"/>
<File RelativePath=".\SRC_MASTER\innode_2c.f"/>
<File RelativePath=".\SRC_MASTER\inpatran.f"/>
<File RelativePath=".\SRC_MASTER\inperm.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\inpres_2d.f"/>
12. Added new subroutine inpress_solubility to read variables associated with isothermal two phase w solubility; added old mass fraction cnlof() (set to cnlf()
<File RelativePath=".\SRC_MASTER\inptrk.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\input_2g.f"/>
13. Nothing new
<File RelativePath=".\SRC_MASTER\inrestart.f"/>
<File RelativePath=".\SRC_MASTER\inrlp_2b.f90"/>
<File RelativePath=".\SRC_MASTER\inrock.f"/>
<File RelativePath=".\SRC_MASTER\insptr.f"/>
<File RelativePath=".\SRC_MASTER\interblock_iso.f"/>
<File RelativePath=".\SRC_MASTER\interp.f"/>
<File RelativePath=".\SRC_MASTER\interp2.f"/>
<File RelativePath=".\SRC_MASTER\interp4.f"/>
<File RelativePath=".\SRC_MASTER\interp8.f"/>
<File RelativePath=".\SRC_MASTER\interp_bilin.f"/>
<File RelativePath=".\SRC_MASTER\interpolate.f90"/>
<File RelativePath=".\SRC_MASTER\interpolate_2a.f90"/>
<File RelativePath=".\SRC_MASTER\interpolate_2b.f90"/>
<File RelativePath=".\SRC_MASTER\interpolate_2c.f90"/>
<File RelativePath=".\SRC_MASTER\interpolate_q2n_3r.f"/>
<File RelativePath=".\SRC_MASTER\intime.f"/>
<File RelativePath=".\SRC_MASTER\inverf.f"/>
<File RelativePath=".\SRC_MASTER\inzeol.f"/>
<File RelativePath=".\SRC_MASTER\iofile.f"/>
<File RelativePath=".\SRC_MASTER\J2Invariant.f"/>
<File RelativePath=".\SRC_MASTER\KozenyCoeff.f90"/>
<File RelativePath=".\SRC_MASTER\linear.f"/>
<File RelativePath=".\SRC_MASTER\load_omr_flux_array.f"/>
<File RelativePath=".\SRC_MASTER\loadpcnsk.f"/>
<File RelativePath=".\SRC_MASTER\lstsq.f"/>
<File RelativePath=".\SRC_MASTER\lubksb.f"/>
<File RelativePath=".\SRC_MASTER\lubksb0.f"/>
<File RelativePath=".\SRC_MASTER\lubsk.f"/>
<File RelativePath=".\SRC_MASTER\ludcmp.f"/>
<File RelativePath=".\SRC_MASTER\ludcmp0.f"/>
<File RelativePath=".\SRC_MASTER\ludsk.f"/>
<File RelativePath=".\SRC_MASTER\mainrip.f"/>
<File RelativePath=".\SRC_MASTER\massflow.f"/>
<File RelativePath=".\SRC_MASTER\material_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\material_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\mdnodes.f"/>
<File RelativePath=".\SRC_MASTER\methane_properties.f"/>
<File RelativePath=".\SRC_MASTER\methh2o_combine.f"/>
<File RelativePath=".\SRC_MASTER\min_max.f"/>
<File RelativePath=".\SRC_MASTER\mod_eqs_ngas.f"/>
<File RelativePath=".\SRC_MASTER\model_setup_2d.f"/>
<File RelativePath=".\SRC_MASTER\modpod.f"/>
<File RelativePath=".\SRC_MASTER\move_particle.f"/>
<File RelativePath=".\SRC_MASTER\namefile1.f"/>
<File RelativePath=".\SRC_MASTER\namefile2.f"/>
<File RelativePath=".\SRC_MASTER\near3.f"/>
<File RelativePath=".\SRC_MASTER\new_loc_omr.f"/>
<File RelativePath=".\SRC_MASTER\node_midedge_box_size.f"/>
<File RelativePath=".\SRC_MASTER\nopcnr.f"/>
<File RelativePath=".\SRC_MASTER\nopcnv.f"/>
<File RelativePath=".\SRC_MASTER\normal.f"/>
<File RelativePath=".\SRC_MASTER\normal_dof_2b.f"/>
<File RelativePath=".\SRC_MASTER\normal_uc.f"/>
<File RelativePath=".\SRC_MASTER\nr_stop_ctr_2f.f"/>
<File RelativePath=".\SRC_MASTER\nrmlz4.f"/>
<File RelativePath=".\SRC_MASTER\null.f"/>
<File RelativePath=".\SRC_MASTER\null_new.f"/>
<File RelativePath=".\SRC_MASTER\open_file.f"/>
<File RelativePath=".\SRC_MASTER\outbnd_2.f"/>
<File RelativePath=".\SRC_MASTER\paractr.f"/>
<File RelativePath=".\SRC_MASTER\parse_string.f"/>
<File RelativePath=".\SRC_MASTER\parse_string2.f"/>
<File RelativePath=".\SRC_MASTER\part_track.f"/>
<File RelativePath=".\SRC_MASTER\pcp_save.f"/>
<File RelativePath=".\SRC_MASTER\pebi.f"/>
<File RelativePath=".\SRC_MASTER\pebi3.f"/>
<File RelativePath=".\SRC_MASTER\peint.f"/>
<File RelativePath=".\SRC_MASTER\permfac_nodal_out.f"/>
<File RelativePath=".\SRC_MASTER\pest_2.f"/>
<File RelativePath=".\SRC_MASTER\phase_balance_ctr.f"/>
<File RelativePath=".\SRC_MASTER\phase_name.f90"/>
<File RelativePath=".\SRC_MASTER\plasticctr.f"/>
<File RelativePath=".\SRC_MASTER\plot.f"/>
<File RelativePath=".\SRC_MASTER\plot_co2.f"/>
<File RelativePath=".\SRC_MASTER\plot_hyd.f"/>
<File RelativePath=".\SRC_MASTER\plot_new_2c.f"/>
<File RelativePath=".\SRC_MASTER\plotc1.f"/>
<File RelativePath=".\SRC_MASTER\pod_derivatives.f"/>
<File RelativePath=".\SRC_MASTER\PorosCoeff.f90"/>
<File RelativePath=".\SRC_MASTER\porosi.f"/>
<File RelativePath=".\SRC_MASTER\porosity_gradient_log.f"/>
<File RelativePath=".\SRC_MASTER\porosity_gradient_omr.f"/>
<File RelativePath=".\SRC_MASTER\porosity_wrt_displacements.f"/>
<File RelativePath=".\SRC_MASTER\psat.f"/>
<File RelativePath=".\SRC_MASTER\psatl.f"/>
<File RelativePath=".\SRC_MASTER\ptrac1.f"/>
<File RelativePath=".\SRC_MASTER\ptrac3.f"/>
<File RelativePath=".\SRC_MASTER\radius.f"/>
<File RelativePath=".\SRC_MASTER\ran_sp.f"/>
<File RelativePath=".\SRC_MASTER\random_walk.f"/>
<File RelativePath=".\SRC_MASTER\random_walk_omr.f"/>
<File RelativePath=".\SRC_MASTER\rarng.f"/>
<File RelativePath=".\SRC_MASTER\rcm.f"/>
<File RelativePath=".\SRC_MASTER\rd1dof.f"/>
<File RelativePath=".\SRC_MASTER\rdcon.f"/>
<File RelativePath=".\SRC_MASTER\rddpdp.f"/>
<File RelativePath=".\SRC_MASTER\rdof_new.f"/>
<File RelativePath=".\SRC_MASTER\rdthick.f"/>
<File RelativePath=".\SRC_MASTER\rdtr.f90"/>
<File RelativePath=".\SRC_MASTER\react.f"/>
<File RelativePath=".\SRC_MASTER\read_avs_io_2a.f"/>
<File RelativePath=".\SRC_MASTER\read_rtd.f"/>
<File RelativePath=".\SRC_MASTER\read_rxn.f"/>
<File RelativePath=".\SRC_MASTER\read_sx.f"/>
<File RelativePath=".\SRC_MASTER\read_tcurve.f"/>
<File RelativePath=".\SRC_MASTER\reduce.f90"/>
<File RelativePath=".\SRC_MASTER\releasemem.f"/>
<File RelativePath=".\SRC_MASTER\renum_1.f"/>
<File RelativePath=".\SRC_MASTER\renumber_array.f"/>
<File RelativePath=".\SRC_MASTER\resettrc.f"/>
<File RelativePath=".\SRC_MASTER\resetv_2e.f"/>
<File RelativePath=".\SRC_MASTER\residual.f"/>
<File RelativePath=".\SRC_MASTER\river_ctr.f"/>
<File RelativePath=".\SRC_MASTER\rlp_cap_2a.f90"/>
<File RelativePath=".\SRC_MASTER\rlp_cap_table.f"/>
<File RelativePath=".\SRC_MASTER\rlp_frac.f"/>
<File RelativePath=".\SRC_MASTER\rlp_frac2.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\rlperm_2h.f"/>
<File RelativePath=".\SRC_MASTER\rlperm_co2.f"/>
<File RelativePath=".\SRC_MASTER\rlperm_hyd.f"/>
<File RelativePath=".\SRC_MASTER\rlperm_wtsi.f"/>
<File RelativePath=".\SRC_MASTER\rootls.f"/>
<File RelativePath=".\SRC_MASTER\rotate.f"/>
<File RelativePath=".\SRC_MASTER\saltctr.f"/>
<File RelativePath=".\SRC_MASTER\scaled_perm.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\scanin_b7.f"/>
14. Modified code so initial liquid mass fraction is not read in ‘air’ keyword but in ‘pres’ keyword
<File RelativePath=".\SRC_MASTER\set_mptr.f"/>
<File RelativePath=".\SRC_MASTER\set_ptrk.f"/>
<File RelativePath=".\SRC_MASTER\setbit.f"/>
<File RelativePath=".\SRC_MASTER\setconnarray_3.f"/>
<File RelativePath=".\SRC_MASTER\setord.f"/>
<File RelativePath=".\SRC_MASTER\setparams_2c.f"/>
<File RelativePath=".\SRC_MASTER\setunits.f"/>
<File RelativePath=".\SRC_MASTER\setup_edgePointers_3D.f"/>
<File RelativePath=".\SRC_MASTER\setup_NodeElems.f"/>
<File RelativePath=".\SRC_MASTER\setzone.f"/>
<File RelativePath=".\SRC_MASTER\sfn2r.f"/>
<File RelativePath=".\SRC_MASTER\sfn3r.f"/>
<File RelativePath=".\SRC_MASTER\shap2r.f"/>
<File RelativePath=".\SRC_MASTER\shap3p.f"/>
<File RelativePath=".\SRC_MASTER\shap3r.f"/>
<File RelativePath=".\SRC_MASTER\sice.f"/>
<File RelativePath=".\SRC_MASTER\simplify_ncon.f"/>
<File RelativePath=".\SRC_MASTER\skip_macro.f"/>
<File RelativePath=".\SRC_MASTER\slvesu.f"/>
<File RelativePath=".\SRC_MASTER\solstore.f"/>
<File RelativePath=".\SRC_MASTER\solve_cubic.f90"/>
<File RelativePath=".\SRC_MASTER\solve_dual.f"/>
<File RelativePath=".\SRC_MASTER\solve_new.f"/>
<File RelativePath=".\SRC_MASTER\solve_rdof.f"/>
<File RelativePath=".\SRC_MASTER\sort_values.f"/>
<File RelativePath=".\SRC_MASTER\sortranl.f"/>
<File RelativePath=".\SRC_MASTER\split.f"/>
<File RelativePath=".\SRC_MASTER\sptr_save.f"/>
<File RelativePath=".\SRC_MASTER\ss_trans.f"/>
<File RelativePath=".\SRC_MASTER\start_macro.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_1\startup_2i.f"/>
15. Added coding to set phase state(ieos) based on initial conditions (including iso 2phase w sol)
<File RelativePath=".\SRC_MASTER\steady_4.f"/>
<File RelativePath=".\SRC_MASTER\sther_2.f"/>
<File RelativePath=".\SRC_MASTER\stone_rlp.f"/>
<File RelativePath=".\SRC_MASTER\storag.f"/>
<File RelativePath=".\SRC_MASTER\storage_derivatives.f"/>
<File RelativePath=".\SRC_MASTER\storsx.f"/>
<File RelativePath=".\SRC_MASTER\storsx_read.f"/>
<File RelativePath=".\SRC_MASTER\storsx_write.f"/>
<File RelativePath=".\SRC_MASTER\stream_tube.f"/>
<File RelativePath=".\SRC_MASTER\stress.f"/>
<File RelativePath=".\SRC_MASTER\stress_2D_post.f"/>
<File RelativePath=".\SRC_MASTER\stress_2D_post_fem.f"/>
<File RelativePath=".\SRC_MASTER\stress_3D_post.f"/>
<File RelativePath=".\SRC_MASTER\stress_3D_post_fem.f"/>
<File RelativePath=".\SRC_MASTER\stress_combine.f"/>
<File RelativePath=".\SRC_MASTER\stress_fluid_mech_props.f"/>
<File RelativePath=".\SRC_MASTER\stress_mech_props.f"/>
<File RelativePath=".\SRC_MASTER\stress_perm.f"/>
<File RelativePath=".\SRC_MASTER\stress_uncoupled.f"/>
<File RelativePath=".\SRC_MASTER\stressctr.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_1.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_11.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_2.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_21.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_22.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_222.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_23.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_24.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_25.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_3.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_31.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_4.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_5.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_6.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_7.f90"/>
<File RelativePath=".\SRC_MASTER\stressperm_8.f90"/>
<File RelativePath=".\SRC_MASTER\stressperm_91.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_allocate.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_allocate_3.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_deallocate_3.f"/>
<File RelativePath=".\SRC_MASTER\stressperm_setup.f90"/>
<File RelativePath=".\SRC_MASTER\structured_2a.f"/>
<File RelativePath=".\SRC_MASTER\sub_av.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgs1.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgs2.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgs3.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgstab1.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgstabn1.f"/>
<File RelativePath=".\SRC_MASTER\sub_bcgstabnn.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub1.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub2.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub3.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub31.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub4.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub5.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub6.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksub62.f"/>
<File RelativePath=".\SRC_MASTER\sub_bksubn.f"/>
<File RelativePath=".\SRC_MASTER\Sub_FEHM_tec_to_vtk_version_2o.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres1.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres2.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres3.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres31.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres4.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres5.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres6.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmres62.f"/>
<File RelativePath=".\SRC_MASTER\sub_gmresn.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu1.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu2.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu3.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu4.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu5.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilu6.f"/>
<File RelativePath=".\SRC_MASTER\sub_ilun.f"/>
<File RelativePath=".\SRC_MASTER\sub_lapack.f"/>
<File RelativePath=".\SRC_MASTER\subdivide_rtd.f"/>
<File RelativePath=".\SRC_MASTER\submodel_bc_2d.f"/>
<File RelativePath=".\SRC_MASTER\subsidence.f"/>
<File RelativePath=".\SRC_MASTER\svdcmp_new.f"/>
<File RelativePath=".\SRC_MASTER\switch.f"/>
<File RelativePath=".\SRC_MASTER\switchb.f"/>
<File RelativePath=".\SRC_MASTER\sx_combine_2da.f"/>
<File RelativePath=".\SRC_MASTER\sx_combine_ani.f"/>
<File RelativePath=".\SRC_MASTER\tec_write_grid.f"/>
<File RelativePath=".\SRC_MASTER\termin_2.f"/>
<File RelativePath=".\SRC_MASTER\termio.f"/>
<File RelativePath=".\SRC_MASTER\ther_co2_h2o.f"/>
<File RelativePath=".\SRC_MASTER\ther_meth_h2o.f"/>
<File RelativePath=".\SRC_MASTER\thermc.f"/>
<File RelativePath=".\SRC_MASTER\thermw_prop_2g.f"/>
<File RelativePath=".\SRC_MASTER\thickness.f"/>
<File RelativePath=".\SRC_MASTER\thomeercap.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\thrair_3q.f"/>
16. Modified boundary condition (kq = -202), corrected but in solubility calculation
<File RelativePath=".\SRC_MASTER\thrmwc_9_prop7eb.f"/>
<File RelativePath=".\SRC_MASTER\thrsznapl_2d.f"/>
<File RelativePath=".\SRC_MASTER\timcrl.f"/>
<File RelativePath=".\SRC_MASTER\time_adjust.f"/>
<File RelativePath=".\SRC_MASTER\time_delay.f"/>
<File RelativePath=".\SRC_MASTER\time_diff.f"/>
<File RelativePath=".\SRC_MASTER\time_disp2.f"/>
<File RelativePath=".\SRC_MASTER\tracout.f"/>
<File RelativePath=".\SRC_MASTER\tree_search.f"/>
<File RelativePath=".\SRC_MASTER\tree_search_porosity.f"/>
<File RelativePath=".\SRC_MASTER\trxnvars.f90"/>
<File RelativePath=".\SRC_MASTER\tyming.f"/>
<File RelativePath=".\SRC_MASTER\update_permfactors.f"/>
<File RelativePath=".\SRC_MASTER\user_ymp1.f"/>
<File RelativePath=".\SRC_MASTER\userc.f"/>
<File RelativePath=".\SRC_MASTER\uz_wtctr.f"/>
<File RelativePath=".\SRC_MASTER\vaporl.f"/>
<File RelativePath=".\SRC_MASTER\varchk_12f.f"/>
<File RelativePath=".\SRC_MASTER\varchk_AWH_a3.f"/>
<File RelativePath=".\SRC_MASTER\varchk_part.f"/>
<File RelativePath=".\SRC_MASTER\varchk_simple_awh_2c.f"/>
<File RelativePath=".\SRC_MASTER\varctr.f"/>
<File RelativePath=".\SRC_MASTER\vboun.f"/>
<File RelativePath=".\SRC_MASTER\vcap-ek.f"/>
<File RelativePath=".\SRC_MASTER\vcon.f"/>
<File RelativePath=".\SRC_MASTER\veloc.f"/>
<File RelativePath=".\SRC_MASTER\velocities.f"/>
<File RelativePath=".\SRC_MASTER\velocity_derivatives.f"/>
<File RelativePath=".\SRC_MASTER\velocity_derivatives_lsq.f"/>
<File RelativePath=".\SRC_MASTER\vfcal.f"/>
<File RelativePath=".\SRC_MASTER\vg_regions.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_2e.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_fit.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_fit3.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_inv_calc_2d.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_inv_test.f"/>
<File RelativePath=".\SRC_MASTER\vgcap_match.f"/>
<File RelativePath=".\SRC_MASTER\vgrlp2.f"/>
<File RelativePath=".\SRC_MASTER\vgrlp_2a.f"/>
<File RelativePath=".\SRC_MASTER\vgrlpa.f"/>
<File RelativePath=".\SRC_MASTER_UPDATE_2\vgrlps_2f.f"/>
17. Modified code so upper cutoff is passed through subroutine header
<File RelativePath=".\SRC_MASTER\vonMises_stiffness.f"/>
<File RelativePath=".\SRC_MASTER\vonMises_stress_update.f"/>
<File RelativePath=".\SRC_MASTER\vrock_ctr.f"/>
<File RelativePath=".\SRC_MASTER\WaterProps.f90"/>
<File RelativePath=".\SRC_MASTER\wellimped_ctr_2d.f"/>
<File RelativePath=".\SRC_MASTER\wellphysicsctr_2d.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_head_s_2g.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_con_2c.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_h.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_hf_3.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_mat_s_2i.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_s_2i.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_node_v_2a.f"/>
<File RelativePath=".\SRC_MASTER\write_avs_ucd_header.f"/>
<File RelativePath=".\SRC_MASTER\write_copyright.f90"/>
<File RelativePath=".\SRC_MASTER\write_rlp_hyd.f"/>
<File RelativePath=".\SRC_MASTER\writeio.f"/>
<File RelativePath=".\SRC_MASTER\wrtcon.f"/>
<File RelativePath=".\SRC_MASTER\wrtout_2d.f"/>
<File RelativePath=".\SRC_MASTER\wrtptrk.f"/>
<File RelativePath=".\SRC_MASTER\wtrise.f"/>
<File RelativePath=".\SRC_MASTER\wtsi_column.f"/>
<File RelativePath=".\SRC_MASTER\wtsictr.f"/>
<File RelativePath=".\SRC_MASTER\zeolites.f"/>
<File RelativePath=".\SRC_MASTER\zone_2g.f"/></Filter></Files>