10. ソースコード

時間積分の箇所 (leap frog スキーム)

       !---------------------- 時間積分 ----------------------------------------------

       w_Zeta =  wa_Zeta + dt * min(2,it) *                                              &
       ( - w_Div_xy_xy( ( xy_Coli + xy_w(wb_Zeta) ) * xy_GradLon_w(w_Chi) / R0,       &
                         ( xy_Coli + xy_w(wb_Zeta) ) * xy_GradLat_w(w_Chi) / R0 ) / R0 &
         + w_Jacobian_w_w( w_xy( xy_Coli + xy_w(wb_Zeta ) ), w_Psi ) / R0**2 )

       w_D = delta *  ( wa_D + dt * min(2,it) *                                          & 
       ( + w_Div_xy_xy( ( xy_Coli + xy_w(wb_Zeta) ) * xy_GradLon_w(w_Psi) / R0,       &
                         ( xy_Coli + xy_w(wb_Zeta) ) * xy_GradLat_w(w_Psi) / R0 ) / R0 &
         + w_Jacobian_w_w( w_xy(  xy_Coli + xy_w(wb_Zeta) ), w_Chi ) / R0**2         &
         - w_Lapla_w( Grav*wb_H + w_E  ) / R0**2 ) ) 

       w_H =  wa_H + dt * min(2,it) * delta *                                           &
       ( - w_Div_xy_xy( xy_w(wb_H) * xy_GradLon_w(w_Chi) / R0,                        &
                        xy_w(wb_H) * xy_GradLat_w(w_Chi) / R0 ) / R0                  &
         + w_Jacobian_w_w( wb_H, w_Psi ) / R0**2 )

       !---------------------- 時間積分おわり ----------------------------------------