fix a typo, add comments

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

Hide file tree

Showing 3 changed files with 5 additions and 3 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
@@ -178,7 +178,8 @@ 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)

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,8 @@ 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)