================================================================= ------------------------------------------ subroutine interpol(x,y,rho_mat,rho_int) subroutine campoel(v,ex,ey) COMMONS: