! ! $Header$ ! SUBROUTINE gr_v_scal(nx,x_v,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, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm IMPLICIT NONE !----------------------------------------------------------------------- ! Declararations: ! --------------- ! Arguments: ! ---------- INTEGER :: nx REAL :: x_v(ip1jm,nx),x_scal(ip1jmp1,nx) ! Local: ! ------ INTEGER :: l,ij !----------------------------------------------------------------------- DO l=1,nx DO ij=iip2,ip1jm x_scal(ij,l)= & (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) & /(airev(ij-iip1)+airev(ij)) ENDDO DO ij=1,iip1 x_scal(ij,l)=0. ENDDO DO ij=ip1jm+1,ip1jmp1 x_scal(ij,l)=0. ENDDO ENDDO RETURN END SUBROUTINE gr_v_scal