Skip to content

Commit

Permalink
correct bug in paleopiezometer
Browse files Browse the repository at this point in the history
- Van de Wal piezometer correction
  • Loading branch information
marcoalopez committed Apr 9, 2020
1 parent c5dc25a commit d7c2fbb
Show file tree
Hide file tree
Showing 4 changed files with 18 additions and 19 deletions.
8 changes: 4 additions & 4 deletions DOCS/brief_tutorial.md
Original file line number Diff line number Diff line change
Expand Up @@ -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 |



Expand Down
8 changes: 4 additions & 4 deletions DOCS/getting_started.md
Original file line number Diff line number Diff line change
Expand Up @@ -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 |



Expand Down
13 changes: 6 additions & 7 deletions grain_size_tools/GrainSizeTools_script.py
Original file line number Diff line number Diff line change
Expand Up @@ -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

Expand Down
8 changes: 4 additions & 4 deletions grain_size_tools/piezometers.py
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand All @@ -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()
Expand Down

0 comments on commit d7c2fbb

Please sign in to comment.