From 84f526c2ec03fcb8905f31f55c56473af2c80d5a Mon Sep 17 00:00:00 2001
From: "Fabian P. Schmidt" <kerel@mailbox.org>
Date: Sat, 14 Jan 2023 13:23:09 +0100
Subject: [PATCH] rffit: only call get_site(rsite_id) if rsite_id is non-zero

---
 rffit.c | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/rffit.c b/rffit.c
index 1cb8b6c..7d0ff6c 100644
--- a/rffit.c
+++ b/rffit.c
@@ -555,14 +555,16 @@ int main(int argc,char *argv[])
 	
 	// Loop over sites for plotting model
 	for (j=0;j<nsite;j++) {
-	  s0=get_site(site_number[j]);
-	  s1=get_site(d.p[0].rsite_id);
+	  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++) {
 	    mjd=xmin+d.mjd0+(xmax-xmin)*(float) i/(float) (NMAX-1);
 	    t=(float) (mjd-d.mjd0);
-	    if (d.p[0].rsite_id!=0) {
+	    if (d.p[0].rsite_id != 0) {
 	      velocity(orb,mjd,s1,&v1,&azi,&alt);
 	      velocity(orb,mjd,s0,&v,&azi,&alt);
 	      f=(float) ((1.0-v/C)*(1.0-v1/C)*d.ffit-d.f0);