ColliderReverser: Reverses a time step to coliision point

Description Usage Arguments Details Value Author(s) Examples

View source: R/ColliderReverser.R

Description

If a collision happens within a time bin establishes the proportion of time elapsed

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
ColliderReverser(
  min_separation,
  continent_1_longitude_t0,
  continent_1_latitude_t0,
  continent_1_longitude_t1,
  continent_1_latitude_t1,
  continent_2_longitude_t0,
  continent_2_latitude_t0,
  continent_2_longitude_t1,
  continent_2_latitude_t1,
  continent_1_euler_longitude,
  continent_1_euler_latitude,
  continent_2_euler_longitude,
  continent_2_euler_latitude,
  continent_1_degrees_per_step,
  continent_2_degrees_per_step,
  EarthRad = 6367.4447,
  Warn = TRUE
)

Arguments

min_separation

The minimum separation between continents in kilometres.

continent_1_longitude_t0

Decimalised longitude of first continent at t0.

continent_1_latitude_t0

Decimalised latitude of first continent at t0.

continent_1_longitude_t1

Decimalised longitude of first continent at t1.

continent_1_latitude_t1

Decimalised latitude of first continent at t1.

continent_2_longitude_t0

Decimalised longitude of second continent at t0.

continent_2_latitude_t0

Decimalised latitude of second continent at t0.

continent_2_longitude_t1

Decimalised longitude of second continent at t1.

continent_2_latitude_t1

Decimalised latitude of second continent at t1.

continent_1_euler_longitude

Decimalised longitude of Euler pole for first continent.

continent_1_euler_latitude

Decimalised latitude of Euler pole for first continent.

continent_2_euler_longitude

Decimalised longitude of Euler pole for second continent.

continent_2_euler_latitude

Decimalised latitude of Euler pole for second continent.

continent_1_degrees_per_step

Degrees per step (speed) of first continent.

continent_2_degrees_per_step

Degrees per step (speed) of second continent.

EarthRad

Radius of the Earth in kilometres.

Warn

Whether or not to print warnings.

Details

Nothing yet.

Value

Proportion (0 to 1) of time step at which the minimum separation distance collision occurs.

Author(s)

Graeme T. Lloyd graemetlloyd@gmail.com

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
min_separation <- 500

continent_1_euler_longitude <- 2
continent_1_euler_latitude <- 89

continent_2_euler_longitude <- 4
continent_2_euler_latitude <- -89

continent_1_degrees_per_step <- 1
continent_2_degrees_per_step <- 1

continent_1_longitude_t0 <- 3
continent_1_latitude_t0 <- 0
continent_2_longitude_t0 <- -2
continent_2_latitude_t0 <- 1

continent_1_start_bearing <- 
  BearingBetweenTwoLongLatPoints(continent_1_euler_longitude,
  continent_1_euler_latitude, continent_1_longitude_t0, continent_1_latitude_t0)
continent_2_start_bearing <-
  BearingBetweenTwoLongLatPoints(continent_2_euler_longitude,
  continent_2_euler_latitude, continent_2_longitude_t0, continent_2_latitude_t0)
continent_1_euler_distance <-
  GreatCircleDistanceFromLongLat(continent_1_euler_longitude, continent_1_euler_latitude,
  continent_1_longitude_t0, continent_1_latitude_t0)
continent_2_euler_distance <-
  GreatCircleDistanceFromLongLat(continent_2_euler_longitude, continent_2_euler_latitude,
  continent_2_longitude_t0, continent_2_latitude_t0)

continent_1_longitude_t1 <- EndPoint(continent_1_euler_longitude, continent_1_euler_latitude,
  continent_1_start_bearing + continent_1_degrees_per_step, continent_1_euler_distance,
  EarthRad = 6367.4447)$longitude
continent_1_latitude_t1 <- EndPoint(continent_1_euler_longitude, continent_1_euler_latitude,
  continent_1_start_bearing + continent_1_degrees_per_step, continent_1_euler_distance,
  EarthRad = 6367.4447)$latitude
continent_2_longitude_t1 <- EndPoint(continent_2_euler_longitude, continent_2_euler_latitude,
  continent_2_start_bearing + continent_2_degrees_per_step, continent_2_euler_distance,
  EarthRad = 6367.4447)$longitude
continent_2_latitude_t1 <- EndPoint(continent_2_euler_longitude, continent_2_euler_latitude,
  continent_2_start_bearing + continent_2_degrees_per_step, continent_2_euler_distance,
  EarthRad = 6367.4447)$latitude

GreatCircleDistanceFromLongLat(continent_1_longitude_t0, continent_1_latitude_t0,
  continent_2_longitude_t0, continent_2_latitude_t0)
GreatCircleDistanceFromLongLat(continent_1_longitude_t1, continent_1_latitude_t1,
  continent_2_longitude_t1, continent_2_latitude_t1)

ColliderReverser(min_separation, continent_1_longitude_t0, continent_1_latitude_t0,
  continent_1_longitude_t1, continent_1_latitude_t1, continent_2_longitude_t0,
  continent_2_latitude_t0, continent_2_longitude_t1, continent_2_latitude_t1,
  continent_1_euler_longitude, continent_1_euler_latitude, continent_2_euler_longitude,
  continent_2_euler_latitude, continent_1_degrees_per_step, continent_2_degrees_per_step)

laurasoul/dispeRse documentation built on May 25, 2021, 4:50 a.m.