#! /bin/bash
# Last edited on 2021-07-31 15:42:39 by jstolfi

# Assumes that the input file has fields
#  "{k} {weight[k]} {xs[k]} {ys[k]} {B[0][k]} ... {B[nb-1][k]}"

dfile="$1"; shift   # File name.
elix="$1"; shift    # Basis element index, from 0.
title="$1"; shift   # Plot title.

gnuplot --persist <<EOF
set term x11
set size square

set pm3d map
set palette rgbformulae 22,13,-31
set title "${title}"

bcol = 5 + ${elix}

splot "${dfile}" using 3:4:(column(bcol)) title "${title}"
 
EOF