gradsで発散・収束、渦度を計算

水平二次元の発散を計算する関数「hdivg」
風の成分u,vを用いて
ga-> d hdivg(u,v)
とすると描ける。

また、収束は
ga-> d -hdivg(u,v)
のようにマイナスをつければ良い。

ついでに、渦度を計算する関数「hcurl」
同様に風の成分u,vを用いて
ga-> d hcurl(u,v)
で描ける。