R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Find Ray Parameter for Distance

P4X.lin

R Documentation

Find Ray Parameter for Distance

Description

Calculates ray parameter and azimuth of a ray that reaches some given point, using a 2-D iterative optimization search.

Usage

P4X.lin(x, y, zs, zr, ATM = CheckAtm.lin(list()), maxerror = 3)

Arguments

`x`	Abscissa of receiver relative to source (m)
`y`	Ordinate of receiver relative to source (m)
`zs`	Elevation of source (m)
`zr`	Elevation of receiver (m)
`ATM`	Linear atmosphere
`maxerror`	Maximum permitted error for solution (m)

Details

'maxerror' should be set with the necessary arrival time precision in mind. Smaller values of maxerror will produce more accurate arrival locations (and arrival times), but will require longer calculation times.

Value

List with following elements:

`p`	ray parameter
`az`	azimuth
`error`	distance between inputs (x,y) and final position of result ray

Author(s)

Jake Anderson

Examples

ATM = CheckAtm.lin()
P4X.lin(100, 100, 100, 00, ATM)
P4X.lin(100, 100, 100, 00, ATM, 0.01)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(AtmRay)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/AtmRay/P4X.lin.Rd_%03d_medium.png", width=480, height=480)
> ### Name: P4X.lin
> ### Title: Find Ray Parameter for Distance
> ### Aliases: P4X.lin
> ### Keywords: misc
> 
> ### ** Examples
> 
> ATM = CheckAtm.lin()
> P4X.lin(100, 100, 100, 00, ATM)
$p
[1] 0.002474232

$az
[1] 45

$error
[1] 0.0006428255

> P4X.lin(100, 100, 100, 00, ATM, 0.01)
$p
[1] 0.002474232

$az
[1] 45

$error
[1] 0.0006428255

> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>