def draw_cube_panel(self, panel, **kargs):
kwds = dict(c='0.5', lw=1., ls='-')
kwds.update(kargs)
for side, (corner1,corner2) in panel_side2corners.items():
xy1, xy2 = panel_corner_xys[corner1], panel_corner_xys[corner2]
latlon1 = xyp2latlon(xy1[0], xy1[1], panel)
latlon2 = xyp2latlon(xy2[0], xy2[1], panel)
draw_line(self.bmap, latlon1, latlon2, **kwds)
ps = PlotSphere(38, 127, figsize=(15,15)) # Korea centered
ps.draw_cube_panel(2)
#ps.draw_cube_panel(3)
ps.draw_cube_panel(6)
#latlon1 = xyp2latlon(0.5, 0.4, 2)
#latlon2 = xyp2latlon(0.4, 0.5, 6)
#latlon1 = xyp2latlon(at, 0.4, 2)
#latlon2 = xyp2latlon(0.4, at, 6)
#latlon1 = xyp2latlon(0.4, 0.3, 2)
#latlon2 = xyp2latlon(-0.4, 0.2, 3)
latlon1 = xyp2latlon(0.5, 0.4, 2)
latlon2 = xyp2latlon(0.4, 0.5, 6)
latlon3 = xyp2latlon(0.35, 0.5, 6)
latlon4 = xyp2latlon(0.5, 0.35, 2)
#latlon1 = xyp2latlon(0.5, 0.5, 2)
#latlon2 = xyp2latlon(0.35, 0.5, 6)
#latlon3 = xyp2latlon(at, at, 6)
#latlon4 = xyp2latlon(0.5, 0.35, 2)
#draw_greatcircle_line(ps.bmap, latlon1, latlon2, c='r')
#draw_greatcircle_line(ps.bmap, latlon3, latlon4, c='b')
latlons = [latlon1, latlon4, latlon3, latlon2]
poly = draw_polygon(ps.bmap, latlons)
poly.update(dict(fc='r'))
