add functions to output dNdy and vn differential in T and tau
fix convention for the boundaries