PROGRAM NEWPSM IMPLICIT REAL*8 (A-H,O-Z) PARAMETER(LDPAR=250) DIMENSION GC(LDPAR),THETAI(LDPAR),WPP(3) LDA=LDPAR CALL RDSIM0(1,LDPAR,NITER,NI) CALL FOPEN('Reweighted simulation matrix',2,.TRUE.) DO 30 ITER=1,NITER CALL RDSIM(1,ITER,WPP,THETAI) CALL CLIENT(NI,THETAI,WPP(2),NC,GC,PRIORC) IF(ITER.GT.1)GO TO 20 CALL PINTI("Client's functions of interest",NC,1,LDPAR) CALL WTSIM0(2,NITER,NC) 20 WPP(1)=WPP(1)+PRIORC-WPP(2) WPP(2)=PRIORC 30 CALL WTSIM(2,ITER,WPP,GC) STOP 110 CALL PEND("FIRST LINE OF INVESTIGATOR'S FILE MISSING/INCOMPLETE") 120 CALL PEND("READ ERROR, FIRST LINE OF INVESTIGATOR'S FILE") END