# NOT RUN {
{
# Coelho et al. (2015):
input <- system.file("extdata",
"uwnd.mon.mean_200hPa_2014JFM.nc",
package = "raytracing")
b <- betaks(u = input)
rt <- ray(betam = b$betam,
u = b$u,
lat = b$lat,
K = 3,
itime = 30,
x0 = -135,
y0 = -30,
dt = 6,
direction = -1)
rp <- ray_path(x = rt$lon, y = rt$lat)
plot(rp, axes = TRUE, graticule = TRUE)
}
# }
Run the code above in your browser using DataLab