diff --git a/examples/EAGLE_ICs/EAGLE_100/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_100/vrconfig_3dfof_subhalos_SO_hydro.cfg
index 8590cbf5bc77e8d7a956d210339cced4bbdc692c..bb01aeed3acd0fc104aaf19472e99d3bef6d671d 100644
--- a/examples/EAGLE_ICs/EAGLE_100/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_100/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -1,5 +1,5 @@
 #Configuration file for analysing Hydro
-#runs 3DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
+#runs 6DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
 #Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
 #To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
 
@@ -26,6 +26,7 @@ Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
 #assumes input is in 1e10 msun, Mpc and km/s and output units are the same
 Gravity=43.0211349 #for 1e10 Msun, km/s and Mpc
 Hubble_unit=100.0 # assuming units are km/s and Mpc, then value of Hubble in km/s/Mpc
+
 #converting hydro quantities
 Stellar_age_input_is_cosmological_scalefactor=1
 Metallicity_input_unit_conversion_to_output_unit=1.0
@@ -39,8 +40,8 @@ Length_unit_to_kpc=1000.0
 Velocity_to_kms=1.0
 #conversion of output mass units to solar masses
 Mass_to_solarmass=1.0e10
-#1 / 0.012
-Metallicity_to_solarmetallicity=83.33
+
+Metallicity_to_solarmetallicity=83.33 #1 / 0.012
 Star_formation_rate_to_solarmassperyear=97.78
 Stellar_age_to_yr=1.0
 #ensures that output is physical and not comoving distances per little h
@@ -112,7 +113,7 @@ Halo_core_num_loops=8 #number of loops to iteratively search for cores
 Halo_core_loop_ellx_fac=0.75 #how much to change the configuration space linking per iteration
 Halo_core_loop_ellv_fac=1.0 #how much to change the velocity space linking per iteration
 Halo_core_loop_elln_fac=1.2 #how much to change the min number of particles per iteration
-Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
+Halo_core_phase_significance=5.0 #how significant a core must be in terms of dispersions (sigma) significance
 
 ################################
 #Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
@@ -134,6 +135,7 @@ Frac_pot_ref=1.0
 Min_npot_ref=20
 #reference frame only meaningful if calculating velocity frame using subset of particles in object. Can use radially sorted fraction of particles about minimum potential or centre of mass
 Kinetic_reference_frame_type=0
+#extra options in new unbinding optimisation
 Unbinding_max_unbound_removal_fraction_per_iteration=0.5
 Unbinding_max_unbound_fraction=0.95
 Unbinding_max_unbound_fraction_allowed=0.005
@@ -166,6 +168,23 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass
 Sort_by_binding_energy=1 #sort particles by binding energy
 Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties 
 
+################################
+#Extra fields
+################################
+
+# Also compute the H abundance as well as the HI, HII and H_2 fractions
+Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
+Gas_internal_property_index_in_file=0,0,1,2,
+Gas_internal_property_input_output_unit_conversion_factors=1.0,1.0,1.0,1.0
+Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
+Gas_internal_property_output_units=unitless,unitless,unitless,unitless,
+
+# Collect the BH subgrid masses and compute the max, min, average and total mass in apertures
+BH_internal_property_names=SubgridMasses,SubgridMasses,SubgridMasses,SubgridMasses,
+BH_internal_property_input_output_unit_conversion_factors=1.0e10,1.0e10,1.0e10,1.0e10,
+BH_internal_property_calculation_type=max,min,average,aperture_total,
+BH_internal_property_output_units=solar_mass,solar_mass,solar_mass,solar_mass,
+
 ################################
 #output related
 ################################
@@ -186,6 +205,4 @@ Snapshot_value=SNAP
 ################################
 #other options
 ################################
-Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
-
-
+Verbose=1 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
diff --git a/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
index ad342fa599468bf459df22f9ba145209db476170..bb01aeed3acd0fc104aaf19472e99d3bef6d671d 100644
--- a/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -1,5 +1,5 @@
 #Configuration file for analysing Hydro
-#runs 3DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
+#runs 6DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
 #Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
 #To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
 
@@ -26,6 +26,7 @@ Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
 #assumes input is in 1e10 msun, Mpc and km/s and output units are the same
 Gravity=43.0211349 #for 1e10 Msun, km/s and Mpc
 Hubble_unit=100.0 # assuming units are km/s and Mpc, then value of Hubble in km/s/Mpc
+
 #converting hydro quantities
 Stellar_age_input_is_cosmological_scalefactor=1
 Metallicity_input_unit_conversion_to_output_unit=1.0
@@ -39,8 +40,8 @@ Length_unit_to_kpc=1000.0
 Velocity_to_kms=1.0
 #conversion of output mass units to solar masses
 Mass_to_solarmass=1.0e10
