from pycmbs.mapping import SingleMap
from pycmbs.data import Data
import matplotlib.pyplot as plt
from pycmbs.polygon import Polygon
import numpy as np
import cartopy.crs as ccrs
plt.close('all')
poly1 = [(-105., 60.), (-168.022, 60.), (-168.022, 72.554), (-105., 72.554)] # ALA 1
P1 = Polygon(1, poly1)
P1.value=100.
poly2 = [(-66.377, -20.), (-79.729, -1.239), (-68.8, 11.439), (-50., 11.439), (-50., -20.)]
P2 = Polygon(7, poly2)
P2.value=200.
#~ AMZ 7 (20.000S, 66.377W) (1.239S, 79.729W) (11.439N, 68.800W) (11.439N, 50.000W) (20.000S, 50.000W)
#~ CAM 6 (11.439N, 68.800W) (1.239S, 79.729W) (28.566N, 118.323W) (28.566N, 90.315W)
tmp = np.ones((180, 360))
d = Data(None, None)
d.data = np.ma.array(tmp, mask=tmp!=tmp)
