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)