clear;
close all;

grid.N=101;
[grid.x, grid.w, grid.A, grid.B] = FiniteDifferenceGrid (grid.N);

param.Pe = 5.0e-5;
param.Film = 14;
param.KcCa0 = 500;
param.Gk = 100;


sim = AdsorptionPackedBed (grid, param);

Nt = size(sim.t)
Nc = size(sim.c,2)