1      subroutine qhop_dsite(ga,rsa,rind)
2c
3c $Id$
4c
5      implicit none
6c
7#include "qhop_common.fh"
8#include "mafdecls.fh"
9c
10      integer ga,rind,rsa
11c
12      call qhop_dsite2(ga,rsa,rind,
13     + int_mb(i_ptseq),int_mb(i_ptarat),int_mb(i_iarat),
14     + dbl_mb(i_racs),dbl_mb(i_racs+mxar),int_mb(i_ptpar),
15     + dbl_mb(i_par),dbl_mb(i_par+7*maxpar),dbl_mb(i_par+14*maxpar),
16     + dbl_mb(i_par+17*maxpar),dbl_mb(i_par+23*maxpar),
17     + dbl_mb(i_par+28*maxpar))
18c
19      return
20      end
21      subroutine qhop_dsite2(ga,rsa,rind,
22     + ptseq,ptarat,arat,deq,aneq,ptpar,tunnel,defe12,zpef,
23     + tdsgl,tstval,e12fxy)
24c
25c $Id$
26c
27      implicit none
28c
29#include "qhop_common.fh"
30c
31      integer ga,rind,rsa
32c
33      real*8 tunnel(maxpar,7),defe12(maxpar,7),zpef(maxpar,3),
34     + tdsgl(maxpar,6),tstval(maxpar,5),e12fxy(maxpar,3)
35      real*8 deq(mxar),aneq(mxar)
36      integer ptseq(maxseq,2),ptpar(maxpar),arat(mxar,4),ptarat(mxseq)
37c
38      rind=arat(ptarat(ptseq(rsa,1))+ga-
39     + ptseq(rsa,2),1)
40c
41      return
42      end
43