      SUBROUTINE GSSTCD(Z00,SATCRI,ISST)
      COMMON IDATE(4),
     * GESHEM( 50 , 42 ),TSEA( 50 , 42 ),GRGRID( 50 , 42 )
      COMMON DPDPHI( 306 ),DPDLAM( 289 ),RELVOR( 7 ),ABSVOR( 7 ),
     1  Q( 289 ), TE( 289 , 7 ), DI( 289 , 7 ), ZE( 289 , 7 ),
     2  Z( 289 ),  Y( 289 , 7 ),  X( 289 , 7 ),  W( 289 , 7 ),
     3 QM( 289 ),TEM( 289 , 7 ),DIM( 289 , 7 ),ZEM( 289 , 7 ),
     5 GZ( 289 ), RQ( 289 , 6  ), RM( 289 , 6  ), RT( 289 , 6  )
      COMPLEX GZ,Q,TE,DI,ZE,Z,Y,X,W,QM,TEM,DIM,ZEM,RELVOR,
     1DPDLAM,DPDPHI,ABSVOR,RQ,RM,RT
      COMMON
     * PS( 50 ), DPHI( 50 ), DLAM( 50 ),
     * QF( 50 ),DPHIF( 50 ),DLAMF( 50 ),
     *  EF( 50 , 7 ), EEF( 50 , 7 ), CG( 50 , 7 ),
     *  TF( 50 , 7 ), TEF( 50 , 7 ),TAU( 50 , 7 ),
     *  AF( 50 , 7 ), ZEF( 50 , 7 ),  A( 50 , 7 ),
     *  BF( 50 , 7 ), DIF( 50 , 7 ),  B( 50 , 7 ),
     *  FF( 50 , 7 ),  UF( 50 , 7 ),  F( 50 , 7 ),
     *  GF( 50 , 7 ),  VF( 50 , 7 ),  G( 50 , 7 ),
     * UQF( 50 , 6 ), UQFF( 50 , 6 ),
     * VQF( 50 , 6 ), VQFF( 50 , 6 ),
     * RQF( 50 , 6 ), RQFF( 50 , 6 ), RTG( 50 , 6 )
      COMMON
     * U1( 50 ),V1( 50 ),T1( 50 ),Q1( 50 ),
     * U2( 50 ),V2( 50 ),T2( 50 ),Q2( 50 )
      COMMON/PLNCOM/PLN( 306 ),DER( 289 ),PLNWCS( 289 ),EPS( 306 ),
     1 COLRAD( 21 ),WGT( 21 ),WGTCS( 21 ),RCS2( 21 )
      DOUBLE PRECISION COLRAD,EPS
      COMMON/VERCOM/AM( 7 , 7 ),HM( 7 , 7 ),TM( 7 , 7 ),
     O              BM( 7 , 7 ),CM( 7 , 7 ),EKIN( 7 ),
     1 SI( 8 ),SL( 7 ),DEL( 7 ),RDEL2( 7 ),RMSDOT( 6 ),
     2 CI( 8 ),CL( 7 ),TOV( 7 ),   SV( 7 ),   RPI( 6 ),
     3 P1( 7 ),P2( 7 ), H1( 7 ),   H2( 7 ),RPIREC( 6 ),
     4 ROTSIN,ROTCOS,SEADRY,SL1KAP,C1,C2,RLRV,SL100K
C ----------------------MEMBER GCOM---------------------------------
C     RPI(K) = (SL(K+1)/SL(K))**RK  FROM SETSIG  K=1... 6
      DO 1 K=1, 6
1     RPIREC(K) = 1./RPI(K)
      DO 2 K=1, 7
2     RDEL2(K)=1./(2.*DEL(K))
      REWIND ISST
      READ(ISST) GRGRID
      GA2= 9.8 /( 6370000. * 6370000. )
      Z00=REAL(GZ(1))
      DO 4 LL=1, 17
      N=LL-2
      DO 4 I=1, 17
      N=N+1
      FACT=N*(N+1.)*GA2
      GZ((LL-1)* 17 +I)=CMPLX(
     1 FACT* REAL(GZ((LL-1)* 17 +I)),
     1 FACT*AIMAG(GZ((LL-1)* 17 +I)))
4     CONTINUE
C     NOW READ TSEA
      MONM1=IDATE(2)-1
      IF(MONM1.EQ.0) GO TO 55
      DO 5 I=1,MONM1
      READ(ISST)
5     CONTINUE
55    CONTINUE
      READ(ISST)TSEA
      REWIND ISST
C-------
      RLRV=-2.51E+10/(1.61*2.87E+6)
      CONST=6.11*EXP(-RLRV/273.16)/10.
      C1=CONST* 0.622
      C2=CONST*(-0.378)
C-------
      HBOUND= 287.05 *280.*DEL(1)/( 9.8 *SL(1))
      SEADRY=(2.*0.07*0.001)/HBOUND
      RK= 287.05 / 1005.
      SL1KAP=1./SL(1)**RK
      SL100K = (SL(1)/100.)**RK
      ROTSIN=SIN(20./57.295779)
      ROTCOS=COS(20./57.295779)
      WRITE(6,9877)HBOUND
9877  FORMAT(1H0,' HBOUND=',E12.3)
      SQR2=SQRT(2.)
      DO 444 LA=1, 42
      DO 444 LO=1, 50
444   GRGRID(LO,LA)=GRGRID(LO,LA)*SQR2/HBOUND
102   FORMAT(1H0,'GSSTCD COMPLETED FOR ',4I4)
      WRITE(6,102) IDATE
C(DEBUGLINE) VERSION CREATED JULY 8 1982 BY J.SELA,NMC.
      RETURN
      END
