Harvey, thanks a lot!
Yes, I noticed the minor issue with the argument of bc_func, but that much I knew how to figure out. As you mentioned, the result is unaltered as the slab is initially at a constant temperature.
Pdesolve is DEFINITELY needed in Prime as numol is not quite as versatile as pdesolve (not that pdesolve is perfect, but considering...).
Anyway, thanks again.
Best,
Josué
