Dear Patrik,
> When I replace the 'if continuous' section with the following section
> it will complain.
>
> if continuous:
> def lead_shape(site):
> (x, y, z) = site.pos
> return (x == -0.7 and y == 0.0 and z <= -a)
> t00=0.0
> Leadham
> ="t00*
Thank you Joe for pointing out the mistake!
I defined the lead_shape function that return True for sites in the 'lat1'.
I have already added the lat1(0, 0, -1) site to the system.
When I replace the 'if continuous' section with the following section it
will complain.
if continuous:
def