source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/convflu.F @ 4327

Last change on this file since 4327 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.7 KB
Line 
1      SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
2c
3c  P. Le Van
4c
5c
6c    *******************************************************************
7c  ... calcule la (convergence horiz. * aire locale)du flux ayant pour
8c      composantes xflu et yflu ,variables extensives .  ......
9c    *******************************************************************
10c      xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
11c      convfl                est  un argument de sortie pour le s-pg .
12c
13c     njxflu  est le nombre de lignes de latitude de xflu,
14c     ( = jjm ou jjp1 )
15c     nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
16c
17      IMPLICIT NONE
18c
19#include "dimensions.h"
20#include "paramet.h"
21      REAL       xflu,yflu,convfl,convpn,convps
22      INTEGER    l,ij,nbniv
23      DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) ,
24     *         convfl( ip1jmp1,nbniv )
25c
26      EXTERNAL   SSUM
27      REAL       SSUM
28c
29c
30#include "comgeom.h"
31c
32      DO 5 l = 1,nbniv
33c
34      DO 2  ij = iip2, ip1jm - 1
35      convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   +
36     *                      yflu(ij +1,l ) - yflu( ij -iim,l )
37   2  CONTINUE
38c
39c
40
41c     ....  correction pour  convfl( 1,j,l)  ......
42c     ....   convfl(1,j,l)= convfl(iip1,j,l) ...
43c
44CDIR$ IVDEP
45      DO 3 ij = iip2,ip1jm,iip1
46      convfl( ij,l ) = convfl( ij + iim,l )
47   3  CONTINUE
48c
49c     ......  calcul aux poles  .......
50c
51      convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
52      convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
53      DO 4 ij = 1,iip1
54      convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
55      convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
56   4  CONTINUE
57c
58   5  CONTINUE
59      RETURN
60      END
Note: See TracBrowser for help on using the repository browser.