subroutine getxsp(j,k) IMPLICIT NONE include 'parameter.par' include 'common.com' integer j, k, kk, jtop, jbot c------------------------------------------------------------------------------ if(j.eq.0)then do kk=k+1,ng if(jv(kk).ne.0)goto 5010 end do stop 'getxsp: bad exit at bot' 5010 jtop = jv(kk) do kk=1,nsp xsp(kk) = xs(jtop,kk) end do elseif(j.eq.n)then do kk=k-1,1,-1 if(jv(kk).ne.n)goto 5020 end do stop 'getxsp: bad exit at top' 5020 jbot = jv(kk) do kk=1,nsp xsp(kk) = xs(jbot,kk) end do else stop 'getxsp: error' endif c------------------------------------------------------------------------------ return end