hi aq1q,
As I have said, their PDE is exactly like mine, their boundary condition at z = 0 is also same as mine, the only difference is their bc at z = infinity, in which they have:
u(infnity, t) = uc (a constant), and I have du/dz = 0 instaed. But let me remark that they did not claim to...