For the third weekend in a row, U.S. Coast Guard officials searched
  the area below the Piscataqua River Bridge Sunday for a person
  suspected to have jumped from the Interstate 95 passageway. 



  Authorities received a report of a white male between 15 and 18
  years old having climbed over the rail of the bridge at about 2
  p.m., according to Bill Abt, day officer at the Coast Guard's
  Portsmouth Harbor station in New Castle. 

xstar
   type 2 dummies 2  end
end

#!/bin/bash



declare -r FNO=0
declare -r LNO=999999
declare -r INFILE=EEL5kL01-00001.sgy
declare -r OUTFILE=EEL5kL01-00001-wbt.sgy
declare -r GAINFILE=default.gains 



OS=`uname -s`
sioseis=${HOME}/x/src/sioseis/sioseis.${OS}
$sioseis << EOF 
procs diskin header wbt3 gains prout diskoa end
!
diskin
  ipath $INFILE end
end
!
header 
 fno $FNO lno $LNO
 l6 = l3 end
end
! ORIGINAL thres 0.000000005 track 0.005 end
wbt3
  thres 0.000000005 track 0.05 end
end
!
gains
 fno $FNO lno $LNO addwb yes type 9 TGP
`cat $GAINFILE`
 end
end
!
! wbt = l50
prout
 indices l1 l2 l3 l4 l6 l50 fno $FNO lno $LNO end
end
!
diskoa
 opath $OUTFILE fno $FNO lno $LNO ftr 0 ltr 1 end
end
!
end
EOF

    0.0000    1.0000
    0.0005    1.0000
    0.0020    5.0000
    0.005     10.0000
    0.0200    15.0000 
    0.0500    15.0000
    0.1       40.000

./wbt.bash | grep -v WARNING | awk '{print $1,$4,$6}' | \
grep '1\.0' > foo

gnuplot
plot 'foo' using 1:3 with l

wbt3
  thres 0.0000000001 track 0.05 end
end

#!/bin/bash
declare -ri VERBOSITY=2 # TRACE



project="Eel_River_Basin"



for file in *.sgy; do
    echo
    echo ========    $file    ========
    name=${file%%.sgy}



    db="${file}.db"



    if [ ! -f $db ]; then
	echo "WARNING: missing database ($db)... generating assuming SEGY REV 1"
	segysql.py -f $db $file --verbosity=$VERBOSITY --warn-only
    fi



    ShotpointName=EnsembleNo
    minShot=`sqlite $db "select min($ShotpointName) from segyTrace"`
    maxShot=`sqlite $db "select max($ShotpointName) from segyTrace"`
    echo "Shot range:      $minShot ... $maxShot"



    minDelay=`sqlite $db 'SELECT min(Delay) FROM segyTrace'`
    maxDelay=`sqlite $db 'SELECT max(Delay) FROM segyTrace'`
    echo "delay range:     $minDelay ... $maxDelay"



    TraceSamples=`sqlite $db 'SELECT min(TraceSamples) FROM segyTrace'`
    SampleInterval=`sqlite  $db 'select max(sampleinterval) from segyTrace'`
    SampleInterval=`python -c "print \"%.7f\" %($SampleInterval / 1000000.)"`
    echo "Sample Interval: $SampleInterval"
    #SampleInterval=0.000083 # Force if reading SampleInterval from xstar file.



    timeStart=`python -c "print \"%0.2f\" % ($minDelay/1000.)"`
    timeEnd=`python -c "print ($maxDelay/1000. + $SampleInterval * $TraceSamples)"`
    echo "Time range:      $timeStart ... $timeEnd"



    pltWidth=350 # The default
    # Screen width / number of shots we want to see
    spa=`python -c "print \"%.2f\" % ($pltWidth/float($maxShot-$minShot))"`
    echo "Spacing:         $spa"



    # SPA of 0.5 showed some
    # FIX: how to calculate a good timescale?
    #      time scale. The scale for the time axis in mm./second.
    echo
    makepltsegy.py --verbosity=$VERBOSITY -o $name.pltsegy -i $name.sgy \
	-s $minShot -S $maxShot \
	--spa=$spa --timescale=400 \
	--tbi=$timeStart --tfi=$timeEnd \
	--title="$project $name $minShot ... $maxShot"
done

#!/bin/bash
# Why am I writing this in bash???
declare -ri VERBOSITY=2
project="Eel_River_Basin"
name=EEL5kL01-00001
db="${name}.xstar.db"
if [ ! -f $db ]; then
    echo "ERROR: missing database ($db)"
    exit
fi
minShot=`sqlite $db 'select min(Shotpoint) from segyTrace'`
maxShot=`sqlite $db 'select max(Shotpoint) from segyTrace'`



echo shotrange = $minShot ... $maxShot



minDelay=`sqlite $db 'SELECT min(Delay) FROM segyTrace'`
maxDelay=`sqlite $db 'SELECT max(Delay) FROM segyTrace'`



# FIX: make sure min/max are the same
TraceSamples=`sqlite $db 'SELECT min(TraceSamples) FROM segyTrace'`



#SampleInterval=`sqlite  $db 'select max(sampleinterval) from segyTrace'`
# 41 that is in the db file is NOT ok.  It is 83
SampleInterval=0.000083



timeStart=`python -c "print \"%0.2f\" % ($minDelay/1000.)"`
timeEnd=`python -c "print ($minDelay/1000. + $SampleInterval * $TraceSamples)"`



echo $minDelay  $SampleInterval $TraceSamples



echo "Time range: $timeStart $timeEnd"



pltWidth=350 # The default



# Screen width / number of shots we want to see
spa=`python -c "print \"%.2f\" % ($pltWidth/float($maxShot-$minShot))"`
echo "Spacing = $spa"



# SPA of 0.5 showed some
makepltsegy.py --verbosity=$VERBOSITY -o EEL5kL01-00001.pltsegy -i $name.sgy \
  -s $minShot -S $maxShot \
  --spa=$spa --timescale=200 \
  --title="$project $name $minShot ... $maxShot"

19705360 EEL5kL05.xstar



18789077 EEL5kL05.xstar-compress.Z
15711363 EEL5kL05.xstar-zip.zip
15704232 EEL5kL05.xstar-gzip.gz
13580310 EEL5kL05.xstar-7z.7z
13549235 EEL5kL05.xstar-rzip.rz
13539313 EEL5kL05.xstar-bzip2.bz2



TESTING: seismic.xstar



163843600 seismic.xstar
161645283 seismic.xstar-compress.Z
134228345 seismic.xstar-zip.zip
134142038 seismic.xstar-gzip.gz
119504153 seismic.xstar-rzip.rz
119409223 seismic.xstar-bzip2.bz2
118475932  seismic.xstar-7z.7z

#!/bin/bash
# Testing for compression programs
infile=test0.txt
infile=EEL5kL05.xstar
infile=seismic.xstar
echo TESTING: $infile



rm -f $infile-*



for type in bzip2 gzip rzip zip compress; do
    cp $infile $infile-$type
done



bzip2 -9 $infile-bzip2
gzip -9 $infile-gzip
rzip -9 $infile-rzip
zip $infile-zip.zip $infile-zip
compress $infile-compress
7za a $infile-7z.7z -mx=9 $infile-7z 



ls -l $infile | awk '{print $5,$9}'
ls -lS $infile-*.* | awk '{print $5,$9}'



rm -f $infile-*

