Coastal & Marine Geology InfoBank

Home FACS Activities Atlas Geology School Related Sites More

USGS InfoBank program -- seisendpointsfolder

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

       
AUTHOR
  Clint Steele    3/2/2008

  Sumy Tom        5/29/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 (SeisEndPointsName, '(1024x)')
      write (SeisEndPointsUnit, '(2x,a)') '<Folde
      write (SeisEndPointsUnit, '(2x,a)') '<name>'//ID/
     *      'Line End Points</nam
      write (SeisEndPointsUnit, '(4x,a)') '<Folde
      write (SeisEndPointsUnit, '(4x,a)')
     *    '<name> Line '//ALine(IndexFromLatLong:ALineLength)//'&l
      write (SeisEndPointsUnit, '(4x,a)') '<Placemar
      write (SeisEndPointsUnit, '(6x,a)')
     *    '<styleUrl>#Triangle</s
      write (SeisEndPointsUnit, '(6x,a)')
     *    '<visibility>0</vis
      write (SeisEndPointsUnit, '(6x,a)') '<name&
     *    ALine(IndexFromLatLong:ALineLength)//
     *    '</nam
      write (SeisEndPointsUnit, '(6x,a)') '<Poin
      write (SeisEndPointsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *    '<coordinates>',Long,',
     *    ',0</coordinate
      write (SeisEndPointsUnit, '(6x,a)') '</Poin
      write (SeisEndPointsUnit, '(4x,a)') '</Placemar
          write (PrevLine, '(1024x)')
              write (SeisEndPointsUnit, '(4x,a)') '<Placemar
              write (SeisEndPointsUnit, '(6x,a)')
     *            '<styleUrl>#Square</s
              write (SeisEndPointsUnit, '(6x,a)')
     *            '<visibility>0</vis
              write (SeisEndPointsUnit, '(6x,a)')
     *            '<name&
     *            PrevLine(IndexFromLatLong:PrevLineLength)//'</nam
              write (SeisEndPointsUnit, '(6x,a)') '<Poin
              write (SeisEndPointsUnit, '(8x,a)') '<coordinate
              write (SeisEndPointsUnit, '(sp,10x,f11.6,a,f10.6,a)')
     *            Long,',',Lat,',0'
              write (SeisEndPointsUnit, '(8x,a)') '</coordinates>'       
              write (SeisEndPointsUnit, '(6x,a)') '</Poin
              write (SeisEndPointsUnit, '(4x,a)') '</Placemar
              write (SeisEndPointsUnit, '(a)') '  <LookA
              write (SeisEndPointsUnit, '(sp,a,f11.6,a)') 
     *             '<longitude>', Long3 , '</lo
              write (SeisEndPointsUnit, '(sp,a,f10.6,a)')  
     *             '<latitude>',Lat,'</l
              write (SeisEndPointsUnit, '(sp,a,f,a)')  
     *             ' <range> ', Range, ' <
              write (SeisEndPointsUnit, '(a)') '    <tilt>0&l
              write (SeisEndPointsUnit, '(a)')
     *            '    <heading>0</
              write (SeisEndPointsUnit, '(a)') '  </LookA
              write (SeisEndPointsUnit, '(4x,a)') '</Folde
              write (SeisEndPointsUnit, '(/,4x,a)') '<Folde
              write (SeisEndPointsUnit, '(4x,a)')
     *            '<name> Line '//ALine(IndexFromLatLong:ALineLen
     *            '</nam
              write (SeisEndPointsUnit, '(4x,a)') '<Placemar
              write (SeisEndPointsUnit, '(6x,a)')
     *            '<styleUrl>#Triangle</s
              write (SeisEndPointsUnit, '(6x,a)')
     *            '<visibility>0</vis
              write (SeisEndPointsUnit, '(6x,a)') '<name&
     *            ALine(IndexFromLatLong:ALineLength)//
     *            '</nam
              write (SeisEndPointsUnit, '(6x,a)') '<Poin
              write (SeisEndPointsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *             '<coordinates>',Long,',
     *             ',0</coordinate
              write (SeisEndPointsUnit, '(6x,a)') '</Poin
              write (SeisEndPointsUnit, '(4x,a)') '</Placemar
      write (SeisEndPointsUnit, '(4x,a)') '<Placemar
      write (SeisEndPointsUnit, '(6x,a)')
     *        '<styleUrl>#Square</s
      write (SeisEndPointsUnit, '(6x,a)')
     *        '<visibility>0</vis
      write (SeisEndPointsUnit, '(6x,a)')
     *        '<name&
     *        PrevLine(IndexFromLatLong:PrevLineLength)//'</nam
      write (SeisEndPointsUnit, '(6x,a)') '<Poin
      write (SeisEndPointsUnit, '(sp,8x,a,f11.6,a,f10.6,a)')
     *        '<coordinates>',Long,',',Lat,',0</coor
      write (SeisEndPointsUnit, '(6x,a)') '</Poin
      write (SeisEndPointsUnit, '(4x,a)') '</Placemar
      write (SeisEndPointsUnit, '(4x,a)') '</Folde
      write (SeisEndPointsUnit, '(1x,a)') '</Folde
Opens
      open (unit = SeisEndPointsUnit,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = SeisEndPointsName(1:SeisEndPointsNameLength))
Calls
      call stringlen (KmlUnitName, KmlUnitNameLength)
      call stringlen (SeisEndPointsName, SeisEndPointsNameLength)
          call system ('/bin/rm '//
     *                   SeisEndPointsName(1:SeisEndPointsNameLength))
      call describeseis (SeisEndPointsUnit, ID2, IDLength2,ID, 
     *               AkaIDLength, ALine, ALineLength, LongID,
     *               ALine(IndexFromLatLong:ALineLength)) 
      call describeseis (SeisEndPointsUnit, ID2, IDLength2,ID, 
     *            AkaIDLength, PrevLine, PrevLineLength, LongID,
     *            PrevLine(IndexFromLatLong:PrevLineLength)) 
              call  getrange( Lat, Long3,
     *          North,South, East ,West, Range)              
      call describeseis (SeisEndPointsUnit, ID2, IDLength2,ID, 
     *             AkaIDLength, ALine, ALineLength, LongID,
     *             ALine(IndexFromLatLong:ALineLength))     
      call describeseis (SeisEndPointsUnit, ID2, IDLength2,ID, 
     *            AkaIDLength, PrevLine, PrevLineLength, LongID,
     *            PrevLine(IndexFromLatLong: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/seisendpointsfolder.doc.html
Page Contact Information: InfoBank staff
Page Last Modified: Mon Sep 16 03:39:04 PDT 2013  (chd)