7. ソースコード例: 時間積分

リープフロッグスキームを用いた場合

    ──────────────────────────────────────────
    w_EtaA =  w_EtaB + 2 * dt * & 						! 渦度方程式
        ( - w_Div_xy_xy( xy_w(w_Eta) * xy_GradLon_w(w_Chi) / R0,			&
                                      xy_w(w_Eta) * xy_GradLat_w(w_Chi) / R0 ) / R0			&
          + w_Jacobian_w_w( w_Eta, w_Psi ) / R0**2 )

    w_DivA = w_DivB + 2 * dt * & 						! 発散方程式
        ( + w_Div_xy_xy( xy_w(w_Eta) * xy_GradLon_w(w_Psi) / R0,			&
                                       xy_w(w_Eta) * xy_GradLat_w(w_Psi) / R0 ) / R0			&
          + w_Jacobian_w_w( w_Eta, w_Chi ) / R0**2					&
          - w_Lapla_w( Grav*w_Hsfc + w_E  ) / R0**2 )  

    w_HsfcA =  w_HsfcB + 2 * dt * & 					! 質量保存の式
        ( - w_Div_xy_xy( xy_w(w_Hsfc) * xy_GradLon_w(w_Chi) / R0,			&
                                      xy_w(w_Hsfc) * xy_GradLat_w(w_Chi) / R0 ) / R0			&
          + w_Jacobian_w_w( w_Hsfc, w_Psi ) / R0**2 )