#include <stdio.h>
#define LG 1000
/* this program reads iteration record from posterior simulation file
   hfile to iter, wtlog, pri, pdat and g[] */
extern FILE *fp1;
int readrec(int nf, int *iter, double *wtlog, double *pri, 
    double *pdat, double *g)
{
         int i, ss;
         ss=fscanf(fp1,"%d %le %le %le", iter, wtlog, pri, pdat); 
         if (ss!=4) {
            printf("error in reading input file iteration %d \n", *iter); 
            return 1;}
         ss=0;
         for (i=0; i<=nf-1; i++) 
             ss=ss+fscanf(fp1,"%le", g+i);
         if (ss!=nf) {
            printf("error in reading input file iteration %d \n", *iter); 
            return 1;}
         else
            return 0;
         
}