source: LMDZ.3.3/trunk/libf/dyn3d/adv_hor.F @ 3431

Last change on this file since 3431 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1      SUBROUTINE adv_hor(h,pbaru,pbarv,dh,lfiltre)
2      IMPLICIT NONE
3
4c=======================================================================
5c
6c   Auteur:  P. Le Van
7c   -------
8c
9c   ********************************************************************
10c   .... calcul du terme de convergence horizontale du flux d'enthalpie
11c        potentielle   ......
12c   ********************************************************************
13c     h,pbaru et pbarv sont des arguments d'entree  pour le s-pg ....
14c     dh               sont des arguments de sortie pour le s-pg ....
15c
16c=======================================================================
17
18
19#include "dimensions.h"
20#include "paramet.h"
21#include "logic.h"
22
23      REAL h( ip1jmp1,llm ),pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm)
24      REAL dh( ip1jmp1,llm )
25      INTEGER   l,ij
26
27      REAL hbyv( ip1jm,llm ), hbxu( ip1jmp1,llm )
28
29      EXTERNAL     convflu
30      EXTERNAL     filtreg
31      LOGICAL lfiltre
32
33
34      DO l = 1,llm
35         DO ij = iip2, ip1jm - 1
36            hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( h(ij,l) + h(ij+1,l) )
37         ENDDO
38         DO ij=1,ip1jm
39            hbyv(ij,l)=pbarv(ij,l)*0.5*(h(ij, l)+h(ij+iip1,l))
40         ENDDO
41      ENDDO
42
43c    .... correction pour  hbxu(iip1,j,l)  .....
44c    ....   hbxu(iip1,j,l)= hbxu(1,j,l) ....
45      CALL SCOPY(jjp1*llm,hbxu(iip1+iip1,1),iip1,hbxu(iip1+1,1),iip1)
46
47      CALL  convflu ( hbxu, hbyv, llm, dh )
48
49      IF (lfiltre) CALL filtreg(dh,jjp1,llm,2,2,.true.,1)
50
51      RETURN
52      END
Note: See TracBrowser for help on using the repository browser.