PMS probes and Nevzorov

PMS probes and Nevzorov#

The following example presents the data collected by the PMS and Nevzorov probes on Polar 5 during the AFLUX and MOSAiC-ACA campaigns. The PMS probes where attached to the wings of the aircraft and operated by the group at DLR e.V. and JGU Mainz namely: C. Voigt, M. Moser, Yvonne Boose, and V. Hahn. If you have questions or if you would like to use the data for a publication, please get in contact with the dataset authors as stated above and in the dataset attributes contact or author. The data has been uploaded to the PANGAEA data base for both campaigns: AFLUX and MOSAiC-ACA.

Data access#

import os
from dotenv import load_dotenv

load_dotenv()

ac3cloud_username = os.environ['AC3_USER']
ac3cloud_password = os.environ['AC3_PASSWORD']

credentials = dict(user=ac3cloud_username, password=ac3cloud_password)

# local caching
kwds = {'simplecache': dict(
    cache_storage=os.environ['INTAKE_CACHE'], 
    same_names=True
)}

---------------------------------------------------------------------------
ModuleNotFoundError                       Traceback (most recent call last)
Cell In[1], line 2
      1 import os
----> 2 from dotenv import load_dotenv
      4 load_dotenv()
      6 ac3cloud_username = os.environ['AC3_USER']

ModuleNotFoundError: No module named 'dotenv'

To analyse the data they first have to be loaded by importing the (AC)³ airborne meta data catalogue. To do so the ac3airborne package has to be installed. More information on how to do that and about the catalog can be found here.

import ac3airborne
cat = ac3airborne.get_intake_catalog()

Load Polar 5 flight phase information#

Polar 5 flights are divided into segments to easily access start and end times of flight patterns. For more information have a look at the respective github repository.

At first we want to load the flight segments of (AC)³airborne

meta = ac3airborne.get_flight_segments()
flight = meta['MOSAiC-ACA']['P5'][flight_id]

Here, we want to have a look on the size distributions measured during different legs of a racetrack_pattern. In order to simplify things we can import the module flightphase from the ac3airborne.tools.

from ac3airborne.tools import flightphase

We can now select only the racetrack_pattern. This kind of pattern is available two times.

flight_query = flightphase.FlightPhaseFile(flight)
queried = flight_query.selectKind(['racetrack_pattern'])
len(queried)

/home/mech/.local/lib/python3.8/site-packages/ac3airborne/tools/flightphase.py:27: UserWarning: the segment MOSAiC-ACA_P5_RF11_rt01 contains following irregularities: first leg begins far away from starting point of the other legs and is therefore longer
  warnings.warn(str)

Plots#

%matplotlib inline
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
plt.style.use("../../mplstyle/book")

fig, (ax1, ax2, ax3, ax4) = plt.subplots(4, 1, sharex=True, gridspec_kw=dict(height_ratios=(1, 1, 1, 1)))
ax1.semilogy(ds_cdp.time, ds_cdp.N, color='k')
ax1.set_title('CDP')
ax2.semilogy(ds_cip.time, ds_cip.N, color='b')
ax2.set_title('CIP')
ax3.semilogy(ds_pip.time, ds_pip.N, color='g')
ax3.set_title('PIP')
ax4.plot(ds_nevzorov.time, ds_nevzorov.lwc, '--', color='g', label='LWC')
ax4.plot(ds_nevzorov.time, ds_nevzorov.twc, '--', color='b', label='TWC')
ax4.set_title('Nevzorov')

ax1.set_ylabel('N [$1/m^3$]')
ax2.set_ylabel('N [$1/m^3$]')
ax3.set_ylabel('N [$1/m^3$]')

ax4.set_ylabel('water content [kg/m^3]')

ax4.legend(frameon=False, loc='upper right')

ax4.xaxis.set_major_formatter(mdates.DateFormatter('%H:%M'))

plt.show()

../../_images/a258f5b2c21100a4f3ac04426b4e60b79b4442d3f15a756687497e50fa206b8a.png

from matplotlib import cm
import numpy as np
def colors(n):
    """Creates set of random colors of length n"""
    
    cmap = cm.get_cmap('gist_rainbow')
    rnd = np.random.uniform(low=0, high=1, size=n)
    cols = cmap(rnd)
    
    return cols

titles = ['ice','ocean']
fig, axs = plt.subplots(1, 2, gridspec_kw=dict())
for i, rtp in enumerate(queried):
    col_segments = colors(len(rtp['parts']))
    for j,part in enumerate(rtp['parts']):
        if 'leg' in part['name']:
            ds_sel = ds_pms_combined.sel(time=slice(part['start'],part['end']))
            axs[i].loglog(ds_sel.bin_mid, ds_sel.dNdD.mean(dim='time'), 'o',color=col_segments[j],label=part['levels'][0])
    axs[i].legend(frameon=False)
    axs[i].set_xlabel('D[$\mu m$]')
    axs[i].set_ylabel('N[$m^{-4}$]')
    axs[i].set_title('Over '+titles[i])

../../_images/43ea5a0afce18f8a6b7399effe2c01a544575ddecfa30f2bb5ff3881e16d9ce0.png

PMS probes and Nevzorov

Contents

PMS probes and Nevzorov#

Data access#

Nevzorov#

PMS#

Load Polar 5 flight phase information#

Plots#