pyrset.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyrset.f

 
C*********************************************************************
 
C...PYRSET
C...Reads a state of the random number generator from a file
C...for subsequent generation from this state onwards.
 
      SUBROUTINE pyrset(LFN,MOVE)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Commonblocks.
      common/pydatr/mrpy(6),rrpy(100)
      SAVE /pydatr/
C...Local character variable.
      CHARACTER cherr*8
 
C...Backspace required number of records (or as many as there are).
      IF(move.LT.0) THEN
        nbck=min(mrpy(6),-move)
        DO 100 ibck=1,nbck
          backspace(lfn,err=120,iostat=ierr)
  100   CONTINUE
        mrpy(6)=mrpy(6)-nbck
      ENDIF
 
C...Unformatted read from unit LFN.
      nfor=1+max(0,move)
      DO 110 ifor=1,nfor
        READ(lfn,err=120,iostat=ierr) (mrpy(i1),i1=1,5),
     &  (rrpy(i2),i2=1,100)
  110 CONTINUE
      mrpy(6)=mrpy(6)+nfor
      RETURN
 
C...Write error.
  120 WRITE(cherr,'(I8)') ierr
      CALL pyerrm(18,'(PYRSET:) error when accessing file, IOSTAT ='//
     &cherr)
 
      RETURN
      END