Steven Decker
Steven Decker
### What went wrong? Using `.metpy.parse_cf()` generates: `Could not find variable corresponding to the value of grid_mapping: LatLon_Projection` multiple times. ### Operating System Linux ### Version 1.2.0 ### Python Version...
Here's example 1: ```python from datetime import datetime from metpy.io import GempakGrid import metpy print(metpy.__version__) gem_file_name = '/ldmdata/gempak/model/gfs/21092206_thin.gem' gem_file = GempakGrid(gem_file_name) plot_time = datetime(2021, 9, 25, 6) ht = gem_file.gdxarray(parameter='HGHT',...
Some MetPy functions, such as gradient and advection, allow the user to supply arrays of deltas rather than a single value, but vorticity (and its relatives) does not. Vorticity is...
### Example Code ``` import matplotlib.pyplot as plt import numpy as np import netCDF4 url = ('http://thredds-test.unidata.ucar.edu/thredds/dodsC/' 'satellite/goes16/GRB16/ABI/CONUS/Channel02/current/' 'OR_ABI-L1b-RadC-M3C02_G16_s20182341637323_e20182341640096_c20182341640138.nc') ds = netCDF4.Dataset(url) values = ds.variables['Rad'][:,:] print(np.count_nonzero(values==-20.289911)) plt.imshow(values) plt.show() ``` ###...
The code: ```python import datetime as dt import numpy as np import xarray as xr from xarray.backends import NetCDF4DataStore from scipy.ndimage import generic_filter from siphon.catalog import TDSCatalog from metpy.units import...
This code ``` from datetime import datetime, timedelta from siphon.catalog import TDSCatalog yesterday = date.utcnow() - timedelta(days=1) start_time = yesterday.replace(hour=12, minute=0, second=0, microsecond=0) end_time = start_time + timedelta(hours=6) cat =...
This is a feature request identifying what I see as a "hole" in the current ecosystem. There appear to be some GitHub issues related to some of the details; this...
I have successfully compiled and tested netcdf 4.8.1, but netcdf-fortran 4.5.4 generates numerous test failures. I am using gcc/gfortran 12.1. Basic log: ``` FAIL: nf_test FAIL: tst_int64 PASS: tst_f77_v2 FAIL:...
I'm not sure if `gpnexr2` is supposed to work these days, but I get: ``` Creating process: xw for queue 152600624 GPMAP PARAMETERS: Device: xw Projection: RAD Graphics area name:...
These images are processed with the follow LDM pattern-action: ``` UNIWISC ^pnga2area Q. (U1) (.*) (.*) (.*) (.*) (........) (....) FILE -close ldmdata/gempak/sat/ANTARCTIC/4km/IR/\6/IR_\6_\7 ``` Attempting to use `gpmap` to plot...