Color cross-section layers by sampled data

[mdtanker]

Description of the desired feature:

It would be good to implement the option to color certain cross-section layers by sampled datasets. For example, the ice layer could be colored according to ice velocity; as in Fig 4b of Tankersley et al. 2022. See code below:

pygmt.makecpt(cmap = 'cool', series = '0/250/1', output = '../data/plotting/vel.cpt', continuous=True)

fig.plot(x=earth_profiles.Distance, y=earth_profiles[earth_profilelist[0]],
        style='b1000u', cmap='../data/plotting/vel.cpt', color=earth_profiles[ice_vel])

Are you willing to help implement and maintain this feature?

Yes