[Distutils] setuptools-0.4a2: Eggs, scripts, and __file__
Ryan Tomayko rtomayko at gmail.com
Tue Jun 14 00:00:56 CEST 2005



    * Previous message: [Distutils] setuptools-0.4a2: Eggs, scripts, and __file__
    * Next message: [Distutils] setuptools-0.4a2: Eggs, scripts, and __file__
    * Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]



On Jun 13, 2005, at 12:21 PM, Phillip J. Eby wrote:
> Are you distributing an application, or a library?  If you're  
> distributing a library, you don't need require() in library code.   
> If it's an application, you can handle your own dependencies by  
> force-installing the eggs in the application script directory using  
> EasyInstall, and then just require() your main package.



Libraries, applications, and libraries with helper scripts.



>> This should give me all of the benefits of
>> eggs when I'm using them and fallback to the old-style manual
>> dependency management otherwise. Does that make sense?
>
> Um, yeah, except I don't think you really need to fall back, just  
> because people have other stuff installed.  The worst that's going  
> to happen is that you're going to force reinstallation of  
> dependencies they already have, just to get them into eggs.  (Or  
> make them create .egg-info directories to tell the system the stuff  
> is already installed.)



That's the problem. I'm trying to figure out a general plan of attack  
that Linux/BSD package maintainers can adopt for python packages that  
want to use eggs / setuptools. Here's some numbers on how many python  
packages are included in a couple different distributions:



On a OS X darwinports box:



     $ port list | grep -e '^py-' | wc -l
     237



On a Fedora 3 box (Core + Extras):



     $ yum list all | grep -e '^py' -e 'python' | wc -l
     78



I don't have any debian or gentoo boxes handy but I imagine they'd  
weigh in somewhere around the darwinports number.



None of these packages are currently provided as eggs or with .egg- 
info directories when they are installed to site-packages and they  
have complex dependency relationships that are managed by the  
distribution's utility (port, yum, apt-get, emerge, etc.) This  
creates a problem for these packages because it means that they can  
not assume dependencies will always be egg managed. If they start  
adding require() calls to their scripts, they will break under these  
environments. require() is an all or nothing proposition for  
distributions and that means there will need to be a planned  
"upgrade" period or something for all packages.



As a more specific example, I contribute to two packages that are  
distributed with Fedora Core: python-urlgrabber and yum. yum depends  
on python-urlgrabber and python-elementtree. Now, if I wanted to move  
yum to be egg based and use require(), I would also need to ensure  
that all yum's dependencies are egg based. When yum (and its  
dependencies) are installed from RPM, they must all be in egg format  
(or at least provide .egg-info dirs). If not, the yum script will fail.



So a single package using require() can cause a snowball effect where  
many other packages would need to be upgraded to egg format as well.  
In time, this may be a good thing because it could accelerate  
adoption of eggs but for the time being it makes it really hard to  
use require().



> Hm.  What if you created .egg directories and symlinked the  
> dependencies into them during your installation process?  Or if  
> there were some way to create the .egg-info directories  
> automatically from packaging system databases, or from inspecting  
> module contents?  Setuptools has some code that can look for the  
> setting of constants or presence of symbols in modules, without  
> importing them.  Perhaps I could extend this somehow so that a  
> transitional package like yours could include additional info in  
> the setup script, that checks for these dependencies and tags them  
> somehow?



Yes, yes. Along those lines.



> Or maybe this could be done by metadata -- you put a legacy.py file  
> in your egg-info, and when processing your egg's dependencies, if  
> pkg_resources can't find a package you need, it would call a  
> function in legacy.py that would check for the dependency using  
> setuptools' inspection facilities, and return a path and a guess at  
> a version number.
>
> How does that sound?



That would solve my problem perfectly.



>>  If that's the case, I might be
>> able to make my scripts as simple as::
>>
>>     from pkg_resources import require, find_distributions
>>     if list(find_distributions('MyPackage')):
>>
>
> Don't do this.  find_distributions() yields distributions found in  
> a directory or zipfile; it doesn't take a package name, it takes a  
> sys.path entry.



Ahhh..



> I'm not sure exactly what you're trying to do here.  If you just  
> want to know if your script is running from a development location  
> (and therefore needs to call require() to set up dependencies),  
> couldn't you just check for 'MyPackage.egg-info' in the sys.path[0]  
> (script) directory?
>
> e.g.:
>
>     import sys, os
>     if os.path.isdir(os.path.join(sys.path[0],"MyPackage.egg-info")):
>         from pkg_resources import require
>         require("MyPackage")  # ensures dependencies get processed
>
> If this is what you want, perhaps we can create a standard recipe  
> in pkg_resources, like maybe 'script_package("MyPackage")', that  
> only does the require if you're a development egg and not being run  
> from run_main().
>
>> The downside to this approach is that I would have to be sure to NOT
>> distribute MyPackage.egg-info with RPMs and other packages, which
>> kind of rules out any phased approach to bringing egg based packaging
>> to Fedora's stock RPMs.
>>
>
> I don't know if .egg-info is a good idea for RPMs.  .egg-info is  
> primarily intended for *development*, not deployment, because you  
> can't easily override a package installed with .egg-info in site- 
> packages.  In fact, the only way you can normally override it is to  
> install an egg alongside the script.
>
> My current idea for how RPMs and other packagers should install  
> eggs is just to dump them in site-packages as egg files or  
> directories, and let people use require() or else use EasyInstall  
> to set the active package.



Right. But that's going to require significant lobbying and effort to  
get all core python packages included in a distribution migrated over  
to eggs. This is the root of my dilemma. The only realistic approach  
seems to be supporting dual egg and non-egg deployment for a little  
while. The legacy.py proposal is one method and I threw out a couple  
earlier. But yea.. I think this is the root of the problem.



> Hey, wait a second...  if you can put install/uninstall scripts in  
> packages, couldn't installing or uninstalling an RPM ask  
> EasyInstall to fix up the easyinstall.pth file?  This would let  
> packagers distribute as eggs, but without breaking users'  
> expectations that the package would be available to "just import".   
> If somebody explicitly wants to support multiversion for a package,  
> they can run 'EasyInstall -m PackageName' to reset it to multi- 
> version after installing a new version.



I think that would be great but it assumes modification to a whole  
lot of existing packages. Ideally, I'd like to be able to start using  
eggs and require() in my packages without having to convince everyone  
else to do so just yet. Not that I won't be trying to convince  
people, I just don't want to rely on other packages being eggified in  
Fedora Core, darwinports, etc for upcoming releases.



> EasyInstall doesn't have everything that's needed to do this yet  
> (no "uninstall" mode), but perhaps it's a good option to add, and  
> then packagers could standardize on this approach.



I'd be happy to advocate to / work with packagers once we get a basic  
set of best practices together. It seems like there are a lot of  
options here - we just need to iron out the details.



Thanks,



Ryan Tomayko

... The Lake Project, published by the esteemed Nazraeli Press. This
oversized, 14"x14" book is comprised of surreal aerial images made in
the vicinity of Owens Lake, California.

#!/bin/bash



# Convert all original xstar type files into 
#  1. chunks with uniform chirp band
#  2. use sioseis to convert to standard sgy with 1 channel

rm -f sioseis_tmp[0-9]*



