SUBROUTINE SPSHUT

#include "com2d.h"
	COMMON/SPSH/ROCKALT(9),ROCK25(9),ROCK35(9),SSEM(9),T4EM(9),
     *  RKHCL25(Z$),RKHCL35(Z$)
	COMMON/ALOX/ALO25(9),ALO35(9),RKALO25(Z$),RKALO35(Z$),
     *  AERAL(L$,Z$)
	COMMON/SSCON/ISSK,IT4K,IT4V,CON25K,CON35K,CON35V,CONAO25K,
     *  CONAO35K,CONAO35V,SSEMAO,T4EMAO

	dimension hslss(31),hslt4(31)
	data hslss/0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,2.,
     *  3.,4.,5.,9.,1.,0.0,2.,5.,6.,6.,8.,7.,7.,7.,7.,7.,7.,7.,7./
	data hslt4/6., 8., 9., 9., 7., 8., 9., 7., 6., 5., 3., 5.,
     *  5.,3.,7.,1.,1.,3., 2.,5.,3.,2.,3.,1.,5.,3.,3.,3.,3.,3.,3./

c  Start launch rate in 1970 and go through 2000
	ilaunch=iyr-34

	if(ilaunch .le. 0)ilaunch = 0

	SSK=ISSK
	T4K=IT4K
	T4V=IT4V
	SSTOTAL=0.0E0
	T4TOTAL=0.0E0
	DO 5100 I=1,9
	   if(ilaunch .eq. 0)then
	      ROCK25(I)=0.0E0
	      ROCK35(I)=0.0E0
	      GO TO 5100
	   endif
	ROCK25(I)=CON25K*(HSLSS(ILAUNCH)*SSEM(I) + 
     *  (2./3.)*HSLT4(ILAUNCH)*T4EM(I))
	ROCK35(I)=CON35K*(HSLSS(ILAUNCH)*SSEM(I) + 
     *  (2./3.)*HSLT4(ILAUNCH)*T4EM(I)) +
     *  CON35V*((1./3.)*HSLT4(ILAUNCH)*T4EM(I))
	SSTOTAL=SSTOTAL+SSEM(I)
	T4TOTAL=T4TOTAL+T4EM(I)
5100	CONTINUE


	DO 5200 I=1,9
	   if(ilaunch .eq. 0)then
	      ALO25(I)=0.0E0
	      ALO35(I)=0.0E0
	      GO TO 5200
	   endif
	ALO25(I)=CONAO25K*(HSLSS(ILAUNCH)*SSEMAO*(SSEM(I)/SSTOTAL) + 
     *  (2./3.)*HSLT4(ILAUNCH)*T4EMAO*(T4EM(I)/T4TOTAL))
	ALO35(I)=CONAO35K*(HSLSS(ILAUNCH)*SSEMAO*(SSEM(I)/SSTOTAL) + 
     *  (2./3.)*HSLT4(ILAUNCH)*T4EMAO*(T4EM(I)/T4TOTAL)) +
     *  CONAO35V*((1./3.)*HSLT4(ILAUNCH)*T4EMAO*(T4EM(I)/T4TOTAL))
5200	CONTINUE


	DO 100 IK=1,Z$
	   RKHCL25(IK)=0.0E0
	   RKALO25(IK)=0.0E0
	   DO 200 IRK=1,8
	   IF(ZKM(12,IK) .GT. ROCKALT(IRK) .AND.
     *      ZKM(12,IK) .LT. ROCKALT(IRK+1))THEN
	      IRKUSE=IRK
	      GO TO 300
	   ENDIF
 200	   CONTINUE
 300	   IF(ZKM(12,IK) .GE. 52.5)GO TO 1100
	   IF(ZKM(12,IK) .LE. 12.5)GO TO 1100
	   DRKA=ROCKALT(IRKUSE+1)-ROCKALT(IRKUSE)
	   FRAB=(ZKM(12,IK)-ROCKALT(IRKUSE))/DRKA
	   FRBE=(ROCKALT(IRKUSE+1)-ZKM(12,IK))/DRKA
	   RKHCL25(IK)=FRBE*ROCK25(IRKUSE) + FRAB*ROCK25(IRKUSE+1)
	   RKALO25(IK)=FRBE*ALO25(IRKUSE) + FRAB*ALO25(IRKUSE+1)
 1100	   continue
	   if(iday360.eq.75)then
	   print *,' AT 25 IK=',IK,' HCL=',RKHCL25(IK),' CM-3 S-1'
	   print *,' AT 25 IK=',IK,' Al2O3=',RKALO25(IK),' CM-3 S-1'
	   endif
 100	CONTINUE

	DO 2100 IK=1,Z$
	   RKHCL35(IK)=0.0E0
	   RKALO35(IK)=0.0E0
	   DO 2200 IRK=1,8
	   IF(ZKM(13,IK) .GT. ROCKALT(IRK) .AND.
     *      ZKM(13,IK) .LT. ROCKALT(IRK+1))THEN
	      IRKUSE=IRK
	      GO TO 2300
	   ENDIF
 2200	   CONTINUE
 2300	   IF(ZKM(13,IK) .GE. 52.5)GO TO 4100
	   IF(ZKM(13,IK) .LE. 12.5)GO TO 4100
	   DRKA=ROCKALT(IRKUSE+1)-ROCKALT(IRKUSE)
	   FRAB=(ZKM(13,IK)-ROCKALT(IRKUSE))/DRKA
	   FRBE=(ROCKALT(IRKUSE+1)-ZKM(13,IK))/DRKA
	   RKHCL35(IK)=FRBE*ROCK35(IRKUSE) + FRAB*ROCK35(IRKUSE+1)
	   RKALO35(IK)=FRBE*ALO35(IRKUSE) + FRAB*ALO35(IRKUSE+1)
 4100	   continue
	   if(iday360.eq.75)then
	   print *,' AT 35 IK=',IK,' HCL=',RKHCL35(IK),' CM-3 S-1'
	   print *,' AT 35 IK=',IK,' Al2O3=',RKALO35(IK),' CM-3 S-1'
	   endif
 2100	CONTINUE
	SAVE
	RETURN
	END