And for those who toil in the CFD vineyard and are looking for a Fortran specific approach for computing flux Jacobians using AD, I can also recommend the following arXiv pre-print and the subsequent Journal of Computational Physics article.
See also JCP, Vol 399, (2019) 108942
https://www.sciencedirect.com/science/article/abs/pii/S0021999119306473