From e53c47246ddce4ae36fc44160125dd72feae21cf Mon Sep 17 00:00:00 2001
From: Jonas Krimmer <jonas.krimmer@kit.edu>
Date: Thu, 27 Mar 2025 14:05:33 +0000
Subject: [PATCH] feat: dump rytov variances

---
 examples/repeated_wave_propagation.jl | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/examples/repeated_wave_propagation.jl b/examples/repeated_wave_propagation.jl
index 57b2a97..3cb44da 100644
--- a/examples/repeated_wave_propagation.jl
+++ b/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
-- 
GitLab