Temperature comparison#
Let’s see how temperature readings at BCO compare with the ERA5 reanalysis:
import intake
import healpy as hp
import matplotlib.pylab as plt
cat = intake.open_catalog("https://tcodata.mpimet.mpg.de/internal.yaml")
wxt = cat.BCO.surfacemet_wxt_v1.to_dask()
era5 = cat.HERA5(time="PT1H").to_dask()
i_bco = hp.ang2pix(2**7, wxt.lon, wxt.lat, nest=True, lonlat=True)
era5["2t"].isel(cell=i_bco).sel(time=slice("2020-01", "2020-02")).plot(label="ERA5")
(wxt["T"].sel(time=slice("2020-01", "2020-02")).resample(time="1H").mean() + 272.15).plot(label="WXT")
plt.legend();
/home/runner/miniconda3/envs/orcestra_book/lib/python3.12/site-packages/xarray/groupers.py:326: FutureWarning: 'H' is deprecated and will be removed in a future version, please use 'h' instead.
self.index_grouper = pd.Grouper(
![_images/d4b6bf657019d3f76a655021b69c142f64d7c3d71589987ced5989e8c5d38636.png](_images/d4b6bf657019d3f76a655021b69c142f64d7c3d71589987ced5989e8c5d38636.png)