-#1 / 0.012
-Metallicity_to_solarmetallicity=83.33
+
+Metallicity_to_solarmetallicity=83.33 #1 / 0.012
 Star_formation_rate_to_solarmassperyear=97.78
 Stellar_age_to_yr=1.0
 #ensures that output is physical and not comoving distances per little h
@@ -112,7 +113,7 @@ Halo_core_num_loops=8 #number of loops to iteratively search for cores
 Halo_core_loop_ellx_fac=0.75 #how much to change the configuration space linking per iteration
 Halo_core_loop_ellv_fac=1.0 #how much to change the velocity space linking per iteration
 Halo_core_loop_elln_fac=1.2 #how much to change the min number of particles per iteration
-Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
+Halo_core_phase_significance=5.0 #how significant a core must be in terms of dispersions (sigma) significance
 
 ################################
 #Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
@@ -134,6 +135,7 @@ Frac_pot_ref=1.0
 Min_npot_ref=20
 #reference frame only meaningful if calculating velocity frame using subset of particles in object. Can use radially sorted fraction of particles about minimum potential or centre of mass
 Kinetic_reference_frame_type=0
+#extra options in new unbinding optimisation
 Unbinding_max_unbound_removal_fraction_per_iteration=0.5
 Unbinding_max_unbound_fraction=0.95
 Unbinding_max_unbound_fraction_allowed=0.005
@@ -203,6 +205,4 @@ Snapshot_value=SNAP
 ################################
 #other options
 ################################
-Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
-
-
+Verbose=1 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
diff --git a/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
index ad342fa599468bf459df22f9ba145209db476170..bb01aeed3acd0fc104aaf19472e99d3bef6d671d 100644
--- a/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -1,5 +1,5 @@
 #Configuration file for analysing Hydro
-#runs 3DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
+#runs 6DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
 #Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
 #To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
 
@@ -26,6 +26,7 @@ Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
 #assumes input is in 1e10 msun, Mpc and km/s and output units are the same
 Gravity=43.0211349 #for 1e10 Msun, km/s and Mpc
 Hubble_unit=100.0 # assuming units are km/s and Mpc, then value of Hubble in km/s/Mpc
+
 #converting hydro quantities
 Stellar_age_input_is_cosmological_scalefactor=1
 Metallicity_input_unit_conversion_to_output_unit=1.0
@@ -39,8 +40,8 @@ Length_unit_to_kpc=1000.0
 Velocity_to_kms=1.0
 #conversion of output mass units to solar masses
 Mass_to_solarmass=1.0e10
-#1 / 0.012
-Metallicity_to_solarmetallicity=83.33
+
+Metallicity_to_solarmetallicity=83.33 #1 / 0.012
 Star_formation_rate_to_solarmassperyear=97.78
 Stellar_age_to_yr=1.0
 #ensures that output is physical and not comoving distances per little h
@@ -112,7 +113,7 @@ Halo_core_num_loops=8 #number of loops to iteratively search for cores
 Halo_core_loop_ellx_fac=0.75 #how much to change the configuration space linking per iteration
 Halo_core_loop_ellv_fac=1.0 #how much to change the velocity space linking per iteration
 Halo_core_loop_elln_fac=1.2 #how much to change the min number of particles per iteration
-Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
+Halo_core_phase_significance=5.0 #how significant a core must be in terms of dispersions (sigma) significance
 
 ################################
 #Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
@@ -134,6 +135,7 @@ Frac_pot_ref=1.0
 Min_npot_ref=20
 #reference frame only meaningful if calculating velocity frame using subset of particles in object. Can use radially sorted fraction of particles about minimum potential or centre of mass
 Kinetic_reference_frame_type=0
+#extra options in new unbinding optimisation
 Unbinding_max_unbound_removal_fraction_per_iteration=0.5
 Unbinding_max_unbound_fraction=0.95
 Unbinding_max_unbound_fraction_allowed=0.005
@@ -203,6 +205,4 @@ Snapshot_value=SNAP
 ################################
 #other options
 ################################
-Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
-
-
+Verbose=1 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
diff --git a/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
index ad342fa599468bf459df22f9ba145209db476170..bb01aeed3acd0fc104aaf19472e99d3bef6d671d 100644
--- a/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -1,5 +1,5 @@
 #Configuration file for analysing Hydro
-#runs 3DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
+#runs 6DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
 #Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
 #To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
 
@@ -26,6 +26,7 @@ Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
 #assumes input is in 1e10 msun, Mpc and km/s and output units are the same
 Gravity=43.0211349 #for 1e10 Msun, km/s and Mpc
 Hubble_unit=100.0 # assuming units are km/s and Mpc, then value of Hubble in km/s/Mpc
