consistent writing and reading of BL lines params during checkpointing

src/BayesLine.c

@@ -704,12 +704,23 @@ void BayesLineLorentzSplineMCMC(struct BayesLineParams *bayesline, double heat,
     logPsx = logprior_gaussian(priors->mean, priors->sigma, Snx, 0, ncut);
   }
 
-  if(logPsx<-1)
+  if(logPsx<0)
   {
     printf("how did PSD come in outside of the prior?\n");
+    printf("%lg\n",logPsx);
     FILE *temp=fopen("failed_psd.dat","w");
     for(i=0; i<ncut; i++) fprintf(temp,"%i %lg %lg %lg\n",i,Snx[i],priors->lower[i],priors->upper[i]);
     
+    for(i=0; i<ncut; i++)
+      {
+	if(Snx[i]>priors->upper[i] || Snf[i]<priors->lower[i])
+	  {
+	    
+	    printf("The PSD is outside the prior. Offending point: %i %lg %lg %lg\n",i,Snx[i],priors->lower[i],priors->upper[i]);                                                                      
+	  }
+      }
+
+
     print_line_model(stdout, bayesline);
 
     exit(1);
@@ -1658,7 +1669,7 @@ void print_line_model(FILE *fptr, struct BayesLineParams *bayesline)
   int j;
   
   fprintf(fptr,"%i ", bayesline->lines_x->n);
-  for(j=0; j< bayesline->lines_x->n; j++) fprintf(fptr,"%e %e %e ", bayesline->lines_x->f[j], bayesline->lines_x->A[j], bayesline->lines_x->Q[j]);
+  for(j=0; j< bayesline->lines_x->n; j++) fprintf(fptr,"%lg %lg %lg ", bayesline->lines_x->f[j], bayesline->lines_x->A[j], bayesline->lines_x->Q[j]);
   fprintf(fptr,"\n");
 
 }