conccalc Subroutine

public subroutine conccalc(itime, weight)

                                                                        *
 Calculation of the concentrations on a regular grid using volume       *
 sampling                                                               *
                                                                        *
 Author: A. Stohl                                                       *
                                                                        *
 24 May 1996                                                            *
                                                                        *
 April 2000: Update to calculate age spectra                            *
             Bug fix to avoid negative conc. at the domain boundaries,  *
             as suggested by Petra Seibert                              *
                                                                        *
 2 July 2002: re-order if-statements in order to optimize CPU time      *
                                                                        *
 2021, LB: OpenMP parallelisation                                       *
                                                                        *

                                                                        *

Variables: * nspeciesdim = nspec for forward runs, 1 for backward runs * *

For special runs, interpolate the air density to the particle position

AF IND_SOURCE switches between different units for concentrations at the source Af NOTE that in backward simulations the release of particles takes place Af at the receptor and the sampling at the source. Af 1="mass" Af 2="mass mixing ratio" Af IND_RECEPTOR switches between different units for concentrations at the receptor Af 1="mass" Af 2="mass mixing ratio" Af switches for the conccalcfile: AF IND_SAMP = 0 : xmass * 1 Af IND_SAMP = -1 : xmass / rho Af ind_samp is defined in readcommand.f

1. Evaluate grid concentrations using a uniform kernel of bandwidths dx, dy

For backward simulations, look from which release point the particle comes from For domain-filling trajectory option, npoint contains a consecutive particle number, not the release point information. Therefore, nrelpointer is set to 1 for the domain-filling option.

Do everything for mother domain

Do everything for the nested domain

Arguments