documented new plot features in SkyTEM gradbased

GeoscienceAustralia · Sep 13, 2023 · cc6d58e · cc6d58e
1 parent 50dcbe1
commit cc6d58e
Show file tree

Hide file tree

Showing 10 changed files with 33 additions and 31 deletions.
diff --git a/examples/SkyTEM1D/Menindee/McMC/02_set_options.jl b/examples/SkyTEM1D/Menindee/McMC/02_set_options.jl
@@ -19,12 +19,12 @@ nmin, nmax = 2, 40
 K = transD_GP.GP.OrstUhn()
 demean = false
 sampledc = true
-fbounds = [-1 3.0]
-sdpos = 0.05
-sdprop = 0.05
-sddc = 0.01
-λ, δ = [2], 0.1
-save_freq = 100
+fbounds = [-1 3.0] # log10 resistivity
+sdpos = 0.05 # fraction of fbounds
+sdprop = 0.05 # fraction of nlayers
+sddc = 0.01 # fraction of fbounds
+λ, δ = [2], 0.1 # GP scale length in layers, and nugget in log10 
+save_freq = 50
 Tmax = 2.50
 nchainsatone = 1
 restart = false

diff --git a/examples/SkyTEM1D/Menindee/gradientbased/02_set_options.jl b/examples/SkyTEM1D/Menindee/gradientbased/02_set_options.jl
@@ -11,8 +11,8 @@ nfreqsperdecade    = 5
 regtype            = :R1
 nstepsmax          = 40
 ntries             = 6
-lo                 = -2.5
-hi                 = 0.5
+lo                 = -2.5 # log10 conductivity
+hi                 = 0.5  # log10 conductivity 
 λ²min              = -0.5
 λ²max              = 8
 β²                 = 0.1

diff --git a/examples/SkyTEM1D/Menindee/gradientbased/04_plot.jl b/examples/SkyTEM1D/Menindee/gradientbased/04_plot.jl
@@ -1,6 +1,8 @@
-transD_GP.plotconvandlast(soundings, 50, 2; zall, 
-    figsize=(15,7), vmin=lo, vmax=hi, postfix="_β²_$(β²)_$(regtype)_bg_$(round(10. ^σ0, sigdigits=4))Spm", prefix=zipsaveprefix,
+dr, dz = 50, 2 # interpolation in m
+lnames = [200613] # array of lines wanted in survey to plot
+idx = [[1, 5]] # array of index arrays with indices per line in lnames 
+transD_GP.plotconvandlast(soundings, dr, dz; zall, 
+    figsize=(15,7), vmin=lo, vmax=hi, postfix="_β²_$(β²)_$(regtype)_bg_$(round(10. ^σ0, sigdigits=4))Spm", 
+    prefix=zipsaveprefix, lnames, idx, plotforward=true, aem_in=aem,
     preferEright=true, showplot=true, logscale=true)
-@info "writing PDF ..."
-run(`./pdfscript.sh`)
 @info "done"
diff --git a/examples/VTEM/DarlingParoo/McMC/02_set_options.jl b/examples/VTEM/DarlingParoo/McMC/02_set_options.jl
@@ -20,12 +20,12 @@ nmin, nmax = 2, 40
 K = transD_GP.GP.OrstUhn()
 demean = false
 sampledc = true
-fbounds = [-1 3.0]
-sdpos = 0.05
-sdprop = 0.05
-sddc = 0.01
-λ, δ = [2], 0.1
-save_freq = 100
+fbounds = [-1 3.0] # log10 resistivity
+sdpos = 0.05 # fraction of fbounds
+sdprop = 0.05 # fraction of nlayers
+sddc = 0.01 # fraction of fbounds
+λ, δ = [2], 0.1 # GP scale length in layers, and nugget in log10
+save_freq = 50
 Tmax = 2.50
 nchainsatone = 1
 ## plot 1 random sounding and a background response

diff --git a/examples/VTEM/DarlingParoo/McMC/04_plot.jl b/examples/VTEM/DarlingParoo/McMC/04_plot.jl
@@ -1,5 +1,5 @@
 dr = 20
-idx = 1:100:length(soundings)
+idx = [1,2]
 burninfrac = 0.5
 vmin, vmax = -2.5, 0.5
 transD_GP.summaryAEMimages(soundings, opt;

diff --git a/examples/VTEM/DarlingParoo/gradientbased/02_set_options.jl b/examples/VTEM/DarlingParoo/gradientbased/02_set_options.jl
@@ -12,8 +12,8 @@ nfreqsperdecade    = 5
 regtype            = :R1
 nstepsmax          = 40
 ntries             = 6
-lo                 = -3.
-hi                 = 1.
+lo                 = -3. # log10 conductivity
+hi                 = 1.  # log10 conductivity
 λ²min              = -0.5
 λ²max              = 8
 β²                 = 0.01

diff --git a/examples/tempest/Menindee/McMC/02_set_options.jl b/examples/tempest/Menindee/McMC/02_set_options.jl
@@ -16,12 +16,12 @@ nmin, nmax = 2, 40
 K = transD_GP.GP.OrstUhn()
 demean = false
 sampledc = true
-fbounds = [-1.0 3.0]
-sddc = 0.01 # fraction
-sdpos = 0.05 # fraction
-sdprop = 0.05 # fraction
-λ, δ = [2.], 0.1
-save_freq = 25
+fbounds = [-1.0 3.0] # log 10 resistivity
+sddc = 0.01 # fraction of fbounds
+sdpos = 0.05 # fraction of nlayers 
+sdprop = 0.05 # fraction of fbounds 
+λ, δ = [2.], 0.1 # GP scale length in layers, and nugget in log10
+save_freq = 50
 Tmax = 2.50
 nchainsatone = 1
 vectorsum = true

diff --git a/examples/tempest/Menindee/McMC/04_plot.jl b/examples/tempest/Menindee/McMC/04_plot.jl
@@ -1,5 +1,5 @@
 
-idxplot = [2]
+idxplot = [1,2]
 burninfrac = 0.5
 nbins = 50
 computeforwards, nforwards = true, 2

diff --git a/examples/tempest/Menindee/gradientbased/01_read_data.jl b/examples/tempest/Menindee/gradientbased/01_read_data.jl
@@ -25,7 +25,7 @@ soundings = transD_GP.TEMPEST1DInversion.read_survey_files(;
 	Z = 45,
 	fid = 46,
 	startfrom = 1,
-	skipevery = 1,
+	skipevery = 50,
 	linenum = 47,
 	makeqcplots=true)
 
diff --git a/examples/tempest/Menindee/gradientbased/02_set_options.jl b/examples/tempest/Menindee/gradientbased/02_set_options.jl
@@ -11,8 +11,8 @@ nfreqsperdecade    = 5
 regtype            = :R1
 nstepsmax          = 40
 ntries             = 10
-lo                 = -2.5
-hi                 = 0.5
+lo                 = -2.5 # log10 conductivity
+hi                 = 0.5  # log10 conductivity
 λ²min              = -0.5
 λ²max              = 8
 β²                 = 0.0001