sconstruct

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

from rsfproj import *

Flow('mod',None,
     '''
     spike n1=100 d1=0.015 n2=20 d2=0.15
     nsp=5 k1=50,83,14,65,30
     k2=7,5,19,6,1 l2=7,5,19,13,20
     mag=-1,1,1,0.3,0.3 p2=0,0,0,2,0
     ''')

Flow('mkn','mod','kirchnew adj=n v0=3')
Flow('kne','mkn','kirchnew adj=y v0=3')

Flow('mki','mod',
     '''
     spray axis=2 o=0 n=1 |
     transp plane=24 |
     halfint inv=y adj=n |
     preconstkirch vel=3 inv=y aal=n |
     window
     ''')
Flow('kir','mki',
     '''
     spray axis=2 o=0 n=1 |
     preconstkirch vel=3 aal=n |
     halfint inv=y adj=y |
     window
     ''')

Flow('mka','mod',
     '''
     spray axis=2 o=0 n=1 |
     transp plane=24 |
     halfint inv=y adj=n |
     preconstkirch vel=3 inv=y |
     window
     ''')
Flow('kaa','mka',
     '''
     spray axis=2 o=0 n=1 |
     preconstkirch vel=3 |
     halfint inv=y adj=y |
     window
     ''')

mod = ['mod','mki','mkn','mka']
mig = ['mod','kir','kne','kaa']

for case in mod+mig:
    Plot(case,
         '''
         dots gaineach=0 dots=0 seemean=1 label1= 
         strings=1 constsep=1 overlap=5.5 transp=1
         ''')

Result('amomod',mod,'TwoRows')
Result('amomig',mig,'TwoRows')

End()