sconstruct

国外免费地震资料处理软件包

from rsfproj import *

tpow = {'02':2.0,'20':1.0} # doubled because there is no tpow in grey
gpow = {'02':0.5,'20':0.8}

for case in ('02','20'):
    data = 'wz.%s.H' % case
    Fetch(data,'wz')
    Result('yc'+case,data,
           '''
           dd form=native | window n1=625 | tpow tpow=%g | grey gpow=%g 
           ''' % (tpow[case],gpow[case]))

###########################################################################
End()