def METHOD_NAME(axs, field, field_i, omega_gi, output_dir, rank):
    """
    Plot the local ang global field DOFs local to the subdomain on the task
    with the given `rank`.
    """
    mesh = field.domain.mesh
    ax = pc.plot_wireframe(axs[0], mesh.cmesh)
    coors = field_i.get_coor()
    econn = field_i.econn
    ax = pd.plot_global_dofs(ax, coors, econn)
    ax.set_title('local DOFs on task %d' % rank)
    ax.figure.savefig(os.path.join(output_dir, 'local_dofs_%02d.png' % rank),
                      bbox_inches='tight')
    ax = pc.plot_wireframe(axs[1], mesh.cmesh)
    coors = field.get_coor()
    econn = field.get_econn('volume', omega_gi, 0)
    ax = pd.plot_global_dofs(ax, coors, econn)
    ax.set_title('global DOFs on task %d' % rank)
    ax.figure.savefig(os.path.join(output_dir, 'local_global_%02d.png' % rank),
                      bbox_inches='tight')