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


CSE Eliminated

t_13 := log(varidx(p, [i, j, k])) * varidx(p, [i, j, k])