SUBROUTINE potmatn (n, rmid, w, g0, g1, g2, dmu, psi, gamma,  &
                    rhoval, xksq, energy, w0, w1, w2, wtemp,  &
                    diag, cor, asym, ndim)
!
USE FileUnits_MODULE
USE dip_MODULE
IMPLICIT NONE
!
!#include <Inc.dip>
!#include <Inc.unitno>
!
!---------------------------------------------------------------------
! written by g. a. parker
! calculates the interaction matrix
!---------------------------------------------------------------------
LOGICAL little, medium, full
INTEGER n, ithcll, ithsub, ndim, i
REAL(Kind=WP_Kind)  rmid, w(n,n), g0(n,n),g1(n,n),g2(n,n),  &
        dmu(n,n), psi(n), gamma(n), xksq(n), energy(n),     &
        w0(n,n), w1(n,n), w2(n,n),                          &
        wtemp(n,n), diag(n,n), cor(n,n), asym(n,n), rhoval(*)
SAVE

DATA ithcll/0/,ithsub/0/
DATA little/.false./,medium/.false./,full/.false./
CALL popt('potmatn ', little, medium, full, ithcll, ithsub)

CALL aphmatn(rmid, n, w, w0, w1, w2, xksq, energy, rhoval, wtemp, ndim, diag, cor, asym)

IF(lphoto.AND.rmid<=enddip)THEN
    CALL phot_mat(rmid, n, g0, g1, g2, dmu, psi, gamma, rhoval)
ENDIF
!

!         medium=.true.
IF(medium)THEN
    WRITE(Out_Unit,170) rmid
    CALL MxOutD(w, n, n, 1)
ENDIF

CALL enlvls(rmid, n, w)

RETURN
!----------------***end-potmatn***--------------------------------------
!
170 FORMAT(/,"the potential at r=",f10.5)
END