allow wpm in md_sea

src/juobs_tools.jl

@@ -388,7 +388,7 @@ m_mdl, m_mds = md_sea(m, [md1, md2], pp, ADerrors.wsg)
 m_shifted = m + 2 * dml * m_mdl + dms * m_mds
 ```
 """
-function md_sea(a::uwreal, md::Vector{Matrix{Float64}}, ow::uwreal, w::Union{uwreal, Nothing}=nothing, ws::ADerrors.wspace=ADerrors.wsg)
+function md_sea(a::uwreal, md::Vector{Matrix{Float64}}, ow::uwreal, w::Union{uwreal, Nothing}=nothing, ws::ADerrors.wspace=ADerrors.wsg, wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)
     nid = neid(a)
     p = findall(t-> t==1, a.prop)
 
@@ -420,7 +420,7 @@ function md_sea(a::uwreal, md::Vector{Matrix{Float64}}, ow::uwreal, w::Union{uwr
     end
     
     nrw = size(md_aux, 1)
-    uwerr(ow)
+    isnothing(wpm) ? uwerr(ow) : uwerr(ow,wpm)
     ow_data = ow.mean .+ mchist(ow, id)
 
     if isnothing(w)