1      subroutine tce_multipole_destroy(d_mo_mp1,l_mo_mp1_offset,
2     &                                 multipole)
3c
4c $Id$
5c
6c Written by Jeff Hammond, January 2008.
7c
8      implicit none
9#include "mafdecls.fh"
10#include "tcgmsg.fh"
11#include "global.fh"
12#include "bas.fh"
13c#include "geom.fh"
14c#include "rtdb.fh"
15c#include "sym.fh"
16#include "util.fh"
17#include "msgids.fh"
18#include "stdio.fh"
19#include "sf.fh"
20#include "inp.fh"
21#include "errquit.fh"
22c#include "tce.fh"
23c#include "tce_main.fh"
24#include "tce_ints.fh"
25c#include "tce_prop.fh"
26      integer maxaxis,multipole
27
28      maxaxis=0
29      if (multipole.eq.1) then
30        maxaxis=3
31      elseif (multipole.eq.1) then
32        maxaxis=9
33      elseif (multipole.eq.1) then
34        maxaxis=19
35      endif
36
37      do axis = maxaxis,1,-1
38        call deletefile(d_mo_mp1(axis))
39        if (.not.ma_pop_stack(l_mo_mp1_offset(axis)))
40     1    call errquit("tce_multipole_destroy: MA problem",
41     2                 l_mo_mp1_offset(axis),MA_ERR)
42      enddo
43
44      return
45      end
46