source: PSPA/parmelaPSPA/trunk/srchfd.f @ 418

Last change on this file since 418 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 682 bytes
Line 
1      subroutine srchfd(m,y,n,x,j,d,f)
2c-----------------------------------------------------------------------
3c     locate vector y(m) in an ordered table x(n) and return
4c     the index locations in j(m). vectors d and f contain
5c     the differences: d = x(j+1)-x(j)
6c     f = (y-x(j))/d
7c-----------------------------------------------------------------------
8      dimension y(*), x(*), j(*), d(*), f(*)
9c-----
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
17c     
18      return
19        end
20c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.