source: LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90 @ 5273

!
! $Header$
!
SUBROUTINE gr_u_scal(nx,x_u,x_scal)
  !%W%    %G%
  !=======================================================================
  !
  !   Author:    Frederic Hourdin      original: 11/11/92
  !   -------
  !
  !   Subject:
  !   ------
  !
  !   Method:
  !   --------
  !
  !   Interface:
  !   ----------
  !
  !  Input:
  !  ------
  !
  !  Output:
  !  -------
  !
  !=======================================================================
  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
IMPLICIT NONE
  !-----------------------------------------------------------------------
  !   Declararations:
  !   ---------------



  INCLUDE "comgeom.h"

  !   Arguments:
  !   ----------

  INTEGER :: nx
  REAL :: x_u(ip1jmp1,nx),x_scal(ip1jmp1,nx)

  !   Local:
  !   ------

  INTEGER :: l,ij

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

  DO l=1,nx
     DO ij=ip1jmp1,2,-1
        x_scal(ij,l)= &
              (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l)) &
              /(aireu(ij)+aireu(ij-1))
     ENDDO
  ENDDO

  CALL SCOPY(nx*jjp1,x_scal(iip1,1),iip1,x_scal(1,1),iip1)

  RETURN
END SUBROUTINE gr_u_scal