Coastal & Marine Geology InfoBank

Home FACS Activities Atlas Geology School Related Sites More

USGS InfoBank program -- seisidsfolder

Skip navigational links
Loading
InfoBank Programs: by Name   by Topic  
Expanded Description
Topic Description
Name
seisidsfolder
Synopsis
/infobank/programs/kml/seis4kml/seisidsfolder.for
Description
  Subroutine to create a .seisids.kmlinc file

       
AUTHOR
  Clint Steele    3/2/2008

  Sumy Tom        6/5/2008   
  Carolyn Degnan, 8/3/2012, Added sp format specification to output positive lat/lons with plus signs.
  Carolyn Degnan, 8/3/2012, Removed blanks before <range> v

    
Reads
      read (StandardUnit, '(q,a)', end = 99) 
     *    ALineLength, ALine(1:ALineLength)
         read (ALine, '(t16,f9.5,t26,f10.5)')
     *       Lat, Long
         read (ALine, '(t17,f10.6,t28,f11.6)')
     *    Lat, Long
          read (StandardUnit, '(q,a)', end = 99)
     *        ALineLength, ALine(1:ALineLength)
             read (ALine, '(t16, f9.5, t26, f10.5)')
     *             Lat, Long
             read (ALine, '(t17, f10.6, t28, f11.6)')
     *             Lat, Long
                 read (PrevLine, '(t16,f9.5,t26,f10.5)')
     *               Lat, Long
                 read (PrevLine, '(t17,f10.6,t28,f11.6)')
     *               Lat, Long
                 read (ALine, '(t16,f9.5,t26,f10.5)')
     *               Lat, Long
                 read (ALine, '(t17,f10.6,t28,f11.6)')
     *               Lat, Long
            read (PrevLine, '(t16,f9.5,t26,f10.5)') Lat, Long
            read (PrevLine, '(t17,f10.6,t28,f11.6)') Lat, Long
Writes
      write (SeisIDsName, '(1024x)')
      write (SeisIDsUnit, '(1x,a)') '<Folde
      write (SeisIDsUnit, '(2x,a)') '<name&
     *      Label//' Tags</nam
      write (SeisIDsUnit, '(2x,a)') '<Folde
      write (SeisIDsUnit, '(4x,a)')
     *    '<name>'//Label//' Li
     *    ALine(37:ALineLength)//'</nam
      write (SeisIDsUnit, '(4x,a)') '<Placemar
      write (SeisIDsUnit, '(6x,a)')
     *    '<styleUrl>#Triangle</s
      write (SeisIDsUnit, '(6x,a)')
     *    '<visibility>0</vis
      write (SeisIDsUnit, '(6x,a)') '<name&
     *    Label//
     *    '</nam
      write (SeisIDsUnit, '(6x,a)') '<Poin
      write (SeisIDsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *    '<coordinates>',Long,',
     *    ',0</coordinate
      write (SeisIDsUnit, '(6x,a)') '</Poin
      write (SeisIDsUnit, '(4x,a)') '</Placemar
          write (PrevLine, '(1024x)')
              write (SeisIDsUnit, '(4x,a)') '<Placemar
              write (SeisIDsUnit, '(6x,a)')
     *            '<styleUrl>#Square</s
              write (SeisIDsUnit, '(6x,a)')
     *            '<visibility>0</vis
              write (SeisIDsUnit, '(6x,a)')
     *            '<name&
     *            Label//'</nam
              write (SeisIDsUnit, '(6x,a)') '<Poin
              write (SeisIDsUnit, '(8x,a)') '<coordinate
              write (SeisIDsUnit, '(sp,10x,f11.6,a,f10.6,a)')
     *            Long,',',Lat,',0'
              write (SeisIDsUnit, '(8x,a)') '</coordinate
              write (SeisIDsUnit, '(6x,a)') '</Poin
              write (SeisIDsUnit, '(4x,a)') '</Placemar
              write (SeisIDsUnit, '(a)') '  <LookA
              write (SeisIDsUnit, '(sp,a,f11.6,a)') 
     *             '<longitude>', Long , '</lo
              write (SeisIDsUnit, '(sp,a,f10.6,a)')  
     *             '<latitude>',Lat,'</l
              write (SeisIDsUnit, '(sp,a,f,a)')  
     *             ' <range> ', Range, ' <
              write (SeisIDsUnit, '(a)') '    <tilt>0&l
              write (SeisIDsUnit, '(a)')
     *            '    <heading>0</
              write (SeisIDsUnit, '(a)') '  </LookA
              write (SeisIDsUnit, '(2x,a)') '</Folde
              write (SeisIDsUnit, '(/,2x,a)') '<Folde
              write (SeisIDsUnit, '(4x,a)')
     *            '<name>'//Label//' Li
     *            ALine(37:ALineLength)//'</nam
              write (SeisIDsUnit, '(4x,a)') '<Placemar
              write (SeisIDsUnit, '(6x,a)')
     *            '<styleUrl>#Triangle</s
              write (SeisIDsUnit, '(6x,a)')
     *            '<visibility>0</vis
              write (SeisIDsUnit, '(6x,a)') '<name&
     *            Label//
     *            '</nam
              write (SeisIDsUnit, '(6x,a)') '<Poin
              write (SeisIDsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *             '<coordinates>',Long,',
     *             ',0</coordinates
              write (SeisIDsUnit, '(6x,a)') '</Poin
              write (SeisIDsUnit, '(4x,a)') '</Placemar
      write (SeisIDsUnit, '(4x,a)') '<Placemar
      write (SeisIDsUnit, '(6x,a)')
     *        '<styleUrl>#Square</s
      write (SeisIDsUnit, '(6x,a)')
     *        '<visibility>0</vis
      write (SeisIDsUnit, '(6x,a)')
     *        '<name&
     *        Label//'</nam
      write (SeisIDsUnit, '(6x,a)') '<Poin
      write (SeisIDsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *        '<coordinates>',Long,',',Lat,',0</coor
      write (SeisIDsUnit, '(6x,a)') '</Poin
      write (SeisIDsUnit, '(4x,a)') '</Placemar
      write (SeisIDsUnit, '(4x,a)') '</Folde
      write (SeisIDsUnit, '(1x,a)') '</Folde
Opens
      open (unit = SeisIDsUnit,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = SeisIDsName(1:SeisIDsNameLength))
Calls
      call stringlen (KmlUnitName, KmlUnitNameLength)
      call stringlen (SeisIDsName, SeisIDsNameLength)
          call system ('/bin/rm '//
     *                   SeisIDsName(1:SeisIDsNameLength))
       call describeseis (SeisIDsUnit, ID2, IDLength2,ID, 
     *          AkaIDLength, ALine, ALineLength, LongID,
     *          ALine(37:ALineLength))     
      call describeseis (SeisIDsUnit, ID2, IDLength2,ID, AkaIDLength, 
     *                    PrevLine, PrevLineLength, LongID,
     *                    PrevLine(37:PrevLineLength))     
              call  getrange( Lat, Long3,
     *          North,South, East ,West, Range)              
      call describeseis (SeisIDsUnit, ID2, IDLength2,ID, AkaIDLength, 
     *                    ALine, ALineLength, LongID,
     *                    ALine(37:ALineLength))         
      call describeseis (SeisIDsUnit, ID2, IDLength2,ID, AkaIDLength, 
     *                    prevLine, PrevLineLength, LongID,
     *                    prevLine(37:PrevLineLength))        

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/seis4kml/seisidsfolder.doc.html
Page Contact Information: InfoBank staff
Page Last Modified: Mon Sep 16 03:39:04 PDT 2013  (chd)