*define model Example-1 # model name *info Ray tracing in model {model} from Rayinvr with Receiver Functions ###################################################################### *setup noapp units km km ### some symbolic definitions *define vcont 6.0 # velocity level contoured in the model fname Example-1.v.in.save # file name to save the model in when "Save" is pressed ## Generate synthetic Receiver Functions for three locations *call genrate 50 10000 -1000 gauss 1.0 0.0 200 hlist x real 30 120 220 hlist p real 0.1 ### generate synthetic RF *call genrf rf-rays Example-1 ray 1.0 p r x p 1 # 1 1 *call print rf.print samples 0 1000 traces 0 100 ### Map the RF into (Z,X) along the S-wave ray *call depcon1 s-p x p bott Example-1 *call print test.print samples 0 1000 traces 0 100 *call table * x p timstrt _tnrec_ lasttr ### plot the synthetic Receiver Functions along S-wave rays *call plot wgva+- up median 1.0 0.15 name rf-section 0 0 gain 5.0e2 setamp peak -10.0 10.0 offset =comp= 1 compxyz 2 0 trange -10000 30000.0 *end # remove all traces ### Rayinvr-like model: *call rayinvr gui ps Example-1 *incl Example-1.igm column vp 100 source 1 Sp1 0 forw source 2 Sp2 300 back raycode full 30 60 60 1 1.1 3u summary # print brief model summary edit run ### Define various graphical elements *call graphic backgr =same= foregr =same= palette pal-vp1 -0.3 red 0 blue 0 white 0.3 red #palette pal-vp buor 5.5 8 palette pal-vp gmtcopper 1 2 # 4.5 8.2 palhue pal-vs color cont 2 5. 4 0.05 .95 palette pal-raycodes gmtjet 0 3 palette pal-arrivals gmtjet 0 4 palette pal-ampl gmtseis -1 1 #palproc edit loop pal-vp # add a reversed copy of itself to the palette #palproc edit repeat 2 pal-vp # repeat the palette palproc edit scale pal-vp # scale the Vp palette by Vp in the model vp Example-1 #palproc edit hbal pal-vp # histogram-balance the Vp palette #vp accrete line solid 1 red l-layers line solid 1 blue l-columns line solid 1 cyan l-contour line solid 1 blue cells-vert line solid 1 yellow cells-diag line solid 1 emerald l-rays1 line solid 1 plum l-rays2 line solid 1 brown l-tau line solid 1 melon line-grcalc # line for calculated gravity anomaly line solid 1 blue ttimes-forw line solid 1 red ttimes-rev line solid 1 yellow ttimes-observed line solid 1 yellow wiggle-line # line for the wiggles fill solid green velocity fill solid blue wet fill solid blue lakes fill solid yellow pos-fill fill solid gray neg-fill marker star 0.2 m-sources axis axis-tt bar left _label T [s] _ticks axis axis-grav bar left _label mGal _ticks axis axis-dist bar top _label X [km] _ticks axis axis-depth bar left _label Z [km] # buttons and the corresponding commands to the ray tracer button button1 Rerun -flow- reset reset button button3 Reload Example-1 load-s {fname} button button4 Save Example-1 save-s {fname} button button5 Print Example-1 ps ### Create cross-section image *call image im-section ##### plot the model and rays 1 1. 6 6 project xy width range-x -20 310 range-y -5 45 yes velmod Example-1 layers l-layers #cells vert cells-vert #cells diag cells-diag #vs shade pal-vs #vp shade pal-vp vp column l-columns vp contour l-contour sources m-sources sources label Shot %d sources traces source-%d # interactive data input rays forward l-rays1 #rays reverse l-rays2 #rays pal-raycodes tau velocity 5.7 tau line l-tau traces rf-section wiggle line wiggle-line wiggle fill-pos pos-fill wiggle fill-neg neg-fill object axis-depth object axis-dist # buttons to control modelling object button3 object button4 object button5 object -spacer- ### Display the interactive images *call gui im-section # add interactive editor for parameters of "genrate" proprty genrate.1 # make axes X coordinated in all sections layout vert x zoom caption Interactive modeling of Receiver Functions #dockers #otree #otree graph prop