#!/bin/bash

config=eprem.cfg
stream_row=( 3 )
stream_column=( 1 3 5 )
stream_face=( 2 3 4 5 )
elevation=90
azimuth=-90
rotation=-90
grid_lims=( -3 3 )
days=( 2 3 4 5 )
savedir=images

python ~/open-eprem/eprem-analysis/scatter3D.py \
    --config "$config" \
    --stream-row ${stream_row[@]} \
    --stream-column ${stream_column[@]} \
    --stream-face ${stream_face[@]} \
    --show-points \
    --elevation "$elevation" --azimuth "$azimuth" --rotation "$rotation" \
    --grid-lims ${grid_lims[@]} \
    --time ${days[@]} --time-unit day \
    --quantity Vr \
    --unit "km/s" \
    --min 300 --max 1200 --nbins 9 \
    --colormap viridis \
    --save-path "${savedir}"/stream-Vr.png \
    --verbose

python ~/open-eprem/eprem-analysis/scatter3D.py \
    --config "$config" \
    --stream-row ${stream_row[@]} \
    --stream-column ${stream_column[@]} \
    --stream-face ${stream_face[@]} \
    --show-points \
    --elevation "$elevation" --azimuth "$azimuth" --rotation "$rotation" \
    --grid-lims ${grid_lims[@]} \
    --time ${days[@]} --time-unit day \
    --quantity Bphi \
    --unit "nT" \
    --min -20 --max 0 --nbins 10 \
    --colormap viridis_r \
    --save-path "${savedir}"/stream-Bphi.png \
    --verbose

python ~/open-eprem/eprem-analysis/scatter3D.py \
    --config "$config" \
    --stream-row ${stream_row[@]} \
    --stream-column ${stream_column[@]} \
    --stream-face ${stream_face[@]} \
    --show-points \
    --elevation "$elevation" --azimuth "$azimuth" --rotation "$rotation" \
    --grid-lims ${grid_lims[@]} \
    --time ${days[@]} --time-unit day \
    --quantity rho \
    --unit "cm^-3" \
    --min 1e-1 --max 1e2 --log \
    --colormap viridis \
    --save-path "${savedir}"/stream-rho.png \
    --verbose

