from rsf.proj import *

 ####################################### 
 # Sigsbee 2A velocity model construct # 
 #######################################

PRFX = 'sigsbee2a_'
SUFX = '.sgy'

for c in ('migration_velocity','stratigraphy','reflection_coefficients'):
    if(c=='migration_velocity'):
        v='vmig2A'
        o=3.055
        d=.01143
        s=.0003048
        l='Migration\ Velocity'
        a='y'
        u='km/s'
    if(c=='stratigraphy'):
        v='vstr2A'
        o=3.048
        d=.00762
        s=.0003048
        l='Stratigraphic\ Velocity'
        a='y'
        u='km/s'
    if(c=='reflection_coefficients'):
        v='reflectionCoefficients'
        o=3.048
        d=.00762
        s=1
        l='Reflection\ Coefficients'
        a='n'
        u=''
    h = c + 'head'
    t = PRFX + c + SUFX
    Fetch(t,'sigsbee')

    Flow([v,h],t,
         '''
         segyread tfile=${TARGETS[1]} |
         put
         o1=0  d1=.00762   label1=Depth    unit1=km
         o2=%f d2=%f       label2=Distance unit2=km
         | scale rscale=%f
         ''' %(o,d,s))

    Result(v,v,'''window f1=1 | grey color=i scalebar=y allpos=%s screenratio=.375 screenht=3 barreverse=y
    wheretitle=t title=%s labelsz=4 titlesz=6 barlabel=Velocity barunit=%s''' % (a,l,u))

End()