dynlev = 0
IF ( rC_dyn(L).GE.dptry_pedge(ntry+1) ) dynlev = L
DO L = ntry+1,MIN(Nrphys,ntry+Nr-dynlev)
pref(L) = rC_dyn(dynlev+L-ntry)
if(rF_pmid(L).ge.pref(ntry)) dynlev = L
IF ( rF_pmid(dynlev).ge.pref(ntry)-25. ) THEN
upperBnd = ntry+Nr-dynlev
upperBnd = ntry+Nr-dynlev-1
if(rF_pmid(dynlev).ge.pref(ntry)-25.) then
pref(L) = rF_pmid(dynlev+L-ntry)
pref(ntry) = rF_pmid(dynlev)
pref(L) = rF_pmid(dynlev+L-ntry+1)
WRITE(ioUnit,'(3(A,I5),A)') ' dynlev=', dynlev,
integer n,L,dynlev