Kurt's Weblog

Version 0.18 - 31 Aug 2005
  * decimate.py: added preserve-begin and preserve-end options
  * segy.py: added verbosity
  * segydump.py: uses segy.py verbosity
  * segysql.py: added --warn-only for some types of bad segy files
  * segysqlgmt.py: --skip-line to avoid write points along each line
  * segysqlgmt.py: guess what --font-size does 
  * segysqlgmt.py: --key allows lines to go off something other than shotpoint
  * segysqlgmt.py: --coord-unit to force any segy coord type number like arc sec
  * segysqlgmt.py: Dividing by 3600.0 is not totally hard coded any more
  * segysqlgmt.py: this file is getting ugly
  * segysqlgmt.py.help2man: added this file to start better documentation

  About my Idea Generation Methods website
  By Martin Leith	



  This website lists and explains every idea generation method I've
  encountered during the past 15 years. It is the result of extensive
  research; my many sources include books, management journals,
  websites, academics, consultants and colleagues. 



  The methods have been drawn not just from the worlds of creative
  problem solving and innovation, but also from other worlds such as
  organisational change, strategic planning, psychotherapy, the new
  sciences and the creative arts. 



  The methods are listed below. Each is linked to a description, and
  in some cases you will find full instructions for using the method
  to generate ideas.

\documentclass[agupp]{aguplus}
\usepackage{graphicx}

\documentclass[agums]{aguplus} % For review - double spaced
\usepackage{graphicx}

\begin{landscape}
\begin{figure}
\includegraphics[width=6in]{FiguresCores/Core4-AMS-water-v5} % FIX: only for draft
%\includegraphics{FiguresCores/Core4-AMS-water-v5}
\caption{Core 4 is blah blah blah}
\label{fig:core4}
\end{figure}
\end{landscape}

  OPeNDAP/DODS is a framework that simplifies all aspects of
  scientific data networking, making local data accessible to remote
  locations regardless of local storage format. 



  This Python module offers both a client and a server
  implementation. The client allows transparent and efficient access
  to datasets available on OPeNDAP/DODS server through the internet,
  using a syntax that follows very close the one used by Scientific
  Python to access NetCDF files. 



  The module also comes with three OPeNDAP/DODS servers: a CGI script,
  a Twisted-based server that accepts multiple connections, and a very
  simple server that only allows single connections (mostly used for
  testing purposes). All share the same architecture, and can be
  extended through plugins to serve different types of data. 

  Author's Note: This column was originally intended to be the final
  disaster in the "Disasters Waiting to Happen" series. As I was
  developing the hypothetical situation depicting a devastating
  hurricane striking New Orleans, Louisiana, the disaster waiting to
  happen threatened to become a reality: Hurricane Ivan, a category 4
  hurricane (with 140 mph winds) fluctuating to a category 5 (up to
  155 mph winds), was slowly moving directly toward New
  Orleans. Forecasters were predicting a one-in-four chance that Ivan
  would remain on this direct path and would be an "extreme storm" at
  landfall. In reality, the storm veered to the north and made
  landfall east of Mobile Bay, Alabama, causing devastation and
  destruction well into the central Gulf shoreline and throughout the
  Southeast and the Mid-Atlantic states.



What if Ivan Had Hit New Orleans?



  New Orleans was spared, this time, but had it not been, Hurricane
  Ivan would have: 



    Pushed a 17-foot storm surge into Lake Pontchartrain;



    Caused the levees between the lake and the city to overtop and
    fill the city "bowl" with water from lake levee to river levee, in
    some places as deep as 20 feet; 



    Flooded the north shore suburbs of Lake Pontchartrain with waters
    pushing as much as seven miles inland; and 



    Inundated inhabited areas south of the Mississippi River.



  Up to 80 percent of the structures in these flooded areas would have
  been severely damaged from wind and water. The potential for such
  extensive flooding and the resulting damage is the result of a levee
  system that is unable to keep up with the increasing flood threats
  from a rapidly eroding coastline and thus unable to protect the
  ever-subsiding landscape. 



