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
|
Rev | Line | |
---|
[2] | 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.