timemanager Subroutine

public subroutine timemanager()


                                                                        *

Handles the computation of trajectories, i.e. determines which * trajectories have to be computed at what time. * Manages dry+wet deposition routines, radioactive decay and the computation * of concentrations. * * Author: A. Stohl * * 20 May 1996 * *


Changes, Bernd C. Krueger, Feb. 2001: * Call of convmix when new windfield is read * ------------------------------------ * Changes Petra Seibert, Sept 2002 * fix wet scavenging problem * Code may not be correct for decay of deposition! * Changes Petra Seibert, Nov 2002 * call convection BEFORE new fields are read in BWD mode * Changes Caroline Forster, Feb 2005 * new interface between flexpart and convection scheme * Emanuel's latest subroutine convect43c.f is used * Changes Stefan Henne, Harald Sodemann, 2013-2014 * added netcdf output code * Changes Espen Sollum 2014 * For compatibility with MPI version, * variables uap,ucp,uzp,us,vs,ws,cbt now in module com_mod * Unified ECMWF and GFS builds * Marian Harustak, 12.5.2017 * Changes L Bakels 2022: - OpenMP parallelisation * - converting input to ETA coordinates * - spawning particles from part_ic.nc *


                                                                        *

Variables: * DEP .true. if either wet or dry deposition is switched on * decay(maxspec) [1/s] decay constant for radioactive decay * DRYDEP .true. if dry deposition is switched on * ideltas [s] modelling period * itime [s] actual temporal position of calculation * ldeltat [s] time since computation of radioact. decay of depositions* loutaver [s] averaging period for concentration calculations * loutend [s] end of averaging for concentration calculations * loutnext [s] next time at which output fields shall be centered * loutsample [s] sampling interval for averaging of concentrations * loutstart [s] start of averaging for concentration calculations * loutstep [s] time interval for which concentrations shall be * calculated * loutrestart [s] time interval for which restart files will be produced * npoint index, which starting point the trajectory has * starting positions of trajectories * nstop serves as indicator for fate of particles * in the particle loop * nstop1 serves as indicator for wind fields (see getfields) * outnum number of samples for each concentration calculation * recoutnum number of samples for each receptor calculation * recoutnumsat number of samples for each satellite calculation * prob probability of absorption at ground due to dry * deposition * WETDEP .true. if wet deposition is switched on * weight weight for each concentration sample (1/2 or 1) * *


openmp change



Loop over the whole modelling period in time steps of mintime seconds






















Various variables for testing reason of CBL scheme, by mc


Determine age class of the particle





Before the particle is moved the calculation of the scavenged mass shall only be done once after release xscav_frac1 was initialised with a negative value




Determine age class of the particle










Arguments

None