@ -556,7 +556,9 @@ int main(int argc,char *argv[])
// Loop over sites for plotting model
for (j=0;j<nsite;j++) {
s0 = get_site(site_number[j]);
if (d.p[0].rsite_id != 0) {
s1 = get_site(d.p[0].rsite_id);
}
color=j+2;
for (i=0;i<NMAX;i++) {