note that boundary conditions were created with wrong offset
for velocities due to line below in mk_run_template1.m

    if ( fld{1}=='U' | fld{1}=='U' )
        bgn=2;
    end

bug 1: the second "U" should be a "V"

bug 2: this should only be applied to normal, not trnsversal velocities