5!0L 2HРvT~vJorddata5pQ>@PPP@ArialdtabkmL)@@Y@i@r@y@@@@@@@@@@@ @@@0@@P@@@@p@@@ @@@@h@0@@@v@@@P@@@@p@8@@vptabkm/sL(\@)@@Q @)\(!@\(!@RQ#@{Gz$@p= #%@Q&@\(&@Qk'@{G'@(@(\)@Gz)@(\*@ ףp=*@+@HzG+@HzG+@333333 @q= ףp @q= ף @Q!@Gz"@(\"@Q#@ffffff#@(\#@Q$@Q$@zG$@ffffff#@RQ&@\(\&@Qk&@&@ ףp=&@Gz&@&@Q&@vstabkm/sL333333@)ffffff@@ffffff@ ףp= @(\@HzG@p= ף@Q@Q@q= ףp@@{Gz@Q@RQ@Gz@Gz@{Gz@ףp= @333333@333333@ @ @ @ @ @ @ @ @rhotabg/cmLRQ @)(\ @ ףp= @(\ @ffffff @ ףp= @Q@Q@RQ@(\@Q@333333@p= ף@Q@Q@(\@\(\@(\@)\(@Q@Q@\("@p= ##@33333#@Q8$@33333$@Q%@Q%@zG%@Q8&@ ףp=&@R&@Gz'@ffffff,@0@0@(\0@ 1@Q1@)\(1@q= ף01@3333331@ Klar - Lkke5%xx?Q2t{ Beregning af seismiske blger i en flad Jord } { Peter Husby 19. marts 1998. } TrigPref:=Degrees bp:=b0 { Startvinkel fo P-blgen i grader. } bs:=b0 { Startvinkel fo S-blgen i grader. } LagNr:=0 { Nulstilling af parameter til bregning. } xp,xs,y:=0,0,0 { Blgernes startpunkt. } Nmax:=50 { Maksimalt antal lag } x1,x2:=-500,15000 { Vrdier til tegning af lagene } For nr:=0 To Nmax Do { Data til tegning af lagene beregnes } xlag[2*nr]:=x1 ylag[2*nr]:=-nr*dL xlag[2*nr+1]:=x2 ylag[2*nr+1]:=-nr*dL x1,x2:=x2,x1 Next {Definition af vP(d) : P-blgehastighed i dybden d. } {vP(d) findes ved liner interpolation af tabelvrdierne: } Func vP(d) If d >= 0 Then nr:=0 while d >= dtab[nr] do nr:+1 {finder den mindste tabelvrdi strre end d} v:=vptab[nr]-(vptab[nr]-vptab[nr-1])*(dtab[nr]-d)/(dtab[nr]-dtab[nr-1]) return v Else return vptab[0] Endif Endfunc {Definition af vS(d) : S-blgehastighed i dybden d. } {vS(d) findes ved liner interpolation af tabelvrdierne: } Func vS(d) If d >= 0 Then nr:=0 while d >= dtab[nr] do nr:+1 {finder den mindste tabelvrdi strre end d} v:=vstab[nr]-(vstab[nr]-vstab[nr-1])*(dtab[nr]-d)/(dtab[nr]-dtab[nr-1]) return v Else return vstab[0] Endif Endfunc LagNr:=LagNr+1 y:=-LagNr*dL dxp:=-y*tan(bp) { Beregning af P-blgen } xp:=xp+dxp dxs:=-y*tan(bs) { Beregning af S-blgen } xs:=xs+dxs vp1:=vP((LagNr-.5)*dL) { P-blgehastighed midt i det aktuelle jordlag.} vp2:=vP((LagNr+.5)*dL) { P-blgehastighed midt i nste lag.} ip:=bp { Indfaldsvinklen ved det nste lag = brydningsvinklen i det aktuelle.} bp:=asin(sin(ip)*vp2/vp1) { Brydningsvinklen i nste lag (fra brydningsloven).} vs1:=vS((LagNr-.5)*dL) { P-blgehastighed midt i det aktuelle jordlag.} vs2:=vS((LagNr+.5)*dL) { P-blgehastighed midt i nste lag.} is:=bs { Indfaldsvinklen ved det nste lag = brydningsvinklen i det aktuelle.} bs:=asin(sin(is)*vs2/vs1) { Brydningsvinklen i nste lag (fra brydningsloven).} If LagNr >= AntalLag Then Stop  Konstanter5P4X>@F@@9@@*@ArialLbedybde i kmLL@@Antal lag i modellenAntalLagL@Lbedybde i hvert lagdL L/AntalLagLy@Startvinkel i graderb0L>@LLLLagene5puL>@PPP@ArialxlagL@eL@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@L@@@L@ylagLeyy@@  @@00PPpp  @@hh00PPpp88XX  xx@@hh00\\$$Tg/Model for udbredelse af P- og S-blger. { Blge-forlbetJorden betragtes som flad.  Konstanter5pVandret afstand fra epicenter i km.@1@Dybde under overfladen i km.L00"Arial"Arial"Arial"Arial xlagylagxpyxsy L@@@x10x20-L5e30]6aC.Gg/7%{%0!5PA