Evacuation Challenges



  Researchers have estimated that prior to a "big one," approximately
  700,000 residents of the greater New Orleans area (out of 1.2
  million) would evacuate. In the case of Hurricane Ivan, officials
  estimate that up to 600,000 evacuated from metropolitan New Orleans
  between daybreak on Monday, September 13 and noon on Wednesday,
  September 15, when the storm turned and major roads finally started
  to clear. 



  To aid in the evacuation, transportation officials instituted
  contraflow evacuation for the first time in the area's history
  whereby both lanes of a 12-mile stretch of Interstate 10 were used
  to facilitate the significantly increased outbound flow of traffic
  toward the northwest and Baton Rouge. The distance of the contraflow
  was limited due to state police concerns about the need for staff to
  close the exits. And, although officials were initially pleased with
  the results, evacuees felt the short distance merely shifted the
  location of the major jams. 



  These feelings were justified by the amount of time it took
  residents to evacuate-up to 11 hours to go the distance usually
  traveled in less than 1.5. For many who evacuated into Texas, total
  evacuation time frequently exceeded 20 hours. Since the storm, a
  consensus has developed that to alleviate this congestion much more
  secondary highway coordination is necessary throughout the state,
  contraflow needs to be considered for much greater distances,
  residents who are able and willing to evacuate early must be doubly
  encouraged to do so, families with multiple cars need to be
  discouraged from taking more than one unless they are needed to
  accommodate evacuees, and all modes of transportation in their
  various configurations must be fully considered for the
  contributions they can make to a safe and effective evacuation. 



  The major challenge to evacuation is the extremely limited number of
  evacuation routes, which is the result of the same topography and
  hydrology responsible for the area's high level of hurricane
  risk. The presence of the Mississippi River, several lakes and bays,
  and associated marshes and swamps necessitates very expensive
  roadway construction techniques that are generally destructive to
  the environment, making the addition of more arteries increasingly
  challenging. This problem of limited evacuation routes also plagues
  the rest of the delta plain of southeast and south central
  Louisiana. 






  The fact that 600,000 residents evacuated means an equal number did
  not. Recent evacuation surveys show that two thirds of nonevacuees
  with the means to evacuate chose not to leave because they felt safe
  in their homes. Other nonevacuees with means relied on a cultural
  tradition of not leaving or were discouraged by negative experiences
  with past evacuations. 



  For those without means, the medically challenged, residents without
  personal transportation, and the homeless, evacuation requires
  significant assistance. The medically challenged often rely on life
  support equipment and are in such fragile states of health that they
  can only be moved short distances to medically equipped
  shelters. While a large storm-resistant structure with appropriate
  equipment has yet to be constructed or retrofitted, the Superdome
  was used to shelter nonevacuees during Ivan. 



  Residents who did not have personal transportation were unable to
  evacuate even if they wanted to. Approximately 120,000 residents
  (51,000 housing units x 2.4 persons/unit) do not have cars. A
  proposal made after the evacuation for Hurricane Georges to use
  public transit buses to assist in their evacuation out of the city
  was not implemented for Ivan. If Ivan had struck New Orleans
  directly it is estimated that 40-60,000 residents of the area would
  have perished. 



  Unwilling to merely accept this reality, emergency managers and
  representatives of nongovernmental disaster organizations, local
  universities, and faith based organizations have formed a working
  group to engage additional faith-based organizations in developing
  ride-sharing programs between congregation members with cars and
  those without. In the wake of Ivan's near miss, this faith-based
  initiative has become a catalyst in the movement to make evacuation
  assistance for marginalized groups (those without means of
  evacuation) a top priority for all levels of government. 



To the Rescue



  If a hurricane of a magnitude similar to Ivan does strike New
  Orleans, the challenges surrounding rescue efforts for those who
  have not evacuated will be different from other coastal
  areas. Rescue teams would have to don special breathing equipment to
  protect themselves from floodwaters contaminated with chemicals and
  toxins released from commercial sources within the city and the
  petrochemical plants that dot the river's edge. Additionally, tank
  cars carrying hazardous materials, which constantly pass through the
  city, would likely be damaged, leaking their contents into the
  floodwater and adding to the "brew." The floodwater could become so
  polluted that the Environmental Protection Agency might consider it
  to be hazardous waste and prohibit it from being pumped out of the
  leveed areas into the lake and marshes until treated. 



  Regional and national rescue resources would have to respond as
  rapidly as possible and would require augmentation by local private
  vessels (assuming some survived). And, even with this help, federal
  and state governments have estimated that it would take 10 days to
  rescue all those stranded within the city. No shelters within the
  city would be free of risk from rising water. Because of this
  threat, the American Red Cross will not open shelters in New Orleans
  during hurricanes greater than category 2; staffing them would put
  employees and volunteers at risk. For Ivan, only the Superdome was
  made available as a refuge of last resort for the medically
  challenged and the homeless. 



