C     Print SUGRA parameters and results
C
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/sslun.inc"
#include "isajet/sugxin.inc"
#include "isajet/sugmg.inc"
#include "isajet/sugpas.inc"
#include "isajet/sspar.inc"
#include "isajet/sugnu.inc"
      REAL PI,GPX,SIN2W,ALEMI,AS
      INTEGER IMODEL,J,K
C
C          Entry
C
      PI=4.*ATAN(1.)
      GPX=SQRT(.6)*GSS(1)
      SIN2W=GPX**2/(GSS(2)**2+GPX**2)
      ALEMI=4*PI/GSS(2)**2/SIN2W
      AS=GSS(3)**2/4./PI
C
C          Print inputs and couplings
C          
      IF (IMODEL.EQ.1) THEN
        WRITE(LOUT,1000) XSUGIN
1000    FORMAT(//' ISASUGRA input:'/
     $  ' M_0,  M_(1/2),  A_0,  tan(beta),  sgn(mu),  M_t ='
     $  /4F10.3,2X,F6.1,F10.3)
        WRITE(LOUT,1001) MGUTSS,GGUTSS,AGUTSS
1001    FORMAT(/' ISASUGRA unification:'/' M_GUT      =',E10.3,
     $  '   g_GUT          =',F5.3,3X,'   alpha_GUT =',F5.3)
      ELSE IF (IMODEL.EQ.2) THEN
        WRITE(LOUT,1002) (XGMIN(J),J=1,7)
1002    FORMAT(//' GMSB model input:'/
     $  ' Lambda,  M_mes,  N_5,  tan(beta),  sgn(mu),  M_t,  C_grav='
     $  /2E10.3,2F10.3,2X,F6.1,F10.3,1X,E10.3)
        WRITE(LOUT,1020) (XGMIN(J),J=8,14)
1020    FORMAT(/' GMSB2 model input:'/
     $  ' Rsl,    dmH_d^2,   dmH_u^2,     d_Y,     N5_1,  N5_2,  N5_3='
     $  /F5.3,1X,E10.3,1X,E10.3,1X,E10.3,2X,3F7.3)
        WRITE(LOUT,1003) AMGVSS
