Commit 8073aaf8 by Francois Gygi

qbox_maxforce.sh replaced by qbox_maxforce.py


git-svn-id: http://qboxcode.org/svn/qb/trunk@1846 cba15fb0-1239-40c8-b417-11db7ca47a34
parent ee8f4405
#!/bin/bash
# get the largest force component in a file of "<force> fx fy fz </force>"
# use: qbox_maxforce nat file.r
grep '<force>' $2 | tail -$1 | \
awk '{
if ( mx*mx < $2*$2 ) mx = $2;
if ( my*my < $3*$3 ) my = $3;
if ( mz*mz < $4*$4 ) mz = $4;
} END {printf("%9.2e %9.2e %9.2e\n", mx, my, mz)}' -
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment