(* enter the following *)
hfile="sim";  (* investigator's simulation file *)
outfile="newpsmout";  (* output file for error messages to be written *)
simout="client1"; (* client's simulation file (output) *)

(* DO NOT MODIFY ANYTHING BELOW *)

Off[General::spell];
Off[General::spell1];
LDPAR=250;
<<readrec.m; <<for.m; << writerec.m; <<client.m

OpenRead::noopen="Cannot open input file";
Check[h=OpenRead[hfile], Quit[1]];
out=OpenWrite[outfile];                               
outs=OpenWrite[simout];                               

SetStreamPosition[hfile,0];
Check[niter=Read[h, Number],
      WriteString[out, "read error in simulation input file, line 1 \n"];
      Quit[1]];
If[niter==EndOfFile,
   WriteString[out,"unexpected end of simulation input file, line 1 \n"];
   Quit[1];
  ];
Check[npar=Read[h,Number],
      WriteString[out,"read error in simulation input file, line 1 \n"];
      Quit[1]];
If[npar==EndOfFile,
   WriteString[out,"unexpected end of simulation input file, line 1 \n"];
   Quit[1];
  ];    
If[niter <= 0 || npar <=0,
   WriteString[out,"illegal first record of simulation input file "];
   WriteString[out,hfile, niter, npar];
   Quit[1];
  ];    
If[npar > LDPAR,
   WriteString[out,"input simulation file "];
   WriteString[out,hfile,npar, " parameters exceed limit of ", LDPAR];
   Quit[1];
  ];    


(* First pass *)
For[iter=1, iter <= niter, iter++, 
   readrec[h,npar]; 
   client[npar,thetai,pri];
   If[iter==1,
       If[nc<1 || nc>LDPAR, 
          WriteString[out,"Client's function of interest ", nc];
          WriteString[out," not in range 1 to ", LDPAR];
          Quit[1];
         ];
       WriteString[outs, niter, "  ", nc, "\n"];
     ];     
   wtlog=wtlog+priorc-pri;
   pri=priorc;
   writerec[];
   ];  (* for iter=1 to niter *)      



Quit[1];
Close[hfile];