replace resolution by size in julia viz

tduretz · Jun 13, 2024 · 3a21e0e · 3a21e0e
1 parent d05944e
commit 3a21e0e
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 7 deletions.
diff --git a/JuliaVisualisation/Main_Visualisation_FaultSmearing_MD7.jl b/JuliaVisualisation/Main_Visualisation_FaultSmearing_MD7.jl
@@ -78,11 +78,13 @@ function main()
         "/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t3/",
         "/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t3_nonconv/",
         "/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t2_nonconv/",
+        "/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t1_nonconv/",
     )
     names = (
         "T3: Fully converged",
         "T3: 1 iteration",
         "T2: 1 iteration",
+        "T1: 1 iteration",
     )
     n_files = length(path)
 
@@ -149,7 +151,7 @@ function main()
     cm_yr = 100.0*3600.0*24.0*365.25
 
     # Time loop
-    f = Figure(resolution = (Lx/Lz*resol*1.2, resol), fontsize=25)
+    f = Figure(size = (Lx/Lz*resol*1.2, resol), fontsize=25)
 
     for istep=file_start:file_step:file_end
 
@@ -161,7 +163,7 @@ function main()
 
         #####################################
         empty!(f)
-        f = Figure(resolution = (Lx/Lz*resol*1.2, resol), fontsize=25)
+        f = Figure(size = (Lx/Lz*resol*1.2, resol), fontsize=25)
 
         if field==:EffectiveFrictionTime
             for i_file=1:n_files

diff --git a/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl b/JuliaVisualisation/Main_Visualisation_Makie_MD7.jl
@@ -45,6 +45,7 @@ function main()
     path ="/Users/tduretz/REPO/MDOODZ7.0/MDLIB/TEST_ROMAN_ANI3_00_MR/"
     path ="/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t3/"
     path ="/Users/tduretz/REPO/MDOODZ7.0/MDLIB/_p10_e18_t1_nonconv/"
+    path ="/Users/tduretz/REPO/MDOODZ7.0/MDLIB/"
 
     # path ="/Users/tduretz/Downloads/"
     # path ="/Users/tduretz/REPO/MDOODZ7.0/MDLIB/DoubleSubduction_OMP16/"
@@ -57,9 +58,9 @@ function main()
     # path ="/Users/tduretz/REPO/MDOODZ7.0/RUNS/1_NR09/"
 
     # File numbers
-    file_start = 000
+    file_start = 130
     file_step  = 100
-    file_end   = 2000
+    file_end   = 130
 
     # Select field to visualise
     field = :Phases
@@ -89,7 +90,7 @@ function main()
     framerate   = 3
     PlotOnTop = (
         ph_contours = false,  # add phase contours
-        T_contours  = false,   # add temperature contours
+        T_contours  = true,   # add temperature contours
         fabric      = false,  # add fabric quiver (normal to director)
         topo        = false,
         σ1_axis     = false,
@@ -116,7 +117,7 @@ function main()
     cm_yr = 100.0*3600.0*24.0*365.25
 
     # Time loop
-    f = Figure(resolution = (Lx/Lz*resol*1.2, resol), fontsize=25)
+    f = Figure(size = (Lx/Lz*resol*1.2, resol), fontsize=25)
 
     for istep=file_start:file_step:file_end
 
@@ -233,7 +234,7 @@ function main()
 
         #####################################
         empty!(f)
-        f = Figure(resolution = (Lx/Lz*resol*1.2, resol), fontsize=25)
+        f = Figure(size = (Lx/Lz*resol*1.2, resol), fontsize=25)
 
         if field==:Phases
             ax1 = Axis(f[1, 1], title = L"Phases at $t$ = %$(tMy) Ma", xlabel = L"$x$ [m]", ylabel = L"$y$ [m]")