You can do this with cgScatter2d from the Coyote Library:
data_1 = cgDemoData(1)+ RandomU(seed, 101) * 10
data_2 = cgDemoData(1)+ RandomU(seed, 101) * 20 + 50
data_3 = cgDemoData(1)
cgLoadCT, 33
cgScatter2d, data_1, data_2, COLOR=BytScl(data_3)
You can learn more about the Coyote Library at Coyote's Guide to IDL Programming:
http://www.idlcoyote.com/
Get the latest version! :-)