change to normal grad plots for fixed_wing

GeoscienceAustralia · Dec 5, 2023 · acab573 · acab573
1 parent 58b1a37
commit acab573
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/AEMnoNuisanceGradientInversionTools.jl b/src/AEMnoNuisanceGradientInversionTools.jl
@@ -54,6 +54,7 @@ function plotconvandlast(soundings, delr, delz;
         yl = nothing,
         plotforward = false,
         aem_in = nothing,
+        nnu = 0,
         dpi=400)
     linestartidx = splitsoundingsbyline(soundings)                    
     nlines = length(linestartidx)
@@ -62,7 +63,7 @@ function plotconvandlast(soundings, delr, delz;
     isfile(fnamecheck) && compress(soundings, zall) # write everything in one file if not done yet
     fzipped = prefix == "" ? "zipped.dat" : prefix*"_zipped.dat"
     A = readdlm(fzipped)
-    σ = A[:,end-nlayers:end-1]
+    σ = A[:,end-nlayers-nnu:end-nnu-1] # so we can plot TEMPEST and SPECTREM similar to heli
     ϕd = A[:,end]
     for i in 1:nlines
         a, b = linestartend(linestartidx, i, nlines, soundings)