connected to the flare.

Methods
=======

Input of the program is a datacube in FITS format where two coordinates are spatial and the third is time. The size in both space directions must be the same. There is also an optional free parameter and that is the multiple $\sigma_M$ of the standard deviation further described in the next paragraph. If it is omitted the program uses default value defined in the code which was set after analysing a number of different values. This process is also explained later in the next paragraphs.

The program reads input and stores it in a 3D array. For each spatial point we find and eliminate long-time trend in intensity. It is achieved by fitting the time-dependent intensity in every spatial point with five degree polynomial. Then we subtract this polynomial from the original values. After that we count the expected value and the standard deviation of intensity for each point in the space domain. The program then stores every point from the 3D array, where intensity is higher than the sum of the expected value and a multiple $\sigma_M$ of the standard deviation.

Now we find groups of points where a WLF could occur. 