From e0f4c757d09cc69e3b2f832e8e49f45390f6d2e5 Mon Sep 17 00:00:00 2001
From: Thibault Duretz <48455309+tduretz@users.noreply.github.com>
Date: Sun, 23 Jun 2024 09:23:48 +0200
Subject: [PATCH] Add necking review setups (#127)

* recovered NeckingReview and added RiftingRoman

* edit

* done with adding setup
---
 JuliaVisualisation/Main_Visualisation_Makie_MD7.jl  | 2 +-
 SETS/NeckinReviewSystematics/NeckingReview_NR01.txt | 8 ++++----
 SETS/NeckinReviewSystematics/NeckingReview_NR03.txt | 8 ++++----
 SETS/NeckinReviewSystematics/NeckingReview_NR07.txt | 6 +++---
 SETS/NeckinReviewSystematics/NeckingReview_NR09.txt | 8 ++++----
 SETS/NeckingReview.c                                | 1 +
 6 files changed, 17 insertions(+), 16 deletions(-)

diff --git a/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl b/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl
index 0d47fe47..f1757560 100644
--- a/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl
+++ b/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl
@@ -85,7 +85,7 @@ function main()
     # field = :EffectiveFrictionTime
 
     # Switches
-    printfig    = false  # print figures to disk
+    printfig    = true  # print figures to disk
     printvid    = false
     framerate   = 3
     PlotOnTop = (
diff --git a/SETS/NeckinReviewSystematics/NeckingReview_NR01.txt b/SETS/NeckinReviewSystematics/NeckingReview_NR01.txt
index 131d2659..8a1f26b0 100755
--- a/SETS/NeckinReviewSystematics/NeckingReview_NR01.txt
+++ b/SETS/NeckinReviewSystematics/NeckingReview_NR01.txt
@@ -73,7 +73,7 @@ Nb_phases = 4
 
 ID   = 0
 rho  = 2800.00 / Upper Crust - Westerly granite
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -97,7 +97,7 @@ eta_vp   = 22.5e20
 
 ID   = 1
 rho  = 2800.00 / Lower crust
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -121,7 +121,7 @@ ID   = 2
 rho  = 3330.00 / lithospheric mantle - Dry Olivine
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
@@ -143,7 +143,7 @@ ID   = 3
 rho  = 3330.00 / asthenosphere 
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
diff --git a/SETS/NeckinReviewSystematics/NeckingReview_NR03.txt b/SETS/NeckinReviewSystematics/NeckingReview_NR03.txt
index 534fe651..76c23390 100755
--- a/SETS/NeckinReviewSystematics/NeckingReview_NR03.txt
+++ b/SETS/NeckinReviewSystematics/NeckingReview_NR03.txt
@@ -73,7 +73,7 @@ Nb_phases = 4
 
 ID   = 0
 rho  = 2800.00 / Upper Crust - Westerly granite
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -97,7 +97,7 @@ eta_vp   = 2.5e20
 
 ID   = 1
 rho  = 2800.00 / Lower crust
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -121,7 +121,7 @@ ID   = 2
 rho  = 3330.00 / lithospheric mantle - Dry Olivine
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
@@ -143,7 +143,7 @@ ID   = 3
 rho  = 3330.00 / asthenosphere 
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
diff --git a/SETS/NeckinReviewSystematics/NeckingReview_NR07.txt b/SETS/NeckinReviewSystematics/NeckingReview_NR07.txt
index a127783b..218f787b 100755
--- a/SETS/NeckinReviewSystematics/NeckingReview_NR07.txt
+++ b/SETS/NeckinReviewSystematics/NeckingReview_NR07.txt
@@ -97,7 +97,7 @@ eta_vp   = 22.5e20
 
 ID   = 1
 rho  = 2800.00 / Lower crust
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -121,7 +121,7 @@ ID   = 2
 rho  = 3330.00 / lithospheric mantle - Dry Olivine
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
@@ -143,7 +143,7 @@ ID   = 3
 rho  = 3330.00 / asthenosphere 
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
diff --git a/SETS/NeckinReviewSystematics/NeckingReview_NR09.txt b/SETS/NeckinReviewSystematics/NeckingReview_NR09.txt
index 1371a7f2..b2554377 100755
--- a/SETS/NeckinReviewSystematics/NeckingReview_NR09.txt
+++ b/SETS/NeckinReviewSystematics/NeckingReview_NR09.txt
@@ -77,7 +77,7 @@ Nb_phases = 4
 
 ID   = 0
 rho  = 2800.00 / Upper Crust - Westerly granite
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -101,7 +101,7 @@ eta_vp   = 2.5e20
 
 ID   = 1
 rho  = 2800.00 / Lower crust
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 2.5
 Qr   = 0.55e-6
@@ -125,7 +125,7 @@ ID   = 2
 rho  = 3330.00 / lithospheric mantle - Dry Olivine
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
@@ -147,7 +147,7 @@ ID   = 3
 rho  = 3330.00 / asthenosphere 
 alp  = 32.0e-6
 bet  = 1.5e-11
-G   = 1.0e10
+G    = 1.0e10
 Cp   = 1050.0
 k    = 3.0
 Qr   = 0.0
diff --git a/SETS/NeckingReview.c b/SETS/NeckingReview.c
index f4145a20..7e30b6e2 100755
--- a/SETS/NeckingReview.c
+++ b/SETS/NeckingReview.c
@@ -145,6 +145,7 @@ int main(int nargs, char *args[]) {
     asprintf(&input_file, "NeckingReview.txt"); // Default
   }
   else {
+
     asprintf(&input_file, "%s", args[1]);     // Custom
   }
   printf("Running MDoodz7.0 using %s\n", input_file);