From 59702394469b653c92d3d07d3c39c55b292d9292 Mon Sep 17 00:00:00 2001 From: Juliane Dannberg Date: Wed, 14 Feb 2024 13:19:45 -0500 Subject: [PATCH] address comments, fix units, add changelog entry --- CHANGELOG.md | 1 + .../3d_cartesian_transform_fault.wb | 18 +++++++++++++++--- .../3d_cartesian_transform_fault/doc/README.md | 8 ++++---- 3 files changed, 20 insertions(+), 7 deletions(-) diff --git a/CHANGELOG.md b/CHANGELOG.md index 0ec3e61ed..d08a2be75 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -12,6 +12,7 @@ This project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.htm ## [Unreleased] ### Added - Added an option to use the plate model as the reference model for the mass conserving temperature of the slab. \[Haoyuan Li; 2024-02-02; [#471](https://github.com/GeodynamicWorldBuilder/WorldBuilder/pull/471) +- Added a cookbook for making a transform fault and using this model in ASPECT. \[Juliane Dannberg; 2024-02-14; [#563](https://github.com/GeodynamicWorldBuilder/WorldBuilder/pull/563) ### Changed diff --git a/cookbooks/3d_cartesian_transform_fault/3d_cartesian_transform_fault.wb b/cookbooks/3d_cartesian_transform_fault/3d_cartesian_transform_fault.wb index 809997b1a..2bbe9eaa6 100644 --- a/cookbooks/3d_cartesian_transform_fault/3d_cartesian_transform_fault.wb +++ b/cookbooks/3d_cartesian_transform_fault/3d_cartesian_transform_fault.wb @@ -6,8 +6,20 @@ "thermal diffusivity":1.06060606e-6, "features": [ - {"model":"oceanic plate", "name":"oceanic plate A", "coordinates":[[-1e3,-1e3],[251e3,-1e3],[251e3,101e3],[-1e3,101e3]], - "temperature models":[{"model":"half space model", "max depth":100e3, "spreading velocity":0.03, "top temperature":273.15, "ridge coordinates":[[[200e3,-1e3],[200e3,50e3]],[[50e3,50e3],[50e3,101e3]]]}] + { + "model":"oceanic plate", + "name":"oceanic plate A", + "coordinates":[[-1e3,-1e3],[251e3,-1e3],[251e3,101e3],[-1e3,101e3]], + "temperature models": + [ + { + "model":"half space model", + "max depth":100e3, + "spreading velocity":0.03, + "top temperature":273.15, + "ridge coordinates":[[[200e3,-1e3],[200e3,50e3]],[[50e3,50e3],[50e3,101e3]]] + } + ] } ] -} +} \ No newline at end of file diff --git a/cookbooks/3d_cartesian_transform_fault/doc/README.md b/cookbooks/3d_cartesian_transform_fault/doc/README.md index 171272d0e..658e7aa4c 100644 --- a/cookbooks/3d_cartesian_transform_fault/doc/README.md +++ b/cookbooks/3d_cartesian_transform_fault/doc/README.md @@ -6,7 +6,7 @@ The goal of this tutorial is to learn how to set up a model of a mid-ocean ridge The first step is to prescribe some global parameters for the GWB. This is important because we use the half-space cooling model to compute the temperature distribution within the oceanic plates, and this model uses the thermal diffusivity. In addition, the GWB needs to know what the adiabatic gradient is to compute the correct temperature, and we need to specify the surface and mantle temperature to be consistent with the study we want to reproduce. -Specifically, we here want to set the `surface temperature` to 0 degrees Celsius (273.25 K) and the `potential mantle temperature` to 1300 degrees Celsius (1573.25 K) as in Behn et al., 2007. Since that study does not include adiabatic heating, we need to set the `thermal expansion coefficient` to zero. To reproduce the temperature profile shown in Figure 2 in Behn et al. (2007), we set the `thermal diffusivity` to 1.06060606e-6 m2/s (which corresponds to a thermal conductivity of 3.5 W/m/K assuming that the density is 3300 kg/m3 and the specific keat is 1000 J/kg/K). Note that we need to make sure that these properties are consistent between the Geodynamic World Builder input file and the input file of the geodynamic modeling code that we are using this file as an initial condition for, in this case ASPECT. +Specifically, we here want to set the `surface temperature` to 0 degrees Celsius (273.25 K) and the `potential mantle temperature` to 1300 degrees Celsius (1573.25 K) as in Behn et al., 2007. Since that study does not include adiabatic heating, we need to set the `thermal expansion coefficient` to zero. To reproduce the temperature profile shown in Figure 2 in Behn et al. (2007), we set the `thermal diffusivity` to 1.06060606e-6 m2/s (which corresponds to a thermal conductivity of 3.5 W/m/K assuming that the density is 3300 kg/m3 and the specific keat is 1000 J/kg/K). Note that we need to make sure that these properties are consistent between the Geodynamic World Builder input file and the input file of the geodynamic modeling code that we are using this file as an initial condition for, in this case ASPECT. ::::::{tab-set} @@ -67,8 +67,8 @@ We also set the [top temperature](open_features_items_oneOf_4_temperature-models ```{literalinclude} ../3d_cartesian_transform_fault.wb :language: json :lineno-start: 7 -:lines: 7-12 -:emphasize-lines: 3,4 +:lines: 7-24 +:emphasize-lines: 4-14 ``` ::::{grid} 2 :::{grid-item-card} 3d_cartesian_transform_fault.wb @@ -87,7 +87,7 @@ We also set the [top temperature](open_features_items_oneOf_4_temperature-models ```{literalinclude} ../3d_cartesian_transform_fault.wb :language: json :lineno-start: 1 -:emphasize-lines: 9,10 +:emphasize-lines: 10-20 ``` ::::{grid} 2