up [pdf]
from rsf.proj import *

# clip percentiles
perc = {'03': 95,
        '30': 85}

for case in perc.keys():
    data = 'wz.%s.H' % case
    plot = 'wz'+case

    Fetch(data,'wz')
    Result(plot,data,'grey pclip=%d' % perc[case])

End()

sfgrey

data/wz/wz.03.H
data/wz/wz.30.H