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

advect_u &
     = -1.0D0 / 3.0D0 * Bar_X(Bar_Y(v)) &
         * ( 4.0D0 * d_2y(u) - d_4y(u) ) 




可読性が飛躍的に向上！！！！