AMSR2 Sea Ice Maps

Here we describe the method to include the AMSR2 sea ice maps as provided by University of Bremen. We grab them and convert the data set to netcdf every morning at 10 a.m. from the day before.

from ac3airborne.tools import get_amsr2_seaice as seaice

sic = seaice.get_amsr2_seaice('20220208')

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import numpy as np

# prepare for plotting
proj   = ccrs.NorthPolarStereo()
#extent = (-5.0, 24.0, 78.0, 83.0)

fig = plt.figure(figsize=(15, 15))
ax = plt.axes(projection=proj)

ax.gridlines()
ax.coastlines()

#ax.scatter(sic.lon.values.ravel(),sic.lat.values.ravel(),c=sic.seaice.values.ravel(), s= 0.1,cmap='Blues_r',transform=ccrs.PlateCarree())
ax.pcolormesh(sic.lon,sic.lat,sic.seaice,cmap='Blues_r',transform=ccrs.PlateCarree())
nya_lat = 78.924444
nya_lon = 11.928611

ax.plot(nya_lon, nya_lat, 'ro', transform=ccrs.PlateCarree())

[<matplotlib.lines.Line2D at 0x7f4c5edb8220>]