ADD: divergence.f90