Writing a differentiable fluid solver in Julia