#!/bin/bash

ni=96 # number of images
l=`wc ../1-pot/pot.txt | awk '{print $1}'` # number of data points (you can use only a portion of the data)

# using generalized wham algorithm (exclude the block number)
awk '{for(i=1;i<=NF;i++) if(i!=2) printf "%s ",$i; printf"\n"}' ../1-pot/pot.txt | \
../../../../codes/npwham.0/npwham -w $ni -l $l -t 310 > density.txt 2> density.err
# use npwham -h for the description of arguments (-w: # of windows, -l: # of data points, -t: temperature)