for orig in xstar/*.xstar; do
    file=${orig##xstar/}
    name=${file%%.xstar}
    #echo $orig $file
    echo ==========    $file   ==========
    echo
    # Do a scary hard link
    ln $orig .




    # Make sure we have only one frequence band per file.
    if [ ! -f ${name}-00001.xstar ]; then
	echo SLICING $name
	segyslice.py -d xstar -k -i $file SweepFreqStart SweepFreqEnd --verbosity=2 -x .xstar -o ${file%%.xstar}
    else
	echo no need to slice.  already done
    fi



    if [ ! -f ${name}-00001.sgy ]; then
	for xstar in ${name}*-[0-9][0-9][0-9][0-9]*.xstar; do
	    #echo ./xstar2segy.bash $xstar
	    ./xstar2segy.bash $xstar
	done
    fi



    if [ ! -f ${file}.db ]; then
	segysql.py -d xstar $file -f $file.db
    fi



    #
    # CLEANUP
    #
    rm -f $file  # Remove the hard link
    echo
    echo
done



for file in *-[0-9]*.xstar; do
    echo checking $file for db
    if [ ! -f $file.db ]; then
	segysql.py -d xstar $file -f $file.db
    fi
done



# Just in case sioseis barfed
ls -l sioseis_tmp*
rm -f sioseis_tmp[0-9]*

  The geolocator library provides a common interface to
  location-related data from various sources and provides means to
  calculate distances between places. Currently the MaxMind country &
  city geoip data sources are supported. There is also some
  preliminary support for using GNS (GEOnet Names Server )data. For
  distance calculation, the Haversine function is used. 



  The library can be useful for example to implement location-based
  features in different www frameworks; for an example of that, see
  the GeoLocation Zope / CMF / Plone package from which this library
  was spun off. Please note, however, that there is nothing
  www-specific in the library, nor does it have any external
  dependencies. 

segyslice.py -d xstar -k -i EEL5kL01.xstar \
  SweepFreqStart SweepFreqEnd \
  --verbosity=2 -M 1500  \
   -x .xstar -o EEL5kL01



FILE:  EEL5kL01-0001.xstar
SweepFreqStart 1000
SweepFreqEnd 6000
  Value change:  15000  ->  6000  for key  SweepFreqEnd
closing file:  EEL5kL01-0001.xstar
opening file:  EEL5kL01-1477.xstar

segysqldump.py -w 6 -d segy.db SweepLen SweepFreqStart SweepFreqEnd            



traceN   X        Y        Min      Shotpo   SweepL   SweepF   SweepF   



500      -44792   146998   36       66744    50       1000     6000    
1000     -44792   146997   37       66994    50       1000     6000
Error occured here 
1500     -44791   146996   38       67244    10       1000     15000   
2000     -44791   146995   39       67497    10       1000     15000   
2500     -44790   146994   40       67747    10       1000     15000   
3000     -44790   146994   41       67997    10       1000     15000   
3500     -44789   146993   42       68247    10       1000     15000   
4000     -44789   146992   43       68497    10       1000     15000   
4500     -44788   146991   44       68747    10       1000     15000   
5000     -44788   146990   45       68997    10       1000     15000   
5500     -44787   146989   47       69247    10       1000     15000   
6000     -44787   146988   48       69497    10       1000     15000   
6500     -44787   146987   49       69747    10       1000     15000   
7000     -44786   146986   50       69997    10       1000     15000   
7500     -44786   146986   51       70247    10       1000     15000   
8000     -44785   146985   52       70498    10       1000     15000   
8500     -44785   146984   53       70756    30       1000     15000   
9000     -44784   146983   54       71006    30       1000     15000   
9500     -44784   146982   55       71267    10       1000     15000   
10000    -44783   146981   56       71517    10       1000     15000   
10500    -44783   146980   57       71767    10       1000     15000   
11000    -44782   146979   58       72017    10       1000     15000   
11500    -44782   146979   59       72267    10       1000     15000   
12000    -44782   146978   0        72517    10       1000     15000   
12500    -44781   146977   1        72767    10       1000     15000   
13000    -44781   146976   2        73017    10       1000     15000   
13500    -44780   146975   3        73267    10       1000     15000   

segyslice.py -i EEL5kL01.xstar -o EEL5kL01-0000-1476.xstar --max=1476
segyslice.py -i EEL5kL01.xstar -o EEL5kL01-1477.xstar --min=1477

• segy.py: Added getEntireRawHdr, getRawTrace to allow duplicating segy
• segyslice.py: New command to chop traces. Need more functionality
• segysqldist.py: can now do UTM, skip zeros, fileKey ranges
• segysqldump.py: CoordTransform, variable width, any keys as args
• segysqlgmt.py: Don't allow scalercoord of 0. Turn zeros into ones.
• .py: lots more --verbosity all around

./xstar2segy.bash yr1999day217-2244z.xstar 
converting yr1999day217-2244z.xstar to yr1999day217-2244z.sgy
./sioseis.Darwin
  SIOSEIS ver 2005.4 (20 May. 2005)(C) Regents of U.C.
  
 procs diskin xstar header diskoa end 
  
 diskin 
     format edgetech 
     ipath yr1999day217-2244z.xstar end 
 end 
  
 xstar 
    type 2 dummies 2 end 
 end 
  
 header 
      fno 0 lno 999999 ftr 0 ltr 999999 
      l3 = l1 l6 = l3 end 
 end 
  
 diskoa 
   opath yr1999day217-2244z.sgy end 
 end 
  
 end 
 ****    0 ERRORS IN THIS JOB   ****
  Dropping trace 2 due to different sample intervals.
  impossible - icount =0 and its the same ping no
Done.

lsd yr1999day217-2244z.xstar  | head
      SHOT   TR       RP    TR ID  RANGE DELAY NSAMPS    SI   YR DAY HR MIN SEC
         0     0       0     0  1      0    90   3976    40 1999 217 22  51  38
         0     1       0     0  1      0    90   3976    40 1999 217 22  51  38
         0     1       0     0  1      0    90   3976    41 1999 217 22  51  38
         0     0       0     0  1      0    90   3976    41 1999 217 22  51  38
         0     0       0     0  1      0    90   3976    41 1999 217 22  51  38
         0     1       0     0  1      0    90   3976    41 1999 217 22  51  39
         0     1       0     0  1      0    90   3976    41 1999 217 22  51  39
         0     0       0     0  1      0    90   3976    41 1999 217 22  51  39
         0     1       0     0  1      0    90   3976    41 1999 217 22  51  39

#!/usr/bin/env python
import sqlite
cx = sqlite.connect('segy.db')
cu = cx.cursor()



cu.execute('select sql from sqlite_master where name = "segyTrace"')
for row in cu.fetchall():
    print row.keys()
    print row
    foo = row['sql'].find(' NMEAyear' )
    print foo
    # Exists, so I get 153 (a positive int)
    bar = row['sql'].find(' nonexistant ')
    print bar
    # Does not exist so -1

suwind < line43.ch2.su key=tracl min=4000 count=3000 > line43.ch2.4000-7000.su
sushift < line43.ch2.4000-7000.su tmin=0.510 tmax=0.866 > line43.ch2.4000-7000.nodelay.su



supsimage x2beg=5500 x2end=6350 d1num=0.01 x1beg=.61 x1end=.68 \
    hbox=10 wbox=18 verbose=1 d1=.000125  \
    <  line43.ch2.4000-7000.nodelay.su > ps2 && gv ps2

  supsimage verbose=1 d1=.000125 wbox=18 < line-4000-7000.nodelay.su  > ps
  gv ps

lrwxr-xr-x    1 schwehr  staff          23 Sep 14 11:19 line.su -> line43.ch1.4000-7000.su

# ~/x/src/density/endian
BIGENDIAN
# segyread tape=line43.sgy conv=1 endian=1 verbose=1 > line43.su



    103772752 Sep  9 13:34 line43.sgy
    201794624 Sep 13 14:18 line43.su

surange < line43.su



23932 traces:
 tracl=(1,11966)  tracr=(1,11966)  fldr=(1,11966)  tracf=(1,2)  trid=1
 sx=(364964,382927)  sy=(4519354,4524110)  gx=(364964,382927)  gy=(4519354,4524110)  delrt=(100,710) 
 ns=2048 dt=125 year=95 day=256 hour=(19,22) 
 minute=(0,59)  sec=(0,59)  timbas=2



suwind < line43.su key=tracf min=1 max=1 > line43.ch1.su
suwind < line43.su key=tracf min=2 max=2 > line43.ch2.su

segysqldist.py -X -124.50211667 -Y 40.8388500 -n 1 -u -z 10 --verbosity=3



BEST fields: fileKey, shotpoint, x, y, distance, EnsembleNo,FieldRecNo
BEST:  1 0 -124.502286138 40.8382818195 0.000592915087008 0 5861



CMP = 5861

suwind < line43.ch1.su key=tracl min=4000 count=3000 > line43.ch1.4000-7000.su



surange < line43.ch1.4000-7000.su 
3000 traces:
 tracl=(4000,6999)  tracr=(4000,6999)  fldr=(4000,6999)  tracf=1 trid=1
 sx=(370650,374952)  sy=(4521460,4522604)  gx=(370650,374952)  gy=(4521460,4522604)  delrt=(510,610) 
 ns=2048 dt=125 year=95 day=256 hour=(20,21) 
 minute=(0,59)  sec=(0,59)  timbas=2

# python
delrt=(510,610)
ns=2048
dt=125/100000.
tmax = delrt[1]/1000. + dt*ns
print tmax



  0.866



tmin = delrt/1000.
print tmin



  0.510



import sys; sys.exit ('see ya')



# sushift < line43.ch1.4000-7000.su tmin=0.510 tmax=0.866 > line43.ch1.4000-7000.nodelay.su 

# supswigb < line43.ch1.4000-7000.nodelay.su key=tracl nbpi=250 interp=1 xbox=1.0 ybox=1.2 hbox=9.0 wbox=428  > foo.ps
sh: line 1: /sw/lib/cwp-su/bin/bin/pswigb: No such file or directory



# export CWPROOT=/sw/lib/cwp-su
# !sups

sudo -s
cd /sw/lib/cwp-su/bin
ln -s psimage _psimage 

# supsimage < line43.ch1.4000-7000.nodelay.su  > foo.ps



# gv foo.ps

  The OpenFOAM (Open Field Operation and Manipulation) CFD Toolbox can
  simulate anything from complex fluid flows involving chemical
  reactions, turbulence and heat transfer, to solid dynamics,
  electromagnetics and the pricing of financial options. 



  The core technology of OpenFOAM is a flexible set of efficient C++
  modules. These are used to build a wealth of: solvers, to simulate
  specific problems in engineering mechanics; utilities, to perform
  pre- and post-processing tasks ranging from simple data
  manipulations to visualisation and mesh processing; libraries, to
  create toolboxes that are accessible to the solvers/utilities, such
  as libraries of physical models.

  TEMPE, Ariz.--(BUSINESS WIRE)--Sept. 12, 2005--Esser Design
  (www.esserdesign.com), a leading branding and marketing firm, and
  Terralever (www.terralever.com), a top web development and
  interactive marketing strategy company, today announced that the two
  companies, in partnership, have been selected by University of
  Arizona's Phoenix Mars Mission to redevelop the program's web
  site. The Phoenix Mission is operated for NASA by the Lunar and
  Planetary Laboratory at the University of Arizona, in partnership
  with the Jet Propulsion Laboratory, Lockheed Martin and the Canadian
  Space Agency. The Mission, scheduled for launch in August 2007 and
  arrival in May 2008, is the first program in NASA's Mars Scout
  Program of competitively selected missions. Scouts are innovative
  and relatively low-cost complements to the core missions of the
  agency's Mars Exploration Program. 
...

• sqlite3.info (3.2.5)
• pysqlite2-py.info (2.0.4)

suspike nt=$nt dt=$dt ntr=1 nspk=1 ix1=1 it1=20  | sugain scale=0.9 > spike1.data
echo -n '#' > spike1.xy
suascii < spike1.data >> spike1.xy



gnuplot
plot 'spike1.xy' with lp

#!/bin/bash



nt=251 
dt=0.001



suspike nt=$nt dt=$dt ntr=1 nspk=1 ix1=1 it1=20  | sugain scale=0.9 > spike1.data
echo -n '#' > spike1.xy
suascii < spike1.data >> spike1.xy



suspike nt=$nt dt=$dt ntr=1 nspk=1 ix1=1 it1=60  | sugain scale=0.6 > spike2.data
suspike nt=$nt dt=$dt ntr=1 nspk=1 ix1=1 it1=120 | sugain scale=0.4 > spike3.data
suspike nt=$nt dt=$dt ntr=1 nspk=1 ix1=1 it1=190 | sugain scale=0.5 > spike4.data



susum spike1.data spike2.data > a.tmp
susum spike3.data spike4.data > b.tmp
susum a.tmp b.tmp > spike.data
rm -f [ab].tmp spike?.data



echo -n '#' > spike.xy
suascii < spike.data >> spike.xy



suvibro dt=$dt tv=0.06 t1=0 t2=0.05 phz=-33 > sweep.data
suascii bare=1 < sweep.data > sweep.z



suconv < spike.data sufile=sweep.data | sushift tmin=0 tmax=0.251 > conv.data
suascii bare=1 < conv.data > conv.z



suxcor < conv.data sufile=conv.data | sushift tmin=0.251 tmax=0.502 > acor.data
suascii bare=1 < acor.data > acor.z



suxwigb < acor.data wbox=300 hbox=500 title="Autocorrelation" xbox=10 ybox=10

ChangeLog
• segy.py: Added a traceKey field to the segyTrace DB table. Unique ID
• segysql.py: More verbosity
• segysqlgmt.py: --fledermaus-pts
• segytextheader.py: --histogram to see diff between EBCDIC and ASCII

python
>>> a = "abc123"
>>> type (a)
<type 'str'>
>>> type (a[1])
<type 'str'>
>>> print a[1],ord(a[1])
b 98
>>> chr (98)
'b'

• segysqlgmt.py: added --fledermaus, rewrote internals. Much happier.
• segysqlgmt.py: UTM to LatLon. New CoordTransform class
• segytextheader.py: New file. Decode ASCII or EBCDIC text headers

Dear LazyWeb,



Does anyone know an easy way to distinguish (or at least guess) if a
block of data (here 3200 characters) is ASCII or EBCDIC.  I suppose I
could do a histogram of character code frequency and it should be
pretty obvious, but is there some simpler trick?



Please email me if you know something I don't (about this particular
problem)!



Thanks,
-kurt

  Thomas Benz wrote a Master's report on the use of Seismic Un*x in
  teaching of seismic processing. We are in the process of converting
  this document for general use, but have been asked by many people to
  make it available as soon as possible. To meet these requests, we
  have taken that report and simply converted it to Adobe Acrobat
  portable document format and are making it available at this
  site. We appreciate your comments on the material included. Seismic
  Un*x is a free software suite of programs that run under Unix on a
  wide range of platforms. It is provided by the Center for Wave
  Phenomena at the Colorado School of Mines, from their web site at
  http://www.cwp.mines.edu/cwpcodes/. 
...
  The main approach used in this report is straightforward: explain
  the principles of 2D seismic data processing and immediately use
  them in the processing of data, whether synthetic or real. The
  report starts with examples of running simple scripts, in exercises
  designed to provide some familiarity with time-series analysis. Then
  it creates a synthetic data set from a structurally-complex model
  (watch out -- this step can use a lot of CPU time), and then
  processes the shot gathers and CMP gathers from that synthetic data
  set. The report uses a proprietary data set made available to
  Michigan Tech for application of these procedures on real data. This
  data set is not currently available for release to others. 



  We are hoping to obtain an appropriate 2D data set for release
  through the Society of Exploration Geophysicists. Ideally, this data
  set would consist of a hundred or so shot points, with sufficient
  fold to provide good resolution (say, 32 at the minimum). It would
  have minimal statics correction required, and exhibit some
  structural feature of interest. 

  In accordance with the present invention, a method and apparatus for
  launching different application functions, based on the length of
  time an application button is pressed, are provided. The method and
  apparatus for launching the application function comprises:
  detecting that an application button has been pressed; determining
  the length of time the application button has been pressed; and
  launching alternative application functions based on the length of
  time the application button was pressed. If an application button is
  pressed for a short period of time, e.g., less than one second, an
  application function, such as displaying a view is launched. For
  example, if the application is a notes application, a short button
  press displays a list of summary information for the stored
  notes. If the same application button is pressed for a long period
  of time, e.g., at least one second, an alternative application
  function is launched. Using the notes example, a new note is opened
  for a long button press. Yet a third application function can be
  launched for a multi-press (e.g., double click) of an application
  button. Using the notes example, a new note is opened with a
  predefined template, for a multi-press of an application button.

The best parts are probably



* Edit VERSION to update to a new version (e.g. 1.9.1 or 1.10 or
  2.0... whatever your pref) 



* the top make files tries to help



    cd src
    make
    # It tells you all the options... hopefully the discriptions are
    # not too cryptic 

* Build with compiler optimizations so things go as fast as possible



     make clean; make all-opt



* Make a tar of the binaries for a platform... e.g. this will make
  pmag-1.9-tiger.tar.bz2 



    make bdist PLATFORM=tiger



* Make a release for fink... this removes all binaries, .o files, and
  fort.xx files and builds a tar.bz2 



   make pmag.info
   scp pmag-1.9.tar.bz2 sorcerer:/Library/WebServer/Documents/software/
   # then email me pmag.info

make



      Welcome to pmag 1.9
            by Lisa Tauxe



Options:



 make all         build everything
 make all-opt     build everything for SPEED
 make install     build everything and put in ../bin



 make clean       remove all build files
 make clean-dist  cleanup in prep for a dist
 make real-clean  nuke the dist tars and do clean/clean-dist



 make sdist       create a SOURCE only distribution
 make bdist       create a BINARY only distribution



 make pmag.info   create a fink info package and tar for Mac OSX

segytextheader.py  ~/Desktop/HuntecLine43/line43.sgy -e






C 1 USGS MARFAC                                                 CREW
NO        
C 2 LINE            AREA                        MAP ID                         
C 3 REEL NO           DAY-START OF REEL     YEAR      OBSERVER                 
C 4 INSTRUMENT: MUDSEIS NT     MODEL            SERIAL NO                      
C 5 DATA TRACES/RECORD        AUXILIARY TRACES/RECORD         CDP FOLD         
C 6 SAMPLE INTERVAL         SAMPLES/TRACE       BITS/IN
BYTES/SAMPLE      
C 7 RECORDING FORMAT        FORMAT THIS REEL        MEASUREMENT SYSTEM         
C 8 SAMPLE CODE:                                                               
C 9 GAIN TYPE:                                                                 
C10 FILTERS: ALIAS     HZ  NOTCH     HZ  BAND     -     HZ  SLOPE    -
DB/OCT
C11 SOURCE:                 NUMBER/POINT         POINT INTERVAL                
C12     PATTERN:                           LENGTH        WIDTH                 
C13 SWEEP:           HZ          HZ  LENGTH      MS  CHANNEL NO
TYPE       
C14 TAPER:                    MS                   MS  TYPE                    
C15 SPREAD:                                                                    
C16 GEOPHONES:                                                                 
C17     PATTERN:                                                               
C18 TRACES SORTED BY: RECORD                                                   
C19 AMPLITUDE RECOVERY:                                                        
C20 MAP PROJECTION                                                             
11
C21 PROCESSING:                                                                
C22 PROCESSING:                                                                
C23                                                                            
C24                                                                            
C25                                                                            
C26                                                                            
C27                                                                            
C28                                                                            
C29                                                                            
C30                                                                            
C31                                                                            
C32                                                                            
C33                                                                            
C34                                                                            
C35                                                                            
C36                                                                            
C37                                                                            
C38                                                                            
C39                                                                            
C40 END EBCDIC                                                                 

• xyzdespike.py: new program to remove easy to chop bad points
• makepltsegy.py: added --timescale to be the same as --timscal
• segysql.py: --verbosity=5 is like -v -v -v -v -v
• segysqlgmt.py: --verbosity, --force-hemisphere

  Mapnik is an toolkit for developing GIS applications. At the core is
  a C++ shared library providing algorithms/patterns for spatial data
  access and visualization. Essentially a collection of geographic
  objects (map, layer, datasource, feature, and geometry), the library
  doesn't rely on "windowing systems" and can be deployed in any
  server environment. It is intended to play fair in a multi-threaded
  environment and is aimed primarily, but not exclusively, at
  Web-based development. High-level Python bindings (boost.python)
  facilitate rapid application development, targeting zope3, django,
  etc.

Iowa City, Iowa  (Hi Geb!)
El Dorado, Arkansas
Barcelona, Cataluna (I really enjoyed visiting Barcelona)
Winfield, Alabama
Copenhagen, Kobenhavn

Open House in Celebration of WHOI's 75th Anniversary



Saturday, September 10, 2005  . 10 a.m. to 4 p.m.  .  Woods Hole
Oceanographic Institution (rain or shine) 






The Woods Hole Oceanographic Institution (WHOI) is celebrating 75
years of ocean research, education, and exploration in 2005. Founded
in 1930 to help advance the role of the United States in worldwide
oceanographic research, WHOI has grown from a summer laboratory to a
thriving, year-round, major research operation. Today the Institution
has a staff and student population of nearly a thousand people working
and studying on two campuses in Woods Hole, Massachusetts, on Cape Cod
and at research locations around the world. It operates three large
research vessels, several smaller boats, the submersible Alvin, and a
fleet of deep-diving remotely operated and autonomous vehicles. 



This will be WHOI's first open house in 25 years. Both the Woods Hole
and Quissett campuses will be included, with shuttle buses running
between them. Parking will be available on the Quissett Campus and in
the village lots of various Woods Hole research institutions. 



Open house features include robotic vehicles, the coastal research
vessel Tioga, scientific exhibits and demonstrations, a just-published
book on the institution's history, and 20 talks on subjects ranging
from right whales and red tides to coastal erosion and archaeological
mapping in Greek waters. An incomplete list of open house components
(more are being added) appears on the web site included below. 



For more information on the Open House and other 75th anniversary
events, visit http://www.whoi.edu/75th/index.html.

  Some of the highest quality images ever taken of the Earth's lower
  crust reveal that the upper and lower crust form in two distinctly
  different ways. 



  "This new way of studying the ocean crust is the equivalent of a new
  telescope in astronomy," said Bruce Malfait, head of the National
  Science Foundation (NSF)'s marine geosciences section, which funded
  the research. "It allows us to look at Earth processes and
  composition at a remarkable new level of detail." 



  To form the images, the researchers, led by a team from Columbia
  University's Lamont-Doherty Earth Observatory (LDEO), analyzed sound
  waves bounced off structures deep in the Earth, a process similar to
  creating an ultrasound image. They published the results of their
  work in the Aug. 25 issue of the journal Nature. 



  The Earth's oceanic crust is formed from molten rock, or magma,
  located beneath mid-ocean ridges. Magma chambers located in the
  mid-crust are known to be responsible for generating the
  upper-oceanic crust. But it has so far been unclear if the lower
  crust is formed from the same magma source, or if it arises
  primarily from "magma lenses"--smaller bodies located at or near the
  base of the crust. 



  The resulting images are the first of their kind to show solidified
  lenses and sills--narrow lateral intrusions of magma--embedded
  within the crust-mantle boundary known as the Moho transition
  zone. The presence of such lenses and sills deep near a mid-ocean
  ridge suggests the lower crust is at least partially formed from
  several smaller sources of magma rather than from a single large
  source located in the middle of the crust. 



  "This demonstrates that the process of crustal formation is more
  complex than believed," said LDEO scientist and lead author of the
  study, Mladen Nedimovic. "It also favors the emerging view that
  volcanoes have a complicated plumbing system consisting of many
  interconnected sills and magma conduits. We still have a lot to
  learn about what goes on beneath the surface of the Earth." 



  Researchers from LDEO as well as the Scripps Institution of
  Oceanography and the Woods Hole Oceanographic Institution collected
  seismic data from the oceanic crust beneath the sea floor off the
  coast of Wash., Ore., and northern Calif. 

Info2: <<



Package: pyproj-py%type_pkg[python]
Version: 1.1
Revision: 1
Description: Projection coversions in python
DescDetail: <<
<<
#http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/python/pyproj.tar.gz
Source: http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/python/pyproj-%v.tar.gz
Source-MD5: f13b2e26389de92a79e5f0efa7ecd0c4
Type: python (2.4)
Depends: python%type_pkg[python], proj
BuildDepend:
CompileScript: <<
 #!/bin/bash
 export PROJ_DIR=%p
 %p/bin/python%type_raw[python] setup.py build
<<
InstallScript: <<
 %p/bin/python%type_raw[python] setup.py install
<<
DocFiles: README
License: BSD
Homepage: http://cheeseshop.python.org/pypi/pyproj/1.0
Maintainer: Kurt Schwehr <goatbar _at_ users.sourceforge.net>
# Info2
<<

#!/usr/bin/env python



from pyproj import Proj



params = {}
params['proj'] = 'utm'
params['zone'] = 10
p = Proj(params)
print p(-120.108, 34.36116666)
#765975.64       3805993.13 0

#include <fenv.h>
#include <stdio.h>



#define FALSE (0==1)
#define TRUE (0==0)
void check_f_except__ (int *result) {
  int found=FALSE;



  if (fetestexcept(FE_INEXACT)) {
    found=TRUE;
    printf("underflow\n");
    feclearexcept(FE_INEXACT);
  }
  if (fetestexcept(FE_DIVBYZERO)) {
    found=TRUE;
    printf("divbyzero\n");
    feclearexcept(FE_DIVBYZERO);
  }
  if (fetestexcept(FE_UNDERFLOW)) {
    found=TRUE;
    printf("underflow\n");
    feclearexcept(FE_UNDERFLOW);
  }
  if (fetestexcept(FE_OVERFLOW)) {
    found=TRUE;
    printf("overflow\n");
    feclearexcept(FE_OVERFLOW);
  }
  if (fetestexcept(FE_INVALID)) {
    found=TRUE;
    printf("invalid\n");
    feclearexcept(FE_INVALID);
  }
  if (fetestexcept(FE_ALL_EXCEPT)) {
    found=TRUE;
    printf("all\n");
    feclearexcept(FE_ALL_EXCEPT);
  }



  *result = found;
}

      double precision d



      i = 999
      call check_f_except(i)
      write (*,*) 'result = ',i



      d = (1.0d+00 - 1.0d+00)
      i = 999
      call check_f_except(i)
      write (*,*) 'result = ',i



      end

underflow
divbyzero
invalid
 result =  1
 result =  0

  It predicted that 200,000 people or more would be unwilling or
  unable to heed evacuation orders and that thousands would die. It
  also predicted that people would be housed in the Superdome, that
  aid workers would find it difficult to gain access to the city as
  roads became impassable, as well as many other of the consequences
  that actually unfolded after Katrina hit this week... 

1. Darwin
2. DarwinSource
3. 10.4.2
4. gdb-384.tar.gz

2005-01-19  Jason Molenda  (jmolenda@apple.com)



	* Makefile (GDB_VERSION): Update to reflect the actual date of the
	last FSF merge.

1. mmalloc
2. readline
3. intl
4. libiberty
5. bfd
6. opcodes
7. binutils
8. gdb

./configure --prefix=${HOME}/darwincvs 



...



long double support in compiler... yes
checking for long double support in printf... yes
checking for long double support in scanf... yes
checking whether  has __NR_tkill... no
checking compiler warning flags...  -Wimplicit -Wreturn-type -Wcomment -Wtrigraphs -Wformat -Wparentheses -Wpointer-arith -Wformat-nonliteral -Wunused-label -Wunused-function -Wuninitialized
checking for cygwin... no
checking for X... libraries /usr/X11R6/lib, headers /usr/X11R6/include
configure: error: "*** Gdb does not support native target powerpc-apple-darwin7.9.0"

powerpc-apple-macosx-gdb sioseis
GNU gdb 6.3.50.20050905
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB.  Type "show warranty" for details.
This GDB was configured as "--host= --target=powerpc-apple-macosx"..."/Users/schwehr/projects/xcore/src/sioseis/sioseis": not in executable format: File format not recognized

  Avahi is a system which facilitates service discovery on a local
  network -- this means that you can plug your laptop or computer into
  a network and instantly be able to view other people who you can
  chat with, find printers to print to or find files being
  shared. This kind of technology is already found in MacOS X (branded
  'Rendezvous', 'Bonjour' and sometimes 'ZeroConf') and is very
  convenient.



  Avahi implements the DNS-SD over Multicast DNS. It uses D-BUS to
  communicate from the application to the system daemon. A system
  daemon is used to coordinate application efforts in caching replies,
  necessary to reduce noise on networks generated by protocols such as
  MDNS. The goals of this project are for implementation into the
  GNOME desktop, however that does not stop it from being used by any
  other applications or desktop environments as long as you have a
  working installation of dbus, which most modern distributions now
  come with by default. 



  Although other implementations exist (such as Howl), none exist that
  use fully open and free GPL-compatible code. The howl library is
  hindered by the APSL among other things, which makes it undesirable
  to some people, especially distributions such as Debian which has
  strict free software guidelines.

gcc -g -O0 -Wall -Wimplicit -pedantic -W -Wstrict-prototypes \
    -Wredundant-decls -Dosx=1   -c -o shifts.o shifts.c
f77 -g -O2 -Wall   -c -o contro.o contro.f

gdb sioseis 
run < yr1999day217-2244z.xstar

bt
#0  0x000eccb0 in s_copy ()
#1  0x0000a514 in getoke_ (cbufo=0xbffff4d0, nchars=0xbffff564, __g77_length_cbufo=1) at getoke.f:94
#2  0x0000610c in MAIN__ () at contro.f:352
#3  0x000e9dec in main ()
#4  0x0000236c in _start (argc=-1073744540, argv=0x3c6c0a0, envp=0x3c6bfcf) at /SourceCache/Csu/Csu-47/crt.c:267
#5  0x8fe1a278 in __dyld__dyld_start ()



break main
Breakpoint 1 at 0xe9dd0
(gdb) r
The program being debugged has been started already.
Start it from the beginning? (y or n) y
Starting program: /Users/schwehr/projects/xcore/src/sioseis/sioseis < yr1999day217-2244z.xstar



Breakpoint 1, 0x000e9dd0 in main ()
(gdb) list
1             SUBROUTINE GETOKE(CBUFO,NCHARS)
2       c  ****   Watch out - getoke writes a null after the last char,
3       c         so cbufo must be bigger than the number of characters
4       C
5       C   GETOKE RETURNS CONSECUTIVE TOKENS (ITEMS BETWEEN A DELIMITER), ONE PER CALL,

./xstar2segy.bash yr1999day217-2244z.xstar 
converting yr1999day217-2244z.xstar to yr1999day217-2244z.sgy
  SIOSEIS ver 2005.4 (20 May. 2005)(C) Regents of U.C. kds                       
  
 procs diskin xstar header diskoa end 
 !procs diskin  header diskoa end 
  
 diskin 
     format edgetech 
     ipath yr1999day217-2244z.xstar end 
 end 
  
 xstar 
    type 2 dummies 2 end 
 end 
  
 header 
      fno 0 lno 999999 ftr 0 ltr 999999 
      l3 = l1 l6 = l3 end 
 end 
  
 diskoa 
   opath yr1999day217-2244z.sgy end 
 end 
  
 end 
 ****    0 ERRORS IN THIS JOB   ****
  Dropping trace 2 due to different sample intervals.
  impossible - icount =0 and its the same ping no
Done.

(gdb) run < zscript.sioseis   
(bt) bt
No stack.

nm xstarex.o | grep xstarex
00000000 t ___g77_masterfun_xstarex
00002844 T _xstarex_



(gdb) break xstarex_

nm core_logger.o | grep T



00001254 T _esc_help
00000144 T _main
00000004 C _origTerm
00000c48 T _process_io
00000044 T _stablk

(gdb) run
(gdb) bt
#0  0x000b9600 in xstarex_ (buf=0x2877f50, lbuf=0x2877f50, ibuf=0x2877f50, nready=0x1398d0) at xstarex.f:1
#1  0x00007b68 in MAIN__ () at contro.f:1542
#2  0x000e9dec in main ()
#3  0x0000236c in _start (argc=1072604, argv=0x4, envp=0x158) at /SourceCache/Csu/Csu-47/crt.c:267
#4  0x8fe1a278 in __dyld__dyld_start ()



(gdb) up
(gdb) list
1537	          CALL getxstartr(buf(index1),buf(index1),buf(index1),nrdyxstar)
1538	          nxstar = nxstar -1
1539	          iout(i) = iout(i) -1
1540	          GOTO 1010
1541	      ELSE
1542	          CALL xstarex( buf(index1), buf(index1), buf(index1),
1543	     &         nrdyxstar)
1544	          iout(i) = nrdyxstar + 1
1545	          nxstar = nrdyxstar
1546	          IF( nrdyxstar .EQ. 0 ) GOTO 1015

  If there isn't really an ap on the system and we are going to use an
  ap simulator, which uses an internal array specified in contro, make
  it very large!

(gdb) info locals
deglat = 0
deglong = 0
dtemp = 0
dlong = 0
dlat = 0
deltad = 0
dbuf = (0)
mlat = 0
mlong = 0
distance = 99999
deglatlast = 0
deglonglast = 0
dbub = (0, 0, 0)
depth_lat__ = 0
depth_long__ = 0
boatlat = 0
boatlong = 0
deltalat = 0
deltalong = 0
itunit = Cannot access memory at address 0x48

      CALL inap( buf(numhdr+1), numdat )
c****
c****  Treat Triton totally separately from Edgetech
c****
      write (*,*) 'type: ',type
      IF( type .EQ. 3 ) THEN






(gdb) n
 type:  2
(gdb) print *type
Cannot access memory at address 0x34
(gdb) print *type
Cannot access memory at address 0x34

gdb --version
GNU gdb 5.3-20030128 (Apple version gdb-309) (Thu Dec  4 15:41:30 GMT 2003)

  A glitch has forced Nasa's Mars Global Surveyor (MGS) spacecraft to
  shut down its science instruments. 



  The spacecraft has switched into a "safe mode", in which the
  instruments and some other systems are turned off. 



  Team members are racing to get the probe out of this mode so it can
  photograph the presumed crash site of a US Mars mission lost in
  1999. 



  The pictures could decide whether design changes are needed on a
  probe due to launch for Mars in 2007. 



  Next week will provide the last opportunity to take images of the
  crash site of Mars Polar Lander (MPL) for another two years. 

  BibDesk is a graphical BibTeX-bibliography manager for Mac OS
  X. BibDesk is designed to help organize and use bibliographic
  databases in BibTeX .bib format. In addition to manual typing,
  BibDesk lets you drag & drop or cut & paste .bib files into the
  bibliographic database and automatically opens files downloaded from
  PubMed. BibDesk also keeps track of electronic copies of literature
  on your computer and allows for searching your database through
  several keys. 

  An ENTRY statement permits a procedure reference to begin with a
  particular executable statement within the function or subroutine
  subprogram in which the ENTRY statement appears. It may appear
  anywhere within a function subprogram after the FUNCTION statement
  or within a subroutine subprogram after the SUBROUTINE statement,
  except that an ENTRY statement must not appear between a block IF
  statement and its corresponding END IF statement, or between a DO
  statement and the terminal statement of its DO-loop. 

  Kurt is a machine with four injection needles and Chinese
  ink. Attached to the upper arm it tattoos an enduring symbol on the
  user's skin. Kurt is an extreme example of the human-machine
  interface. The intervention is aimed at rising the awareness of the
  close connection between the individual and the machine, in a bond
  as long-lasting as the drawing itself.

  Katrina Information Map - This map is intended for the use of people
  affected by Hurricane Katrina who have or are trying to find
  information about the status of specific locations affected by the
  storm and its aftermath. If you have information about the status of
  an area that is not yet on the map, please contribute by following
  the instructions below so that others may get that much needed
  information. 

NOS/OCS SITREP 1230 Sept 2



Subject: NOS/OCS SITREP 1230 Sept 2



1.  NRT1 has completed surveys Pensacola.  Based upon NOAA surveys, USCG
is opening Pensacola to 31 foot draft commercial and 34 foot draft to
Naval Air Station.  A gasoline tanker is due Sunday, can load deeper
based upon NOAA surveys.



2.  NRT2 is side scanning Gulfport to enable USCG cutter to enter on law
enforcement mission.  Initial report basin is clear with some ruins
attached to bulkhead.  At 1300 making first side scan sweep from basin
out channel to GIWW.



3. NRT4 is conducting survey operations Port Fourchon.  A Washington
Post reporter is scheduled to ride with them tomorrow.



4. NRT2 side scan of Pascagoula yesterday discovered 6 obstructions.
Information was relayed to USACE.  They are planning to put divers down
to investigate.  NRT2 will complete side scan of Pascagoula tomorrow.
Pascagoula was opened to 12 foot draft today by USCG.  Gasoline
inventory reported there that can be removed by barge.  Aids to
navigation need to be restored to open to deeper draft.



5.  NRT1 will side scan Biloxi tomorrow.  Power plant needs coal.  Power
plant projects they could be back online in a week.  Also report of
leaking diesel tank, spill reported contained, needs to be pumped out.



6. NANCY FOSTER ETA Mobile Monday Morning to survey outer approaches.



7. THOMAS JEFFERSON  ETA depart Norfolk tomorrow 1700.



8. DAVIDSON proceeding Galveston for resupply.

Subj: Aerial recon flight of storm damaged area... 
Date: 8/31/2005 1:43:00 AM Eastern Standard Time 



Flight, 
This was forwarded to the Vietnam Helicopter Flight Crewmember Network. 
Pretty amazing report. 
Rick Lester 
JOKER 94 



I just got back from a 3.5 hour flight over the Mississippi Gulf Coast 
from Bay St. Louis to Pascagoula.  Please pass this report along to 
others who have evacuated, and need to know. 



I have no words--I'm not sure there ARE words--for what I just saw.  I 
am in shock.  The destruction of the Gulf Coast was of a level that 
Hurricane Camille only began to begin to approach. 



Personal.  Janet's and my house is gone.  An oak tree from the east side 
is lying on top of what looks like a section of the roof.  Whatever else 
is left is in a pile of debris at the CSX railroad tracks to the north, 
along with every other house in the neighborhood, save the one house 
next door (Archie & Margaret's). Theirs is a spanish-southwestern style 
built in 1925, and it used to stand out against all the other homes. 
Now it IS the only other home. The windows are all out, and it looks 
pretty bad.  Grasslawn, an antibellum home on the other side, is gone. 



Going down the coast from east to west... 



Ocean Springs has heavy damage along the beach.  I don't know how far 
back it goes. 



Biloxi.  Point Cadet is gone.  The Hwy 90 bridge is gone.  The old hwy 
bridge is gone.  The railroad bridge is gone.  All spans are down--every 
last one of them--except for the drawbridge in the middle. The Biloxi 
Casinos are mostly standing, but the barge from the Grand is north of 
hwy 90 now.  Mary Mahoney's Restaurant is standing, but most of the 
surrounding buildings are gone. 



The Seashore Manor Methodist retirement home is the only structure 
standing on that stretch of beach.  The windows are even unbroken.  The 
Broadwater Marina is gone.  What's left of the President Casino barge is 
north and 
west, and across Hwy 90.   Edgewater Mall is standing, but it has large 
holes in the structure, high, all along the south side.  Edgewater 
Village strip mall is gone.  Everything south of the railroad tracks has 
severe damage, if it's still standing.  I didn't see the Back Bay area. 



Gulfport.  The Legacy Condominiums are standing, and apparently did 
quite well.  Nearly everything around them is gone.  The VA hospital 
complex did well.  Most of the homes on either side of it are gone.  The 
water went beyond the railroad tracks in multiple places near there. 
There is a continuous debris pile pushed up against the tracks all along 
what used to be Second Street.  Almost every house there is gone. 
Downtown is mostly standing, but there has been water up to the tracks. 
The roofs of most of the buildings have visible damage.  The Hancock 
Bank building and the MS Power Co building seem to have fared the best. 
The Baptist and Presbyterian churches looked good.  I couldn't see the 
GPT Methodist church. The new courthouse appears to be undamaged, 
including the windows.  The harbor is mostly gone.  The roof of Marine 
Life is down.  The shipping terminal warehouses are there, but they have 
debris on the roofs, indicating that they were underwater.  The Grand 
Casino hotel (and Oasis) did very well, but the casino barge is blocking 
Hwy 90 to the west. The homes in West Gulfport did better than on the 
east side, but only for a space.  Approaching Long Beach it gets worse 
quickly. 



Keesler AFB and the Gulfport CBC.  One hangar at Keesler looks bad.  The 
rest of those structures look intact from where we were flying. 



Long Beach.  Everything from about 1000' south of the RR tracks is gone. 
The Baptist Church is standing, but gutted.  St. Thomas has parts 
standing, but it is mostly gone.  USM Gulf park looks pretty good, and 
it is the only structure on that section of Beach.  The Methodist Church 
is standing, and the New Life Center looks good.  The roofs look good. 
Water came up as high as the house next door to the Associate Pastor's 
parsonage.  Everything 
south of Magnolia Street is gone.   North of Magnolia Street, there is 
heavy 
damage up to about 1000' of the tracks.  On my old street of Buena Vista 
Dr., my old house is the first one standing.  It appears from the air to 
have had water up to the second story.  The High School, Middle School, 
and Reeves Middle School appear undamaged.  Harper McCaughan is hurt, 
but I couldn't see how bad.  It probably had water, and it is about at 
where the survival line forms.  North LB did much better, but there was 
lots of roof damage.  Pecan Park seems to have fared the best.  The Wolf 
River is flooded so badly I couldn't see any structures above water.  It 
is up to within a couple of feet of the I-10 bridge--enough that I worry 
for the integrity of the bridge.  The harbor is gone. Every structure in 
the vicinity of the harbor is gone. 



Pass Christian.  The Wal Mart is standing, but it has been  underwater. 
Some downtown buildings are standing.  St Anne's is there, but it has 
had water  over it.  From the LB line westward, everything south of 
Second Street is gone.  The only exception was Scenic Drive, where about 
half of the houses are standing.  From Downtown west I saw nothing 
standing anywhere in Pass Christian.  The only houses I saw standing in 
West Pass were the newer Ones on the North Side, near the Bay of St. 
Louis, and those two new millionaire  acres-looking houses on the beach 
about halfway to Pass Marianne.  The DuPont site is there, but we didn't 
fly north far enough to assess. 



The Hwy 90 bridge over the Bay of St. Louis is down.  Every span.  The 
railroad bridge is down.  Every span. 



Bay St. Louis appears to have had water in the downtown area.  Most of 
the buildings are standing, but most have roof damage.   All the homes 
southwest of the downtown area, along the beach road, are gone.  The 
homes west of the beach road had the heaviest roof damage I saw anywhere 
except in East Biloxi. 



Waveland.  We didn't fly all the way to Waveland, but I could see down 
the beach road quite a ways, and I saw nothing standing. 



Other.  Hwy 90 is washed out in sections all along the coast, from
Biloxi to Bay St. Louis.  Huge sand drifts obscure many more sections
of pavement.  There is a lot more natural sand along the beach, maybe
another 500-1000'.  Ship Island is mostly gone.  There is a clump of
trees on a sand bar on the east, and there is a sand bar on the
west. Fort Massachusetts is there, but has sand berms about 15' high
along the southwest side.  Cat Island has a new shape.  It doesn't
have "fingers" any more.  It is just a continuous stretch, like Ship
Island used to be, only with trees.



My estimates.  Based on the debris lines and what I know of the 
geography of the MS Gulf Coast, the tides were indeed in the 30-foot 
range, all the way from Bay St. Louis to Biloxi.  Based on the damage I 
saw to roofs, and to structures unaffected by water, the winds would 
have had to be greater than 100 knots, possibly as great as 140 knots. 
I couldn't see looting/looters, but I could see the opportunity for it. 
The casualty numbers, when they come in, will be staggering. 



In closing.  I have never seen any manner of destruction like this.  It 
is bad in Biloxi and Gulfport...as bad as anything I saw after Hurricane 
Camille in 1969.  The devastation of Pass Christian is indescribable. 
The news reporting only scratches the surface.  In my opinion, based on 
my considerable life experience, the federal government should declare 
martial law, and rush every possible measure of aid to this area. 



Janet and I are praying for the survivors on the ground there.  Please 
join us.  I will report more later.