4.fortran90でコーディングする

「座標」を表すために「配列添字」を用いる




さっきの差分式をfortran90で書くと

u_a(0:nx) = u_b(0:nx) &
            - grav * ( h(1:nx+1) - h(0:nx) ) / dx * 2.0 * dt

h_a(1:nx) = h_b(1:nx) &
           - depth * ( u(1:nx) - u(0:nx-1) ) / dx * 2.0 * dt