md_val TSM prop=4 case 1st try

src/juobs_tools.jl

@@ -559,6 +559,33 @@ function md_val(a::uwreal, obs::Corr, derm::Vector{Corr}; new_version::Bool=fals
             end
             der = length(a.der) >= maximum(n_ow1) && length(a.der) >= maximum(n_ow2) ? a.der[n_ow1] * w_mean : 0.0
         end
+    elseif all(count.(prop) .== 4)
+        if TSM == false
+            return nothing
+        else
+            corr_der = getfield.(corr, :der)
+            n = findall.(t-> t==1, prop)
+            n = vcat(n'...)
+            if all(n[:, 1] .== n[1, 1]) # find ow and w
+                n_w = n[1, 1]
+                n_ow1 = n[:, 2]
+                n_ow2 = n[:, 3]
+            elseif all(n[:, 2] .== n[1, 2])
+                n_w = n[1, 2]
+                n_ow1 = n[:, 1]
+                n_ow2 = n[:, 3]
+            elseif all(n[:, 3] .== n[1, 3]) && all(n[:, 4] .== n[1, 4])
+                n_w1 = n[1, 4]
+                n_w2 = n[1, 3]
+                n_ow1 = n[:, 1]
+                n_ow2 = n[:, 2]
+            end
+            w_mean = 1 / getindex.(corr_der, n_ow1)[1]
+            if 1 / getindex.(corr_der, n_ow1)[1] != 1 / getindex.(corr_der, n_ow2)[1]
+                error("w_mean with n_ow1 and with n_ow2 different")
+            end
+            der = length(a.der) >= maximum(n_ow1) && length(a.der) >= maximum(n_ow2) ? a.der[n_ow1] * w_mean : 0.0
+        end
     else
         return nothing
     end