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	SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
2
3	! 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