def get_meshes_vec(step, var):
if step.geom.twod_xz:<EOL><INDENT>xmesh, ymesh = step.geom.x_mesh[:, <NUM_LIT:0>, :], step.geom.z_mesh[:, <NUM_LIT:0>, :]<EOL>vec1 = step.fields[var + '<STR_LIT:1>'][:, <NUM_LIT:0>, :, <NUM_LIT:0>]<EOL>vec2 = step.fields[var + '<STR_LIT:3>'][:, <NUM_LIT:0>, :, <NUM_LIT:0>]<EOL><DEDENT>elif step.geom.cartesian and step.geom.twod_yz:<EOL><INDENT>xmesh, ymesh = step.geom.y_mesh[<NUM_LIT:0>, :, :], step.geom.z_mesh[<NUM_LIT:0>, :, :]<EOL>vec1 = step.fields[var + '<STR_LIT:2>'][<NUM_LIT:0>, :, :, <NUM_LIT:0>]<EOL>vec2 = step.fields[var + '<STR_LIT:3>'][<NUM_LIT:0>, :, :, <NUM_LIT:0>]<EOL><DEDENT>else:  <EOL><INDENT>xmesh, ymesh = step.geom.x_mesh[<NUM_LIT:0>, :, :], step.geom.y_mesh[<NUM_LIT:0>, :, :]<EOL>pmesh = step.geom.p_mesh[<NUM_LIT:0>, :, :]<EOL>vec_phi = step.fields[var + '<STR_LIT:2>'][<NUM_LIT:0>, :, :, <NUM_LIT:0>]<EOL>vec_r = step.fields[var + '<STR_LIT:3>'][<NUM_LIT:0>, :, :, <NUM_LIT:0>]<EOL>vec1 = vec_r * np.cos(pmesh) - vec_phi * np.sin(pmesh)<EOL>vec2 = vec_phi * np.cos(pmesh) + vec_r * np.sin(pmesh)<EOL><DEDENT>return xmesh, ymesh, vec1, vec2<EOL>
Return vector field components along with coordinates meshes.

    Only works properly in 2D geometry.

    Args:
        step (:class:`~stagpy.stagyydata._Step`): a step of a StagyyData
            instance.
        var (str): vector field name.
    Returns:
        tuple of :class:`numpy.array`: xmesh, ymesh, fldx, fldy
            2D arrays containing respectively the x position, y position, x
            component and y component of the requested vector field.