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$)

	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
c	   PRINT *,' frab=',frab,' frbe=',frbe,' drka=',drka
c	   print *,' zkm(12,ik)=',zkm(12,ik),' ik=',ik
c	   print *,' rock25 irkuse=',rock25(irkuse),' rk+1=',
c     *   rockalt(irkuse+1)
c	   print *,' AT 25 IK=',IK,' HCL=',RKHCL25(IK),' CM-3 S-1'
 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
c	   PRINT *,' AT 35 IK=',IK,' HCL=',RKHCL35(IK),' CM-3 S-1'
 2100	CONTINUE
	SAVE
	RETURN
	END