source: LMDZ4/trunk/libf/dyn3dpar/tourpot_p.F @ 1940

Last change on this file since 1940 was 764, checked in by Laurent Fairhead, 18 years ago

Merge entre la version V3_conv et le HEAD
YM, JG, LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.3 KB
Line 
1      SUBROUTINE tourpot_p ( vcov, ucov, massebxy, vorpot )
2      USE parallel
3      IMPLICIT NONE
4
5c=======================================================================
6c
7c   Auteur:  P. Le Van
8c   -------
9c
10c   Objet:
11c   ------
12c
13c    *******************************************************************
14c    .........      calcul du tourbillon potentiel             .........
15c    *******************************************************************
16c
17c     vcov,ucov,fext et pbarxyfl sont des argum. d'entree pour le s-pg .
18c             vorpot            est  un argum.de sortie pour le s-pg .
19c
20c=======================================================================
21
22#include "dimensions.h"
23#include "paramet.h"
24#include "comgeom.h"
25#include "logic.h"
26
27      REAL  rot( ip1jm,llm )
28      REAL vcov( ip1jm,llm ),ucov( ip1jmp1,llm )
29      REAL massebxy( ip1jm,llm ),vorpot( ip1jm,llm )
30
31      INTEGER l, ij ,ije,ijb,jje,jjb
32
33
34      ijb=ij_begin-iip1
35      ije=ij_end
36     
37      if (pole_nord) ijb=ij_begin
38     
39     
40c  ... vorpot = ( Filtre( d(vcov)/dx - d(ucov)/dy ) + fext ) /psbarxy ..
41
42
43
44c    ........  Calcul du rotationnel du vent V  puis filtrage  ........
45c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
46      DO 5 l = 1,llm
47
48      if (pole_sud)  ije=ij_end-iip1-1
49      DO 2 ij = ijb, ije
50      rot( ij,l ) = vcov(ij+1,l)-vcov(ij,l)+ucov(ij+iip1,l)-ucov(ij,l)
51   2  CONTINUE
52
53c    ....  correction pour  rot( iip1,j,l )  .....
54c    ....     rot(iip1,j,l) = rot(1,j,l)    .....
55
56CDIR$ IVDEP
57
58      if (pole_sud)  ije=ij_end-iip1
59     
60      DO 3 ij = ijb+iip1-1, ije, iip1
61      rot( ij,l ) = rot( ij -iim, l )
62   3  CONTINUE
63
64   5  CONTINUE
65c$OMP END DO NOWAIT
66      jjb=jj_begin-1
67      jje=jj_end
68     
69      if (pole_nord) jjb=jjb+1
70      if (pole_sud)  jje=jje-1
71      CALL  filtreg_p( rot, jjb,jje,jjm, llm, 2, 1, .FALSE., 1 )
72
73c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
74      DO 10 l = 1, llm
75     
76      if (pole_sud)  ije=ij_end-iip1-1 
77     
78      DO 6 ij = ijb, ije
79      vorpot( ij,l ) = ( rot(ij,l) + fext(ij) ) / massebxy(ij,l)
80   6  CONTINUE
81
82c    ..... correction pour  vorpot( iip1,j,l)  .....
83c    ....   vorpot(iip1,j,l)= vorpot(1,j,l) ....
84CDIR$ IVDEP
85      if (pole_sud)  ije=ij_end-iip1
86      DO 8 ij = ijb+iip1-1, ije, iip1
87      vorpot( ij,l ) = vorpot( ij -iim,l )
88   8  CONTINUE
89
90  10  CONTINUE
91c$OMP END DO NOWAIT
92      RETURN
93      END
Note: See TracBrowser for help on using the repository browser.