vertical cross sections via GMG Vote Maps

[amrelsharkawy]

I’ve created vote maps following the instructions from GMG' GitHub, and I’ve attached a Paraview screenshot for your reference. However, when I try to create a slice along the profile through the 3D vote map, I encounter an error (pasted below). Do you have any suggestions for resolving this?

+++ The error message I got +++

---

********* Profile 1 *******

---

lon start 9.0
lat start 43.0
lon end 15.0
lat end 49.0
ERROR: LoadError: InexactError: Int64(NaN)
Stacktrace:
  [1] Int64
    @ ./float.jl:994 [inlined]
  [2] convert
    @ ./number.jl:7 [inlined]
  [3] _round_convert
    @ ./rounding.jl:480 [inlined]
  [4] round
    @ ./rounding.jl:479 [inlined]
  [5] round
    @ ./rounding.jl:477 [inlined]
  [6] _broadcast_getindex_evalf
    @ ./broadcast.jl:678 [inlined]
  [7] _broadcast_getindex
    @ ./broadcast.jl:661 [inlined]
  [8] getindex
    @ ./broadcast.jl:610 [inlined]
  [9] macro expansion
    @ ./broadcast.jl:973 [inlined]
 [10] macro expansion
    @ ./simdloop.jl:77 [inlined]
 [11] copyto!
    @ ./broadcast.jl:972 [inlined]
 [12] copyto!
    @ ./broadcast.jl:925 [inlined]
 [13] copy
    @ ./broadcast.jl:897 [inlined]
 [14] materialize(bc::Base.Broadcast.Broadcasted{Base.Broadcast.DefaultArrayStyle{3}, Nothing, typeof(round), Tuple{Base.RefValue{Type{Int64}}, Array{Float64, 3}}})
    @ Base.Broadcast ./broadcast.jl:872
 [15] interpolate_datafields(V::GeoData, Lon::Array{Float64, 3}, Lat::Array{Float64, 3}, Depth::Array{Float64, 3})
    @ GeophysicalModelGenerator ~/.julia/packages/GeophysicalModelGenerator/Q1va2/src/utils.jl:1067
 [16] combine_vol_data(VolData::@NamedTuple{Zhu2015::GeoData, Giacomuzzi2023::GeoData, UU07_Amaru2007::GeoData, MIT08_Li2008::GeoData, Piromallo2003::GeoData, Koulakov2009::GeoData, Paffrath2021::GeoData, Zhao2016::GeoData, Timko2023_Mantle::GeoData, Lippitsch2003::GeoData, Mitterbauer2011::GeoData, Plomerova2022::GeoData, VoteMap_3D::GeoData}; lat::Nothing, lon::Nothing, depth::Nothing, dims::Tuple{Int64, Int64, Int64}, dataset_preferred::Int64)
    @ GeophysicalModelGenerator ~/.julia/packages/GeophysicalModelGenerator/Q1va2/src/ProfileProcessing.jl:245
 [17] combine_vol_data(VolData::@NamedTuple{Zhu2015::GeoData, Giacomuzzi2023::GeoData, UU07_Amaru2007::GeoData, MIT08_Li2008::GeoData, Piromallo2003::GeoData, Koulakov2009::GeoData, Paffrath2021::GeoData, Zhao2016::GeoData, Timko2023_Mantle::GeoData, Lippitsch2003::GeoData, Mitterbauer2011::GeoData, Plomerova2022::GeoData, VoteMap_3D::GeoData})
    @ GeophysicalModelGenerator ~/.julia/packages/GeophysicalModelGenerator/Q1va2/src/ProfileProcessing.jl:218
 [18] extract_ProfileData(ProfileCoordFile::String, ProfileNumber::Int64, DataSetFile::String; DimsVolCross::Tuple{Int64, Int64}, DepthVol::Nothing, DimsSurfCross::Tuple{Int64}, WidthPointProfile::Unitful.Quantity{Int64, 𝐋, Unitful.FreeUnits{(km,), 𝐋, nothing}})
    @ GeophysicalModelGenerator ~/.julia/packages/GeophysicalModelGenerator/Q1va2/src/ProfileProcessing.jl:422
 [19] extract_ProfileData
    @ ~/.julia/packages/GeophysicalModelGenerator/Q1va2/src/ProfileProcessing.jl:404 [inlined]
 [20] top-level scope
    @ /meissnerdata/seismologie/models/GMG_Datasets/PickingData/MediT/Process_easternAlps.jl:67
in expression starting at /meissnerdata/seismologie/models/GMG_Datasets/PickingData/MediT/Process_easternAlps.jl:59

[albert-de-montserrat]

Could you share the exact lines of code to reproduce this error? Some NaNs where produced at some point.