procedure(t0e,t_min,t_max,c0e,kde,kdmin,kdmax,fe,fmin,fmax,leng);

/*

create parameters table to perform regression analysis to function b6_l_sf conservative tracer

*/

call newrow(parameters, 'parameter','value','minimum','maximum','std dev','student t','lower','upper');

call newrow(parameters, 'kd',kde,kdmin,kdmax);

call newrow(parameters, 'f',fe,fmin,fmax);

call newrow(parameters, 'Disp',1,.2,1000);

call newrow(parameters, 'theta',.2,.2,.2);

call newrow(parameters, 'rho',1.8,1.8,1.8);

call newrow(parameters, 't0',t0e,t_min,t_max);

call newrow(parameters, 't1',1,.001,50);

call newrow(parameters, 'c0',c0e,.005,8000);

call newrow(parameters, 'length',leng,leng,leng);

call newrow(parameters, 'n',1,1,1);

call newrow(parameters, 'sigma',0,0,0);

del row 0 of table(parameters);

/*

perform regression analysis by first converting the the date time information that currently

is a number back into time format

*/

if tableexists('xy') then delete table 'xy';

if tableexists('parameter') then delete table 'parameter';

rename table(parameters) to 'parameter';

display parameter;

call sleep(10);

/*

make sure there is 'data' to regress

*/

l_row = lastrow('data');

allocate table(xy) l_row rows by 2 columns;

do ii=1 to l_row;

xy[ii,1]=1.0;

end;

set col 1 of table(xy) to col 7 of table('data');

set col 2 of table(xy) to col 6 of table ('data');

rename table(xy) to 'xy';

/*

perform inversion to convective dispersive equation

*/

lsf('xy',2,3,5000,1000,'parameter',0,11);

end;