/*b6_l_sf.rpl*/

PROCEDURE(koc,f,Disp,theta,rho,t0,t1,c0,Length,X,n,sigma);

/*COMMENTS B6_L_SF 6/98 cge

ACTUAL BREAKTHROUGH CALCULATED FOR A PULSE APPLIED TO A LAYERED SOIL

A SCRATCH INTERNAL DATA TABLE CALLED DATA_TABLE IS USED FOR THE SOILS DATA

THE HYDRAULIC CONDUCTIVITY IS ASSUMED TO BE LOG-NORMALLY DISTRIBUTED

T0 = TIME OF CHEMICAL PULSE IN PORE VOLUMES

SIGMA = STANDARD DEVIATION OF HYDRAULIC CONCUCTIVITY

N NUNBER OF LAYERS

KD = PARTITION COEFFICIENT

THETA = WATER FILLED POROSITY

RHO IS BULK DENSITY LENGTH UNIT CM TO BE CONSISTANT WITH B6_PV

X = TIME T IN PORE VOLUMES

A DATA TABLE IS CONSTRUCTED WITH HYDRAULIC CONDUCTIVITIES BASED

ON THE MEAN AND STANDARD DEVIATION TRANSFERED TO THE PROCEDURE

THE DATA STRUCTURE FOR DATA_TABLE IS AS FOLLOWS:

v interstitial flux IS IN COL 1

ff fraction of flow is in col 2

establish a temporary array for parameters to be passed

through named param there include:

Koc = PARTITION COEFFICIENT TO organic carbon;

f = fraction organic carbon

Disp = dispersivity

THETA IS THE FLUID FILLED POROSITY;

RHO = BULK DENSITY OF SOIL

t0 = time for chemical pulse in pore volumes

t1 = time for one pore volume

c0 = concentration injected

LENGTH IS THE LENGTH OF THE COLUMN

X = ELAPSED TIME FROM BEGINNING OF CHEMICAL APPLICATION in pore volumes

*/

DELT = 1/N;

INI = DELT/2;

TOTAL = 0.;

DO I = 1 TO N;

VALUE = I*DELT - INI;

travl_t = EXP(($PROBNORM_INV(VALUE))*SIGMA);

C = B6(Koc,f,Disp,theta,rho,t0,travl_t*t1,c0,length,x);

TOTAL = TOTAL + C;

END;

RETURN TOTAL/n;

END;