The Aftermath



  In this hypothetical storm scenario, it is estimated that it would
  take nine weeks to pump the water out of the city, and only then
  could assessments begin to determine what buildings were habitable
  or salvageable. Sewer, water, and the extensive forced drainage
  pumping systems would be damaged. National authorities would be
  scrambling to build tent cities to house the hundreds of thousands
  of refugees unable to return to their homes and without other
  relocation options. In the aftermath of such a disaster, New Orleans
  would be dramatically different, and likely extremely diminished,
  from what it is today. Unlike the posthurricane development surges
  that have occurred in coastal beach communities, the cost of
  rebuilding the city of New Orleans' dramatically damaged
  infrastructure would reduce the likelihood of a similar economic
  recovery. And, the unique culture of this American original that
  contributed jazz and so much more to the American culture would be
  lost. 



Accepting the Reality



  Should this disaster become a reality, it would undoubtedly be one
  of the greatest disasters, if not the greatest, to hit the United
  States, with estimated costs exceeding 100 billion
  dollars. According to the American Red Cross, such an event could be
  even more devastating than a major earthquake in
  California. Survivors would have to endure conditions never before
  experienced in a North American disaster. 



  Loss of the coastal marshes that dampened earlier storm surges puts
  the city at increasing risk to hurricanes. Eighty years of
  substantial river leveeing has prevented spring flood deposition of
  new layers of sediment into the marshes, and a similarly lengthy
  period of marsh excavation activities related to oil and gas
  exploration and transportation canals for the petrochemical industry
  have threatened marsh integrity. Sea level rise is expected to
  further accelerate the loss of these valuable coastal wetlands, the
  loss of which jeopardizes the fabric of Louisiana communities by
  threatening the harvesting of natural resources, an integral part of
  coastal culture. Concerted efforts by state and federal agencies are
  underway to develop appropriate restoration technologies and
  adequate funding to implement them. 



The Future is Now



  These solutions may not be able to overtake the speed of coastal
  loss. Strong storms not only threaten human lives, but also the
  physical coast itself. National hurricane experts predict more
  active and powerful hurricane seasons in the Atlantic basin for the
  next 10-40 years. The hurricane scenario for New Orleans that these
  converging risks portend is almost unimaginable. Hurricane Ivan had
  the potential to make the unthinkable a reality. Next time New
  Orleans may not be so fortunate. 



Shirley Laska
Center for Hazards Assessment, Response and Technology
University of New Orleans

  BRL-CAD is a powerful constructive solid geometry solid modeling
  system that includes an interactive geometry editor, ray tracing
  support for rendering and geometric analysis, network distributed
  framebuffer support, image and signal-processing tools.

  Industrial Open Source alternative to proprietary 3D modeling
  kernels.  Open CASCADE is a powerful CAD/CAM/CAE kernel and
  development platform for 3D modeling applications. It consists of
  reusable C++ object libraries and a set of development tools that
  are available in Open Source. 

segysql.py -v -v -v -v -D line45.sgy
FILENAME:  line45.sgy
fileKey =  1
WARNING: seems you have variable length traces that will not be handled!!!
Traceback (most recent call last):
  File "/Users/schwehr/x/src/segy-py/segysql.py", line 120, in ?
    for i in range(1,xstar.getNumberOfTraces()+1):
  File "/Users/schwehr/projects/xcore/src/segy-py/segy.py", line 237, in getNumberOfTraces
    assert(False)
AssertionError

Data file: ieee.segy
Sort file: sort.default.pltsegy.sort
 Data format  5
rgmax = rgmin : x scale not defined

      rgmin (r), rgmax (r) -
           minimum/maximum range of plot. In general the minimum
           and maximum ranges will be determined from the ranges
           of the traces to be plotted. However, it is possible
           to plot only travel time curves in which case dtfile
           is not specified and rgmin, rgmax must be used to
           specify ranges on a plot.  rgmin/rgmax are ignored if
           a dtfile is given.

Hi Kurt.



This seems to be a common issue caused by tidy or by
the tidy.sh plugin command line arguments in mac osx.
I never investigated it fully, but as a workaround you
can disable most of the command line args in tidy.sh
by setting TIDY_ARGS=" " (blank space).



