SUBROUTINE WETDEPO ( ITIME, LTSAMPLE, LOUTNEXT )

Argument Definitions (+ indicates altered content)
INTEGER            ITIME,      LOUTNEXT,   LTSAMPLE
Description
                         i      i        i 
 
 
Calculation of wet deposition using the concept of scavenging coeffici 
For lack of detailed information, washout and rainout are jointly trea 
It is assumed that precipitation does not occur uniformly within the w 
grid cell, but that only a fraction of the grid cell experiences rainf 
This fraction is parameterized from total cloud cover and rates of lar 
scale and convective precipitation. 
 
   Author: A. Stohl 
 
   1 December 1996 
 
Correction by Petra Seibert, Sept 2002: 
use centred precipitation data for integration 
Code may not be correct for decay of deposition! 
 
 
 
Variables: 
cc [0-1]           total cloud cover 
convp [mm/h]       convective precipitation rate 
fraction [0-1]     fraction of grid, for which precipitation occurs 
ix,jy              indices of output grid cell for each particle 
itime [s]          actual simulation time [s] 
jpart              particle index 
ldeltat [s]        interval since radioactive decay was computed 
lfr, cfr           area fraction covered by precipitation for large sc 
                   and convective precipitation (dependent on prec. ra 
loutnext [s]       time for which gridded deposition is next output 
loutstep [s]       interval at which gridded deposition is output 
lsp [mm/h]         large scale precipitation rate 
ltsample [s]       interval over which mass is deposited 
prec [mm/h]        precipitation rate in subgrid, where precipitation 
wetdeposit         mass that is wet deposited 
wetgrid            accumulated deposited mass on output grid 
wetscav            scavenging coefficient 
 
Constants:
Source file:wetdepo.f
Intrinsic Functions Called
INTEGER            INT,        NINT
REAL               SNGL
GENERIC            ABS,        EXP,        MAX
External Functions and Subroutines Called
SUBROUTINE         INTERPOL_RAIN,          INTERPOL_RAIN_NESTS
SUBROUTINE         WETDEPOKERNEL,          WETDEPOKERNEL_NEST
Parameter Variables Used
INTEGER            NUMPATH            (NUMPATH = 4)
INTEGER            IDIFFNORM          (IDIFFNORM = 10800)
INTEGER            NXMAX              (NXMAX = 361)
INTEGER            NYMAX              (NYMAX = 181)
INTEGER            NUVZMAX            (NUVZMAX = 61)
INTEGER            NWZMAX             (NWZMAX = 61)
INTEGER            NZMAX              (NZMAX = 61)
INTEGER            NCONVLEVMAX        (NCONVLEVMAX = NUVZMAX-1)
INTEGER            MAXNESTS           (MAXNESTS = 1)
INTEGER            NXMAXN             (NXMAXN = 1)
INTEGER            NYMAXN             (NYMAXN = 1)
INTEGER            JPACK              (JPACK = 4*NXMAX*NYMAX)
INTEGER            MAXXGRID           (MAXXGRID = 180)
INTEGER            MAXYGRID           (MAXYGRID = 90)
INTEGER            MAXXGRIDN          (MAXXGRIDN = 1)
INTEGER            MAXYGRIDN          (MAXYGRIDN = 1)
INTEGER            MAXAGECLASS        (MAXAGECLASS = 1)
INTEGER            MAXZGRID           (MAXZGRID = 14)
INTEGER            NCLASSUNC          (NCLASSUNC = 10)
INTEGER            MAXRECEPTOR        (MAXRECEPTOR = 200)
INTEGER            MAXPART            (MAXPART = 1000000)
INTEGER            MAXPOINT           (MAXPOINT = 10000)
INTEGER            MAXSPEC            (MAXSPEC = 1)
INTEGER            MAXPOINTSPEC       (MAXPOINTSPEC = MAXSPEC)
INTEGER            MAXWF              (MAXWF = 50000)
INTEGER            MAXTABLE           (MAXTABLE = 1000)
INTEGER            NUMCLASS           (NUMCLASS = 9)
INTEGER            NI                 (NI = 11)
INTEGER            MAXCOLUMN          (MAXCOLUMN = 3000)
INTEGER            MAXRAND            (MAXRAND = 2000000)
REAL               PI                 (PI = 3.14159265)
Local Variables (+ indicates altered content)
INTEGER           +I,         +ITAGE,     +IX,        +J
INTEGER           +JPART,     +JY,        +K,         +LDELTAT
INTEGER           +NAGE,      +NGRID
REAL               CC,        +CFR(5),     CONVP,     +FRACTION
REAL              +LFR(5),     LSP,       +PREC,      +RESTMASS
REAL              +SMALLNUM,  +WETDEPOSIT(MAXSPEC),   +WETSCAV
REAL              +XTN,       +YTN
Referenced Common Block Variables (+ indicates altered content)
GLOBAL4            INTEGER            NUMBNESTS
GLOBAL4            INTEGER            NXN(MAXNESTS)
GLOBAL4            INTEGER            NYN(MAXNESTS)
GLOBAL4            REAL               CONVPRECN(0:NXMAXN-1,0:NYMAXN-1,1,2,MAXNESTS)
GLOBAL4            REAL               LSPRECN(0:NXMAXN-1,0:NYMAXN-1,1,2,MAXNESTS)
GLOBAL4            REAL               TCCN(0:NXMAXN-1,0:NYMAXN-1,1,2,MAXNESTS)
GLOBAL4            REAL               XLN(MAXNESTS)
GLOBAL4            REAL               XRESOLN(0:MAXNESTS)
GLOBAL4            REAL               XRN(MAXNESTS)
GLOBAL4            REAL               YLN(MAXNESTS)
GLOBAL4            REAL               YRESOLN(0:MAXNESTS)
GLOBAL4            REAL               YRN(MAXNESTS)
GLOBALI            INTEGER            ITRA1(MAXPART)
GLOBALI            INTEGER            ITRAMEM(MAXPART)
GLOBALI            INTEGER            LAGE(MAXAGECLASS)
GLOBALI            INTEGER            LDIRECT,   LOUTSTEP
GLOBALI            INTEGER            MEMIND(2)
GLOBALI            INTEGER            MEMTIME(2)
GLOBALI            INTEGER            NAGECLASS
GLOBALI            INTEGER            NCLASS(MAXPART)
GLOBALI            INTEGER            NESTED_OUTPUT
GLOBALI            INTEGER            NSPEC,     NUMPART
GLOBALI            INTEGER            NX,        NY
GLOBALR            DOUBLE PRECISION   XTRA1(MAXPART)
GLOBALR            DOUBLE PRECISION   YTRA1(MAXPART)
GLOBALR            REAL               CONVPREC(0:NXMAX-1,0:NYMAX-1,1,2)
GLOBALR            REAL               DECAY(MAXSPEC)
GLOBALR            REAL               LSPREC(0:NXMAX-1,0:NYMAX-1,1,2)
GLOBALR            REAL               TCC(0:NXMAX-1,0:NYMAX-1,1,2)
GLOBALR            REAL               WETA(MAXSPEC)
GLOBALR            REAL               WETB(MAXSPEC)
GLOBALR            REAL              +XMASS1(MAXPART,MAXSPEC)