It would be helpful if someone could provide an example of using pydda.constraints.get_iem_obs(). Specifically, I’m noticing that I’m getting similar results with or without specifying points=st when retrieving wind data. Could someone please provide guidance on the correct usage and whether I’m missing something?
code:
st=pydda.constraints.get_iem_obs(grid_ktlx, window=60.0)
grid_kict = pydda.initialization.make_constant_wind_field(grid_kict1, (0.0, 0.0, 0.0), vel_field=‘V’)
grids_out, _ = pydda.retrieval.get_dd_wind_field(
[grid_kict, grid_ktlx],
frz=4500.0, Co=1.0, Cm=1500.0,
max_wind_mag=100.0,points=st,
refl_field=‘dBZ’, vel_name=‘V’, filter_window=fw,filter_order=fo, engine=‘scipy’)