Kevin

  magpy is a Python wrapper for the mg full-text search engine, which
  has support for fast full-text indexing/retrieval and ranked/boolean
  queries.

  JabRef is a graphical frontend to manage BibTeX databases, the
  standard LaTeX bibliography reference format. JabRef is build to be
  platform independant (requires Java 1.4.2 or newer). It merges and
  extends the functionalities of BibKeeper (Morten O. Alver) and
  JBibtexManager (Nizar Batada).

#!/bin/bash



fileIn=$1
fileOut=$2



sioseis <<EOF
  procs diskin diskoa end
  diskin lno 99999 ipath $fileIn end end
  diskoa opath $fileOut ofmt 5 end end end
EOF

#!/bin/bash
if [ 0 == 1 ]; then
  echo This does not get run
fi
if [ 0 == 0 ]; then
  echo This DOES RUN!
fi

#!/bin/bash
if [ false ]; then
  echo This does not get run
fi
if [ true ]; then
  echo This DOES RUN!
fi

  # grd2cpt -L-800/0      
  Bus error

  SourcePlease removes some of the insanity of SourceForge file
  listings by linking directly to a specific mirror's copy of a file,
  removing the need to wade through five pages of ads.

  Seashore is an open source image editor for Cocoa. It features
  gradients, textures and anti-aliasing for both text and brush
  strokes. It supports multiple layers and alpha channel editing. It
  is based around the GIMP's technology and uses the same native file
  format.  

THIS CABLE SYSTEM IS 
CONDUCTING A REQUIRED TEST 
OF THE EAS SYSTEM

  When the tide is out, the polyps are contracted and quiet. As the
  tide covers the colonies, "scouts" move out into the border to look
  for empty space to occupy. Larger, well-armed "warriors" inflate
  their stinging arms and swing them around. Towards the center of the
  colony, poorly armed "reproductive" anemones stay out of the fray
  and conduct the clone's business of breeding.

• http://plone.org/
• SOAPpy and ZSI
• twisted (which I wrote the original fink info file for)

  GeoTypes is a Python library that implements both the
  OpenGIS/PostGIS and standard PostgreSQL geometry types. It
  integrates with the psycopg Python/PostgreSQL interface. It provides
  implementations of all of the OpenGIS/PostGIS classes, except
  (x,y,m) and (x,y,z,m). It currently supports the EWKB, HEXEWKB, WKB,
  and WKT formats.

cat blend.file2



mbari-1c.grd -R-124.83333/-124.5/40.75/41 10
eureka.grd -R-124.6/-124.1/40.7/41 5
3sec.grd -R-125.5/-123/39/42 1

export R=-125/-124/40.5/41
grdblend blend.file2 -Geel-combined.grd -R$R -I1c

2004-10-10 gmt_grdio.c PATCH Subroutine opening grdfile for row-by-row reading (for grdblend.c) failed to initialize scale/offset.

cat blend.file2
mbari-1c.grd -R-124.83333/-124.5/40.75/41 3
eureka.grd -R-124.6/-124.1/40.7/41 2
3sec.grd -R-125.5/-123/39/42 1



grdblend blend.file2 -Gfoo.grd -R-124.6/-124.1/40.7/41 -I1c



grdinfo foo.grd 
foo.grd: Title: 
foo.grd: Command: grdblend blend.file2 -Gfoo.grd -R-124.6/-124.1/40.7/41 -I1c
foo.grd: Remark: 
foo.grd: Normal node registration used
foo.grd: grdfile format # 0
foo.grd: x_min: -124.6 x_max: -124.1 x_inc: 0.000277778 units: user_x_unit nx: 1801
foo.grd: y_min: 40.7 y_max: 41 y_inc: 0.000277778 units: user_y_unit ny: 1081
foo.grd: z_min: -750.13 z_max: 242.963 units: user_z_unit
foo.grd: scale_factor: 1 add_offset: 0

There are a couple of questions burried below...



I am having a bit of trouble merging some higher
resolution grids with a lower res 3sec coastal relief
dataset.  Maybe some on here can point me in the right
direction.



* I have been trying to merge the two better grids as
xyz with mbm_grid.  The first question is that mb_grid
seems to ignore -E.0002 or -E1c for grid cell size by
always making 100x100 grids, but -D4000/1800 seems to
make a bigger grid.  What is a proper -E option?



* When I used -Klarge.grd with mb_grid I get a silent
seg fault with no indication of what died, where or
why.  I did a grdcut to my local region and that went
away.  I presume this just a memory allocation
exceeded?  So I just need to do grdcuts in my scripts?



