I will be cutting a new release, PyDDA 1.3.0 that has three major bugfixes:
The weights from the final radar in the specified list were being set to zero by default. Now the coverage for the last radar is calculated correctly (Addressing #89 ).
The yz-section plotting functions were not plotting the contours onto the specified axis, but rather whatever was the matplotlib current axis. This is now fixed.
Smoothness constraints between the three engines are now consistent with each other (numpy engine was not dividing by grid size). I recommend setting Cx, Cy, and Cz to 1e-4 for best results.
New features:
You can now specify the input initial state as u, v, and w fields in the first input Grid to get_dd_wind_field. This is useful if you want to save the state of the retrieval after many iterations, then apply smoothing to reduce propogation of high-frequency noise, or change the cost function weights mid-retrieval.
get_dd_wind_field now returns 2 argments, the list of grids and the parameters used for the retrieval. This allows the user to see the weights that PyDDA pre-calculates for each grid.
I will be cutting a new release, PyDDA 1.3.0 that has three major bugfixes:
New features: