SUBROUTINE COULMOM ( PSI, F4, WGT2, D1, D2, D3, D4, D5, D6, D7, D8, WK1, WK2 )

Argument Definitions (+ indicates altered content)
REAL*8            +D1(NNI,MXNMU),         +D2(NNI,MXNMU)
REAL*8            +D3(NNI,MXNMU),         +D4(NNI,MXNMU)
REAL*8            +D5(NNI,MXNMU),         +D6(NNI,MXNMU)
REAL*8            +D7(NNI,MXNMU),         +D8(NNI,MXNMU)
REAL*8             F4(*),      PSI(*),     WGT2(*),    WK1(*)
REAL*8             WK2(*)
Description
                                                                        * 
  Copyright (C) 1996 Leif Laaksonen, Dage Sundholm                      * 
  Copyright (C) 1996-2010 Jacek Kobus               * 
                                                                        * 
  This program is free software; you can redistribute it and/or modify  * 
  it under the terms of the GNU General Public License version 2 as     * 
  published by the Free Software Foundation.                            * 
                                                                        * 
 
### coulMom ### 
 
    Calculates multipole moments up to order k_max=7.
Source file:coulMom.f
I/O Operations:
Unit ID  Unit No       Access  Form   Operation
      *                   SEQ  FMTD           W  

Operation codes A=rewind,B=backspace,C=close,E=endfile
                I=inquire,O=open,R=read,W=write
Intrinsic Functions Called
DOUBLE PRECISION   DBLE
GENERIC*8          ABS,        SQRT
External Functions and Subroutines Called
REAL*8             DOT
SUBROUTINE         PMTX,       PROD,       PROD2
Local Variables (+ indicates altered content)
INTEGER           +I,         +IBEG,      +IORB,      +MU
INTEGER           +N,         +NGRID,     +NI
REAL*8            +COSTH,     +DOME(10),  +RR,        +XR
REAL*8            +XW
Referenced Common Block Variables (+ indicates altered content)
ADDRESS1           INTEGER            I1B(60)
ADDRESS3           INTEGER            I1SI(60)
CONFIG             INTEGER            NORB
CONSFP             REAL*8             PRECIS
CONSINT            INTEGER            IONE
DEBUG              INTEGER            INCRMU,    INCRNI
DEBUG              INTEGER            IPRINT(1000)
GRIDI              INTEGER            MXNMU,     NNI
GRIDR              REAL*8             VETA(1500)
GRIDR              REAL*8             VETASQ(1500)
GRIDR              REAL*8             VXI(2500)
GRIDR              REAL*8             VXISQ(2500)
INPUTD2            REAL*8             R2
MOMENT             REAL*8            +CMULTI(1200)
SCFI               INTEGER            IHOMON
SORI               INTEGER            ITOUCH(60)
SORI               INTEGER            MPOLE