e.g. 
NaN values used to flag regions with no data
MBIO parameters:
  Ping averaging:       1
  Longitude flipping:   0
  Speed minimum:       0.0 km/hr
1753756 data points processed in Orig/Eureka.xyz
6145312 data points processed in Orig/mbari.eela.xyz



Extracting background from grid file
coastal-relief-3sec.grd...
Segmentation fault
All done!






* This last one is probably a bug with how I packaged
mbsystem in fink.  Any suggestion for how to get
mbsystem to find grdraster.info?



  Speed minimum:       0.0 km/hr
1753756 data points processed in Orig/Eureka.xyz
6145312 data points processed in Orig/mbari.eela.xyz



Extracting background from grdraster dataset 3...
grdraster: ERROR cannot open r
/sw/share/dbase/grdraster.info
grdraster:  ERROR reading grdraster.info file.
grdinfo:  File tmpgrdraster18266.grd not found
grdsample: No such file or directory
[tmpgrdraster18266.grd]
grd2xyz: No such file or directory
[tmpgrdsample18266.grd]
rm: tmpgrd*18266.grd: No such file or directory



When I go looking for grdraster, it is in
/sw/share/gmt/dbase, which is what is required of fink
packages, so I am not allowed to modify the gmt
install.  I tried setting export
GMT_GRIDDIR=/sw/share/gmt/dbase/  but having no luck. 
Then I get this error message:



grdraster ERROR opening
/sw/share/gmt/dbase//haxby5m_faa.i2 for read.
grdinfo:  File tmpgrdraster18294.grd not found
grdsample: No such file or directory
[tmpgrdraster18294.grd]
grd2xyz: No such file or directory
[tmpgrdsample18294.grd]



I am not turning up anything with haxby on my system.






* Lastly, if anyone can suggest good tutorials on
bringing together different resoluion data, I would
really appreciate it.  I started with Nathan Becker's
word document, but seem to be getting tripped up on...
that's why I gave mbsystem a try.

#! /bin/csh -f



set REGION          = -125/-124/40:30/41:00
set INPUT_FILE      = datalist.dst
set INPUT_FORMAT    = -1
set ROOT       	= datalist.dst



mbgrid -I$INPUT_FILE \
	-R$REGION \
	-O$ROOT \
	-A1 -N \
	-K3sec.grd \
	-D4000/1800



#	-Kcoastal-relief-3sec.grd \
#	-E0.0002 \
echo All done!

% cat datalist.dst



Orig/Eureka.xyz      162 10
Orig/mbari.eela.xyz  162 5

awk '{print $1,$2}' eel-strange-block.xyz  > block.xy



  521  bzcat 3sec-sub.txt.bz2 | awk '{print $1,$2,$3,.4}' |gmtconvert -bo > 3sec.xyzw.bin
  530  bzcat Orig/Eureka.xyz.bz2 | awk '{print $1,$2,$3,.8}' |gmtconvert -bo > eureka.xyzw.bin
  531  bzcat Orig/mbari.eela.xyz.bz2 | awk '{print $1,$2,$3,.8}' |gmtconvert -bo > mbari.xyzw.bin
  537  less mbari.xyzw.bin 
  544  blockmean -bi -W {3sec,eureka,mbar}.xyzw.bin -R$R -I1c -bo > all.blockmean.xyzw.bin
  545  blockmean -bi -W {3sec,eureka,mbari}.xyzw.bin -R$R -I1c -bo > all.blockmean.xyzw.bin
  549  blockmean -bi -W {3sec,eureka,mbari}.xyzw.bin -R$R -I1c -bo > all.blockmean.xyzw.bin
  551  nearneighbor -bi all.blockmean.xyzw.bin -R$R -I1c -N4 -S4c -Gall.nearneighbor.grd
  552  man mbdatalist



  576  mbm_grid -Idatalist.dst -R$R -E1c -Kcoastal-relief-3sec.grd 
  604  grdcut coastal-relief-3sec.grd -G3sec.grd -R$R -V -fg

  SWAP is a collaborative project to provide a working set of
  hardware, software and networking configurations to facilitate
  ad-hoc mesh networking between ships, and between ships and shore
  facilities, within the UNOLS oceanographic research fleet. SWAP is
  also being developed to facilitate connections between other kinds
  of nodes, including instrumented buoys, light houses and cabled
  observatories.

  Almost certainly it is 32 bit floating point IEEE, which was not a
  'standard' SEGY rev. 0 format back when these data were processed.

  I went through the segy files for W9605B and made a list of the ones
  with Format=6.  We ran a SIOSEIS script which changed it to 5.   The
  modified files are now available through the website.  At the end of
  this e-mail is the list that I came up with. 



