Compute optimal fields for the instantaneous production of enstrophy in 3D flows