10-05-2011
That always seems to be the first column of that group in your data, though I presume it's not always so.
But that's why it always becomes zero here:
[code]$ awk '{ MIN=0;
for(N=6; N<=8; N++) if($N < MIN) MIN=$N
for(N=6; N<=8; N++) $N -= MIN }' < data
HETATM 1 C UNK 0 0 33.504 34.109 0.00 0.00 C+0
HETATM 2 C UNK 0 0 33.316 32.133 0.00 0.00 C+0
HETATM 3 C UNK 0 0 31.589 30.717 0.00 0.00 C+0
HETATM 4 C UNK 0 0 31.406 28.736 0.00 0.00 C+0
HETATM 5 C UNK 0 0 29.676 27.329 0.00 0.00 C+0
HETATM 6 C UNK 0 0 29.478 25.34 0.00 0.00 C+0
HETATM 7 C UNK 0 0 27.745 23.945 0.00 0.00 C+0
HETATM 8 C UNK 0 0 27.545 21.957 0.00 0.00 C+0
HETATM 9 C UNK 0 0 25.814 20.563 0.00 0.00 C+0
HETATM 10 C UNK 0 0 25.602 18.572 0.00 0.00 C+0
HETATM 11 C UNK 0 0 23.861 17.196 0.00 0.00 C+0
HETATM 12 C UNK 0 0 23.634 15.2 0.00 0.00 C+0
HETATM 13 C UNK 0 0 21.906 13.789 0.00 0.00 C+0
HETATM 14 C UNK 0 0 21.758 11.82 0.00 0.00 C+0
HETATM 15 C UNK 0 0 20.051 10.34 0.00 0.00 C+0
HETATM 16 C UNK 0 0 19.948 8.395 0.00 0.00 C+0
HETATM 17 O UNK 0 0 18.448 7.194 0.00 0.00 O+0
HETATM 18 P UNK 0 0 17.805 5.177 0.00 0.00 P+0
HETATM 19 O UNK 0 0 16.082 4.405 0.00 0.00 O+0
HETATM 20 O UNK 0 0 19.303 5.032 0.00 0.00 O+0
HETATM 21 O UNK 0 0 17.351 4.065 0.00 0.00 O+0
HETATM 22 H UNK 0 0 19.82 7.866 0.00 0.00 H+0
HETATM 23 H UNK 0 0 21.257 8.578 0.00 0.00 H+0
HETATM 24 H UNK 0 0 18.735 10.148 0.00 0.00 H+0
HETATM 25 H UNK 0 0 20.163 10.855 0.00 0.00 H+0
HETATM 26 H UNK 0 0 23.064 11.989 0.00 0.00 H+0
HETATM 27 H UNK 0 0 21.636 11.282 0.00 0.00 H+0
HETATM 28 H UNK 0 0 22.037 14.337 0.00 0.00 H+0
HETATM 29 H UNK 0 0 20.608 13.629 0.00 0.00 H+0
HETATM 30 H UNK 0 0 24.915 15.331 0.00 0.00 H+0
HETATM 31 H UNK 0 0 23.487 14.623 0.00 0.00 H+0
HETATM 32 H UNK 0 0 24.015 17.783 0.00 0.00 H+0
HETATM 33 H UNK 0 0 22.587 17.076 0.00 0.00 H+0
HETATM 34 H UNK 0 0 26.883 18.7 0.00 0.00 H+0
HETATM 35 H UNK 0 0 25.456 17.993 0.00 0.00 H+0
HETATM 36 H UNK 0 0 24.534 20.435 0.00 0.00 H+0
HETATM 37 H UNK 0 0 25.962 21.143 0.00 0.00 H+0
HETATM 38 H UNK 0 0 28.83 22.093 0.00 0.00 H+0
HETATM 39 H UNK 0 0 27.4 21.385 0.00 0.00 H+0
HETATM 40 H UNK 0 0 27.891 24.522 0.00 0.00 H+0
HETATM 41 H UNK 0 0 26.462 23.814 0.00 0.00 H+0
HETATM 42 H UNK 0 0 30.762 25.474 0.00 0.00 H+0
HETATM 43 H UNK 0 0 29.332 24.766 0.00 0.00 H+0
HETATM 44 H UNK 0 0 29.818 27.897 0.00 0.00 H+0
HETATM 45 H UNK 0 0 28.389 27.189 0.00 0.00 H+0
HETATM 46 H UNK 0 0 31.264 28.168 0.00 0.00 H+0
HETATM 47 H UNK 0 0 32.692 28.876 0.00 0.00 H+0
HETATM 48 H UNK 0 0 30.298 30.572 0.00 0.00 H+0
HETATM 49 H UNK 0 0 31.728 31.281 0.00 0.00 H+0
HETATM 50 H UNK 0 0 34.611 32.285 0.00 0.00 H+0
HETATM 51 H UNK 0 0 33.183 31.578 0.00 0.00 H+0
HETATM 52 H UNK 0 0 34.736 35.051 0.00 0.00 H+0
HETATM 53 H UNK 0 0 32.247 34.011 0.00 0.00 H+0
HETATM 54 H UNK 0 0 33.681 34.721 0.00 0.00 H+0
$[/code] Doesn't do what you want.
Last edited by Corona688; 10-05-2011 at 01:59 PM..
Reason: columns off by one
10 More Discussions You Might Find Interesting
1. Shell Programming and Scripting
Hi Gurus,
In my file I have an amount field from position 74 to 87, which contains values starting with '+' as well as '-'. I want to add all positive values in a varible called "CREDIT" and all negative values in a variable "DEBIT". I know, we can use grep to identify values with positive and... (4 Replies)
Discussion started by: berlin_germany
4 Replies
2. Programming
got 1 problem..
can someone help me wit the logic?
Money Money ::operator+(const Money &rhs)const
{
Money temp;
temp.a = a+rhs.a;
temp.b = b+rhs.b;
return temp;
}//end
i got 2 number e.g 6.2 and 3.8 (1 Reply)
Discussion started by: xiaojesus
1 Replies
3. Shell Programming and Scripting
Hello,
For my weather station I have made a little perl script to put the data into cacti. The next problem I have.
I can only get positive numbers or negative numbers.
What do I do:
Though a shell scrip I call the perl script.
Shell script:
#!/bin/sh
cat data.txt | stats.pl
Perl... (4 Replies)
Discussion started by: rbl-blacklight
4 Replies
4. Shell Programming and Scripting
Hi,
I have a file that looks like this:
Jake 2 3 4 6 4 3 -2 -1
Jerry 1 2 3 2 1 7 -6 -1
Timmy -1 -4 -5 -8 9 3 1
I want to find the most positive and negative value for each row and also define its position (based on column #)
So the output would look... (7 Replies)
Discussion started by: gisele_l
7 Replies
5. Shell Programming and Scripting
Hello,
I have a list like this :
1
2
-4
0
-3
-7
5
6 etc.
Is there a way to remove all the positive values and print only the negative values, without using grep, sed or awk?
Thanks,
Prasanna (4 Replies)
Discussion started by: prasanna1157
4 Replies
6. Shell Programming and Scripting
Let, I have three numbers
+00123.25
-00256.54
+00489.23
I need to sum up all those three numbers, after storing them in three variables (say var1, var2, var3).
I used both expr and BC, but they didn't work for me.
But, I am not able to sum up them, as I don't have any idea how to... (13 Replies)
Discussion started by: mady135
13 Replies
7. Shell Programming and Scripting
Dear All,
I have to split a tab delimited file in two files based on the presence of a positive or negative in column number 9 , for example
file:
A 1 5 erg + 6766 0.9889 0.9817 9.01882 erg inside upstream
B 1 8 erg2 + 6766 0.9889 0.9817 -9.22 erg2 inside... (3 Replies)
Discussion started by: paolo.kunder
3 Replies
8. Shell Programming and Scripting
considering the following table:
ID col1 col2 col3 col4
1 -16.06801249 13.49785832 -56.57087607 -27.00500526
2 -1.53315720 0.71731735 -42.03602078 -39.78554623
3 -1.53315190 0.71731587 -42.03601548 ... (3 Replies)
Discussion started by: Birda
3 Replies
9. UNIX for Beginners Questions & Answers
Hi ALL,
I am having semi column separated file as below. I am having negative values for the records starting with 11095. How can I convert that positive number
I tried this below seems not working
sed 's/ \(*\)$/ -\1/;t;s/\(.*\)-/\1/ myfile
myfile... (6 Replies)
Discussion started by: arunkumar_mca
6 Replies
10. UNIX for Beginners Questions & Answers
I have a file that is pipe delimited and in Column F they have number values, both positive and negative. I need to take the one file I am starting with and split it into two separate files based on negative and positive numbers. What is the command to do so? And then I need to also transfer... (4 Replies)
Discussion started by: cckaiser15
4 Replies
LEARN ABOUT DEBIAN
map_slippy
map::slippy(3tcl) Mapping utilities map::slippy(3tcl)
__________________________________________________________________________________________________________________________________________________
NAME
map::slippy - Common code for slippy based map packages
SYNOPSIS
package require Tcl 8.4
package require Tk 8.4
package require map::slippy ?0.4?
::map::slippy length level
::map::slippy tiles level
::map::slippy tile size
::map::slippy tile valid tile levels ?msgvar?
::map::slippy geo 2tile geo
::map::slippy geo 2tile.float geo
::map::slippy geo 2point geo
::map::slippy tile 2geo tile
::map::slippy tile 2point tile
::map::slippy point 2geo point
::map::slippy point 2tile point
_________________________________________________________________
DESCRIPTION
This package provides a number of methods doing things needed by all types of slippy-based map packages.
API
::map::slippy length level
This method returns the width/height of a slippy-based map at the specified zoom level, in pixels. This is, in essence, the result
of
expr { [tiles $level] * [tile size] }
::map::slippy tiles level
This method returns the width/height of a slippy-based map at the specified zoom level, in tiles.
::map::slippy tile size
This method returns the width/height of a tile in a slippy-based map, in pixels.
::map::slippy tile valid tile levels ?msgvar?
This method checks whether tile described a valid tile in a slippy-based map containing that many zoom levels. The result is a bool-
ean value, true if the tile is valid, and false otherwise. For the latter a message is left in the variable named by msgvar, should
it be specified.
A tile identifier as stored in tile is a list containing zoom level, tile row, and tile column, in this order. The command essen-
tially checks this, i.e. the syntax, that the zoom level is between 0 and "levels-1", and that the row/col information is within the
boundaries for the zoom level, i.e. 0 ... "[tiles $zoom]-1".
::map::slippy geo 2tile geo
Converts a geographical location at a zoom level (geo, a list containing zoom level, latitude, and longitude, in this order) to a
tile identifier (list containing zoom level, row, and column) at that level. The tile identifier uses pure integer numbers for the
tile coordinates, for all geographic coordinates mapping to that tile.
::map::slippy geo 2tile.float geo
Converts a geographical location at a zoom level (geo, a list containing zoom level, latitude, and longitude, in this order) to a
tile identifier (list containing zoom level, row, and column) at that level. The tile identifier uses floating point numbers for the
tile coordinates, representing not only the tile the geographic coordinates map to, but also the fractional location inside of that
tile.
::map::slippy geo 2point geo
Converts a geographical location at a zoom level (geo, a list containing zoom level, latitude, and longitude, in this order) to a
pixel position (list containing zoom level, y, and x) at that level.
::map::slippy tile 2geo tile
Converts a tile identifier at a zoom level (tile, list containing zoom level, row, and column) to a geographical location (list con-
taining zoom level, latitude, and longitude, in this order) at that level.
::map::slippy tile 2point tile
Converts a tile identifier at a zoom level (tile, a list containing zoom level, row, and column, in this order) to a pixel position
(list containing zoom level, y, and x) at that level.
::map::slippy point 2geo point
Converts a pixel position at a zoom level (point, list containing zoom level, y, and x) to a geographical location (list containing
zoom level, latitude, and longitude, in this order) at that level.
::map::slippy point 2tile point
Converts a pixel position at a zoom level (point, a list containing zoom level, y, and x, in this order) to a tile identifier (list
containing zoom level, row, and column) at that level.
COORDINATE SYSTEMS
The commands of this package operate on three distinct coordinate systems, which are explained below.
GEOGRAPHIC
Geographical coordinates are represented by Latitude and Longitude, each of which is measured in degrees, as they are essentially angles.
Zero longitude is the Greenwich meridian, with positive values going east, and negative values going west, for a total range of +/- 180
degrees. Note that +180 and -180 longitude are the same meridian, opposite to greenwich.
zero latitude the Equator, with positive values going north and negative values going south. While the true range is +/- 90 degrees the
projection used by the package requires us to cap the range at +/- 85.05112877983284 degrees. This means that north and south pole are not
representable and not part of any map.
TILES
While Geographical coordinates of the previous section are independent of zoom level the tile coordinates are not.
Generally the integer part of tile coordinates represent the row and column number of the tile in question, wheras the fractional parts
signal how far inside the tile the location in question is, with pure integer coordinates (no fractional part) representing the upper left
corner of the tile.
The zero point of the map is at the upper left corner, regardless of zoom level, with larger coordinates going right (east) and down
(south), and smaller coordinates going left (west) and up (north). Again regardless of zxoom level.
Negative tile coordinates are not allowed.
At zoom level 0 the whole map is represented by a single, putting the geographic zero at 1/2, 1/2 of tile coordinates, and the range of
tile coordinates as [0...1].
To go from a zoom level N to the next deeper level N+1 each tile of level N is split into its four quadrants, which then are the tiles of
level N+1.
This means that at zoom level N the map is sliced (horizontally and vertically) into 2^N stripes, for a total of 4^N tiles, with tile coor-
dinates ranging from 0 to 2^N+1.
PIXELS/POINTS
pixel coordinates, also called point coordinates are in essence tile coordinates scaled by the size of the image representing a tile. This
tile size currently has a fixed value, 256.
REFERENCES
[1] http://wiki.openstreetmap.org/wiki/Main_Page
KEYWORDS
geodesy, geography, latitute, location, longitude, map, slippy, zoom
map 0.4 map::slippy(3tcl)