Cheers,
Lisa



ar54.4173.w9605b.eel003.migration.segy
ar54.4176.w9605b.eel007.migration.segy
ar54.4178.w9605b.eel010.migration.segy
ar54.4180.w9605b.eel011.migration.segy
ar54.4181.w9605b.eel014.migration.segy
ar54.4183.w9605b.eel014.migration.segy
ar54.4185.w9605b.eel015.migration.segy
ar54.4186.w9605b.eel018.migration.segy
ar54.4188.w9605b.eel018.migration.segy
ar54.4190.w9605b.eel019.migration.segy
ar54.4192.w9605b.eel022.migration.segy
ar54.4194.w9605b.eel023.migration.segy
ar54.4196.w9605b.eel027.migration.segy
ar54.4198.w9605b.eel027.migration.segy
ar54.4199.w9605b.eel028.migration.segy
ar54.4202.w9605b.eel031.migration.segy
ar54.4204.w9605b.eel035.migration.segy
ar54.4207.w9605b.eel039.migration.segy
ar54.4210.w9605b.eel042.migration.segy
ar54.4212.w9605b.eel043.migration.segy
ar54.4215.w9605b.eel045.migration.segy
ar54.4217.w9605b.eel046.migration.segy
ar54.4219.w9605b.eel046.migration.segy
ar54.4220.w9605b.eel050.migration.segy
ar54.4223.w9605b.eel051.migration.segy
ar54.4225.w9605b.eel053.migration.segy
ar54.4227.w9605b.eel054.migration.segy
ar54.4229.w9605b.eel058.migration.segy
ar54.4231.w9605b.eel059.migration.segy
ar54.4233.w9605b.eel062.migration.segy
ar54.4235.w9605b.eel063.migration.segy
ar54.4237.w9605b.eel063.migration.segy
ar54.4239.w9605b.eel066.migration.segy
ar54.4241.w9605b.eel067.migration.segy
ar54.4243.w9605b.eel068.migration.segy
ar54.4245.w9605b.eel070.migration.segy
ar54.4247.w9605b.eel071.migration.segy
ar54.4249.w9605b.eel075.migration.segy
ar54.4251.w9605b.eel075.migration.segy
ar54.4253.w9605b.eel079.migration.segy
ar54.4255.w9605b.eel083.migration.segy
ar54.4257.w9605b.eel087.migration.segy
ar54.4258.w9605b.eel091.migration.segy
ar54.4261.w9605b.eel095.migration.segy
ar54.4263.w9605b.eel099.migration.segy
ar54.4265.w9605b.eel103.migration.segy
ar54.4267.w9605b.eel107.migration.segy
ar54.4269.w9605b.eel111.migration.segy

In window 1:
   cd ~/Desktop
   watch -n 10 -d ls -l
In window 2:
   cd ~/Desktop
   touch foo; sleep 11; touch bar; sleep 11; rm -f foo bar

...
  Prior to MESSENGER, JHU/APL missions used a raw storage model of
  storing data, but new mission and operational requirements meant
  that MESSENGER would have to incorporate a file system of data
  storage into its spacecraft software architecture. A reliable method
  of downlinking files to the ground had to be found and CFDP was
  chosen by mission planners to get the job done. 
...
  CFDP allows an instrument to record an observation in a file and
  transmit the file to Earth without having to consider whether or not
  physical transmission is possible at that time. Sequestering
  outbound data management and transmission planning functions within
  CFDP can simplify flight and ground software, which reduces mission
  costs - an important benefit to lower cost missions like
  MESSENGER.
...

  Even wonder what people have written about an area? 



  The Blogdigger cognoscenti have added native Google Earth publishing
  to their spatial blog search capabilities. (Fast work!) Visit
  Blogdigger, click on Local, specify your search (what, where,
  click), and then notice the box to the right, and the link at the
  bottom reading, "See it in Google Earth." Guess what that does? Now
  click on one of the icons.  

  The key to all of this is Location Messaging: a symmetric messaging
  technology where every user is both a content provider and a content
  user. Every message, whether by email or by SMS, has a location tag
  that associates it with a place. 

segydump -t 1 ar54.4178.w9605b.eel010.migration.segy



filename="ar54.4178.w9605b.eel010.migration.segy"



