Complex Wave Real and Imaginary Parts

How do we access the real and imaginary parts of a complex wave?
For example,

function cWaveMean()
make/C/o/n=1000 cWave=cmplx(1+gnoise(1),2+gnoise(10))
print sum(cWave)/numpnts(cWave)
end

outputs only the average of the real part of cWave.
Defining new waves for the real and imaginary parts is time intensive.
An inbuilt function to access the real and imaginary parts would be much faster?

I am trying to get the sum of the imaginary part of cWave without calling wavestats.

make/C junk=cmplx(p, 10*p)
print/C sum(junk)

Variable/C cvar
cvar = sum(junk)
print cvar