Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Please help to fix awk script
Top Forums UNIX for Dummies Questions & Answers Please help to fix awk script Post 302794993 by hanson44 on Wednesday 17th of April 2013 05:23:09 AM
Old 04-17-2013
This may be close to what you want, with some modifications to better handle the initial part of the data in each timestep section.
Code:
$ cat atoms.awk
BEGIN { printf ("CRYST1 200.000 200.000 200.000 90.00 90.00 90.00 P 1 1\n")
  maxatoms = 1000;
  natom = 0
  found_struct = 0
  found_bond = 0
  }

{ if (NF == 5) {
  foundff = 0
  natom++
  fftype[natom] = "UNKNOWN"
  if ($1 ~ /CT/) {
    foundff=1; fftype[natom] = "C"
    }
  else if ($1 ~ /OH/) {
    foundff=1; fftype[natom] = "O"
    }
  else if ($1 ~ /HC/) {
    foundff=1; fftype[natom] = "H"
    }
  else if ($1 ~ /N/) {
    foundff=1; fftype[natom] = "N"
    }
  else if ($1 ~ /H1/) {
    foundff=1; fftype[natom] = "H"
    }
  else if ($1 ~ /HO/) {
    foundff=1; fftype[natom] = "H"
    }
  else if ($1 = "C") {
    foundff=1; fftype[natom] = "C"
    }
  else if ($1 = "O") {
    foundff=1; fftype[natom] = "O"
    }

  getline

  x[natom] = $1
  y[natom] = $2
  z[natom] = $3

  if (foundff == 0)
    printf("PROBLEM : Atom ff type %s not known\n", $6)
    }
  }

function report_frame (iatom) {
  for (iatom=1; iatom <= natom; iatom++) {
    printf ("HETATM %d %2s %d %14.9f %14.9f %14.9f\n",
    iatom, fftype[iatom], iatom, x[iatom], y[iatom], z[iatom])
    }
  printf ("END of frame\n")
  natom = 0
  }

END { report_frame() }

/timestep/ { report_frame() }

Code:
$ awk -f atoms.awk atoms.dat
CRYST1 200.000 200.000 200.000 90.00 90.00 90.00 P 1 1
HETATM 1  C 1    0.000000000  500.000000000   25.000000000
END of frame
HETATM 1  C 1    0.098499632   -3.415112416    1.363374178
HETATM 2  H 2   -0.212166256   -3.863455405    0.308072192
HETATM 3  H 3   -1.497545961   -4.274924074    2.443930026
HETATM 4  H 4   -0.020467827   -3.884045345    3.393900210
HETATM 5  H 5    0.066450838   -5.274859543    2.207872210
HETATM 6  H 6    2.012043588   -4.255526008    1.282247861
HETATM 7  H 7    2.030435317   -2.972730305    2.435165857
HETATM 8  H 8    1.986543784   -2.614081904    0.665157533
HETATM 9  H 9   -4.516000060   -0.639058444    0.182144985
HETATM 10  H 10   -3.847187217   -2.017469231   -0.687331929
HETATM 11  H 11   -4.558244417   -0.651631383   -1.645346738
HETATM 12  C 12   -4.027556611   -0.930595682   -0.756787954
HETATM 13  C 13    1.660000184   -3.278731061    1.464017002
HETATM 14  C 14   -0.373640371   -4.330624723    2.450383850
HETATM 15  C 15   -0.536691403   -1.980928569    1.574219882
HETATM 16  C 16   -0.427892225   -1.357753029    2.570890083
HETATM 17  N 17   -1.208805366   -1.532235679    0.567346363
HETATM 18  C 18   -1.136311398   -2.139234863   -0.278615463
HETATM 19  C 19   -1.903616044   -0.273908404    0.490796116
HETATM 20  C 20   -2.713660887   -0.194578984   -0.823379185
HETATM 21  H 21   -2.144431000   -0.619396144   -1.696199396
HETATM 22  H 22   -1.113333313    0.459747980    0.487359101
HETATM 23  H 23   -2.605232776   -0.253914802    1.355851660
HETATM 24  O 24   -2.883721417    1.170430901   -1.063007617
HETATM 25  H 25   -3.755698183    1.390412071   -0.726998504
END of frame
HETATM 1  C 1    0.593081343   -3.059647383    0.909317053
HETATM 2  H 2    0.432526110   -2.863254539    1.968758844
HETATM 3  H 3    1.943889135   -1.400824972    0.444681557
HETATM 4  H 4    1.997587984   -2.723903724   -0.633048117
HETATM 5  H 5    2.776704671   -2.790314949    0.997173445
HETATM 6  H 6    1.303777680   -5.015186286    1.216356852
HETATM 7  H 7    0.544448266   -4.839385307   -0.392673991
HETATM 8  H 8   -0.335159064   -4.945216690    1.087335425
HETATM 9  H 9   -2.335162001    1.703747125    1.766118970
HETATM 10  H 10   -1.485614704    0.523233248    2.667235397
HETATM 11  H 11   -3.119713679    1.071097349    3.299766915
HETATM 12  C 12   -2.493930014    0.859930348    2.445446594
HETATM 13  C 13    0.538271151   -4.518100324    0.726813102
HETATM 14  C 14    1.873711712   -2.498659898    0.436192443
HETATM 15  C 15   -0.608192609   -2.279925921    0.186396952
HETATM 16  C 16   -0.777784481   -2.458871627   -1.043373240
HETATM 17  N 17   -1.381080530   -1.599443295    0.993979052
HETATM 18  C 18   -1.248085227   -1.605021545    2.010060279
HETATM 19  C 19   -2.357481513   -0.717738092    0.491640491
HETATM 20  C 20   -3.243082218   -0.162632033    1.574424959
HETATM 21  H 21   -3.471572515   -0.988386471    2.248890966
END of frame