+
 #converting hydro quantities
 Stellar_age_input_is_cosmological_scalefactor=1
 Metallicity_input_unit_conversion_to_output_unit=1.0
@@ -39,8 +40,8 @@ Length_unit_to_kpc=1000.0
 Velocity_to_kms=1.0
 #conversion of output mass units to solar masses
 Mass_to_solarmass=1.0e10
-#1 / 0.012
-Metallicity_to_solarmetallicity=83.33
+
+Metallicity_to_solarmetallicity=83.33 #1 / 0.012
 Star_formation_rate_to_solarmassperyear=97.78
 Stellar_age_to_yr=1.0
 #ensures that output is physical and not comoving distances per little h
@@ -112,7 +113,7 @@ Halo_core_num_loops=8 #number of loops to iteratively search for cores
 Halo_core_loop_ellx_fac=0.75 #how much to change the configuration space linking per iteration
 Halo_core_loop_ellv_fac=1.0 #how much to change the velocity space linking per iteration
 Halo_core_loop_elln_fac=1.2 #how much to change the min number of particles per iteration
-Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
+Halo_core_phase_significance=5.0 #how significant a core must be in terms of dispersions (sigma) significance
 
 ################################
 #Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
@@ -134,6 +135,7 @@ Frac_pot_ref=1.0
 Min_npot_ref=20
 #reference frame only meaningful if calculating velocity frame using subset of particles in object. Can use radially sorted fraction of particles about minimum potential or centre of mass
 Kinetic_reference_frame_type=0
+#extra options in new unbinding optimisation
 Unbinding_max_unbound_removal_fraction_per_iteration=0.5
 Unbinding_max_unbound_fraction=0.95
 Unbinding_max_unbound_fraction_allowed=0.005
@@ -203,6 +205,4 @@ Snapshot_value=SNAP
 ################################
 #other options
 ################################
-Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
-
-
+Verbose=1 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
diff --git a/examples/EAGLE_ICs/EAGLE_6/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_6/vrconfig_3dfof_subhalos_SO_hydro.cfg
index ad342fa599468bf459df22f9ba145209db476170..bb01aeed3acd0fc104aaf19472e99d3bef6d671d 100644
--- a/examples/EAGLE_ICs/EAGLE_6/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_6/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -1,5 +1,5 @@
 #Configuration file for analysing Hydro
-#runs 3DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
+#runs 6DFOF + substructure algorithm, demands subhalos and FOF halos be self-bound, calculates many properties
 #Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
 #To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
 
@@ -26,6 +26,7 @@ Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
 #assumes input is in 1e10 msun, Mpc and km/s and output units are the same
 Gravity=43.0211349 #for 1e10 Msun, km/s and Mpc
 Hubble_unit=100.0 # assuming units are km/s and Mpc, then value of Hubble in km/s/Mpc
+
 #converting hydro quantities
 Stellar_age_input_is_cosmological_scalefactor=1
 Metallicity_input_unit_conversion_to_output_unit=1.0
@@ -39,8 +40,8 @@ Length_unit_to_kpc=1000.0
 Velocity_to_kms=1.0
 #conversion of output mass units to solar masses
 Mass_to_solarmass=1.0e10
-#1 / 0.012
-Metallicity_to_solarmetallicity=83.33
+
+Metallicity_to_solarmetallicity=83.33 #1 / 0.012
 Star_formation_rate_to_solarmassperyear=97.78
 Stellar_age_to_yr=1.0
 #ensures that output is physical and not comoving distances per little h
@@ -112,7 +113,7 @@ Halo_core_num_loops=8 #number of loops to iteratively search for cores
 Halo_core_loop_ellx_fac=0.75 #how much to change the configuration space linking per iteration
 Halo_core_loop_ellv_fac=1.0 #how much to change the velocity space linking per iteration
 Halo_core_loop_elln_fac=1.2 #how much to change the min number of particles per iteration
-Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
+Halo_core_phase_significance=5.0 #how significant a core must be in terms of dispersions (sigma) significance
 
 ################################
 #Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
@@ -134,6 +135,7 @@ Frac_pot_ref=1.0
 Min_npot_ref=20
 #reference frame only meaningful if calculating velocity frame using subset of particles in object. Can use radially sorted fraction of particles about minimum potential or centre of mass
 Kinetic_reference_frame_type=0
+#extra options in new unbinding optimisation
 Unbinding_max_unbound_removal_fraction_per_iteration=0.5
 Unbinding_max_unbound_fraction=0.95
 Unbinding_max_unbound_fraction_allowed=0.005
@@ -203,6 +205,4 @@ Snapshot_value=SNAP
 ################################
 #other options
 ################################
-Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox
-
-
+Verbose=1 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox