In the context of the study of the integrated Sachs Wolfe effect (ISW), we revisit the angular cross correlation of WMAP cosmic microwave background (CMB) data with the NRAO Very Large Array Sky Survey (NVSS). We compute 2-point cross functions between the two surveys, both in real and in Fourier space, paying particular attention on the dependence of results on the flux of NVSS radio sources, the angular scales where correlations arise and the comparison with theoretical expectations. We reproduce previous results that claim an excess of correlation in the angular correlation function (ACF), and we also find some (low significance, ∼2-σ) similarity between the CMB and radio galaxy data in the multipole range l ∈ [10, 25]. However, the signal to noise ratio (S /N) in the ACFs increases with higher flux thresholds for NVSS sources, but drops a ∼30−50% in separations of the order of a pixel size, suggesting some residual point source contribution. When restricting our analyses to multipoles l < 60, we fail to find any evidence for cross correlation in the range l ∈ [2, 10], where according to the model predictions and our simulations ∼50% of the S /N is supposed to arise. Also, the accumulated S /N for l < 60 is below 1, far from the theoretical expectation of S /N ∼ 5. Part of this disagreement may be caused by an inaccurate modeling of the NVSS source population: as in previous works, we find a level of large scale (l < 70) clustering in the NVSS catalog that seems incompatible with a high redshift population. This large scale clustering excess is unlikely to be caused by contaminants or systematics, since it is independent of flux threshold, and hence present for the brightest, most clearly detected (>30σ) NVSS sources. Either our NVSS catalogs are not probing the high redshift, large scale gravitational potential wells, or there is a clear mismatch between the ISW component present in WMAP data and theoretical expectations.