SUBROUTINE SPDR(IJC,IKC)

#include "com2d.h"
#include "comphot.h"

C     SOLAR PHOTODISSOCIATION 
C     RATES FOR CHEMISTRY MODULE IN POINT CALL FORMAT.

C     NO CALCULATION OF J COEFFICIENTS OR REDUCED FLUXES IF THE SUN IS
C     NOT UP AT THE CURRENT POINT.
      IF(ZGRZ(IJC).LE.0.0)THEN
        DO 200 JL=1,IL$
200	RFLUX(JL,IJC,IKC)=0.0
        DO 220 JP=1,PH$
220	J(JP,IJC,IKC)=0.0
	SAVE
        RETURN
      END IF

C     SOLAR FLUX AND ABSORPTION CROSS-SECTIONS (PHOTODISSOCIATION DATA)
C     FOR SPECIES IN OPTICAL DEPTH CALCULATION ARE DEFINED AT IL$
C     WAVELENGTHS. 

C     INITIALIZE ARRAYS FOR J-COEFFICIENTS AND 
	DO 320 JP=1,PH$
	J(JP,IJC,IKC)=0.0
320 	CONTINUE


C     J COEFFICIENTS FOR EACH PROCESS WHEN GRAZING HEIGHT > 0KM.
C     COMPUTE EACH J-COEFFICIENT OVER ITS DEFINED WAVELENGTH RANGE.
C     I IS THE WAVELENGTH INDEX; J, THE CROSS-SECTION INDEX. DO NOT
C     COMPUTE J-COEFFICIENTS CORRESPONDING TO OMITTED PROCESSES.
	DO 400 JL=1,IL$
	DO 400 JP=1,PH$
	J(JP,IJC,IKC)=J(JP,IJC,IKC) + XSECT(JP,JL)*RFLUX(JL,IJC,IKC)
c	if(ijc.eq.9 .and. jp.eq.30 .and. ikc.eq.4)then
c	print *,' jl=',jl,' j=',j(jp,ijc,ikc),
c     *  ' xsect(1,jl)=',xsect(jp,jl),' rflux=',rflux(jl,ijc,ikc)
c	endif
400	CONTINUE

	j(1,ijc,ikc)=o2jint(ijc,ikc)
c J(O2) from look-up table - 4/5/95

c        if(ijc.eq.9 .and. ikc.eq.46)then
c           write(6,1101)(j(30,ijc,ikout),ikout=1,z$)
c 1101      format(' j1301=',/,1p6e12.6,/,1p6e12.6,/,1p6e12.6,/,
c     *  1p6e12.6,/,1p6e12.6,/,1p6e12.6,/,1p6e12.6)
c           endif
	SAVE
	RETURN
	END