PROGRAM SFCCONV
c
c     $Revision: 1.1 $
c
csfcconv  - program to test the Geographic to/from Geomagnetic Coordinate
c             subprogram to be run at AFSFC
c
c
c Comments: This program and the SFC subroutine have been developed on a
c            DEC VAX 7000 using VMS FORTRAN and trying to maintain the
c         SFC style of coding.  The calling sequence is identical to the
c            previous version of the subprogram.
c
c    This version has been ported to the unix environment.  The original
c     program used a loop and ran through a set of altitudes from
c     0 km to 2000 km.  It opened the file SFCOUT.DAT for each height
c     and made use of the fact the VMS keeps multiple versions of
c     files to allow you to have a separate file for each altitude.
c
c     For the UNIX implementation the automatic loop through the
c     altitudes has been replaced by a request for the user to
c     input an arbitrary altitude to be used.
c
c     $Log: not supported by cvs2svn $
c     Revision 1.1  94/10/14  10:59:07  10:59:07  baker (Kile Baker S1G)
c     Initial revision
c
c
C     BEGIN PROGRAM
C     RUN FOR ALTITUDES FROM 0 KM TO 2000 KM
C      DO 200 JH = 0, 20
C     .. Local Scalars ..
      DOUBLE PRECISION CGLAT,CGLON,GLAT,GLON,GMLAT_,GMLON_,R_HEIGHT_IN
      INTEGER I_ERROR,LAT,LON
C     ..
C     .. External Subroutines ..
      EXTERNAL CONVERT_GEO_COORD
C     ..
C     .. Intrinsic Functions ..
      INTRINSIC DBLE
C     ..
      R_HEIGHT_IN = 0.0D0
C      type *,'enter the height to be used'
C      accept *,r_height_in
C
C
      OPEN (1,FILE='SFCOUT.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
CRUN FOR ALL GEOGRAPHICS LATITUDES (STEPPING BY 2 DEGREES) AND LONGITUDES
C     (STEPPING BY TEN DEGREES)
      DO 20 LAT = 67,69,2
         GLAT = DBLE(LAT)
         DO 10 LON = 309,319,10
            GLON = DBLE(LON)
C           Convert from Geograhic (Geocentric) to Geomagnetic
            CALL CONVERT_GEO_COORD(GLAT,GLON,R_HEIGHT_IN,CGLAT,CGLON,1,
     *                             I_ERROR)
            IF (I_ERROR.NE.0) THEN
               PRINT *,'err: glat, glong, i_err',GLAT,GLON,I_ERROR
            END IF
C           Convert back to Geographics as a sanity check
            CALL CONVERT_GEO_COORD(CGLAT,CGLON,R_HEIGHT_IN,GMLAT_,
     *                             GMLON_,2,I_ERROR)
            IF (I_ERROR.NE.0) THEN
               PRINT *,'err: mlat, mlong, i_err',CGLAT,CGLON,I_ERROR
            END IF
C
            WRITE (1,FMT=9000) R_HEIGHT_IN,GLAT,GLON,CGLAT,CGLON,
     *        GMLAT_,GMLON_
   10    CONTINUE
C
   20 CONTINUE
      CLOSE (1)
   30 CONTINUE
 9000 FORMAT (3f8.2,4F13.5)
      STOP
      END