shell scripting on linux

ibabhelix · 10-23-2009, 01:39 AM

@ghostdog74 - thank you so much

the code u have written what is it called ? how do i learn it ?

acid_kewpie · 10-23-2009, 02:02 AM

it's called shell script, like you asked for.

ibabhelix · 10-23-2009, 02:44 AM

oh okay ,

like in other progrmming languages is there a command to search for a particular string ?

like if suppose i have two lines

Code:

ATOM    902  C4'   C A  29      
ATOM      3  C4'   G A   1

like this and they have C4' in the same position, when such lines come continuously i want to add a empty line between them, how do i do that ?

i want it like this

Code:

ATOM    902  C4'   C A  29      
 
ATOM      3  C4'   G A   1

please dont think i'am asking everything here and not trying myself, i am trying to learn, but i need help.

could anyone please tel me how to do this ?
please

ibabhelix · 11-01-2009, 11:28 PM

pixellany · 11-02-2009, 12:42 AM

Earlier in this thread you were given some reading suggestions. I asked if you had read them, but no answer.

So now please tell us what you have read......

Quote:

like in other progrmming languages is there a command to search for a particular string ?

MANY....You are using 2 of them in the code above.

ibabhelix · 11-02-2009, 03:27 AM

hey i have a doubt, suppose in the script i include an exe type of file for ex. ./a.out, which takes input as a pdb file namely 1sn3.pdb, now i want my output text file to contain the name 1sn3.txt , how do i do that using the shell script ???
is there any method ???? or any command ???

please let me know

acid_kewpie · 11-02-2009, 03:29 AM

well that depends what on earth the program is doing with the input doesn't it?? If you mean that you need to provide an output filename on the command line, then as everm there are dozens of ways, like using "basename" to remove the main suffix from the variable and adding your own to the shorter string.

ibabhelix · 11-02-2009, 03:33 AM

ok its like this :

Code:

./a.out
clear
cut -c 13-28 *.pdb > prt1t
exec 6<"prt1t" # assign file descriptor 6 to file1

while read  a <&6 # read a line from file1
do
    read  b <&6 # read the 2nd line from file1
    echo "$a" >> prt1  # print both lines from file1
    echo "$b" >> prt1
    echo "$b" >> prt1

done
exec >&6

sed '1,2d' prt1 > prt2
paste prt1 prt2 a6.txt > ans1t
awk '{ if($9 != "") print $0 }' ans1t > ans1
 
rm prt1t prt1 prt2 a6.txt b6.txt c6.txt ans1t

here for the ./a.out , will ask for a file name and if i give tat as 1sn3.pdb, then i want ans1 to be 1sn3.

the input file which is 1sn3.pdb in this case can be any pdb file.like 1ds4.pdb or 4hgf.pdb