7. ソースコード例: 時間積分 (つづき)

リープフロッグ + 台形 semi-implicit の場合 (質量保存の式)


    ──────────────────────────────────────────
    call get_dHsfc	
    w_HsfcA =  ((1 + dt**2 * Grav * H0* n(:,1))w_HsfcB + 2 * dt * w_dHsfc) / (1 - dt**2 * Grav * H0 *n(:,1))
    ...

    end
    contains 
    soubroutine get_dHsfc

    ...
       w_dHsfc = - (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 )				&
                           - H0 * w_DivB								&
                           - H0 * 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( w_E  ) / R0**2 )