advection 2d plane: convert the compute-intensive part to Fortran90