移流項の計算(gms使用前)


v_BarX_BarY(0:nx,1:ny,1:nz) = .....

advect_u(0:nx,1:ny,1:nz) &
     =  -1.0 / 3.0 * v_BarX_BarY(0:nx,1:ny,1:nz) &
         * ( 4.0 &
            * (   &
                 u(0:nx,2:ny+1,1:nz,1)&
               - u(0:nx,0:ny-1,1:nz,1)  &
               ) / 2.0 / dy &
            - (  &
                 u(0:nx,3:ny+2,1:nz,1) &
               - u(0:nx,-1:ny-2,1:nz,1)  &
               ) / 4.0 / dy &
            )