SUBROUTINE WCI_NEWSYS( NAME, EQNX, EPOCH, ID, STATUS )
*+
*  Name:
*     WCI_NEWSYS

*  Purpose:
*     Create a new coordinate system object

*  Language:
*     Starlink Fortran

*  Invocation:
*     CALL WCI_NEWSYS( NAME, EQNX, EPOCH, ID, STATUS )

*  Description:
*     Creates a coordinate system object which fully defines an absolute
*     celestial coordinate system. 

*  Arguments:
*     NAME = CHARACTER*(*) (given)
*        Name of the coordinate system. Can take values FK4, FK5, ECLIPTIC,
*        GALACTIC or SUPERGALACTIC, or any abbreviation thereof.
*     EQNX = REAL (given)
*        The epoch of the mean equator and equinox in years.
*     EPOCH = DOUBLE (given)
*        Epoch of the observation in years. If the value WCI__FLAG is 
*        supplied then the year of the EQNX is used.
*     ID = INTEGER (returned)
*        The ADI identifier of the new CoordSystem object
*     STATUS = INTEGER (given and returned)
*        The global status.

*  Examples:
*     CALL WCI_NEWSYS( 'FK4', 1950.0, WCI__FLAG, ID, STATUS )
*        {routine_example_description}

*  Pitfalls:
*     {pitfall_description}...

*  Notes:
*     {routine_notes}...

*  Prior Requirements:
*     {routine_prior_requirements}...

*  Side Effects:
*     {routine_side_effects}...

*  Algorithm:
*     {algorithm_description}...

*  Accuracy:
*     {routine_accuracy}

*  Implementation Deficiencies:
*     {routine_deficiencies}...

*  References:
*     WCI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/wci.html

*  Keywords:
*     package:wci, usage:public

*  Copyright:
*     Copyright (C) University of Birmingham, 1995

*  Authors:
*     DJA: David J. Allan (Jet-X, University of Birmingham)
*     {enter_new_authors_here}

*  History:
*     4 Jan 1995 (DJA):
*        Original version.
*     {enter_changes_here}

*  Bugs:
*     {note_any_bugs_here}

*-
      
*  Type Definitions:
      IMPLICIT NONE              ! No implicit typing

*  Global Constants:
      INCLUDE 'SAE_PAR'          	! Standard SAE constants
      INCLUDE 'WCI_PAR'          	! ASTERIX WCI constants
      INCLUDE 'AST_PKG'

*  Arguments Given:
      CHARACTER*(*)		NAME 
      REAL			EQNX
      DOUBLE PRECISION  	EPOCH

*  Arguments Returned:
      INTEGER			ID

*  Status:
      INTEGER 			STATUS             	! Global status

*  External References:
      EXTERNAL			AST_QPKGI
        LOGICAL			AST_QPKGI

*  Local variables:
      CHARACTER*1		EFORM			! Epoch code
      CHARACTER*3		LNAME			! Coord system name

      DOUBLE PRECISION		LEPOCH			! Local epoch
*.

*  Check inherited global status.
      IF ( STATUS .NE. SAI__OK ) RETURN

*  Check initialised
      IF ( .NOT. AST_QPKGI( WCI__PKG ) ) CALL WCI1_INIT( STATUS )

*  Check system name
      CALL WCI1_CHKSNM( NAME, STATUS )

*  Create new instance of CoordSystem
      CALL ADI_NEW0( 'CoordSystem', ID, STATUS )

*  Store attributes
      LNAME = NAME
      CALL CHR_UCASE( LNAME )
      CALL ADI_CPUT0C( ID, 'NAME', LNAME, STATUS )
      CALL ADI_CPUT0R( ID, 'EQUINOX', EQNX, STATUS )

*  Use EQNX for the epoch if EPOCH has the flag value
      IF ( EPOCH .EQ. WCI__FLAG ) THEN
        LEPOCH = EQNX
      ELSE
        LEPOCH = EPOCH
      END IF

*  Choose the epoch type, Besselian or Julian
      IF ( LNAME .EQ. 'FK4' ) THEN
        EFORM = 'B'
      ELSE IF ( LNAME .EQ. 'FK5' ) THEN
        EFORM = 'J'
      ELSE
        IF ( LEPOCH .LT. 1984D0 ) THEN
          EFORM = 'J'
        ELSE
          EFORM = 'J'
        END IF
      END IF
      CALL ADI_CPUT0C( ID, 'EFORM', EFORM, STATUS )

*  Convert epoch to years if flag value not used
      CALL ADI_CPUT0D( ID, 'EPOCH', LEPOCH, STATUS )

*  Report any errors
      IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'WCI_NEWSYS', STATUS )

      END