feat: dump rytov variances

examples/repeated_wave_propagation.jl

@@ -21,9 +21,9 @@ function repeated_wave_propagation(utx::AbstractArray{Complex{T}}, midx, domain:
     # estimate Fried parameters for partial propagations
     r0 = partial_r0_pw(2π / λ, domain.z, Cn²)
 
+    σr²eff = rytov_pw(2π / λ, domain.Dz, channel.Cn²)
     # compute Rytov variance per step
-    σr² = rytov_pw(2π / λ, domain.dz, channel.Cn²)
-    (σr² >= 0.1) && @warn "Insufficient number of planes, not within weak fluctuation regime."
+    σr² = rytov_pw.(2π / λ, diff(domain.z), channel.Cn²)
 
     # estimated beam waist radius at receiver
     zr = 2π / λ * w0^2 / 2
@@ -84,6 +84,8 @@ function repeated_wave_propagation(utx::AbstractArray{Complex{T}}, midx, domain:
             attrs(fid)["D_r0"] = D_r0
             attrs(fid)["r0eff"] = r0eff
             attrs(fid)["r0"] = r0
+            attrs(fid)["sigmaR2eff"] = σr²eff
+            attrs(fid)["sigmaR2"] = σr²
             attrs(fid)["Cn2"] = channel.Cn²
             attrs(fid)["l0"] = channel.l0
             attrs(fid)["L0"] = channel.L0