USGS InfoBank program -- section41

section41

  Subroutine to create HORIZONTAL COORDINATE SYSTEM DEFINITION of formal metadata.


       

          write (OutUnit, '(1x)')

      call head (2, OutUnit,
     *            'Horizontal_Coordinate_System_Definition:')
      call head (3, OutUnit,
     *            'Geographic:')
      call value (4, OutUnit,
     *            'Latitude_Resolution:            ', 
     *            '0.00001')
      call value (4, OutUnit,
     *            'Longitude_Resolution:           ', 
     *            '0.00001')
      call value (4, OutUnit,
     *            'Geographic_Coordinate_Units:    ', 
     *            'Decimal degrees')
          call head (3, OutUnit,
     *            'Geodetic_Model:')
              call value (4, OutUnit,
     *            'Horizontal_Datum_Name:              ', 
     *            'North American Datum of 1927')
              call value (4, OutUnit,
     *            'Ellipsoid_Name:                     ', 
     *            'Clarke 1866')
              call value (4, OutUnit,
     *            'Semi-major_Axis:                    ', 
     *            '6378206.4')
              call value (4, OutUnit,
     *            'Denominator_of_Flattening_Ratio:    ', 
     *            '294.98')
              call value (4, OutUnit,
     *            'Horizontal_Datum_Name:              ', 
     *            'North American Datum of 1983')
              call value (4, OutUnit,
     *            'Ellipsoid_Name:                     ', 
     *            'Geodetic Reference System 80')
              call value (4, OutUnit,
     *            'Semi-major_Axis:                    ', 
     *            '6378137')
              call value (4, OutUnit,
     *            'Denominator_of_Flattening_Ratio:    ', 
     *            '298.257')