#! /bin/csh
set reduceatom=`grep -v REMARK $1|grep -v CRYS|awk 'NR==1 {print $2-1}'`
set reduceresi=`grep -v REMARK $1|grep -v CRYS|awk 'NR==1 {print $6-1}'`

echo "$reduceatom $reduceresi"

if ( $reduceatom != 0 && $reduceresi != 0 ) then
        if ( $2 == write ) then
                grep -v REMARK $1|grep -v CRYS|awk -v var=$reduceatom -v varb=$reduceresi 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2-var,$3,$4,$6-varb,$7,$8,$9)}'|sed 's/HOH /TIP3/' >tmp_$1 ; mv tmp_$1 $1
                echo TER >> $1
        else
                grep -v REMARK $1|grep -v CRYS|awk -v var=$reduceatom -v varb=$reduceresi 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2-var,$3,$4,$6-varb,$7,$8,$9)}' $1
                echo TER
        endif
endif

if ( $reduceatom == 0 && $reduceresi != 0 ) then
        if ( $2 == write ) then
                grep -v REMARK $1|grep -v CRYS|awk -v varb=$reduceresi 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2,$3,$4,$6-varb,$7,$8,$9)}'|sed 's/HOH /TIP3/' >tmp_$1 ; mv tmp_$1 $1
                echo TER >> $1
        else
                grep -v REMARK $1|grep -v CRYS|awk -v varb=$reduceresi 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2,$3,$4,$6-varb,$7,$8,$9)}'
                echo TER
        endif
endif

if ( $reduceatom != 0 && $reduceresi == 0 ) then
        if ( $2 == write ) then
                grep -v REMARK $1|grep -v CRYS|awk -v var=$reduceatom 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2-var,$3,$4,$6,$7,$8,$9)}'|sed 's/HOH /TIP3/' >tmp_$1 ; mv tmp_$1 $1
                echo TER >> $1
        else
                grep -v REMARK $1|grep -v CRYS|awk -v var=$reduceatom 'NF>5 {printf("ATOM%7d  %-4s%3s%6s%12s%8s%8s        1.00      Z\n",$2-var,$3,$4,$6,$7,$8,$9)}'
                echo TER
        endif
endif