This User Gave Thanks to hanson44 For This Post:
snysmumrik
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

awk / shell - Fix broken lines and data

Gurus, I am struggling with a issue and thought I could use some of your expertise. Need Help with this I have a flat file that has millions of records 24|john|account ~ info |56| 25|kuo|account ~ journal |58| 27|kim|account ~ journal |59| 28|San|account ~ journal |60|... (3 Replies)
Discussion started by: rimss
3 Replies

2. Shell Programming and Scripting

fix a problem in this script

z=9 i=0 h=02 min=55 while do cat /home/barmecha/test | grep $h:$min >> /home/barmecha/file1 min=`expr $min + 1` if ; then h=`expr $h + 1` fi i=`expr $i + 1` done i have a log file with time wise log in it, this script help me to pull out logs of the give time interval...but the... (8 Replies)
Discussion started by: abhishek27
8 Replies

3. Shell Programming and Scripting

AWK record length fix

Hi Friends, Need some help in AWK. Working on AIX 5 Have been trying the following functionality to make the record length fixed: if( length(record) < 300 ) { printf("%-300s\n", record); } In my opinion it will apply some fillers in the end. Its is not making any... (4 Replies)
Discussion started by: kanu_pathak
4 Replies

4. Shell Programming and Scripting

how to fix this awk script?

i have a log file while looks like this ++ user_a blabla blabla nas_b blabla user_d this is a user_a junk line another junk line user_c nas_m blabla ++ basically most of the lines contain a "user" keywords, and the rest of the lines do not have "user" at all. So I have the... (17 Replies)
Discussion started by: fedora
17 Replies

5. Shell Programming and Scripting

how to fix the column length in a file using Awk Prog

Hi I use the following code to read the file and to fix the length of the column of the record in the file 'Sample.txt' ls Samp* | awk ' { a=$1 } END{ FS="n" for(i=1;i<=NR;i++) { while( getline < a ) { f1=$0; print("Line::",f1); f2=substr(f1,1,10) print("Field1::",f2);... (10 Replies)
Discussion started by: meva
10 Replies

6. Shell Programming and Scripting

How to fix this awk

I have a script which will mask the 9th and 15th column in a record starting with BPR. The record looks like below before my script BPR*C*160860.04*C*ACH*CTX*01*072000326*DA*1548843*3006968523**01*071000013*DA*5529085*100323*VEN The record will be masked after my script parses this... (19 Replies)
Discussion started by: Muthuraj K
19 Replies

7. Shell Programming and Scripting

Fix timestamp with Sed or Awk

Hi I am dealing with the following string: Date: Thur, 13 March 2011 01:01:10 +0000 I asked for help in another topic that converted a similar string: Date: Thur, 13 March 2011 9:50 AM To a 24 hr standard. The problem is that it comes out as: Date: Thur, 13 March 2011 9:50:00 +0000... (4 Replies)
Discussion started by: duonut
4 Replies

8. Shell Programming and Scripting

My script failed and can't fix it ?

Hi , I'd like to give you a little bit idea about my script which is used to get any generated file from remote server using ftp session then organized those file into directories based on their date ( at the end I supposed to have 1 months directories 20130401 20130402 ....20130430 ,... (27 Replies)
Discussion started by: arm
27 Replies

9. Shell Programming and Scripting

awk if condition match and fix print decimal place

Hi All, I have problem in the middle of implementing to users, whereby the complaint is all about the decimal place which is too long. I need two decimal places only, but the outcome from command is always fixed to 6. See the sample : before: Sort Total Site Sort SortName Parts ... (3 Replies)
Discussion started by: horsepower
3 Replies

10. UNIX for Dummies Questions & Answers

awk file subtracted by a fix value - conditioned

Hi all... i have been trying to make this work but I have been failing for 6 hours .. I know it should be something simple that I am missing to it would be great if you can help me ... I want to subtract a fixed value (lets set 1) from any value >=1 from the whole file my file looks like ... (4 Replies)
Discussion started by: A-V
4 Replies
All times are GMT -4. The time now is 12:53 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy