USGS InfoBank program -- oceans4kml

oceans4kml

  Program to create virtual globe with oceans lit up

       

  Alex Ma     6/21/2006

  based on atlaskmlsuite
  Clint Steele     2/22/2008     upgrade

    

          read (89, '(q,a)', end = 98)
     *        ALineLength, ALine(1:ALineLength)
          read (90, '(q,a)', end=199) Record1Length, 
     *                   Record1(1:Record1Length)
      read (Record1(RecTab(1)+1:RecTab(2)-1), '(f16)') North
      read (Record1(RecTab(2)+1:RecTab(3)-1), '(f16)') South
      read (Record1(RecTab(3)+1:RecTab(4)-1), '(f16)') East
      read (Record1(RecTab(4)+1:RecTab(5)-1), '(f16)') West

      write (Blank, '(80x)')
          write (88, '(a)') ALine(Tab(8)+1:Tab(9)-1)//
     *        char(9)//ALine(1:Tab(8)-1)//
     *        ALine(Tab(9):ALineLength)//
     *        char(9)//char(9)//char(9)//char(9)
      write (91, '(a)') '<?xml version="1.0" encoding="UTF-8"
      write (91, '(a)') 
     *     '<kml xmlns="http://earth.google.com/kml/2.0
      write (91, '(a)') Blank(1:Level)//'<Documen
      write (91, '(a)') Blank(1:Level)//'<name>Highlight Oceans&l
      write (91, '(a)') Blank(1:Level)//'<LookA
      write (91, '(a)') Blank(1:Level)//'<longitude>-113.33</lo
      write (91, '(a)') Blank(1:Level)//'<latitude>41.1</l
      write (91, '(a)') Blank(1:Level)//'<range>7600000<
      write (91, '(a)') Blank(1:Level)//'</LookA
      write (91, '(a)') Blank(1:Level)//'<visibility>0</vis
          write (91, '(a)') Blank(1:Level)//'</Folde
          write (91, '(a)') Blank(1:Level)//'<Folde
          write (91, '(a)') Blank(1:Level+2)//
     *        '<name>'//Current//'&l
      write (91, '(a)') Blank(1:Level)//'<GroundOverla
      write (91, '(a)') Blank(1:Level)//
     *    '<name>'//AtlasName(1:AtlasNameLength)//'&l
      write (91, '(a)') Blank(1:Level)//'<visibility>1</vis
      write (91, '(a)') Blank(1:Level)//'<Ico
      write (91, '(a)') Blank(1:Level)//
     *    '<href>http://walrus.wr.usgs.gov/infobank/progr
     *    '/html/kml/regions/'//AtlasID(1:AtlasIDLength)//
     *    '/'//AtlasID(1:AtlasIDLength)//'.jpg</hre
      write (91, '(a)') Blank(1:Level)//'</Ico
      write (91, '(a)') Blank(1:Level)//'<Regio
      write (91, '(a)') Blank(1:Level)//'<LatLonAltBo
      write (91, '(a,f9.5,a)') Blank(1:Level)//
     *        '<north>',North,'<
      write (91, '(a,f9.5,a)') Blank(1:Level)//
     *        '<south>',South,'<
      write (91, '(a,f10.5,a)') Blank(1:Level)//
     *        '<east>',East,'&l
      write (91, '(a,f10.5,a)') Blank(1:Level)//
     *        '<west>',West,'&l
      write (91, '(a)') Blank(1:Level)//
     *           '<altitudeMode>relativeToGround</altit
      write (91, '(a)') Blank(1:Level)//'<minAltitude>0</minA
      write (91, '(a)') Blank(1:Level)//'<maxAltitude>0</maxA
      write (91, '(a)') Blank(1:Level)//'</LatLonAltBo
      write (91, '(a)') Blank(1:Level)//'<Lo
      write (91, '(a)') Blank(1:Level)//
     *    '<minLodPixels>750</minLo
      write (91, '(a)') Blank(1:Level)//
     *    '<maxLodPixels>1080</maxLo
      write (91, '(a)') Blank(1:Level)//
     *    '<minFadeExtent>0</minFad
      write (91, '(a)') Blank(1:Level)//
     *    '<maxFadeExtent>0</maxFad
      write (91, '(a)') Blank(1:Level)//'</Lo
      write (91, '(a)') Blank(1:Level)//'</Regio
      write (91, '(a)') Blank(1:Level)//'<LatLonBo
      write (91, '(a,f9.5,a)') Blank(1:Level)//
     *    '<north>',North,'<
      write (91, '(a,f9.5,a)') Blank(1:Level)//
     *    '<south>',South,'<
      write (91, '(a,f10.5,a)') Blank(1:Level)//
     *    '<east>',East,'&l
      write (91, '(a,f10.5,a)') Blank(1:Level)//
     *    '<west>',West,'&l
      write (91, '(a)') Blank(1:Level)//'<rotation>0</r
      write (91, '(a)') Blank(1:Level)//'</LatLonBo
      write (91, '(a)') Blank(1:Level)//'</GroundOverla
      write (91, '(a)') Blank(1:Level)//'</Folde
      write (91, '(a)') Blank(1:Level)//'</Documen
      write (91, '(a)') '</km

      open (unit = 89,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'old',
     *      name = '/infobank/programs/keep/gazette.keep')
      open (unit = 88,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/oceans4kml.gazette.temp')
      open (unit = 90,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'old',
     *      name = '/infobank/tmp/oceans4kml.temp')
      open (unit = 91, 
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/programs/kml'//
     *             '/oceans4kml/oceans4kml.kml')

      call system ('/bin/sort '//
     *             '/infobank/tmp/oceans4kml.gazette.temp '//
     *             '> /infobank/tmp/oceans4kml.tem
      call system ('/bin/rm '//
     *             '/infobank/tmp/oceans4kml.gazette.temp')
      call includeit (91,
     *    '/infobank/programs/keep/infobank_style.keep')
      call changecase ('SMALLS', SMallAtlasID)
      call changecase ('SMALLS', AtlasID(1:AtlasIDLength))
      call changecase ('CAPITALS', CapAtlasID)
      call kmlzip 
     *  ('/infobank/programs/kml/oceans4kml/oceans4kml.kml')