readreleases Subroutine

public subroutine readreleases()


                                                                        *
 This routine reads the release point specifications for the current    *
 model run. Several release points can be used at the same time.        *
                                                                        *
 Author: A. Stohl                                                       *
                                                                        *
 18 May 1996                                                            *
                                                                        *
 Update: 29 January 2001                                                *
 Release altitude can be either in magl or masl                         *
 HSO, 12 August 2013
 Added optional namelist input
                                                                        *

                                                                        *

Variables: * decay decay constant of species * dquer [um] mean particle diameters * dsigma e.g. dsigma=10 or dsigma=0.1 means that 68% of the mass are between 0.1dquer and 10dquer * ireleasestart, ireleaseend [s] starting time and ending time of each * release * kindz 1: zpoint is in m agl, 2: zpoint is in m asl, 3: zpoint is in hPa * npart number of particles to be released * nspec number of species to be released * density [kg/m3] density of the particles * rm [s/m] Mesophyll resistance * species name of species * xmass total mass of each species * xpoint1,ypoint1 geograf. coordinates of lower left corner of release * area * xpoint2,ypoint2 geograf. coordinates of upper right corner of release * area * weta_gas, wetb_gas parameters for below-cloud scavenging (gas) * crain_aero, csnow_aero parameters for below-cloud scavenging (aerosol) * ccn_aero, in_aero parameters for in-cloud scavenging (aerosol) * zpoint1,zpoint2 height range, over which release takes place * num_min_discrete if less, release cannot be randomized and happens at * time mid-point of release interval * lroot true if serial version, or if MPI and root process * *




















Arguments

None

Namelists

Namelist releases_ctrl


Variables

Name Type Default Description
nspec integer None
specnum_rel integer None

Namelist release


Variables

Name Type Default Description
idate1 integer None
itime1 integer None
idate2 integer None
itime2 integer None
lon1 real None
lon2 real None
lat1 real None
lat2 real None
z1 real None
z2 real None
zkind integer(kind=2) None
mass real None
parts integer None
comment character(len=40) None