more fix

Author: szy21

src/saturation_adjustment.jl

@@ -594,6 +594,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,
@@ -624,6 +628,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,
@@ -654,6 +662,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,
@@ -662,6 +674,7 @@ function saturation_adjustment(
     h,
     q_tot;
     maxiter::Int = 2,
+    tol = nothing,
 )
     sa_result = saturation_adjustment_fixed_iters(param_set, ph(), p, h, q_tot, maxiter)
     return (; sa_result.T, sa_result.q_liq, sa_result.q_ice)
@@ -684,6 +697,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,
@@ -692,6 +709,7 @@ function saturation_adjustment(
     θ_li,
     q_tot;
     maxiter::Int = 2,
+    tol = nothing,
 )
     sa_result =
         saturation_adjustment_fixed_iters(param_set, pθ_li(), p, θ_li, q_tot, maxiter)
@@ -715,6 +733,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,
@@ -746,6 +768,10 @@ branch-free execution on GPUs. For typical atmospheric conditions (T < 320 K),
 this achieves better than 0.1 K accuracy.
 
 For more control over solver parameters, use the full signature with explicit method type.
+
+# Returns
+- `NamedTuple` `(; T, q_liq, q_ice)` — note: `converged` is not included; fixed iterations
+  always succeed by construction.
 """
 function saturation_adjustment(
     param_set::APS,