Coastal & Marine Geology InfoBank

Home FACS Activities Atlas Geology School Related Sites More

USGS InfoBank program -- getseis

Skip navigational links
Loading
InfoBank Programs: by Name   by Topic  
Expanded Description
Topic Description
Name
getseis
Synopsis
/infobank/programs/kml/getseis/getseis.for
Description
  Program to display seis by equipment type in google earth

       
AUTHOR
  Sumy Tom        4/25/2008

  Carolyn Degnan, 8/2/2012, Added sp format specification to output positive lat/lons with plus signs.
  Carolyn Degnan, 8/3/2012, Removed blanks before <range> and <altitude

    
Reads
 1      read(88, '(q,a)', end = 199)
     *     EQIDLength, EQIDType
        read(89, '(q,a)', end = 199)
     *     EQEPLength, EQEPType
        read(87, '(q,a)', end = 199)
     *     EQNAVLength, EQNAVType     
      read(EQIDType(FirstDot +1 : LastDot - 1), '(q,a)') 
     * EQNameLength, EQName
Writes
      write (94, '(4x,a)') '<Folde
      write (94, '(6x,a)') 
     * '<name>'//EQName(1:EQNameLength) //'&l
      write (90, '(1x,a)') '<Folde
      write (90, '(2x,a)') 
     * '<name>'// 'Activity IDs'//'<
      write (90, '(2x,a)') 
     * '<visibility>1</visi
      write (90, '(a)') '  <LookA
      write (90, '(sp,a,f11.6,a)') 
     *    '    <longitude>-100</lo
      write (90, '(sp,a,f10.6,a)') '    <latitude>+40</l
      write (90, '(a)') '    <altitude>6437300.</a
      write (90, '(a)') 
     * '    <altitudeMode>relativeToGround</altit
      write (90, '(a)') '    <range>6437300.<
      write (90, '(a)') '    <tilt>0&l
      write (90, '(a)') '    <heading>0</
      write (90, '(a)') '  </LookA
      write (91, '(1x,a)') '<Folde
      write (91, '(2x,a)') 
     * '<name>'// 'EndPoints'//'<
      write (91, '(a)') '  <LookA
      write (91, '(sp,a,f11.6,a)') 
     *    '    <longitude>-100.0</lo
      write (91, '(sp,a,f10.6,a)') '    <latitude>+40</l
      write (91, '(a)') '    <altitude>6437300.</a
      write (91, '(a)') 
     * '    <altitudeMode>relativeToGround</altit
      write (91, '(a)') '    <range>6437300.<
      write (91, '(a)') '    <tilt>0&l
      write (91, '(a)') '    <heading>0</
      write (91, '(a)') '  </LookA
      write (92, '(1x,a)') '<Folde
      write (92, '(2x,a)') 
     * '<name>'// 'Line IDs'//'<
      write (92, '(a)') '  <LookA
      write (92, '(sp,a,f11.6,a)') 
     *    '    <longitude>-100.0</lo
      write (92, '(sp,a,f10.6,a)') '    <latitude>+40</l
      write (92, '(a)') '    <altitude>6437300.</a
      write (92, '(a)') 
     * '    <altitudeMode>relativeToGround</altit
      write (92, '(a)') '    <range>6437300.<
      write (92, '(a)') '    <tilt>0&l
      write (92, '(a)') '    <heading>0</
      write (92, '(a)') '  </LookA
      write (93, '(1x,a)') '</Folde
      write (Header, '(a)') 
     *    '<?xml version="1.0" encoding="UTF-8"
      write (Header, '(a)') 
     *    '<kml xmlns="http://earth.google.com/kml/2.0
      write (Header, '(a)') '<Documen
      write (Header, '(a)') 
     *      '  <name>Seismic&l
      write (Header, '(4x,a)') '<open>0&l
      write (Header, '(a)') '   '
      write (Header, '(a)') '  <LookA
      write (Header, '(sp,a,f11.6,a)') 
     *    '    <longitude>-100.0</lo
      write (Header, '(sp,a,f10.6,a)') '    <latitude>+40</l
      write (Header, '(a)') '    <altitude>6437300.</a
      write (Header, '(a)') 
     * '    <altitudeMode>relativeToGround</altit
      write (Header, '(a)') '    <range>6437300.<
      write (Header, '(a)') '    <tilt>0&l
      write (Header, '(a)') '    <heading>0</
      write (Header, '(a)') '  </LookA
      write (Footer, '(a)') '</Documen
      write (Footer, '(a)') '</km
