added visc4 estimation - seems to work pretty well
Added some checks for PV fluxes in compare_plot