# FILE HEADER
CorrelatedTraces=0
TraceSorting=4
SweepTaperLenStart=0
SweepLen=0
MeasurementSystem=1
SampleFormat=6
SweepTaperLenEnd=0
JobId=0
SegYRev=0
LineNo=10
TracesPerEnsemble=0
SweepType=0
SamplesPerTrace=4200
TaperType=0
VibPolarity=0
BinGainRecovered=0
AmpRecovMeth=0
OrigSampleInterval=500
FixedLenTraceFlag=0
EnsembleFold=0
ImpulsePolarity=0
ReelNo=0
VerticalSum=0
NumExtendedTextHeaders=0
TraceSweepChannel=0
SweepStart=0
OrigSamplesPerTrace=0
SweepEnd=0
SampleInterval=500
AuxTracesPerEnsemble=0
SweepType= Unknown value - '0'
TraceSorting='Horizontally stacked'
SampleFormat='Not currently used'
1:  SweepFreqStart=0
1:  InstrInitGain=0
1:  RecvGrpElev=0
1:  SweepLen=0
1:  OverTravel=0
1:  GapSize=0
1:  SrcEDirTenths=0
1:  LowCutFreq=0
1:  ShotpointScaler=0
1:  ElevRcvGrp=0
1:  Day=206
1:  MuteStart=0
1:  TraceSamples=4200
1:  InstrGainConst=0
1:  AliasFiltFreq=0
1:  X=0
1:  SrcDepth=0
1:  SubWxVel=0
1:  SrcMeasExp=0
1:  Yensemble=0
1:  TransdConstMant=0
1:  HorzStackedNo=0
1:  LineSeqNo=0
1:  GrpY=0
1:  GrpX=0
1:  TimeScaler=0
1:  Unassigned233=4294967296
1:  TimeBasis=0
1:  TraceNoEnsemble=1
1:  FieldRecNo=1
1:  HighCutSlope=0
1:  DataUse=1
1:  HighCutFreq=0
1:  SweepType=0
1:  TaperType=0
1:  SrcMeasMantisa=0
1:  Hour=22
1:  GrpCorr=0
1:  PostInLineNo=0
1:  TotStatic=0
1:  SrcEDir=0
1:  DistanceCenter=0
1:  WaterDepthSrc=0
1:  SweepFreqEnd=0
1:  SampleInterval=500
1:  TraceWeight=0
1:  SrcMeasUnits=0
1:  TransdConstExp=0
1:  GeopGrpTr1=0
1:  TraceNo=1
1:  Min=28
1:  ESrcPtNo=0
1:  SweepTaperLenEnd=0
1:  GainType=0
1:  UpholeTSrc=0
1:  WxVel=0
1:  AliasFiltSlope=0
1:  LagTA=0
1:  Correlated=0
1:  FileSeqNo=1
1:  PostCrossLineNo=0
1:  LowCutSlope=0
1:  TransdUnits=0
1:  UpholeTGrp=0
1:  GeopGrpRlSw1=0
1:  Sec=0
1:  TraceUnits=0
1:  WaterDepthGrp=0
1:  SweepTaperLenStart=0
1:  Delay=0
1:  LagTB=0
1:  Xensemble=0
1:  SurfElevSrc=0
1:  SrcTypeOrient=0
1:  ElevSrc=0
1:  EnsembleNo=200
1:  DevTrId=0
1:  TraceId=1
1:  NotchFiltFreq=0
1:  GeopGrpNoLstTr=0
1:  Y=0
1:  SrcCorr=0
1:  MuteEnd=0
1:  ScalerCoord=-1000
1:  CoordUnits=2
1:  NotchFiltSlope=0
1:  ScalerElevDepth=0
1:  Shotpoint=0
1:  Year=1996
1:  VertSumNo=0
1:  CoordUnits_ascii='Seconds of arc'
1:  TraceId_ascii='Seismic data'
1:  TransdUnits_ascii='Unknown'
1:  SrcTypeOrient_ascii='Unknown'
1:  DataUse_ascii='Production'
1:  SrcMeasUnits_ascii='Unknown'
1:  OverTravel_ascii='__UNKOWN__'
1:  GainType_ascii='__UNKOWN__'
1:  SweepType_ascii='__UNKOWN__'
1:  TimeBasis_ascii='__UNKOWN__'
1:  TaperType_ascii='__UNKOWN__'
1:  TraceUnits_ascii='Unknown'

