goad 1.1.9

Compute the single scattering properties of particles much larger than the wavelength of light with geometric optics and aperture diffraction theory.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47

"""Run GOAD backscattering simulation on the roughened hexagonal plate."""

from pathlib import Path

from goad import Convergence, Param, Settings

output_dir = Path(__file__).parent

# Path to the generated geometry
geometry_file = (
    output_dir
    / "roughened_edge15.0_sigma0p1_merge1p0_hexagonal_column_l5.0_r25.0_6904c8.obj"
)

# Set up simulation settings
# - wavelength 0.532 microns (green laser)
# - ice refractive index 1.31
# - empty zones for faster computation (just forward and backward scattering)
# - high max_tir for backscattering
settings = Settings(
    geom_path=str(geometry_file),
    wavelength=0.532,
    particle_refr_index_re=1.31,
    particle_refr_index_im=0.0,
    zones=[],
    max_tir=20,
    directory=str(output_dir),
    beam_area_threshold_fac=0.01,
    beam_power_threshold=0.01,
    cutoff=0.999,
)

# Set up convergence solver
convergence = Convergence(settings)

# Converge on lidar ratio and depolarisation ratio within 10% error
convergence.add_target(Param.LidarRatio, 0.1)
convergence.add_target(Param.DepolarizationRatio, 0.1)

# Run the simulation
convergence.solve()

# Save results
convergence.save()

print("Simulation complete!")
print(f"Results saved to: {output_dir}")