!	SUBROUTINE MAPINI (PARMAP,LCODE)GC_TITLE  MAPINI  Query for a map projection; and also initiate  MAPROJ.C_ARGUMENTS<	REAL*4 PARMAP(8)   ! [O] Array of map definition parameters(	INTEGER*4 LCODE(2) ! [O] Line type codeC_DESCRIPTION LC The subroutine MAPINI prompts for parameters necessary to initialize a mapMC projection.  After calls are made to the user interface MAPINI passes them KC back to the main routine. MAPINI uses TAE as the user interface and will GC prompt if one of the critical variables has been left blank.  MAPINI OC initializes the MAP projection routine MAPROJ.  See module MAPROJ for further
C details.C_VARIABLES	INCLUDE 'RAD$INC:XXTAE.INC'"C					Uses: LOUT.  Modifies: none."        INCLUDE 'RAD$INC:GRID.INC'"C					Uses: KPROJ. Modifies: none.%        INCLUDE 'RAD$INC:TAEBAS.INC'	"C 					TAE interface link library.C_CALLS  Tae routines3C_HISTORY  86May30  R.Gurule  USGS, Flagstaff 860014C	   87NOV17  R.Gurule, create includes for commons.C_END
	LOGICAL PFLG	REAL*4 DFAULT(8),DUM,ATMP(3)	CHARACTER*4 JECT,PROJ(7)?	DATA DFAULT /3393.4, .101929, .017453, 2.54E-5  !MARS DEFAULTS-     &	,25.0E6, 3*0./	 				!MERCATOR DEFAULTSC        DATA PROJ/'MERC','LAMB','POLA','SIMP','SINU','VIEW','DEGR'/      8        CALL XRSTR(BLOCK,'KPROJ',1,JECT,ICNT,NSTR,ISTAT)8        CALL XRREAL(BLOCK,'MPINIX',4,PARMAP,ICNT,ISTAT) 8        CALL XRREAL(BLOCK,'MPINIY',3,ATMP(1),ICNT,ISTAT)	PARMAP(5) = ATMP(1)	PARMAP(7) = ATMP(2)	PARMAP(8) = ATMP(3)         DO J=1,7          IP = J&	  IF(JECT(1:4) .EQ. PROJ(J)(1:4))THEN	    KPROJ=J  	    GOTO 50	  ENDIF
        ENDDO6   50   CALL XRREAL(BLOCK,'CLON',1,ATMP(1),ICNT,STAT)		PARMAP(6)=ATMP(1)	N!-------IF projection = view then rearrange the order or parmap so that x(6)=xH!       x(7) = y, x(8) = x, instead of X(7) = X, X(8) = Y, and X(6) = Z.	IF(KPROJ.EQ.6)THEN	  TMPCOR = PARMAP(6)	  PARMAP(6) = PARMAP(7)	  PARMAP(7) = PARMAP(8)	  PARMAP(8) = TMPCOR	ENDIF	/        IF(PARMAP(6).NE.0.0.and.KPROJ.NE.6)THEN!        WRITE(LOUT,1003)PARMAP(6), 1003   FORMAT(' Central Meridian = ',F12.3)
        ENDIF/        IF(PARMAP(7).NE.0.0.AND.KPROJ.EQ.2)THEN-          WRITE(LOUT,1004)PARMAP(3),PARMAP(4)C 1004     FORMAT(' THE 1st and 2nd Standard parallels are ',2F12.3)
        ENDIF1        IF(KPROJ.EQ.6 .AND. PARMAP(8).NE.0.0)THEN7          WRITE(LOUT,1005)PARMAP(6),PARMAP(7),PARMAP(8)B 1005     FORMAT(' The coordinates of the viewer''s eye =',3F12.1)
        ENDIF        IF(KPROJ.NE.0)THEN?          CALL MAPROJ(-KPROJ,ALIN1,ALIN2,PARMAP(1),PROJ(KPROJ))3          CALL MAPROJ(0,ALIN1,ALIN2,PARMAP(1),PROJ)           LCODE(1)= JIFIX(ALIN1)           LCODE(2)= JIFIX(ALIN2)"          WRITE(LOUT,1000)PROJ(IP)( 1000	  FORMAT(' The projection is ',A4)+          WRITE(LOUT,1001)(PARMAP(I),I=1,4): 1001     FORMAT(' Maproj initialization parameters are ',&     1  /,' Semi major axis = ',F12.4,"     2  /,' Eccentricity = ',F9.8,?     3  /,' Angular units of long and lat in radians = ',F10.2,#     4  /,' Output units = ',F10.2)
        ENDIF        RETURN
      	END