GitLab

Because boundaries for alaska, arctic ocean, and pacific WWIII area boundaries intersect with the 180th meridian, the area boundary is defined as spanning the entire range of longitudes. Since the earth is spherical and not rectangular this causes problems when checking for area intersection.

class Boundary():
    def __init__(self, south, north, west, east, fetchvar=''):

wwiii_regions = [ 
        Boundary( 15,  47,  -99,  -60, 'at_4m'),  # atlantic, also available : at_10m for 10 minute resolution
        Boundary( 15,  50, -165, -116, 'wc_4m'),  # US west, also available : wc_10m for 10 minute resolution
        Boundary( 48,  74, -180,  180, 'ak_4m'),  # alaska, also available : ak_10m for 10 minute resolution
        Boundary( 65,  90, -180,  180, 'ao_30m'), # arctic ocean
        Boundary(-20,  30, -180,  180, 'ep_10m'), # pacific
        wwiii_global
    ]