Doubt about plotting seismicity in cross-section

Dear gurus,

I did my cross-section to plot earthquakes using “project”, basically something like this (sorry about the variables, it is just an example that works using a color “gray” to fill in the symbols):

awk ‘{ if ($4 <= 50) print $3, $2, $4, $5*.075 }’ FS=“,” $usgs | gmtselect $Fgmtselect $fgmtselect | project -Ccat tmp1 | tr ' ' '/' -Ecat tmp2 | tr ' ' '/' -Q -W-85/85 -V | awk ‘{ print $5, $3*-1, $4 }’ | psxy -J -R -Sc -Ggray -W -O >> $ps

My doubt is about a figure I found here:

This figure has a warm-to-cold color gradient to plot seismicity with depth (besides its size is based on magnitude). I am wondering if there is a way to use a single psxy call to create this behavior or should I use several psxy calls once seismicity is filtered out by depth (e.g., awk)?

Any support is welcomed,

Best regards,


Add one more depth output column and use -Ccpt to look up colors from those depths.

Simple like that, many thanks and sorry for the noise, I needed an extra cup of coffee to keep awake yesterday in this part of the world, I did not spot that.