diff --git a/contrib/world_builder/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc b/contrib/world_builder/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc index 4794a54f5da..de1ab11921e 100644 --- a/contrib/world_builder/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc +++ b/contrib/world_builder/source/world_builder/features/oceanic_plate_models/composition/tian2019_water_content.cc @@ -178,8 +178,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); - - partition_coefficient = std::min(max_water_content, partition_coefficient); + // 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/contrib/world_builder/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc b/contrib/world_builder/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc index 0835b4c4c58..f0407b08bb1 100644 --- a/contrib/world_builder/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc +++ b/contrib/world_builder/source/world_builder/features/subducting_plate_models/composition/tian2019_water_content.cc @@ -171,8 +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); - - partition_coefficient = std::min(max_water_content, partition_coefficient); + // 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) {