[524] | 1 | ! $Header$ |
---|
[5099] | 2 | |
---|
[5106] | 3 | SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv) |
---|
[5103] | 4 | IMPLICIT NONE |
---|
| 5 | ! |
---|
| 6 | !----------------------------------------------------------------------- |
---|
| 7 | ! |
---|
| 8 | ! Auteur: P. Le Van |
---|
| 9 | ! ------- |
---|
| 10 | ! |
---|
| 11 | ! Objet: |
---|
| 12 | ! ------ |
---|
| 13 | ! calcul du terme de rotation |
---|
| 14 | ! ce terme est ajoute a d(ucov)/dt et a d(vcov)/dt .. |
---|
| 15 | ! vorpot, pbaru et pbarv sont des arguments d'entree pour le s-pg .. |
---|
| 16 | ! du et dv sont des arguments de sortie pour le s-pg .. |
---|
| 17 | ! |
---|
| 18 | !----------------------------------------------------------------------- |
---|
[524] | 19 | |
---|
[5134] | 20 | INCLUDE "dimensions.h" |
---|
| 21 | INCLUDE "paramet.h" |
---|
[524] | 22 | |
---|
[5103] | 23 | REAL :: vorpot(ip1jm, llm), pbaru(ip1jmp1, llm), & |
---|
| 24 | pbarv(ip1jm, llm), du(ip1jmp1, llm), dv(ip1jm, llm) |
---|
| 25 | INTEGER :: l, ij |
---|
| 26 | ! |
---|
| 27 | ! |
---|
| 28 | DO l = 1, llm |
---|
| 29 | ! |
---|
| 30 | DO ij = iip2, ip1jm - 1 |
---|
| 31 | du(ij, l) = 0.125 * (vorpot(ij - iip1, l) + vorpot(ij, l)) * & |
---|
| 32 | (pbarv(ij - iip1, l) + pbarv(ij - iim, l) + & |
---|
| 33 | pbarv(ij, l) + pbarv(ij + 1, l)) |
---|
| 34 | END DO |
---|
| 35 | ! |
---|
| 36 | DO ij = 1, ip1jm - 1 |
---|
| 37 | dv(ij + 1, l) = - 0.125 * (vorpot(ij, l) + vorpot(ij + 1, l)) * & |
---|
| 38 | (pbaru(ij, l) + pbaru(ij + 1, l) + & |
---|
| 39 | pbaru(ij + iip1, l) + pbaru(ij + iip2, l)) |
---|
| 40 | END DO |
---|
| 41 | ! |
---|
| 42 | ! .... correction pour dv( 1,j,l ) ..... |
---|
| 43 | ! .... dv(1,j,l)= dv(iip1,j,l) .... |
---|
| 44 | ! |
---|
| 45 | !DIR$ IVDEP |
---|
| 46 | DO ij = 1, ip1jm, iip1 |
---|
| 47 | dv(ij, l) = dv(ij + iim, l) |
---|
| 48 | END DO |
---|
| 49 | ! |
---|
| 50 | END DO |
---|
[5105] | 51 | |
---|
[5103] | 52 | END SUBROUTINE dudv1 |
---|