Orientation analysis

Archived notebook using Python and R

Code
using DrWatson
@quickactivate
include(srcdir("main.jl"))
log10 (generic function with 17 methods)
Code
combine(dfs) = reduce(vcat, dfs, cols=:intersect)

wind_df = load_wind();
jno_df = load_jno();

df = combine([jno_df, wind_df]);
┌ Warning: automatically converting Arrow.Timestamp with precision = MICROSECOND to `DateTime` which only supports millisecond precision; conversion may be lossy; to avoid converting, pass `Arrow.Table(source; convert=false)
└ @ Arrow /Users/zijin/.julia/packages/Arrow/5pHqZ/src/eltypes.jl:273
┌ Warning: automatically converting Arrow.Timestamp with precision = NANOSECOND to `DateTime` which only supports millisecond precision; conversion may be lossy; to avoid converting, pass `Arrow.Table(source; convert=false)
└ @ Arrow /Users/zijin/.julia/packages/Arrow/5pHqZ/src/eltypes.jl:273
Code
axis = (yscale=identity,);
layer = data(df) * mapping(color=year_map, layout=:dataset);
Code
maps = [:θ_vn :θ_vk :θ_nk]
plot_dist(layer; maps = maps, axis = axis)
Code
plot_dist(layer; maps = [:Vn_x :Vn_y :Vn_z], axis = axis, datalimits = datalimits_f)
Code
plot_dist(layer; maps = [:v_x :v_y :v_z], axis = axis, datalimits = datalimits_f)