| Summary: | Optimization of dry deposition velocity calculation has been of great interest. Every time, determining the value of the concentration boundary layer (CBL) thickness led to a waste of numerical calculation time, which appears as a huge time in large-scale climate models. The goal of this study is to optimize the numerical calculation time in the three-layer deposition model for smooth surfaces through the development of a MATLAB code that can parameterize the appropriate concentration boundary layer height (<i>y<sup>+</sup><sub>cbl</sub></i>) and internal integral calculation intervals for each particle diameter <i>D<sub>p</sub></i> (0.01–100 µm) and friction velocity <i>u*</i> (0.01–100 m/s). The particle concentration, as a solution to the particle flux equation, is obtained and modeled numerically by performing the left Riemann sum using MATLAB software. On the other hand, the number of subdivisions <i>N</i> of the Riemann sum was also parameterized for each <i>D<sub>p</sub></i> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msup><mi>u</mi><mo>*</mo></msup></mrow></semantics></math></inline-formula> in order to lessen the numerical calculation time. From a numerical point of view, the new parameterizations were tested by several computers; about 78% on the average of the computation time was saved when compared with the original algorithm. In other words, on average, about 1.2 s/calculation was gained, which is valuable in climate models simulations when millions of dry deposition calculations are needed.
|
|---|