integer ierr, nff,nlf, indf
nff = IFNBLNK(fname)
& fname(nff:nlf), ''' for read/write access'
integer dlen, id, xtype, nat, nff,nlf, iv
mnc_f_names(indf)(1:(nlf-nff+1)) = fname(nff:nlf)
& ' in file ''', fname(nff:nlf), ''''
& fname(nff:nlf), ''''
INTEGER nff,idum
PARAMETER(nff=55)
INTEGER mj,mk,ma(nff)
real*8 mj,mk,ma(nff)
ma(nff)=mj
DO i=1,nff-1
ii=mod(21*i,nff)
DO i=1,nff
ma(i)=ma(i)-ma(1+mod(i+30,nff))
inext=mod(inext,nff)+1
inextp=mod(inextp,nff)+1