fix a typo, add comments

GeodynamicWorldBuilder · Nov 12, 2024 · d4c29b7 · d4c29b7
1 parent 8a2b88a
commit d4c29b7
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 2 deletions.
diff --git a/doc/JOSS/1.0/paper.md b/doc/JOSS/1.0/paper.md
@@ -48,7 +48,7 @@ affiliations:
    index: 2
  - name: GEOMAR Helmholtz-Zentrum für Ozeanforschung, Kiel, Germany
    index: 3
- - name: Universiy of Florida, USA
+ - name: University of Florida, USA
    index: 4
  - name: Clemson University, USA
    index: 5

diff --git a/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc b/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc
@@ -179,6 +179,8 @@ namespace WorldBuilder
                   double partition_coefficient = calculate_water_content(lithostatic_pressure,
                                                                          slab_temperature);
 
+                  // The partition_coefficient is output as a percentage, but geodynamic modeling software
+                  // typically deal with fractions, so we divide by 100 below
                   partition_coefficient = std::min(max_water_content, partition_coefficient) / 100;
 
                   for (unsigned int i = 0; i < compositions.size(); ++i)

diff --git a/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc b/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc
@@ -171,7 +171,9 @@ namespace WorldBuilder
               const double slab_temperature = world->properties(position_in_cartesian_coordinates.get_array(), depth, {{{1,0,0}}})[0];
               double partition_coefficient = calculate_water_content(lithostatic_pressure,
                                                                      slab_temperature);
-
+
+              // The partition_coefficient is output as a percentage, but geodynamic modeling software
+              // typically deal with fractions, so we divide by 100 below
               partition_coefficient = std::min(max_water_content, partition_coefficient) / 100;
 
               for (unsigned int i = 0; i < compositions.size(); ++i)