ISPACK-F90/N2PACK サンプルプログラム

線形モデル

• 拡散方程式 diffuse.f90 ∂ζ/∂t = ν(∂_xx + ∂_yy)ζ,
  

• 移流方程式 advect.f90
  ∂ζ/∂t = -(c_x∂_x + c_y∂_y)ζ,
  

• KdV方程式 kdv.f90
  ∂ζ/∂t + ζ∂_xζ + ∂_xxxζ=0,
  

• β面順圧系 b-plane.f90 ∂ζ/∂t + β∂ψ/∂x = 0,
  ζ=(∂_xx + ∂_yy)ψ,

• f面浅水系 shall_l.f90 ∂u/∂t - fv = - g ∂h/∂x,
  ∂v/∂t + fu = - g ∂h/∂y,
  ∂h/∂t + H(∂u/∂x+ ∂v/∂y) = 0.

• f面浅水系 : 地衡流調節問題
  • 速度場を初期値として与える場合 gsadj_u.f90
  • 表面変位を初期値として与える場合 gsadj_h.f90

非線形モデル

• 非回転順圧系 baro.f90
  ∂ζ/∂t + J(ψ,ζ) = (-1)^n-1ν(∂_xx + ∂_yy)ⁿζ,
  ζ=(∂_xx + ∂_yy)ψ,

• β面順圧系 baro.f90
  ∂ζ/∂t + J(ψ,ζ) + β∂ψ/∂x = ν(∂_xx + ∂_yy)ζ,
  ζ=(∂_xx + ∂_yy)ψ,

• f面浅水系 shall_n.f90 Du/Dt - fv = - g ∂h/∂x,
  Dv/Dt + fu = - g ∂h/∂y,
  ∂h/∂t + (∂u(H+h)/∂x+ ∂v(H+h)/∂y) = 0.