#include <misc.h>
#include <params.h>
subroutine engy_tdif(cwava ,w ,t ,tm1 ,pdel , & 2,2
difft ,nlon )
!-----------------------------------------------------------------------
!
! Purpose:
! Calculate contribution of current latitude to del-T integral
!
! Method:
!
! Author: J. Olson
!
!-----------------------------------------------------------------------
!
! $Id: engy_tdif.F90,v 1.1.2.3 2003/06/13 15:44:42 hender Exp $
! $Author: hender $
!
!-----------------------------------------------------------------------
!
use shr_kind_mod
, only: r8 => shr_kind_r8
use pmgrid
implicit none
!
!------------------------------Arguments--------------------------------
!
integer , intent(in) :: nlon ! longitude dimension
real(r8), intent(in) :: cwava ! normalization factor l/(g*plon)
real(r8), intent(in) :: w ! gaussian weight this latitude
real(r8), intent(in) :: t (plond,plev) ! temperature
real(r8), intent(in) :: tm1 (plond,plev) ! temperature (previous timestep)
real(r8), intent(in) :: pdel(plond,plev) ! pressure diff between interfaces
real(r8), intent(out) :: difft ! accumulator
!
!---------------------------Local variables-----------------------------
!
integer i,k ! longitude, level indices
real(r8) const ! temporary constant
!
!-----------------------------------------------------------------------
!
! Integration factor (the 0.5 factor arises because gaussian weights sum to 2)
!
const = cwava*w*0.5
difft = 0.
!
! Compute mass integral
!
do k=1,plev
do i=1,nlon
difft = difft + pdel(i,k)
end do
end do
difft = difft*const
return
end subroutine engy_tdif