New issue 1266: ds.add_gradient_fields not working properly on 2D data https://bitbucket.org/yt_analysis/yt/issues/1266/dsadd_gradient_fields-not-w... Dominik Derigs: I want to plot the divergence of the magnetic field. For this, I think, the right approach is to use the function add_gradient_field. However, I have been unable to successfully generate any gradient field on my data (2D, FLASH). There is no suitable public dataset available that could be used to show my problem (2D, periodic, magnetic fields), so I uploaded an exemplary data file [here](https://hera.ph1.uni-koeln.de/~derigs/rotor_mhd_2d_hdf5_chk_0003_cleaning). MWE: ``` #!python import yt ds = yt.load("rotor_mhd_2d_hdf5_chk_0003_cleaning") slc = yt.SlicePlot(ds, 'z', ["pressure", "pressure_gradient_x", "pressure_gradient_y"]) slc.save() ``` **While the pressure plot is correct, the gradient plots are just empty.** pressure:  pressure_gradient_x:  pressure_gradient_y: