source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.f90 @ 5226

! $Header$

SUBROUTINE convflu( xflu,yflu,nbniv,convfl )

  !  P. Le Van


  !    *******************************************************************
  !  ... calcule la (convergence horiz. * aire locale)du flux ayant pour
  !  composantes xflu et yflu ,variables extensives .  ......
  !    *******************************************************************
  !  xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
  !  convfl                est  un argument de sortie pour le s-pg .

  ! njxflu  est le nombre de lignes de latitude de xflu,
  ! ( = jjm ou jjp1 )
  ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu .

  USE lmdz_ssum_scopy, ONLY: ssum
  USE lmdz_comgeom

USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
  USE lmdz_paramet
  IMPLICIT NONE
  !


  REAL :: xflu,yflu,convfl,convpn,convps
  INTEGER :: l,ij,nbniv
  DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) , &
        convfl( ip1jmp1,nbniv )


  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) ...

  !DIR$ IVDEP
  DO ij = iip2,ip1jm,iip1
  convfl( ij,l ) = convfl( ij + iim,l )
  END DO

  ! ......  calcul aux poles  .......

  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
  RETURN
END SUBROUTINE convflu