1czora...compute the matrix elements for the zora corrections
2c
3      subroutine zora_getv(rtdb, g_dens, g_zora, nexc)
4c
5C$Id$
6c
7      implicit none
8c
9#include "mafdecls.fh"
10#include "rtdb.fh"
11#include "stdio.fh"
12#include "global.fh"
13c
14      integer g_dens(2), g_zora(2)
15      integer rtdb
16      integer nexc
17      double precision rho_n   !integrated electron number
18      double precision tmat
19      double precision dummy(2)
20c
21      dummy(1) = 0.d0
22      dummy(2) = 0.d0
23c
24      call grid_quadv0(rtdb, g_dens, g_zora, nexc,rho_n, dummy, tmat)
25c
26      return
27      end
28