{ "cells": [ { "cell_type": "markdown", "id": "a317474d-36fa-4051-92df-1c8915a9fea6", "metadata": {}, "source": [ "## Dust Concentration in Miocene Simulations --- Composite Aerosols patterns\n", "#### 1. Detemining Dust indices\n", "#### 2. Defining High/Low Dust months\n", "#### 3. Examining Composite Aerosols (Black Carbon & Sulfates) Patterns based on the definition\n", "====================================================================================================================" ] }, { "cell_type": "code", "execution_count": 1, "id": "4da43274-62e4-4a7e-b6b6-98d426966714", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import matplotlib.pyplot as plt\n", "from matplotlib import ticker\n", "import matplotlib.ticker as mticker\n", "import matplotlib.patches as patches\n", "import numpy as np\n", "\n", "import cartopy.crs as ccrs\n", "import cartopy.mpl.ticker as cticker\n", "from cartopy.util import add_cyclic_point\n", "\n", "import warnings\n", "warnings.simplefilter(action='ignore', category=FutureWarning)\n", "\n", "def label_latlon(ax,lons,lats):\n", " \"\"\" Add tick labels \"\"\"\n", " # Define the xticks for longtitude\n", " ax.set_xticks(lons,crs=ccrs.PlateCarree())\n", " lon_formatter=cticker.LongitudeFormatter()\n", " ax.xaxis.set_major_formatter(lon_formatter)\n", "\n", " # Define ytick for latitude\n", " ax.set_yticks(lats,crs=ccrs.PlateCarree())\n", " lat_formatter=cticker.LatitudeFormatter()\n", " ax.yaxis.set_major_formatter(lat_formatter)\n", "\n", " return" ] }, { "cell_type": "code", "execution_count": 2, "id": "8fdc1306-873f-498c-af14-9e17b904c2fc", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (time: 1200, lat: 96, lon: 144, nbnd: 2)\n", "Coordinates:\n", " * lat (lat) float64 -90.0 -88.11 -86.21 -84.32 ... 86.21 88.11 90.0\n", " * lon (lon) float64 0.0 2.5 5.0 7.5 10.0 ... 350.0 352.5 355.0 357.5\n", " * time (time) object 1900-02-01 00:00:00 ... 2000-01-01 00:00:00\n", "Dimensions without coordinates: nbnd\n", "Data variables:\n", " BURDENDUST (time, lat, lon) float32 dask.array<chunksize=(1200, 96, 144), meta=np.ndarray>\n", " time_bnds (time, nbnd) object dask.array<chunksize=(1200, 2), meta=np.ndarray>\n", "Attributes:\n", " Conventions: CF-1.0\n", " source: CAM\n", " case: B.MMIOx2_C5_280_WISOon\n", " title: UNSET\n", " logname: pacosta\n", " host: r8i0n31\n", " Version: $Name$\n", " revision_Id: $Id$\n", " initial_file: /glade/scratch/pacosta/B.MMIOx2_C5_280_WISOon/run/B.MMI...\n", " topography_file: /glade/work/pacosta/PaleoBC/heroldn/gx1vMIO/cami_bnd_to...\n", " history: Sun Nov 12 20:18:00 2023: ncrcat ./B.MMIOx2_C5_280_WISO...\n", " NCO: netCDF Operators version 5.1.4 (Homepage = http://nco.s...