#! /bin/bash
#name of the pdb, psf and trajectory file
nom=helix_5_zwit
# Below is a list of all charged residue of the protein to loop through (to change depending on the PROT)
for i in 2 4 15 20;
do
echo i
export i
export name

awk '{print $ 3 }' ${nom}_GOUY_$i.dat > ${nom}_GOUY_$i.tmp
./Egc.o ${nom}_GOUY_$i.tmp 



done
