do while (rms.gt.rmscutfrag(1,j,1) .and.
& ishiff.lt.n_shift(1,j,1))
ishiff=ishiff+1
- rms=rmscalc_frag(-ishiff,1,j,jcon,lprn)
+ rms=rmscalc_frag(-ishiff,1,j,jcon,ipermmin,lprn)
c write(iout,*)"jcon,i,j,ishiff",jcon,i,j,-ishiff,
c & " rms",rms," rmscut",rmscutfrag(1,j,1)
if (lprn) write (iout,*) "rms",rmsfrag(j,1)
if (rms.gt.rmscutfrag(1,j,1)) then
- rms=rmscalc_frag(ishiff,1,j,jcon,lprn)
+ rms=rmscalc_frag(ishiff,1,j,jcon,ipermmin,lprn)
c write (iout,*) "jcon,1,j,ishiff",jcon,1,j,ishiff,
c & " rms",rms
endif
c print *,"jcon,i,j,ishif",jcon,i,j,-ishif," rms",rms
if (lprn) write (iout,*) "rms",rmsfrag(j,i)
if (rms.gt.rmscutfrag(1,j,i)) then
- rms=rmscalc_frag(ishif,i,j,jcon,lprn)
+ rms=rmscalc_frag(ishif,i,j,jcon,ipermmin,lprn)
c print *,"jcon,i,j,ishif",jcon,i,j,ishif," rms",rms
endif
if (lprn) write (iout,*) "rms",rms