diff --git a/DOCS/brief_tutorial.md b/DOCS/brief_tutorial.md index 7785136..8544be9 100644 --- a/DOCS/brief_tutorial.md +++ b/DOCS/brief_tutorial.md @@ -438,10 +438,10 @@ It is key to note that different piezometers require entering **different appare **Table 3**. Parameters relating the average apparent size of dynamically recrystallized grains and differential stress in olivine -| Reference | phase | DRX | A†,‡ | p† | B†,‡ | m† | -| :-----------------------: | :----------: | :--: | :---: | :--: | :-----: | :--: | -| Jung and Karato (2001) | olivine, wet | BLG | 25704 | 1.18 | 5461.03 | 0.85 | -| Van der Wal et al. (1993) | olivine, wet | | 14800 | 1.33 | 4250 | 0.75 | +| Reference | phase | DRX | A†,‡ | p† | B†,‡ | m† | +| :-----------------------: | :--------------: | :--: | :---: | :--: | :-----: | :--: | +| Jung and Karato (2001) | olivine, wet | BLG | 25704 | 1.18 | 5461.03 | 0.85 | +| Van der Wal et al. (1993) | olivine, dry/wet | | 14800 | 1.33 | 1355.4 | 0.75 | diff --git a/DOCS/getting_started.md b/DOCS/getting_started.md index 7f595fe..3d3e2a7 100644 --- a/DOCS/getting_started.md +++ b/DOCS/getting_started.md @@ -552,10 +552,10 @@ It is key to note that different piezometers require entering **different appare **Table 3**. Parameters relating the average apparent size of dynamically recrystallized grains and differential stress in olivine -| Reference | phase | DRX | A†,‡ | p† | B†,‡ | m† | -| :-----------------------: | :----------: | :--: | :---: | :--: | :-----: | :--: | -| Jung and Karato (2001) | olivine, wet | BLG | 25704 | 1.18 | 5461.03 | 0.85 | -| Van der Wal et al. (1993) | olivine, wet | | 14800 | 1.33 | 4250 | 0.75 | +| Reference | phase | DRX | A†,‡ | p† | B†,‡ | m† | +| :-----------------------: | :--------------: | :--: | :---: | :--: | :-----: | :--: | +| Jung and Karato (2001) | olivine, wet | BLG | 25704 | 1.18 | 5461.03 | 0.85 | +| Van der Wal et al. (1993) | olivine, dry/wet | | 15000 | 1.33 | 1355.4 | 0.75 | diff --git a/grain_size_tools/GrainSizeTools_script.py b/grain_size_tools/GrainSizeTools_script.py index 9c58c6a..5eb3f29 100644 --- a/grain_size_tools/GrainSizeTools_script.py +++ b/grain_size_tools/GrainSizeTools_script.py @@ -317,17 +317,16 @@ def calc_diffstress(grain_size, phase, piezometer, correction=False): diff_stress = B * grain_size**(-m) * np.exp(698 / (T + 273.15)) if correction is True: diff_stress = diff_stress * 2 / np.sqrt(3) - print(' ') - print('differential stress = {:0.2f} MPa' .format(diff_stress)) - print(warn) + else: diff_stress = B * grain_size**-m if correction is True: diff_stress = diff_stress * 2 / np.sqrt(3) - print(' ') - print('differential stress = {:0.2f} MPa' .format(diff_stress)) - print(warn) - print(' ') + + print(' ') + print('differential stress = {:0.2f} MPa' .format(diff_stress)) + print(warn) + print(' ') return None diff --git a/grain_size_tools/piezometers.py b/grain_size_tools/piezometers.py index b83a8fb..4258433 100644 --- a/grain_size_tools/piezometers.py +++ b/grain_size_tools/piezometers.py @@ -222,8 +222,8 @@ def olivine(piezometer=None): grain size in microns calculated from equivalent circular diameters (ECD) with no stereological correction. The function will convert automatically this value to linear intercept (LI) grain size using the De Hoff and Rhines (1968) correction. - Since LI was originally multiplied by 1.2 (correction factor), the final relation - is: LI = (1.2 / sqrt(4/pi)) * ECD + It is assumed that LI was multiplied by 1.5 (correction factor), the final relation + is: LI = (1.5 / sqrt(4/pi)) * ECD - The piezometer of Jung and Karato (2001) requires entering the linear mean apparent grain size in microns calculated from equivalent circular diameters @@ -246,10 +246,10 @@ def olivine(piezometer=None): correction_factor = 1.5 elif piezometer == 'VanderWal_wet': - B, m = 4250, 0.75 + B, m = 1355.4, 0.75 warn = 'Ensure that you entered the apparent grain size as the arithmetic mean in linear scale!' linear_interceps = True - correction_factor = 1.2 + correction_factor = 1.5 else: olivine()