Calc_Div

input: Element *E, double *ux, double *uy, double *uz, double *res, double *wk

output:

Description: This function calculates the divergence of flux.