Opens
       open (unit = 88,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'old',
     * name = '/infobank/tmp/getseis.groupids.txt')
       open (unit = 89,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'old',
     * name = '/infobank/tmp/getseis.groupeps.txt') 
       open (unit = 87,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'old',
     * name = '/infobank/tmp/getseis.groupnavs.txt')      
       open (unit = 90,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.openfolderids.txt')
      open (unit = 91,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.openfoldereps.txt')
      open (unit = 92,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.openfoldernavs.txt')     
      open (unit = 93,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.closefolder.txt')
       open (unit = 94,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.eqfolder.txt')
      open (unit = Header,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.lookat.temp')
      open (unit = Footer,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/getseis.close.temp')
Calls
      call system ('/bin/ls -1 '//
     * '/infobank/programs/kml/kmlinc2incgroups/seisids/*.kmlinc
     *  >> /infobank/tmp/getseis.groupids
        call system ('/bin/ls -1 '//
     * '/infobank/programs/kml/kmlinc2incgroups/seiseps/*.kmlinc
     *  >> /infobank/tmp/getseis.groupeps.t
        call system ('/bin/ls -1 '//
     * '/infobank/programs/kml/kmlinc2incgroups/seisnavs/*.kmlinc
     *  >> /infobank/tmp/getseis.groupnavs.txt')
      call includeit (Header,
     *    '/infobank/programs/keep/infobank_style.keep')
      call system('/bin/cat '//
     * '/infobank/tmp/getseis.eqfolder.txt '//
     *  EQIDType(1:EQIDLength)//' '//
     * '/infobank/tmp/getseis.closefolder.txt >&g
     * '/infobank/tmp/getseis.tempids.txt')      
      call system('/bin/cat '//
     * '/infobank/tmp/getseis.eqfolder.txt '//
     *  EQEPType(1:EQEPLength)//' '//
     * '/infobank/tmp/getseis.closefolder.txt >&g
     * '/infobank/tmp/getseis.tempeps.txt') 
      call system('/bin/cat '//
     * '/infobank/tmp/getseis.eqfolder.txt '//
     *  EQNAVType(1:EQNAVLength)//' '//
     * '/infobank/tmp/getseis.closefolder.txt >&g
     * '/infobank/tmp/getseis.tempnavs.txt')      
       call system ('/bin/cat '//
     * '/infobank/tmp/getseis.lookat.temp ' //
      call kmlzip('/infobank/programs/kml'//
     * '/getseis/kml/'//
     *  EQName(1:EQNameLength)//'.kml')
      call system ('/bin/rm /infobank/tmp/getseis.*.temp') 
      call system ('/bin/rm /infobank/tmp/getseis.*.txt')

Skip footer navigational links

Coastal and Marine Science Centers:  Pacific   St. Petersburg   Woods Hole  
InfoBank   Coastal and Marine Geology Program   Geologic Information   Ask-A-Geologist   USGS Disclaimer  


Accessibility FOIA Privacy Policies and Notices

Take Pride in America logo USA.gov logo U.S. Department of the Interior | U.S. Geological Survey
URL: http://walrus.wr.usgs.gov/infobank/programs/kml/getseis/getseis.doc.html
Page Contact Information: InfoBank staff
Page Last Modified: Mon Sep 16 03:38:58 PDT 2013  (chd)