Calculate DOHC, based on GFDL basins, split into N/S and E/W
Add annotations to keep track of basin keys