SHORT DESCRIPTION OF THE SOURCE CODE FILES USED BY THE FLEXPART MODEL
_____________________________________________________________________
FLEXPART.f        main program, manages reading of run specifications, calls
advance.f         advection of the trajectory using a first order scheme for one
assignland.f      assigns the fractions of 8 landuse classes to each ECMWF grid
                  point
calcfluxes.f      calculates the gross N-S, E-W and up-down mass fluxes caused
                  by the motion of a single particle
calcpar.f         calculates boundary layer parameters: friction velocity,
                  convective velocity scale, boundary layer height, Obukhov length
calcpar_nests.f   same as calcpar.f, but for the nested grids
caldate.f         compute calendar date from julian date
cmapf1.0.f        Program package for coordinate transformations between
                  latitude/longitude grid and polar stereographic projection
coordtrafo.f      transformation of release point coordinates from geografical
                  to grid coordinates
dist2.f           calculation of distance (m) between two points given in 
                  geografical coordinates
drydepokernel.f   calculates gridded dry deposition field using a uniform kernel
                  of bandwidths dx, dy
erf.f             error function
ew.f              calculation of saturation vapor pressure at given temperature
                  synchronisation interval; uses mean wind and Langevin equations
                  for parameterizing turbulence and mesoscale wind fluctuations
fixparticles.f    determine heights, release times and masses of all particles
                  timemanager to do the actual calculations
fluxout.f         writes mass flux data to output file and reinitializes fields
getfields.f       management of wind fields in memory
getrb.f           computation of quasilaminar sublayer resistance to dry
                  deposition of gases
getrc.f           calculation of surface resistance to dry deposition of gases
getvdep.f         calculation of dry deposition velocities
gridcheck.f       determines grid specifications (grid area, etc.) from GRIB file
gridcheck_nests.f same as gridcheck.f, but for the nested grids
gridconc.f        calculation of concentrations on a grid using volume sampling
                  and at receptor points using kernel method
gridout.f         output of concentrations on grid and for receptor points
hanna.f           parameterization of the turbulent velocity variance and
                  the respective timescales; Langevin equation for w'/sigw
hanna1.f          alternative model to hanna.f: Langevin equation for w'
                  drift correction velocity after Wilson
hanna_short.f     same as hanna.f, but only for the vertical wind
includecom        include file containing all global variables
includeinterpol   include file containing variables used for interpolation
includehanna      include file used for turbulence parameterization
includepar        include file containing all globally used parameter statements
initialize.f      initialization of the turbulence statistics for each particle
interpol_all.f    interpolation of all data needed by advance for two levels within
                  the boundary layer and for the surface
interpol_all_nests.f  same as interpol_all.f, but for the nested grids
interpol_misslev.f interpolates all levels not yet interpolated by interpol_all.f
interpol_misslev_nests.f same as interpol_misslev.f, but for the nested grids
interpol_rain.f   interpolation of large-scale and convective precipitation and
                  total cloud cover
interpol_rain_nests.f same as interpol_rain.f, but for the nested grids
interpol_vdep.f   interpolation of the deposition velocity
interpol_vdep_nests.f  same as interpol_vdep.f, but for the nested grids
interpol_wind.f   interpolation of wind data above the PBL
interpol_wind_nests.f   same as interpol_wind.f, but for the nested grids
juldate.f         compute julian date from calendar date
mean.f            calculation of mean and standard deviation
obukhov.f         calculation of Obukhov length from surface heat flux
openoutput.f      opens an output file and writes a header; is used for several
                  output files, depending on how it is called
part0.f           calculation of time independent factors for the dry deposition
                  of particles
partdep.f         calculation of dry deposition velocities for particles
partout.f         dumps the particle positions
pbl_profile.f     calculation of surface stress and sensible heat flux using
                  the profile method
psih.f            stability correction term for heat
psim.f            stability correction term for momentum
qvsat.f           saturation specif. humidity at given pressure and temperature
raerod.f          calculation of the aerodynamical resistance to dry deposition
random.f          subroutines for generation of normally distriubted
                  random numbers
readageclasses.f  reads number and time limits of age classes to be used
readavailable.f   read, which wind fields are available within the modelling period
readcommand.f     read user commands, e.g. start of release, duration of release
readdepo.f        read parameters needed for dry deposition of gases
readlanduse.f     read landuse inventory and the respective roughness lenghts
readoro.f         read ecmwf model orography and land sea mask
readoutgrid.f     read output grid specifications
readpartpositions.f reads particle end positions from a previous run's dump file as
                  particle starting positions for the new run
readpaths.f       read pathnames, where input/output files are stored
readreceptors.f   read names and coordinates of receptor points; for receptor
                  points, concentrations are calculated with the kernel method 
readreleases.f    read release point specifications (e.g., coordinates, type of
                  trajectories to be used, number of particles, mass, etc.)
readspecies.f     read physical and chemical properties of species
readwind.f        reads the ECMWF meteorological data fields; GRIB decoding
readwind_nests.f  same as readwind.f, but for the nested grids
richardson.f      calculation of convective velocity scale and boundary layer
                  height using critical Richardson number concept
scalev.f          computation of friction velocity from surface stress
sigma.f           calculation of standard deviation of a 1-D field
skplin.f          skips n lines of a given input file
timemanager.f     time management of FLEXPART; determines when a integration time
                  step of a trajectory is due; determines when a concentration
                  calculation is due, etc, and calls the respective subprograms
verttransform.f   transformation of three-dimensional fields from eta coordinates
                  to meter above ground coordinates
verttransform_nests.f  same as verttransform.f, but for the nested grids
wetdepo.f         calculation of wet deposition
wetdepokernel.f   calculates gridded wet deposition field using a uniform kernel
                  of bandwidths dx, dy
windalign.f       transformation of random velocities from along- and cross-wind
                  components to Cartesian u and v components
