t_61 :: (real ~ "1" on
         (i = 1 .. nlon + 2 by
          j = 1 .. nlat + 1 by
          k = 1 .. nlev + 1
         )
        )


CSE Eliminated

t_61 := - ( 1 / 2
          * scan(( varidx(t_47, [i - 1, j, k])
                 * varidx(u, [i - 1, j, k])
                 ), k = k .. nlev, +)
          )