Code Structure

This page is part of the UV Exposure Model (UVEM) documentation.
It contains a documentation about the realization of the program code of the exposure model.
For an overview of all UVEM-related pages, see the UV Exposure Model main page.

TracNav

  • TOC "doc/tec/uvemtoc" is empty!



Code realization of exposure model

To calculate the human exposure the 3 main exposure model parameters (biologically weighted radiance, human geometry and obstructions) and an array of integration factors are used. The latter one is important for the integration over all solid angles e.g. over all directions of the upper hemisphere. This spatial integration can then easily be done by summation of all final weighted radiance values.

The calculation can be illustrated by using the polar plots of the exposure model parameters (see sections basic model and obstacles) and arrange them as an exploded view (see Figure on the right). The process of weighting the radiance of each solid angle with the human geometry and the existing obstructions can be easily envisioned by multiplying each value of the top polar plot with the corresponding values below. This is indicated in the Figure by the red asterisks and the black arrows.



For the calculation of the exposure Exweighted, the weighted radiances (see framed polar plot in the Figure above) are integrated over all solid angles (directions). For this, the area of each solid angle dΩ on a unit sphere must be calculated. These dimensionless areas have the unit [sr] and can be understood as weighting factors since not all considered solid angles from the different directions have the same size (see also schematic diagram below).

By multiplying the areas of the solid angles with the weighted radiances of each corresponding direction, the final integral can be carried out by summation of the products. The equation can be expressed as:

\begin{equation*}
Ex_{\text{weighted}} \quad=\quad \int_{\Omega} L_{\text{weighted}}(\varepsilon,\varphi)\;\; \text{d}\Omega \quad=\quad \sum_{\Omega}\Big(\:\Omega(\varepsilon,\varphi) \cdot L_{\text{weighted}}(\varepsilon,\varphi)\:\Big) \;.
\end{equation*}



Integration factors





The Figure on the right shows a schematic diagram of an exemplary segmented upper hemisphere, that visualizes the different segments. The numbers and proportions of the segments in the schematic diagram are for demonstration purposes only and do not exactly match the solid angles used in the exposure model.



Storage of obstruction information

To calculate the exposure of a human in an urban (obstructed) environment, information about the present obstructions must be derived. This is described in ''obstacles''. If the obstruction information is derived for an city area, it is stored in a 3-dimesional array. The first two dimensions represent the x and y grid points of the model domain and the third dimension contains the obstruction information of the analyzed directions. Since the obstruction information needs quite a lot of disk space, the size of the obstruction information input file is reduced by using the FORTRAN bit function ibset. This function can set/modify each of the eight values of a Byte (1 Byte = 8 Bit) to 0 or 1 or in other words to true or false. Since the obstruction information only contains true or false information (i.e. is the sky visible or not) this bit function enables the storage of 8 different directions in one byte. The analyzed obstruction information with a resolution of 10° corresponds to 360 directions (36 azimuth directions * 10 zenith directions) and can be stored with the bit function in an array with the z dimension being 45.

In the exposure model the obstruction information are then restored from this array and converted back to integer numbers (zero for obstructed directions and 1 for directions from which radiation can reach the human).

References:

  • Seckmeyer, G., Schrempf, M.Wieczorek, A., Riechelmann, S., Graw, K., Seckmeyer, S., and Zankl, M. 2013. A Novel Method to Calculate Solar UV Exposure Relevant to Vitamin D Production in Humans, Photochem. Photobiol., 89(4), 974-983, DOI: 10.1111/php.12074.
  • Schrempf, M., Thuns, N., Lange, K., and Seckmeyer, G. 2017a. Einuss der Verschattung auf die Vitamin-D-gewichtete UV-Exposition eines Menschen, Aktuelle Derm, DOI: 10.1055/s-0043-105258.
  • Schrempf, M., Thuns, N., Lange, K., and Seckmeyer, G. 2017b. Impact of Orientation on the Vitamin D Weighted Exposure of a Human in an Urban Environment, Int. J. Environ. Res. Public Health, 14(8), 920, DOI: 10.3390/ijerph14080920.
Last modified 5 years ago Last modified on Sep 30, 2019 10:59:10 AM

Attachments (4)

Download all attachments as: .zip