def METHOD_NAME(nside=1, bighp=1):
    clf()
    xlabel('Healpix x')
    ylabel('Healpix y')
    (ralo,rahi,declo,dechi) = healpix_radec_bounds(bighp, 1)
    ragrid = 5
    decgrid = 5
    rastep = 0.5
    decstep = 0.5
    ralo = floor(ralo/ragrid) *ragrid
    rahi = (1+floor(rahi/ragrid)) *ragrid
    declo = floor(declo/decgrid) *decgrid
    dechi = (1+floor(dechi/decgrid)) *decgrid
    for ra in arange(ralo,rahi,ragrid):
        xy=[]
        lastok = True
        for dec in arange(declo, dechi, decstep):
            (bhp, hpx, hpy) = radectohealpixf(ra, dec, nside)
            #print 'ra,dec', ra,dec, 'bighp', bhp, 'x,y', hpx,hpy
            if bhp != bighp:
                if lastok and len(xy):
                    plot([x for x,y in xy], [y for x,y in xy], 'r-')
                    xy = []
                lastok = False
                continue
            lastok = True
            xy.append((hpx,hpy))
        if len(xy):
            plot([x for x,y in xy], [y for x,y in xy], 'r-')
    for dec in arange(declo, dechi, decgrid):
        xy=[]
        lastok = True
        for ra in arange(ralo,rahi,rastep):
            (bhp, hpx, hpy) = radectohealpixf(ra, dec, nside)
            if bhp != bighp:
                if lastok and len(xy):
                    plot([x for x,y in xy], [y for x,y in xy], 'r-')
                    xy = []
                lastok = False
                continue
            lastok = True
            xy.append((hpx,hpy))
        if len(xy):
            plot([x for x,y in xy], [y for x,y in xy], 'r-')
    axhline(y=0, lw=0.5, ls='-', color='0.8')
    axhline(y=nside, lw=0.5, ls='-', color='0.8')
    axvline(x=0, lw=0.5, ls='-', color='0.8')
    axvline(x=nside, lw=0.5, ls='-', color='0.8')
    axis('scaled')
    axis([-0.1, nside+0.1, -0.1, nside+0.1])
    title('RA,Dec lines in healpix space, for big healpix %i' % bighp)
    savefig('hp-radec-%i.png' % bighp)