source: LMDZ4/trunk/libf/dyn3dpar/dteta1_p.F @ 4020

Last change on this file since 4020 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.2 KB
RevLine 
[630]1      SUBROUTINE dteta1_p ( teta, pbaru, pbarv, dteta)
2      USE parallel
3      USE write_field_p
4      IMPLICIT NONE
5
6c=======================================================================
7c
8c   Auteur:  P. Le Van
9c   -------
10c Modif F.Forget 03/94 (on retire q et dq  pour construire dteta1)
11c
12c   ********************************************************************
13c   ... calcul du terme de convergence horizontale du flux d'enthalpie
14c        potentielle   ......
15c   ********************************************************************
16c  .. teta,pbaru et pbarv sont des arguments d'entree  pour le s-pg ....
17c     dteta               sont des arguments de sortie pour le s-pg ....
18c
19c=======================================================================
20
21
22#include "dimensions.h"
23#include "paramet.h"
24#include "logic.h"
25
26      REAL teta( ip1jmp1,llm ),pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm)
27      REAL dteta( ip1jmp1,llm )
28      INTEGER   l,ij
29
30      REAL hbyv( ip1jm,llm ), hbxu( ip1jmp1,llm )
31
32c
33      INTEGER ijb,ije,jjb,jje
34
35     
36      jjb=jj_begin
37      jje=jj_end
[764]38
39c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
[630]40      DO 5 l = 1,llm
41     
42      ijb=ij_begin
43      ije=ij_end
44     
45      if (pole_nord) ijb=ij_begin+iip1
46      if (pole_sud)  ije=ij_end-iip1
47     
48      DO 1  ij = ijb, ije - 1
49        hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( teta(ij,l)+teta(ij+1,l) )
50   1  CONTINUE
51
52c    .... correction pour  hbxu(iip1,j,l)  .....
53c    ....   hbxu(iip1,j,l)= hbxu(1,j,l) ....
54
55CDIR$ IVDEP
56      DO 2 ij = ijb+iip1-1, ije, iip1
57        hbxu( ij, l ) = hbxu( ij - iim, l )
58   2  CONTINUE
59
60      ijb=ij_begin-iip1
61      if (pole_nord) ijb=ij_begin
62     
63      DO 3 ij = ijb,ije
64        hbyv(ij,l)= pbarv(ij, l)* 0.5 * ( teta(ij, l)+teta(ij+iip1,l) )
65   3  CONTINUE
66
67       if (.not. pole_sud) then
68          hbxu(ije+1:ije+iip1,l) = 0
69          hbyv(ije+1:ije+iip1,l) = 0
70        endif
71       
72   5  CONTINUE
[764]73c$OMP END DO NOWAIT
[630]74       
75       
76        CALL  convflu_p ( hbxu, hbyv, llm, dteta )
77
78
79c    stockage dans  dh de la convergence horizont. filtree' du  flux
80c                  ....                           ...........
81c           d'enthalpie potentielle .
82     
83     
84      CALL filtreg_p( dteta,jjb,jje,jjp1, llm, 2, 2, .true., 1)
85     
86     
87      RETURN
88      END
Note: See TracBrowser for help on using the repository browser.