correct bug in paleopiezometer

- Van de Wal piezometer correction
marcoalopez · Apr 9, 2020 · d7c2fbb · d7c2fbb
1 parent c5dc25a
commit d7c2fbb
Show file tree

Hide file tree

Showing 4 changed files with 18 additions and 19 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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()