libf/dyn3d/convflu.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convflu.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
!
      SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
c
c  P. Le Van
c
c
c    *******************************************************************
c  ... calcule la (convergence horiz. * aire locale)du flux ayant pour
c      composantes xflu et yflu ,variables extensives .  ......
c    *******************************************************************
c      xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
c      convfl                est  un argument de sortie pour le s-pg .
c
c     njxflu  est le nombre de lignes de latitude de xflu, 
c     ( = jjm ou jjp1 )
c     nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
c
      use dimens_m
      use paramet_m
      use comgeom
      IMPLICIT NONE
c
      REAL       xflu,yflu,convfl,convpn,convps
      INTEGER    l,ij,nbniv
      DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) ,
     *         convfl( ip1jmp1,nbniv )
c
      REAL       SSUM
c
c
c
      DO 5 l = 1,nbniv
c
      DO 2  ij = iip2, ip1jm - 1
      convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   +
     *                      yflu(ij +1,l ) - yflu( ij -iim,l )
   2  CONTINUE
c
c

c     ....  correction pour  convfl( 1,j,l)  ......
c     ....   convfl(1,j,l)= convfl(iip1,j,l) ...
c
CDIR$ IVDEP
      DO 3 ij = iip2,ip1jm,iip1
      convfl( ij,l ) = convfl( ij + iim,l )
   3  CONTINUE
c
c     ......  calcul aux poles  .......
c
      convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
      convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
      DO 4 ij = 1,iip1
      convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
      convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
   4  CONTINUE
c
   5  CONTINUE
      RETURN
      END
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convflu.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3	!
4	SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
5	c
6	c P. Le Van
7	c
8	c
9	c *******************************************************************
10	c ... calcule la (convergence horiz. * aire locale)du flux ayant pour
11	c composantes xflu et yflu ,variables extensives . ......
12	c *******************************************************************
13	c xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
14	c convfl est un argument de sortie pour le s-pg .
15	c
16	c njxflu est le nombre de lignes de latitude de xflu,
17	c ( = jjm ou jjp1 )
18	c nbniv est le nombre de niveaux vert. de xflu et de yflu .
19	c
20	use dimens_m
21	use paramet_m
22	use comgeom
23	IMPLICIT NONE
24	c
25	REAL xflu,yflu,convfl,convpn,convps
26	INTEGER l,ij,nbniv
27	DIMENSION xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) ,
28	* convfl( ip1jmp1,nbniv )
29	c
30	REAL SSUM
31	c
32	c
33	c
34	DO 5 l = 1,nbniv
35	c
36	DO 2 ij = iip2, ip1jm - 1
37	convfl( ij + 1,l ) = xflu( ij,l ) - xflu( ij + 1,l ) +
38	* yflu(ij +1,l ) - yflu( ij -iim,l )
39	2 CONTINUE
40	c
41	c
42
43	c .... correction pour convfl( 1,j,l) ......
44	c .... convfl(1,j,l)= convfl(iip1,j,l) ...
45	c
46	CDIR$ IVDEP
47	DO 3 ij = iip2,ip1jm,iip1
48	convfl( ij,l ) = convfl( ij + iim,l )
49	3 CONTINUE
50	c
51	c ...... calcul aux poles .......
52	c
53	convpn = SSUM( iim, yflu( 1 ,l ), 1 )
54	convps = - SSUM( iim, yflu( ip1jm-iim,l ), 1 )
55	DO 4 ij = 1,iip1
56	convfl( ij ,l ) = convpn * aire( ij ) / apoln
57	convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
58	4 CONTINUE
59	c
60	5 CONTINUE
61	RETURN
62	END