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