:!0L 77R7V7\7`7h7l7r7v7dUdbredelse af P-blger i Jorden picenter   x akse   y akse 00"Arial"Arial"Arial"ArialR0*xJR0*yJRflyd*xJRflyd*yJRfast*xJRfast*yJxy  XX@XX@-7.000e37.000e3-7.000e37.000e3UUUUUU7@XGb3:Pjo Epicenter^+Ή.mc @"ArialJorddata:p)@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@LLLL Jordradier:p_>@1@1@@@L@ArialJordens gennemsnitsradiusR0kmL@Radius af den flydende kerneRflydkmL"@Radius af den faste kerneRfastkmL@ Klar - Lkke:Q$2{ Beregning af P-blgers udbredelse i (den runde) Jord } { Peter Husby 22. februar 1998. } {I sjletabellen 'Jorddata' er } {dtab : dybde i km (mlt fra overfladen). vptab : P-blgehastighed i km/s } {vstab: S-blgehastighed i km/s . rhotab : densitet i g/cm } {De trigonometriske beregninger indstilles til 'grader': } TrigPref:=Degrees Delta:=0 {Begyndelsesstedets polvinkel } b:=b0 {Blgens startvinkel i grader } dL:=10 {Lbelngde i hvert "lag" i km} r1:=R0 {Startstedets radius} tp:= 0 {Lbetiden 0-stilles} f:=1 {Fortegn i cos-relation: + indad, - udad} N:=-1 {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 specialudgave af asin (arcus-sinus = sin = invers sinus):} Func omvsin(z) If abs(z)<1 Then return asin(z) If z=1 Then f:=-1 {Blgen skal nu udad} return 90 {Symmetripunktet!} Endif If z>1 Then {"Totalreflektion"} f:=-1 {Blgen skal udad} return i {b = i} Endif If z<=-1 Then return -90 Endfunc {Definition af specialudgave af acos (arcus-cosinus = cos = invers cosinus:} Func omvcos(z) If abs(z)<1-1E-10 Then return acos(z) If z>=1-1E-10 Then return 0 If z<=-1 OR r2 <= dL Then {Blgen gr gennem (0,0)} f:=-1 {Blgen skal udad} return 90 Endif Endfunc z{(x,y) er blgefrontens aktuelle sted i sdvanlige koordinater} x:=r1*sin(Delta) {Vinklen Delta er} y:=r1*cos(Delta) {nul p y-aksen! } d1:=R0-r1 {Dybden i 1. lag } vp1:=vp(d1) {P-blgehastighed i det aktuelle jordlag} r2:=sqr(r1^2+dL^2-f*2*r1*dL*cos(b)) {Ny radius fundet fra cos-relationen} dDelta:=omvcos((r1^2+r2^2-dL^2)/(2*r1*r2)) {Vinkelforgelse fundet fra cos-rel.} If b < 0 Then dDelta:=-dDelta {acos giver kun positive vinkler! } If not(b=0 and r2 <= dL) Then Delta:=Delta+dDelta {Ny vinkel} Else Delta:=180 {Blgen er i centrum} f:=-1 {og skal ud igen } Endif dtp := dL/vp1 {Lbetiden i det aktuelle lag i sekunder} tp := tp + dtp {Den samlede lbetid i sekunder} d2:=R0-r2 {Dybden i 2. lag} If d2 <= -dL Then r2:=r1 dDelta:=2*dDelta N:=N+1 If N=M Then STOP {efter M spejlinger stoppes} f:=1 {blgen spejles} Endif vp2:=vp(d2) {P-blgehastighed i nste lag} if d2 > 5121 AND S=1 Then f:=-1 Endif i:=b+f*dDelta {Indfaldsvinklen ved det nste lag = b + dDelta} If abs(i) > 90 Then f:=-1 {Blgen skal udad } i:=180-i {i = 180-(b+dDelta)} Endif b:=omvsin(sin(i)*vp2/vp1) {Brydningsvinklen i nste lag (fra brydningsloven)} r1:=r2 {Gr klar til beregning i nste "lag"}  Enhedscirkel:p BQ@PPPP@Arialss~+2LxJkmcos(s)L?yJkmsin(s)L Resultater:P>@@@+@@4@Arial Lbetid i mintp/60L ;*@ Vinkelafstand fra epicentret i DeltaLcʿ;M@Afstand langs overfladen i km2*pi*R0*Delta/360L}nLe@ Startvinkel:P{>@1@@@@+@ArialStartvinkel i b0LI@ antal refleksioner ved overfladeM4@"S=1 giver spejling p indre kerne S4