trunk/dyn3d/convflu.f

SUBROUTINE convflu( xflu,yflu,nbniv,convfl )

  ! From LMDZ4/libf/dyn3d/convflu.F, version 1.1.1.1 2004/05/19 12:53:05

  !  P. Le Van

  ! Calcule la (convergence horiz. * aire locale) du flux ayant pour
  ! composantes xflu et yflu ,variables extensives .

  ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu.

  use dimens_m
  use paramet_m
  use comgeom

  IMPLICIT NONE

  integer, intent(in):: nbniv
  REAL, intent(in):: xflu( ip1jmp1,nbniv ), yflu( ip1jm,nbniv )
  real, intent(out):: convfl( ip1jmp1,nbniv )

  real convpn,convps
  INTEGER    l,ij
  REAL       SSUM

  !------------------------------------------------------------------

  DO  l = 1,nbniv
     DO  ij = iip2, ip1jm - 1
        convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   + &
             yflu(ij +1,l ) - yflu( ij -iim,l )
     end DO

     ! correction pour  convfl( 1,j,l)
     ! convfl(1,j,l)= convfl(iip1,j,l)

     DO  ij = iip2,ip1jm,iip1
        convfl( ij,l ) = convfl( ij + iim,l )
     end DO

     ! calcul aux pôles

     convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
     convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
     DO  ij = 1,iip1
        convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
        convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
     end DO
  end DO

END SUBROUTINE convflu
1	guez	62	SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
2	guez	3
3	guez	62	! From LMDZ4/libf/dyn3d/convflu.F, version 1.1.1.1 2004/05/19 12:53:05
4
5			! P. Le Van
6
7			! Calcule la (convergence horiz. * aire locale) du flux ayant pour
8			! composantes xflu et yflu ,variables extensives .
9
10			! nbniv est le nombre de niveaux vert. de xflu et de yflu.
11
12			use dimens_m
13			use paramet_m
14			use comgeom
15
16			IMPLICIT NONE
17
18			integer, intent(in):: nbniv
19			REAL, intent(in):: xflu( ip1jmp1,nbniv ), yflu( ip1jm,nbniv )
20			real, intent(out):: convfl( ip1jmp1,nbniv )
21
22			real convpn,convps
23			INTEGER l,ij
24			REAL SSUM
25
26			!------------------------------------------------------------------
27
28			DO l = 1,nbniv
29			DO ij = iip2, ip1jm - 1
30			convfl( ij + 1,l ) = xflu( ij,l ) - xflu( ij + 1,l ) + &
31			yflu(ij +1,l ) - yflu( ij -iim,l )
32			end DO
33
34			! correction pour convfl( 1,j,l)
35			! convfl(1,j,l)= convfl(iip1,j,l)
36
37			DO ij = iip2,ip1jm,iip1
38			convfl( ij,l ) = convfl( ij + iim,l )
39			end DO
40
41			! calcul aux pôles
42
43			convpn = SSUM( iim, yflu( 1 ,l ), 1 )
44			convps = - SSUM( iim, yflu( ip1jm-iim,l ), 1 )
45			DO ij = 1,iip1
46			convfl( ij ,l ) = convpn * aire( ij ) / apoln
47			convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
48			end DO
49			end DO
50
51			END SUBROUTINE convflu