@@ -247,7 +247,7 @@ function `sky()` (this last step requires the 3D scene as input in order to plac
247247the light sources adequately).
248248
249249``` julia
250- function create_sky (; mesh, lat = 52.0 * π/ 180.0 , DOY = 182 )
250+ function create_sky (;mesh, lat = 52.0 * π/ 180.0 , DOY = 182 )
251251 # Fraction of the day and day length
252252 fs = collect (0.1 : 0.1 : 0.9 )
253253 dec = declination (DOY)
@@ -264,14 +264,14 @@ function create_sky(; mesh, lat = 52.0*π/180.0, DOY = 182)
264264 Idir_PAR = f_dir.* Idir
265265 Idif_PAR = f_dif.* Idif
266266 # Create the dome of diffuse light
267- dome = sky (mesh;
268- Idir = 0.0 , # # No direct solar radiation
269- Idif = sum (Idir_PAR)/ 10 * DL, # # Daily Diffuse solar radiation
270- nrays_dif = 1_000_000 , # # Total number of rays for diffuse solar radiation
271- sky_model = StandardSky, # # Angular distribution of solar radiation
272- dome_method = equal_solid_angles, # Discretization of the sky dome
273- ntheta = 9 , # # Number of discretization steps in the zenith angle
274- nphi = 12 ) # # Number of discretization steps in the azimuth angle
267+ dome = sky (mesh,
268+ Idir = 0.0 , # # No direct solar radiation
269+ Idif = sum (Idir_PAR)/ 10 * DL, # # Daily Diffuse solar radiation
270+ nrays_dif = 1_000_000 , # # Total number of rays for diffuse solar radiation
271+ sky_model = StandardSky, # # Angular distribution of solar radiation
272+ dome_method = equal_solid_angles, # Discretization of the sky dome
273+ ntheta = 9 , # # Number of discretization steps in the zenith angle
274+ nphi = 12 ) # # Number of discretization steps in the azimuth angle
275275 # Add direct sources for different times of the day
276276 for I in Idir_PAR
277277 push! (dome, sky (mesh, Idir = I/ 10 * DL, nrays_dir = 100_000 , Idif = 0.0 )[1 ])
@@ -314,7 +314,7 @@ the `Material` objects (see `feed!()` above):
314314function run_raytracer! (forest; DOY = 182 )
315315 mesh = create_scene (forest)
316316 acc_mesh = accelerate (mesh, acceleration = BVH, rule = SAH {3} (5 , 10 ))
317- sources = create_sky (; mesh = acc_mesh, DOY = DOY)
317+ sources = create_sky (mesh = acc_mesh, DOY = DOY)
318318 rtobj = create_raytracer (acc_mesh, sources)
319319 trace! (rtobj)
320320 return nothing
0 commit comments