Puffin — RadiativeTransfer.jl DISORT ocean radiative transfer

Mode?

Water (uniform)? Water type / absorption model Chl (mg/m³)

Phytoplankton absorption model

a_CDOM = a_CDOM(λ_ref)·exp[−S(λ−λ_ref)]

a_CDOM(λ_ref) (m⁻¹)

slope S (nm⁻¹)

λ_ref (nm)

a_NAP = a_NAP(λ_ref)·exp[−S(λ−λ_ref)]

a_NAP(λ_ref) (m⁻¹)

slope S (nm⁻¹)

λ_ref (nm)

Scattering?

Particle phase function

asymmetry g

Particle b_b spectral slope η b_b,p(λ)=b_b,p(550)(550/λ)^η

Inelastic processes?

Chlorophyll fluorescence

fluorescence quantum yield Φ

Raman scattering

Raman coefficient b_R(488) — source

b_R(λ)=b_R(488)(488/λ)⁵ (energy). Bray et al. 2013's 1.84×10⁻⁴ is photon-to-photon at 532 nm → 2.42×10⁻⁴ energy@488 (×(532/488)^5.3×488/585).

Geometry?

Sun zenith (°)

Wind (m/s)

Perfectly flat surface (idealized, σ=0)

Unchecked, wind=0 is a calm sea that still keeps the Cox-Munk residual capillary roughness (σ²=0.003). Check for a truly flat specular interface.

Sky radiance

Bottom depth (m)

Bottom refl.

Pure-water IOPs?

Salinity (PSU)

Temp (°C)

Seawater scattering Pure-water absorption

Default: Mason (≤550 nm) → Pope & Fry (>550 nm) → Smith & Baker (NIR). Sources floored above 310 nm (Pope & Fry, Lee) block the run — grid starts at 310 nm.

Röttgers T/S absorption correction

Spectral range?

λ start (nm)

λ end (nm)

λ step (nm)

The solver runs only [λ start, λ end] (hard cut). Narrower = faster, but Raman/fluorescence near the blue edge lose their sub-λ-start excitation and are under-estimated (shaded on the plots). Raising λ start to a source's floor (Lee ≥350, Pope & Fry ≥380) also unblocks that source.

Display?

Radiance λ (nm) — Kd & polar

Overlay last N runs (R, K_d & R_rs)

Run name — CSV filename + header; blank ⇒ Puffin_<date>_<n>

aCDOM = aCDOM(λref)·exp[−S(λ−λref)]

aNAP = aNAP(λref)·exp[−S(λ−λref)]

Help

Feedback, bug & feature request

a_CDOM = a_CDOM(λ_ref)·exp[−S(λ−λ_ref)]

a_NAP = a_NAP(λ_ref)·exp[−S(λ−λ_ref)]