In line 198, replace meshplt = plot_mesh(mesh_ob,boundary=bound,ax=axs,meshcol='C1',pub=pub,phys=phys,quiet=quiet) with meshplt = plot_mesh(mesh_ob,filename=filename,boundary=bound,ax=axs,meshcol='C1',pub=pub,phys=phys,quiet=quiet)
In line 198, replace
meshplt = plot_mesh(mesh_ob,boundary=bound,ax=axs,meshcol='C1',pub=pub,phys=phys,quiet=quiet)
with
meshplt = plot_mesh(mesh_ob,filename=filename,boundary=bound,ax=axs,meshcol='C1',pub=pub,phys=phys,quiet=quiet)