! ! $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 comgeom_mod_h USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h IMPLICIT NONE !----------------------------------------------------------------------- ! Declararations: ! --------------- ! 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