SEG-Y ASCII HEADER 39 LINES                                                    **********************************************
************************          INPUT STREAM STARTED FROM DISK FILE:                                            /home/mcspad
/spad/data/w9605b/seismic_segy/10.v.1.stolt                          LINE NUMBER FORCED TO    10                              
                       ACQUISTION YEAR FORCED TO 1996                                                 ************************
**********************************************          **********************************************************************
           INPUT SAMPLE INTERVAL is   500 MICROSECONDS                                      INPUT STREAM MODIFIED AS FOLLOWS: 
                                              CDP   TRACES USED 0001 THROUGH 0001                                             
MAX TAPE SECS=  3    MAX STAR MSECS=   3000                                     1ST MSECS TO USE=     0    MAX MSECS TO USE=  
3000                                                                                                             PROGRAM:    S
TAR                   11/15/00  12:17:53                           DIRECTORY:  /home/mcsusers/w9605b                          
                     USER:       jp                                                                  PARAMETER FILE om.stk.cop
y                                                       VELOCITY FILE: /data/mcsdb/vels/ih/w9605b/10.v.1                      
            0 DIFFERENT VELOCITIES. SEE NMO FOR VELOCITY CHANGES                           FIRST CDP   200, LAST CDP  4189, TO
TAL CDPS    24, # OF LAYERS 40             LINE   10  SAMPLE INTERVAL   500 MICROSECONDS OUTPUT                             AL
L TRACES USED                                                               OUTPUT STREAM TO DISK FILE:                       
                              /home/mcspad/ldata/w9605b/stolt/10.v.1.stolt                                    ****************
******************************************************          INPUT STREAM STARTED FROM DISK FILE:                          
                  /data/mcs/w9605b/stolt/10.v.1.stolt                                              ACQUISTION YEAR FORCED TO 1
996                                                  INPUT SAMPLE INTERVAL IS   500 MICROSECONDS                              
       PROGRAM:  TAPECOPY                  11/16/00  10:26:29                          OUTPUT SAMPLE INTERVAL IS   500 MICROSE
CONDS                                     ALL TRACES USED                                                                ACQUI
STION YEAR FORCED TO 1996                                                  INPUT SAMPLE INTERVAL IS   500 MICROSECONDS        
                             PROGRAM:  TAPECOPY                  11/16/00  10:40:09                          OUTPUT SAMPLE INT
ERVAL IS   500 MICROSECONDS                                     ALL TRACES USED                                               
                OUTPUT STREAM TO DISK FILE:                                                     /data/mstore/ldeo/users/mcsdat
a/w9605b/10.v.1.stolt

    sqlite eel-utexas-migration.db 'select x/3600.0/ScalerCoord, y/3600.0/ScalerCoord from segyTrace' \
      | tr '|' ' ' | egrep -v '^-0 -0$' > test.xy



head test.xy 



124.522895 -40.7113669444444
124.522828888889 -40.7113433333333
124.522762777778 -40.7113197222222
124.522696944444 -40.7112963888889
124.522630833333 -40.7112727777778
124.522564722222 -40.7112491666667
124.522498611111 -40.7112255555556
124.5224325 -40.7112019444444
124.522366388889 -40.7111783333333
124.522300555556 -40.7111547222222

  Thank you for your order for Verizon Online DSL service! Your DSL
  order is still being processed, however, we are experiencing a
  further delay in activating the service on your telephone number
  xxx-xxx-xxxx. Please know that we are working hard to provide you
  with service as soon as possible. 



  Your new estimated Service Ready date is now 08-30-2005*.



  We will continue to keep you updated on the status of your order
  through e-mail and will confirm by telephone once your DSL service
  is ready. You can also check the status of your order at
  http://www.verizon.net/ or call us at 1-800-xxx-xxxx if you have any
  questions.

• McMillan LaTeX example source
• McMillan LaTeX typeset (PDF)

  M-files to read and write SEG-Y files from Matlab and Octave,
  implemented using the syntax of the SEG-Y format : SEG-Y revsision 0
  (1975) and revision 1 (May 2002)

ZADAR (O.N. 0322417)
Ship | General Statistics | Engine | Builder
Owner(s) | Authorized Representative | All
Ship
Official Number	0322417	Year Built	1964
Ship Name	ZADAR	Year Rebuilt	-
Former Name	-	Port of Registry	VANCOUVER
IMO Number	-	Registry Date	1964 06 10
 	 	Certificate Expires	2007 04 30
8.3 tons