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