Skip to content

Skip Monte Carlo computations for g-points with high gasous optical depths#68

Open
MennoVeerman wants to merge 5 commits intomicrohh:developfrom
MennoVeerman:develop
Open

Skip Monte Carlo computations for g-points with high gasous optical depths#68
MennoVeerman wants to merge 5 commits intomicrohh:developfrom
MennoVeerman:develop

Conversation

@MennoVeerman
Copy link
Collaborator

G-points with high gaseous optical depth barely show any impact of clouds, but introduce a significant amount of noise. User can define a threshold ratio between the maximum gas and cloud optical depths (per-cell, not vertically integrated). If the ratio exceeds this threshold, use 1D fluxes instead of MC ray tracing. This threshold defaults to 0.1, in which case less than half of the g-points are run in 3D

…ntroduced tau_frac_threshold will be used in next commit) and simplify command line code
…h g-point. If this ratio is above a provided threshold value (command line argumen, with a default of 0.1), compute surface and tod fluxes and flux absorption rates based on the 1D computations
@Chiil
Copy link
Member

Chiil commented Feb 13, 2026

Would it not make more sense in terms of physics to also include the ratio of mean free path to grid spacing? Whether you need 3D also depends on that.

@MennoVeerman
Copy link
Collaborator Author

image Good idea, I agree that also makes sense and is closer to the physics. The attached image visualizes the potential for using mean free path. For each g-point I took the highest gaseous and cloud optical depth and plot the ratio between the two (f_tau) against the mean free path based on both optical depths. Dotted line shows the current grid spacing (200m), solid line the "optimal" f_tau based on visual inspection. The gas mean free path nicely meets the grid spacing at this optimal f_tau

@MennoVeerman
Copy link
Collaborator Author

MennoVeerman commented Feb 13, 2026

Instead of using f_tau, we can switch to using the ratio between the gas mean free path and the grid spacing to decide whether to do monte carlo or not. This ratio can be provided by the user and defaults to ~1?

…st mean free path and the small horizontal grid dimension to decide whether to use 3D monte carlo or simply the 1D solution
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants