CMIP6_export_light_results_NOAA.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "outputs": [],
   "source": [
    "import xarray as xr\n",
    "import numpy as np\n",
    "import warnings\n",
    "import matplotlib.pyplot as plt\n",
    "import os\n",
    "import seaborn as sns\n",
    "from shapely.geometry import box, mapping\n",
    "import geopandas as gpd\n",
    "import global_land_mask\n",
    "import logging\n",
    "import CMIP6_area_calculations\n",
    "\n",
    "logging.getLogger('xclim').setLevel('ERROR')\n",
    "logging.getLogger('root').setLevel('ERROR')"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "outputs": [],
   "source": [
    "def convert_to_180(ds):\n",
    " #   ds=ds.assign_coords(lat=ds.y)\n",
    "    return (ds.assign_coords(lon=(((ds.lon + 180) % 360) - 180))).sortby('lon')\n",
    "\n",
    "def convert_time(ds):\n",
    "    if not ds.indexes[\"time\"].dtype in [\"datetime64[ns]\"]:\n",
    "\n",
    "        time_objects = ds.indexes['time'].to_datetimeindex() \n",
    "        ds=ds.assign_coords({\"time\": time_objects})                   \n",
    "        ds = xr.decode_cf(ds)\n",
    "        \n",
    "    return ds\n",
    "\n",
    "def create_land_ocean_mask(cmip6_grid: xr.Dataset) -> xr.DataArray:\n",
    "    print(\"[create_land_ocean_mask] Running create_land_ocean_mask\")\n",
    "    lon = cmip6_grid.lon.values\n",
    "    lat = cmip6_grid.lat.values\n",
    "    lon_180 = xr.where(lon > 180, lon - 360, lon)\n",
    "    lon_grid, lat_grid = np.meshgrid(lon_180, lat)\n",
    "\n",
    "    mask_data = global_land_mask.globe.is_ocean(lat_grid, lon_grid).astype(int)\n",
    "    #  mask_data[np.isnan(mask_data)] = np.nan\n",
    "\n",
    "    return xr.DataArray(mask_data, coords={'lat': lat, 'lon': lon},\n",
    "                        dims=['lat', 'lon'])\n",
    "\n",
    "def resample_to_monthly(ds:xr.Dataset, var_name):\n",
    "    ds = xr.where(((ds < 1.e-20) | (ds > 1e20)), np.nan, ds)\n",
    "    ds = convert_to_180(ds)\n",
    "    ds = ds.sel(time=slice(start_time,end_time))\n",
    "\n",
    "    ds_lme = convert_time(ds)\n",
    "    # Add land mask\n",
    "    ds_lme[\"mask\"] = create_land_ocean_mask(ds_lme)\n",
    "    ds_lme = ds_lme.where(ds_lme.mask == 1)\n",
    "    ds_lme = get_data_within_LME(ds_lme, var_name, LME, False)\n",
    "    return ds_lme.resample(time=\"M\").mean()\n",
    "\n",
    "def get_area_averaged_ds(fname_perc, fname_stats,  scenario, var_name, LME):\n",
    "\n",
    "    var_name_mean = f\"{var_name}_mean\"\n",
    "\n",
    "    if os.path.exists(fname_stats) and os.path.exists(fname_perc):\n",
    "        ds_stats = xr.open_dataset(fname_stats)\n",
    "        ds_perc = xr.open_dataset(fname_perc)\n",
    "\n",
    "        ds_stats = resample_to_monthly(ds_stats, var_name_mean)\n",
    "        ds_perc = resample_to_monthly(ds_perc, var_name)\n",
    "\n",
    "        # First we calculate the total area covered by the LME. The total areacello is broadcasted along time\n",
    "        # axis and used later.\n",
    "        ds_perc[\"areacello\"] = CMIP6_area_calculations.calculate_areacello(ds_perc.sel(percentiles=2.5), var_name)\n",
    "        total_area_lme = np.nansum(ds_perc[\"areacello\"])\n",
    "\n",
    "        print(\"[calculate_areacello] Average total area {:,.2f} km2\".format(total_area_lme))\n",
    "\n",
    "        # We assess where light more than light_threshold (W/m2) is changing. For Polar cod, light_threshold is\n",
    "        # the minimum light level required to feed. If the area where minimum light increases, this can be indicative\n",
    "        # of areas opening up for feeding given that food levels are increasing too.\n",
    "        light_threshold=0.1\n",
    "        if var_name in [\"uvb_mean\"]:\n",
    "            light_threshold=0.001\n",
    "\n",
    "        ds_perc = get_data_within_LME(ds_perc, var_name, LME, False)\n",
    "        ds_stats = get_data_within_LME(ds_stats, var_name_mean, LME, False)\n",
    "\n",
    "        # Average data spatially within the LME as a function of time\n",
    "        ds = ds_perc.mean({\"lat\",\"lon\"}, skipna=True)\n",
    "        ds_s = ds_stats.mean({\"lat\",\"lon\"}, skipna=True)\n",
    "        ds_s = ds_s.assign(std = ds_stats[f\"{var_name}_stdev\"].mean(dim={\"lat\",\"lon\"}))\n",
    "\n",
    "        # Convert datasets to Pandas dataframes\n",
    "        df = ds.to_dataframe().dropna()\n",
    "        df = df.reset_index()\n",
    "\n",
    "        df_s = ds_s.to_dataframe().dropna()\n",
    "        df_s = df_s.reset_index()\n",
    "        df_s[\"model_scenario\"]=scenario\n",
    "\n",
    " # Save to CSV file\n",
    "        df_s.to_csv(f\"Light_data_stats_{LME}_{scenario}.csv\")\n",
    "        df.to_csv(f\"Light_data_perc_{LME}_{scenario}.csv\")\n",
    "        ds_stats.close()\n",
    "        ds_perc.close()\n",
    "        print(f\"Finished with CSV {LME}_{scenario}\")\n",
    "        return"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "outputs": [],
   "source": [
    "def get_LME_records():\n",
    "    lme_file='/Users/trondkr/DropBox/NIVA/oceanography/Shapefiles/LME66/LMEs66.shp'\n",
    "    return gpd.read_file(lme_file)\n",
    "\n",
    "def get_LME_records_plot():\n",
    "    lme_file='/Users/trondkr/DropBox/NIVA/oceanography/Shapefiles/LME66_180/LME66_180.shp'\n",
    "    return gpd.read_file(lme_file)\n",
    "\n",
    "def get_data_within_LME(ds,var_name,LME,create_maps):\n",
    "    print(\"Working on LME: {}\".format(LME))\n",
    "\n",
    "    # Extract the polygon defining the boundaries of the LME\n",
    "    shdf = get_LME_records()\n",
    "    #for name in shdf['LME_NAME']:\n",
    "    #    print(name)\n",
    "    shdf_sel = shdf[ shdf['LME_NAME']==LME ]\n",
    "\n",
    "    # Rioxarray requires x and y dimensions - we convert these back to lon and lat later.\n",
    "    # We also add the projection (lat-lon) so that rioxarray can do the clipping of the data according to the\n",
    "    # shapefile.\n",
    "    tos=ds.rename({'lon': 'x','lat': 'y'})\n",
    "    tos=tos.rio.write_crs(4326)\n",
    "\n",
    "    # Clip the data within the LME. We have to convert the polygon geometry to a geodataframe using\n",
    "    # `shapely.geometry`. The clipping of data within the polygon is done using rioxarray.clip function\n",
    "    clipped = tos.rio.clip(geometries=shdf_sel.geometry.apply(mapping), crs=tos.rio.crs)\n",
    "    clipped=clipped.rename({'x': 'lon','y': 'lat'}) #.to_dataset()\n",
    "\n",
    "    return clipped"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 2,996,447.57 km2\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Gulf of Alaska_ssp245\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 2,996,447.57 km2\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Gulf of Alaska\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Gulf of Alaska_ssp585\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 27,260,122.62 km2\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Northern Bering - Chukchi Seas_ssp245\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 27,260,122.62 km2\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Northern Bering - Chukchi Seas\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Northern Bering - Chukchi Seas_ssp585\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 24,185,131.43 km2\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV East Bering Sea_ssp245\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 24,185,131.43 km2\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: East Bering Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV East Bering Sea_ssp585\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 2,254,970.66 km2\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Beaufort Sea_ssp245\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "[create_land_ocean_mask] Running create_land_ocean_mask\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/trondkr/anaconda3/envs/cmip6/lib/python3.9/site-packages/iris/analysis/cartography.py:394: UserWarning: Using DEFAULT_SPHERICAL_EARTH_RADIUS.\n",
      "  warnings.warn(\"Using DEFAULT_SPHERICAL_EARTH_RADIUS.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[calculate_areacello] Average total area 2,254,970.66 km2\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Working on LME: Beaufort Sea\n",
      "Baltic Sea\n",
      "Sea of Okhotsk\n",
      "North Sea\n",
      "Gulf of Alaska\n",
      "Labrador - Newfoundland\n",
      "Celtic-Biscay Shelf\n",
      "Sea of Japan\n",
      "Oyashio Current\n",
      "Scotian Shelf\n",
      "California Current\n",
      "Black Sea\n",
      "Northeast U.S. Continental Shelf\n",
      "Mediterranean Sea\n",
      "Iberian Coastal\n",
      "Kuroshio Current\n",
      "Yellow Sea\n",
      "Canary Current\n",
      "East China Sea\n",
      "Southeast U.S. Continental Shelf\n",
      "Gulf of California\n",
      "Gulf of Mexico\n",
      "Arabian Sea\n",
      "Red Sea\n",
      "Insular Pacific-Hawaiian\n",
      "Caribbean Sea\n",
      "South China Sea\n",
      "Bay of Bengal\n",
      "Pacific Central-American Coastal\n",
      "Sulu-Celebes Sea\n",
      "Gulf of Thailand\n",
      "North Brazil Shelf\n",
      "Guinea Current\n",
      "Somali Coastal Current\n",
      "Indonesian Sea\n",
      "East Brazil Shelf\n",
      "Humboldt Current\n",
      "Benguela Current\n",
      "North Australian Shelf\n",
      "Northeast Australian Shelf\n",
      "Agulhas Current\n",
      "Northwest Australian Shelf\n",
      "East Central Australian Shelf\n",
      "South Brazil Shelf\n",
      "New Zealand Shelf\n",
      "Patagonian Shelf\n",
      "Southeast Australian Shelf\n",
      "Antarctica\n",
      "South West Australian Shelf\n",
      "West Central Australian Shelf\n",
      "Faroe Plateau\n",
      "Iceland Shelf and Sea\n",
      "Greenland Sea\n",
      "Norwegian Sea\n",
      "Barents Sea\n",
      "Kara Sea\n",
      "Laptev Sea\n",
      "East Siberian Sea\n",
      "West Bering Sea\n",
      "Aleutian Islands\n",
      "East Bering Sea\n",
      "Canadian Eastern Arctic - West Greenland\n",
      "Hudson Bay Complex\n",
      "Beaufort Sea\n",
      "Canadian High Arctic - North Greenland\n",
      "Central Arctic\n",
      "Northern Bering - Chukchi Seas\n",
      "Finished with CSV Beaufort Sea_ssp585\n"
     ]
    }
   ],
   "source": [
    "scenarios = [\"ssp245\", \"ssp585\"]\n",
    "member_range = 12\n",
    "frequency = \"A\"\n",
    "period = \"1979-01-01-2099-12-16\"\n",
    "start_time = \"1979-01-01\"\n",
    "end_time = \"2099-12-16\"\n",
    "\n",
    "models = [\"ensemble\"]\n",
    "ds_var_names = [\"par\"] #, \"uvb\", \"uv\"]\n",
    "#ds_var_names = [\"uv\",\"uvb\"]  #,\"uv\"]\n",
    "write_stats_to_file = False\n",
    "\n",
    "\n",
    "# These are the regions we extract to fulfill the ACTEA NOAA contract\n",
    "# 31.08.2022\n",
    "LMES = ['Gulf of Alaska', 'Northern Bering - Chukchi Seas', 'East Bering Sea', 'Beaufort Sea']\n",
    "\n",
    "for var_name in ds_var_names:\n",
    "    for LME in LMES:\n",
    "        for scenario in scenarios:\n",
    "\n",
    "            fname_perc = f\"/Users/trondkr/DropBox/NIVA/oceanography/light/ncfiles/ensemble/{var_name}_ensemble_perc_\" \\\n",
    "                         f\"{scenario}.nc\"\n",
    "            fname_stats = f\"/Users/trondkr/DropBox/NIVA/oceanography/light/ncfiles/ensemble/{var_name}_ensemble_stats_{scenario}.nc\"\n",
    "\n",
    "            get_area_averaged_ds(fname_perc, fname_stats, scenario, var_name, LME)\n",
    "\n"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  }
 ],
 "metadata": {
  "kernelspec": {
   "name": "python3",
   "language": "python",
   "display_name": "Python 3 (ipykernel)"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}