My Project
 All Classes Files Functions Variables Macros
tourpot.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE tourpot ( vcov, ucov, massebxy, vorpot )
5  IMPLICIT NONE
6 
7 c=======================================================================
8 c
9 c Auteur: P. Le Van
10 c -------
11 c
12 c Objet:
13 c ------
14 c
15 c *******************************************************************
16 c ......... calcul du tourbillon potentiel .........
17 c *******************************************************************
18 c
19 c vcov,ucov,fext et pbarxyfl sont des argum. d'entree pour le s-pg .
20 c vorpot est un argum.de sortie pour le s-pg .
21 c
22 c=======================================================================
23 
24 #include "dimensions.h"
25 #include "paramet.h"
26 #include "comgeom.h"
27 #include "logic.h"
28 
29  REAL rot( ip1jm,llm )
30  REAL vcov( ip1jm,llm ),ucov( ip1jmp1,llm )
31  REAL massebxy( ip1jm,llm ),vorpot( ip1jm,llm )
32 
33  INTEGER l, ij
34 
35 
36 
37 
38 c ... vorpot = ( Filtre( d(vcov)/dx - d(ucov)/dy ) + fext ) /psbarxy ..
39 
40 
41 
42 c ........ Calcul du rotationnel du vent V puis filtrage ........
43 
44  DO 5 l = 1,llm
45 
46  DO 2 ij = 1, ip1jm - 1
47  rot( ij,l ) = vcov(ij+1,l)-vcov(ij,l)+ucov(ij+iip1,l)-ucov(ij,l)
48  2 CONTINUE
49 
50 c .... correction pour rot( iip1,j,l ) .....
51 c .... rot(iip1,j,l) = rot(1,j,l) .....
52 
53 CDIR$ IVDEP
54 
55  DO 3 ij = iip1, ip1jm, iip1
56  rot( ij,l ) = rot( ij -iim, l )
57  3 CONTINUE
58 
59  5 CONTINUE
60 
61 
62  CALL filtreg( rot, jjm, llm, 2, 1, .false., 1 )
63 
64 
65  DO 10 l = 1, llm
66 
67  DO 6 ij = 1, ip1jm - 1
68  vorpot( ij,l ) = ( rot(ij,l) + fext(ij) ) / massebxy(ij,l)
69  6 CONTINUE
70 
71 c ..... correction pour vorpot( iip1,j,l) .....
72 c .... vorpot(iip1,j,l)= vorpot(1,j,l) ....
73 CDIR$ IVDEP
74  DO 8 ij = iip1, ip1jm, iip1
75  vorpot( ij,l ) = vorpot( ij -iim,l )
76  8 CONTINUE
77 
78  10 CONTINUE
79 
80  RETURN
81  END