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



CSE Eliminated

(t_103 :=
 ( varidx(t_102, [i, j, k])
 * ( varidx(T, [i, j, k]) * varidx(q, [i, j, k])
   + varidx(T, [i - 1, j, k]) * varidx(q, [i - 1, j, k])
   )
 ))