1003    FORMAT(//' M(gravitino)=',E10.3)
      END IF
C
      WRITE(LOUT,1004) ALEMI,SIN2W,AS
1004  FORMAT(' 1/alpha_em =',F8.2,2X,
     $'   sin**2(thetaw) =',F6.4,2X,'   alpha_s   =',F5.3)
      WRITE(LOUT,1005) GSS(7),GSS(8),GSS(9)
1005  FORMAT(' M_1        =',F8.2,2X,
     $'   M_2            =',F8.2,'   M_3       =',F8.2)
      WRITE(LOUT,1006) MU,B,HIGFRZ
1006  FORMAT(' mu(Q)      =',F8.2,2X,
     $'   B(Q)           =',F8.2,'   Q         =',F8.2)
      WRITE(LOUT,1007) GSS(13),GSS(14)
1007  FORMAT(' M_H1^2     =',E10.3,'   M_H2^2         =',E10.3)
C     
C          Write out non-universal GUT scale parameters
C
      IF (XNUSUG(1).LT.1.E19.OR.XNUSUG(2).LT.1.E19.OR.XNUSUG(3)
     $.LT.1.E19) THEN
        WRITE(LOUT,1010) XNUSUG(1),XNUSUG(2),XNUSUG(3)
1010    FORMAT(/' M_1(GUT)= ',F8.2,' M_2(GUT)= ',F8.2,' M_3(GUT)= ',
     $  F8.2) 
      END IF
      IF (XNUSUG(4).LT.1.E19.OR.XNUSUG(5).LT.1.E19.OR.XNUSUG(6)
     $.LT.1.E19) THEN
        WRITE(LOUT,1011) XNUSUG(4),XNUSUG(5),XNUSUG(6)
1011    FORMAT(/' A_tau(GUT)= ',F8.2,' A_b(GUT)= ',F8.2,' A_t(GUT)= ',
     $  F8.2)
      END IF
      IF (XNUSUG(7).LT.1.E19.OR.XNUSUG(8).LT.1.E19) THEN
        WRITE(LOUT,1012) XNUSUG(7),XNUSUG(8)
1012    FORMAT(/' m_Hd(GUT)= ',F8.2,'    m_Hu(GUT)= ',F8.2)
      END IF
      IF (XNUSUG(9).LT.1.E19.OR.XNUSUG(10).LT.1.E19) THEN
        WRITE(LOUT,1013) XNUSUG(9),XNUSUG(10)
1013    FORMAT(/' m_eR(GUT)= ',F8.2,' m_eL(GUT)= ',F8.2)
      END IF
      IF (XNUSUG(11).LT.1.E19.OR.XNUSUG(12).LT.1.E19.OR.XNUSUG(13)
     $.LT.1.E19) THEN
        WRITE(LOUT,1014) XNUSUG(11),XNUSUG(12),XNUSUG(13)
1014    FORMAT(' m_dR(GUT)= ',F8.2,' m_uR(GUT)= ',F8.2,' m_uL(GUT)=',
     $  F8.2)
      END IF
      IF (XNUSUG(14).LT.1.E19.OR.XNUSUG(15).LT.1.E19) THEN
        WRITE(LOUT,1015) XNUSUG(14),XNUSUG(15)
1015    FORMAT(/' m_tauR(GUT)= ',F8.2,' m_tauL(GUT)= ',F8.2)
      END IF
      IF (XNUSUG(16).LT.1.E19.OR.XNUSUG(17).LT.1.E19.OR.XNUSUG(18)
     $.LT.1.E19) THEN
        WRITE(LOUT,1016) XNUSUG(16),XNUSUG(17),XNUSUG(18)
1016    FORMAT('   m_bR(GUT)= ',F8.2,' m_tR(GUT)= ',F8.2,' m_tL(GUT)=',
     $  F8.2)
      END IF
C
C          Print mass spectrum from ISASUGRA 
C
      WRITE(LOUT,2000) MSS(1),MSS(2),MSS(3),MSS(4),MSS(5),MSS(10),
     $MSS(11),MSS(12),MSS(13),MSS(14),MSS(17),MSS(18),MSS(16),
     $MSS(21),MSS(22),MSS(23),MSS(24),MSS(25),MSS(26),MSS(27),
     $MSS(28),MSS(29),MSS(30),MSS(31),MSS(32)
2000  FORMAT(/' ISAJET masses (with signs):'/
     $' M(GL)  =',F9.2/
     $' M(UL)  =',F9.2,'   M(UR)  =',F9.2,'   M(DL)  =',F9.2,
     $'   M(DR) =',F9.2/
     $' M(B1)  =',F9.2,'   M(B2)  =',F9.2,'   M(T1)  =',F9.2,
     $'   M(T2) =',F9.2/
     $' M(SN)  =',F9.2,'   M(EL)  =',F9.2,'   M(ER)  =',F9.2/
     $' M(NTAU)=',F9.2,'   M(TAU1)=',F9.2,'   M(TAU2)=',F9.2/
     $' M(Z1)  =',F9.2,'   M(Z2)  =',F9.2,'   M(Z3)  =',F9.2,
     $'   M(Z4) =',F9.2/
     $' M(W1)  =',F9.2,'   M(W2)  =',F9.2/
     $' M(HL)  =',F9.2,'   M(HH)  =',F9.2,'   M(HA)  =',F9.2,
     $'   M(H+) =',F9.2)
      WRITE(LOUT,2001) THETAT,THETAB,THETAL,ALFAH
2001  FORMAT(/,' theta_t=',F9.4,'   theta_b=',F9.4,
     $'   theta_l=',F9.4,'   alpha_h=',F9.4)
C
C     Write out chargino /neutralino masses/eigenvectors
C
      WRITE(LOUT,3100) AMZ1SS,AMZ2SS,AMZ3SS,AMZ4SS
3100  FORMAT(/' NEUTRALINO MASSES (SIGNED) =',4F10.3)
      DO 100 J=1,4
        WRITE(LOUT,3200) J,(ZMIXSS(K,J),K=1,4)
3200    FORMAT(' EIGENVECTOR ',I1,'       =',4F10.5)
100   CONTINUE
      WRITE(LOUT,3300) AMW1SS,AMW2SS
3300  FORMAT(/' CHARGINO MASSES (SIGNED)  =',2F10.3)
      WRITE(LOUT,3400) GAMMAL,GAMMAR
3400  FORMAT(' GAMMAL, GAMMAR             =',2F10.5/)

C
C          Print ISAJET MSSMi equivalent input
C
      WRITE(LOUT,3000)
3000  FORMAT(/' ISAJET equivalent input:')
      WRITE(LOUT,3001) MSS(1),MU,MSS(31),XSUGIN(4)
3001  FORMAT(' MSSMA: ',4F8.2)
      WRITE(LOUT,3002) SQRT(GSS(19)),SQRT(GSS(17)),SQRT(GSS(18)),
     $SQRT(GSS(16)),SQRT(GSS(15))
3002  FORMAT(' MSSMB: ',5F8.2)
      WRITE(LOUT,3003) SQRT(GSS(24)),SQRT(GSS(22)),SQRT(GSS(23)),
     $SQRT(GSS(21)),SQRT(GSS(20)),GSS(12),GSS(11),GSS(10)
3003  FORMAT(' MSSMC: ',8F8.2)
      RETURN
      END