SUBROUTINE BOUNDC

#include "com2d.h"


C SET UP BOUNDARY CONDITIONS FOR TIME-DEPENDENT SPECIES
        DO 200 II=1,31
        IUSE=II
        IF(IYR .GE. IYRBC(II) .AND. IYR .LT. IYRBC(II+1))GO TO 250
200     CONTINUE

250     DO 300 JJ=1,12

	IS=ISPBC(JJ)

        DELTABC=(BCTDINPUT(IS,IUSE+1)-BCTDINPUT(IS,IUSE))/5.
c	if(iyr .gt. 90)then
c	if(is.eq.34)print *,' is=',is,' iuse=',iuse,' deltabc=',
c     *  deltabc
c	if(is.eq.34)print *,' bciuse+1=',bctdinput(is,iuse+1),
c     *  ' bciuse=',bctdinput(is,iuse)
c	endif
        TEMP1USE=IYR + (DAY360/360.)
        TEMP2USE=IYRBC(IUSE)

	IF(LBCMRTD(IS))THEN
C THIS SECTION IS FOR MIXING RATIO BOUNDARY CONDITIONS

	MIXRATBC=(TEMP1USE-TEMP2USE)*DELTABC + 
     *  BCTDINPUT(IS,IUSE)
c	if(iyr .gt. 90)then
c	if(is.eq.34)print *,' mixrat=',mixratbc,' iyr=',iyr,
c     *  ' day360=',day360,' temp1use=',temp1use,' temp2use=',
c     *  temp2use,' iuse=',iuse
c	endif
	DO 500 IL=1,L$
	BVAL(IS,IL)=TDLAT(IS,IL)*MIXRATBC*1.E-9
c	if(iyr.gt.90 .and. il.eq.9)then
c	if(is.eq.34)print *,' bval=',bval(is,il),' il=',il,
c     *  ' tdlat=',tdlat(is,il),' mixrat=',mixratbc
c	endif
500	CONTINUE

	ENDIF


	IF(.NOT.LBCMRTD(IS))THEN
C THIS SECTION IS FOR FLUX BOUNDARY CONDITIONS

	FLUXBC=(TEMP1USE-TEMP2USE)*DELTABC + 
     *  BCTDINPUT(IS,IUSE)

C NORMALIZE BY MOLECULAR WEIGHT OF MOLECULE
	FLUXBC=FLUXBC*(FRMOLWT(1)/FRMOLWT(JJ))

	DO 1500 IL=1,L$
C USE LATITUDE DEPENDENCE READ IN FOR CFC11 IN ALL CFCs
	BVAL(IS,IL)=TDLAT(ISPBC(1),IL)*FLUXBC
1500	CONTINUE

	ENDIF

300	CONTINUE

	SAVE
	RETURN
	END