Issue with the computation of advection in Metpy

Question

I downloaded ERA5 netcdf data covering the region [30°W-30°E; 0-20°N]. This data downloaded contain horizontal wind components (u, v) and specific humidity (q). I need to calcul