|
Last change
on this file since 68 was
2,
checked in by post_it, 17 years ago
|
|
Initial import from ~/POST_IT/
|
|
File size:
367 bytes
|
| Line | |
|---|
| 1 | ;+ |
|---|
| 2 | ;; |
|---|
| 3 | ;; Calcul du rotationnel d'un champ de vecteur (u,v) [jpi,jpj] |
|---|
| 4 | ;; |
|---|
| 5 | ;; Creation : printemps 98 G. Roullet |
|---|
| 6 | ;; |
|---|
| 7 | ;- |
|---|
| 8 | FUNCTION curl2, u, v |
|---|
| 9 | @common |
|---|
| 10 | @com_eg |
|---|
| 11 | |
|---|
| 12 | zu = u*e1u*umaskr(*, *, 0) |
|---|
| 13 | zv = v*e2v*vmaskr(*, *, 0) |
|---|
| 14 | |
|---|
| 15 | zcurl = shift(zv, -1, 0)-zv+zu-shift(zu, 0, -1) |
|---|
| 16 | zcurl = zcurl*fmaskr(*, *, 0)/(e1f*e2f) |
|---|
| 17 | |
|---|
| 18 | zcurl = boundperio(zcurl, /ff1) |
|---|
| 19 | |
|---|
| 20 | return, zcurl |
|---|
| 21 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
