Skip to content

Instantly share code, notes, and snippets.

@vmirly
Created October 14, 2013 15:05
Show Gist options
  • Select an option

  • Save vmirly/6977133 to your computer and use it in GitHub Desktop.

Select an option

Save vmirly/6977133 to your computer and use it in GitHub Desktop.
Generate scaling factors for each atom of a PDB based on atom type
awk '{
printf "%s\n", substr($3,1,1)
}' amber/heat7.chm.pdb | awk '{
if ($1=="H") {print "Sfj(",NR,")=0.85"}
else if ($1 == "C") {print "Sfj(",NR,")=0.72"}
else if($1 == "N") {print "Sfj(",NR,")=0.79"}
else if($1 == "O"){print "Sfj(",NR,")=0.85"}
else if($1=="P"){print "Sfj(",NR,")=0.86"}
else if($1=="S"){print "Sfj(",NR,")=0.96"}
}'
@vmirly
Copy link
Author

vmirly commented Oct 14, 2013

Amber scaling factors:
H 0.85
C 0.72
N 0.79
O 0.85
P 0.86
S 0.96

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment