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


CSE Eliminated

t_20 := scan(( varidx(T, [i, j, k])
             * varidx(t_18, [i, j, k])
             * (( 0.607824693422519
                * varidx(q, [i, j, k])
                ) + 1)
             ), k = k + 1 .. nlev, +)