Update TPM.jl (#90)

Change the order of the functions to calculate the energy flux. `mutual_shadowing!`, which changes the solar flux, should be placed between `update_flux_sun!` and `update_flux_scat_single!`.
Astroshaper · Oct 14, 2023 · 2f2c511 · 2f2c511
1 parent 08b6c92
commit 2f2c511
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/TPM.jl b/src/TPM.jl
@@ -641,10 +641,9 @@ function run_TPM!(btpm::BinaryTPM, ephem, time_begin::Real, time_end::Real, face
 
         ## Update enegey flux
         update_flux_sun!(btpm, r☉₁, r☉₂)
+        mutual_shadowing!(btpm, r☉₁, rₛ, R₂₁)  # Mutual-shadowing (eclipse)
         update_flux_scat_single!(btpm)
         update_flux_rad_single!(btpm)
-
-        mutual_shadowing!(btpm, r☉₁, rₛ, R₂₁)  # Mutual-shadowing (eclipse)
         mutual_heating!(btpm, rₛ, R₂₁)         # Mutual-heating
 
         update_thermal_force!(btpm)