Note
Go to the end to download the full example code
Simple plot of 2D dataΒΆ
This is an example of how to download and plot ceiliometer data from the SGP site over Oklahoma.
[DOWNLOADING] sgpceilC1.b1.20220114.000001.nc
[DOWNLOADING] sgpceilC1.b1.20220115.000015.nc
[DOWNLOADING] sgpceilC1.b1.20220116.000013.nc
[DOWNLOADING] sgpceilC1.b1.20220117.000010.nc
[DOWNLOADING] sgpceilC1.b1.20220118.000008.nc
If you use these data to prepare a publication, please cite:
Zhang, D., Ermold, B., & Morris, V. Ceilometer (CEIL). Atmospheric Radiation
Measurement (ARM) User Facility. https://doi.org/10.5439/1181954
/home/runner/micromamba/envs/act-docs/lib/python3.11/site-packages/dask/core.py:127: RuntimeWarning: invalid value encountered in log10
return func(*(_execute_task(a, cache) for a in args))
import os
from arm_test_data import DATASETS
import matplotlib.pyplot as plt
import act
# Place your username and token here
username = os.getenv('ARM_USERNAME')
token = os.getenv('ARM_PASSWORD')
# If the username and token are not set, use the existing sample file
if username is None or token is None or len(username) == 0 or len(token) == 0:
filename_ceil = DATASETS.fetch('sgpceilC1.b1.20190101.000000.nc')
ceil_ds = act.io.arm.read_arm_netcdf(filename_ceil, engine='netcdf4')
else:
# Example to show how easy it is to download ARM data if a username/token are set
results = act.discovery.download_arm_data(
username, token, 'sgpceilC1.b1', '2022-01-14', '2022-01-19'
)
ceil_ds = act.io.arm.read_arm_netcdf(results)
# Adjust ceilometer data for plotting
ceil_ds = act.corrections.ceil.correct_ceil(ceil_ds, -9999.0)
# Plot up ceilometer backscatter using HomeyerRainbow CVD friendly colormap
# The same could be done with keyword 'cmap='HomeyerRainbow'
display = act.plotting.TimeSeriesDisplay(ceil_ds, subplot_shape=(1,), figsize=(15, 5))
display.plot('backscatter', subplot_index=(0,), cvd_friendly=True)
plt.show()
Total running time of the script: (0 minutes 14.828 seconds)