From acab5732b6f8daeb2895330bf6e20f6387aef325 Mon Sep 17 00:00:00 2001 From: Anandaroop Ray Date: Tue, 5 Dec 2023 18:09:43 +1100 Subject: [PATCH] change to normal grad plots for fixed_wing --- src/AEMnoNuisanceGradientInversionTools.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/AEMnoNuisanceGradientInversionTools.jl b/src/AEMnoNuisanceGradientInversionTools.jl index 8cc422ef..c8560f2a 100644 --- 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)