=======================
 filterpeaks.awk 
=======================
#!/software/bin/gawk -f
# 
#@  filter peaklist for diagonal,water,lowest/highest intensities 
# Usage:  
#   filterpeaks.awk diagHskip=0.8 \
#           waterHskip=1 waterskip=0.02 hH_2.peaks > hH_2-f.peaks
#
# Command line options:
#   diagHskip=0.03    - skip proton diagonal
#   diagCskip=0.3     - skip carbon diagonal
#
#   diagHoffsetskip=0.24,0.38
#                     - skip peaks that have 0.24<|ppmx-ppmy|<0.38
#                       for removal of sidebands
#
#   waterHskip=0.3    - skip protons closer than 0.3 to water ppm in acquisition
#   waterskip=0.02    - skip protons closer than 0.03 to water ppm in any proton dimension.
#
#    The script uses 
#   waterppm=4.7          : set the ppm-value of water
#   keepass=1    [def=1]  : do not filter out fully assigned peaks 
#   keepass=2    [def=1]  : do not filter out partially assigned peaks 
#   keepass=0    [def=1]  : ignore assignments 
#
# Experimental:
#   skiplo=10%  : skip 10% lowest peaks
#   skiphi=10%  : skip 10% lowest peaks
#   skiplo=10   : skip peaks with int/height <= 10
#   skiphi=10   : skip peaks with int/height >= 10
#  
#