#include <misc.h>#include <params.h>subroutine flxint (w       ,flx     ,flxlat  ,nlon    )!----------------------------------------------------------------------- ! ! Purpose: Calculate contribution of current latitude to energy flux integral! ! Method: ! ! Author: Jerry Olson! !-----------------------------------------------------------------------   use precision   use pmgrid!-----------------------------------------------------------------------   implicit none!-----------------------------------------------------------------------!! Arguments!   real(r8), intent(in) :: w          ! gaussian weight this latitude   real(r8), intent(in) :: flx(plond) ! energy field   integer, intent(in) :: nlon        ! number of longitudes   real(r8), intent(out) :: flxlat    ! accumulator for given latitude!! Local variables!   integer :: i                       ! longitude index!!-----------------------------------------------------------------------!   flxlat = 0.!   do i=1,nlon      flxlat = flxlat + flx(i)   end do!! Integration factor (the 0.5 factor arises because gaussian weights! sum to 2)!   flxlat = flxlat*w*0.5/float(nlon)!   returnend subroutine flxint