proc readrec;
   jter=in[k]; k=k+1;
   wtlog=in[k]; k=k+1;
   pri=in[k]; k=k+1;
   pdat=in[k]; k=k+1;
   for ij (1,npar,1);
      thetai[ij]=in[k]; k=k+1;
   endfor;    
   retp("");
endp;