Rev | Line | |
---|
[101] | 1 | subroutine srchfd(m,y,n,x,j,d,f) |
---|
| 2 | c----------------------------------------------------------------------- |
---|
| 3 | c locate vector y(m) in an ordered table x(n) and return |
---|
| 4 | c the index locations in j(m). vectors d and f contain |
---|
| 5 | c the differences: d = x(j+1)-x(j) |
---|
| 6 | c f = (y-x(j))/d |
---|
| 7 | c----------------------------------------------------------------------- |
---|
| 8 | dimension y(*), x(*), j(*), d(*), f(*) |
---|
| 9 | c----- |
---|
| 10 | do 100 i=1,m |
---|
| 11 | ju=luf(y(i),x,n) |
---|
| 12 | jl = ju-1 |
---|
| 13 | j(i) = jl |
---|
| 14 | d(i) = x(ju)-x(jl) |
---|
| 15 | f(i) = (y(i)-x(jl))/d(i) |
---|
| 16 | 100 continue |
---|
| 17 | c |
---|
| 18 | return |
---|
| 19 | end |
---|
| 20 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|
Note: See
TracBrowser
for help on using the repository browser.