! **************************************************************************************
!   This computer software was developed by Dr. Yang Zhang and her research group      *
!   at North Carolina State University (NCSU) with support  from the NSF Career Award  *
!   No. Atm-0348819, and the Memorandum of Understanding between the                   *
!   U.S. Environmental Protection Agency (EPA) and the U.S. Department of              *
!   Commerce's National Oceanic and Atmospheric Administration (NOAA)		       *
!   and under agreement number DW13921548, and the U.S. EPA/Office of                  *
!   Air Quality Planning & Standards via RTI International contract #4-321-0210288.    *
!										       *
!   NEITHER ANY COSPONSORS, NCSU, NOR ANY PERSON ACTING ON BEHALF                      *
!   OF ANY OF THEM MAKES ANY WARRANTY OR REPRESENTATION                                *
!   WHATSOEVER, EXPRESS OR IMPLIED, OR ASSUMES ANY LIABILITY FOR                       *
!   THE USE OF THIS SOFTWARE.  THIS SOFTWARE OR PART OF IT MAY BE                      *
!   COPYRIGHTED AND IS PERMITTED BY ORIGINAL CODE DEVELOPERS FOR                       *
!   NONPROFIT USE AND SUBJECTED TO RESTRICTIONS                       		       *
!                                                                                      *
!  Contact information: 							       *
!  Dr. Yang Zhang		                                                       * 
!  Principal Investigator						               *
!  Department of Marine, Earth, and Atmospheric Sciences 			       *
!  North Carolina State University						       *
!  Campus Box 8208 								       *
!  Room 5151, Jordan Hall, 2800 Faucette Drive 					       *
!  Raleigh, NC 27695-8208, USA 							       *
!  Tel:  (919) 515-9688 (Office) 						       *
!  Fax:  (919) 515-7802								       *
!  E-Mail:  yang_zhang@ncsu.edu 						       *
!                                                                                      *
! **************************************************************************************
!  FUNCTION: PROVIDE BOUNDARY VALUES FOR CB05 MADRID extension GAS-PHASE SPECIES       *
!  PRECONDITION REQUIRED: use for CB05 MADRID extension gas-phase mechanism            *
!  RETURN VALUES:								       *			         
!  KEY SUBROUTINES AND FUNCTIONS CALLED: None                                          *
!  REVISION HISTORY:                                                                   *
!         This code was based on module_cb05_initmixrats.F,                            *
!               developed by J.P. Huang and Y. Zhang, NCSU, Mar.-Oct. 2006             *
!         Revised by Ying Pan and Yang Zhang, NCSU, Nov. 2007-Nov. 2008	 	       *
!               to couple MADRID with CB05 MADRID extension gas-phase mechanism        * 
!         Revised by Yang Zhang, Xiao-Ming Hu, and Ying Pan, NCSU, Sept.-Nov., 2008    *
!         	Code cleaned up for NOAA WRF/Chem repository checkin     	       *               
!         Revised by Ying Pan and Yang Zhang, NCSU, Sep. 2009                          *
!               to transfer the code to WRF/Chem v3.1.1                                *
!***************************************************************************************


!-----------------------------------------------------------------------
      subroutine bdy_chem_value_cb05 (id_bdy, chem_bv,kk, nch,config_flags,numgas )

        use module_configure, only:  grid_config_rec_type
	use module_input_chem_data, only:  bdy_chem_value

	implicit none

! arguments
	REAL,    INTENT(OUT)  :: chem_bv    ! boundary value for chem(-,-,-,nch)
!	REAL,    INTENT(IN)   :: z          ! height
	INTEGER, INTENT(IN)   :: nch        ! index number of chemical species
        INTEGER, INTENT(IN)   :: numgas     ! index number of last gas species
        INTEGER, INTENT(IN)   :: kk         ! level
        INTEGER, INTENT(IN)   :: id_bdy     ! id of bdy: 1 = YS, 2 = YE, 3=XS, 4 =XE
        TYPE(grid_config_rec_type), INTENT(IN) :: config_flags
! local variables
	real chem_bv_ald, chem_bv_hc3, chem_bv_hc5,   &
	     chem_bv_hc8, chem_bv_ket, chem_bv_oli,   &
	     chem_bv_olt, chem_bv_ora2
	real, parameter :: chem_bv_def = 1.0e-20
      character (len=100)    :: record 
      character (len=20)     ::  chem_var
      real                  :: cmaq_bc9(98,4,34)
      integer               :: iv, i,j,k

! The following boundary values were taken from CMAQ, July, 2001 episode
!
! for NO2  
      data (cmaq_bc9(1,1,k),k=1,34)&
         / 9.9999997E-06,   9.9999997E-06,   9.9999997E-06,   &
       9.9999997E-06,  9.9999997E-06,   9.6428539E-06,   9.2857072E-06,&
       7.4999953E-06,  5.7142838E-06,   5.7142838E-06, 24*0./
      data (cmaq_bc9(1,2,k),k=1,34)&
         / 9.9999997E-06,   9.9999997E-06,   9.9999997E-06,   &
       9.9999997E-06, 9.9999997E-06,   9.6428539E-06,   9.2857072E-06,& 
       7.4999953E-06, 5.7142838E-06,   5.7142838E-06, 24*0./
      data (cmaq_bc9(1,3,k),k=1,34)&
          /1.6700001E-04,  1.6700001E-04,  1.6700001E-04,&
          1.6700001E-04,  1.6700001E-04,  1.6403571E-04, 1.6107140E-04,&
          1.4625001E-04,  1.3142861E-04,  1.3142861E-04, 8.3999999E-05,&
          8.3999999E-05,  8.3999999E-05,  5.6000019E-05, 5.6000019E-05,&
          2.5454530E-05,  2.5454530E-05,  2.5454530E-05, 7.4999953E-06,&
          5.7142838E-06,  5.7142838E-06, 13*0./
      data (cmaq_bc9(1,4,k),k=1,34)&
          /1.6700001E-04,  1.6700001E-04,  1.6700001E-04,& 
          1.6700001E-04,  1.6700001E-04,  1.6403571E-04, 1.6107140E-04,&
          1.4625001E-04,  1.3142861E-04,  1.3142861E-04, 8.3999999E-05,&
          8.3999999E-05,  8.3999999E-05,  5.6000019E-05, 5.6000019E-05,&
          2.5454530E-05,  2.5454530E-05,  2.5454530E-05, 7.4999953E-06,&
          5.7142838E-06,  5.7142838E-06, 13*0./

! for NO
      data (cmaq_bc9(2,1,k),k=1,34)&
          / 9.9999997E-06,9.9999997E-06,9.9999997E-06, &
          9.9999997E-06, 9.9999997E-06,9.6428539E-06, 9.2857072E-06,&
          7.4999953E-06, 5.7142838E-06,5.7142838E-06,24*0./ 
      data (cmaq_bc9(2,2,k),k=1,34)&
         /34*0./
      data (cmaq_bc9(2,3,k),k=1,34)&
         / 5*8.2999999E-05, 8.1535698E-05, 8.0071397E-05,&
           7.2749986E-05, 6.5428569E-05, 6.5428569E-05, 4.2000000E-05,&
           4.2000000E-05, 4.2000000E-05, 2.8000009E-05, 2.8000009E-05,&
           1.2727260E-05, 1.2727260E-05, 1.2727260E-05, 16*0./
      data (cmaq_bc9(2,4,k),k=1,34)&
         / 5*8.2999999E-05, 8.1535698E-05, 8.0071397E-05,& 
           7.2749986E-05,   2*6.5428569E-05,3*4.2000000E-05,& 
           2*2.8000009E-05, 3*1.2727260E-05, 16*0./

! for O
      data ((cmaq_bc9(3,i,k),k=1,34),i=1,4)&
      /136*0./

! for O3 (CAM profile corporated)
      data (cmaq_bc9(4,1,k),k=1,34)&
      /2*2.9999999E-02, 3.0714281E-02,&
      3.1785712E-02, 3.2857139E-02, 3.4107145E-02, 0.034610,&
      0.034933, 0.035268, 0.035629, 0.035993,&
      0.036361, 0.036940, 0.037728, 0.038633,&
      0.039955, 0.041303, 0.042877, 0.044563,&
      0.046555, 0.048089, 0.049924, 0.052937,&
      0.056136, 0.059363, 0.062810, &
      0.064494, 0.065790, 0.069328, 0.074218, &
      0.086511, 0.100520, 0.117374, 0.182958  /
      data (cmaq_bc9(4,2,k),k=1,34)&
      /2*2.9999999E-02, 3.0714281E-02,&
      3.1785712E-02, 3.2857139E-02, 3.4107145E-02, 0.034610,&
      0.034933, 0.035268, 0.035629, 0.035993,&
      0.036361, 0.036940, 0.037728, 0.038633,&
      0.039955, 0.041303, 0.042877, 0.044563,&
      0.046555, 0.048089, 0.049924, 0.052937,&
      0.056136, 0.059363, 0.062810, &
      0.064494, 0.065790, 0.069328, 0.074218, &
      0.086511, 0.100520, 0.117374, 0.182958  /
      data (cmaq_bc9(4,3,k),k=1,34)&
      /5*3.5000000E-02,3.5178576E-02,0.034610,&
       0.034933,0.035268,0.035629,0.035993,&
       0.036361,0.036940,0.037728,0.038633,&
       0.039955,0.041303,0.042877,0.044563,&
       0.046555,0.048089,0.049924,0.052937,&
       0.056136, 0.059363, 0.062810, &
       0.064494, 0.065790, 0.069328, 0.074218, &
       0.086511, 0.100520, 0.117374, 0.182958 /
      data (cmaq_bc9(4,4,k),k=1,34)&
      /5*3.5000000E-02,3.5178576E-02,0.034610,&
       0.034933,0.035268,0.035629,0.035993,&
       0.036361,0.036940,0.037728,0.038633,&
       0.039955,0.041303,0.042877,0.044563,&
       0.046555,0.048089,0.049924,0.052937,&
       0.056136, 0.059363, 0.062810, &
       0.064494, 0.065790, 0.069328, 0.074218, &
       0.086511, 0.100520, 0.117374, 0.182958 /

! for no3
      data ((cmaq_bc9(5,i,k),k=1,34),i=1,4)&
      /136*0./

! for o1d
      data ((cmaq_bc9(6,i,k),k=1,34),i=1,4)&
      /136*0./

! for oh
      data ((cmaq_bc9(7,i,k),k=1,34),i=1,4)&
      /136*0./ 

! for ho2
      data ((cmaq_bc9(8,i,k),k=1,34),i=1,4)&
      /136*0./

! for n2o5
      data ((cmaq_bc9(9,i,k),k=1,34),i=1,4)&
      /136*0./

! for hno3
      data (cmaq_bc9(10,1,k),k=1,34)&
      /25*4.9999999E-05,3*9.9999997E-05,1.2500001E-04,5*1.5000001E-04/
      data (cmaq_bc9(10,2,k),k=1,34)&
      /25*4.9999999E-05,3*9.9999997E-05,1.2500001E-04,5*1.5000001E-04/
      data (cmaq_bc9(10,3,k),k=1,34)&
      /18*4.9999999E-05,2*5.1851850E-05, 5.5370372E-05,&
        5.8888891E-05, 6.3518521E-05, 2*6.8148147E-05,&
        3*8.5000000E-05,9.2499999E-05,5*9.9999997E-05/
      data (cmaq_bc9(10,4,k),k=1,34)&
      /18*4.9999999E-05,2*5.1851850E-05, 5.5370372E-05,&
        5.8888891E-05, 6.3518521E-05, 2*6.8148147E-05,&
        3*8.5000000E-05,9.2499999E-05,5*9.9999997E-05/

! for hono
     data ((cmaq_bc9(11,i,k),k=1,34),i=1,4)&
      /136*9.9999997E-10/

! for pna
     data ((cmaq_bc9(12,i,k),k=1,34),i=1,4)&
      /136*1.9999999E-09/

! for h2o2
     data (cmaq_bc9(13,1,k),k=1,34)&
      /18*2.0000001E-03,2*1.9537040E-03, 1.8657411E-03,&
       1.7777780E-03, 1.6620370E-03, 1.5462959E-03,    &
       1.5462959E-03, 1.2500000E-03, 1.2500000E-03,    &
       1.2500000E-03, 1.1250000E-03, 5*1.0000000E-03/
     data (cmaq_bc9(13,2,k),k=1,34)&
      /18*2.0000001E-03,2*1.9537040E-03, 1.8657411E-03,&
       1.7777780E-03, 1.6620370E-03, 1.5462959E-03,    &
       1.5462959E-03, 1.2500000E-03, 1.2500000E-03,    &
       1.2500000E-03, 1.1250000E-03, 5*1.0000000E-03/
     data (cmaq_bc9(13,3,k),k=1,34)&
      /5*1.0000000E-03, 1.0178576E-03,1.0357150E-03,&
      1.1250004E-03,1.2142860E-03,1.2142860E-03,3*1.5000000E-03,&
      2*1.3333330E-03,3*1.1515151E-03,2*9.8148151E-04,&
      9.4629632E-04,9.1111113E-04,8.6481485E-04,2*8.1851851E-04,&
      3*5.0000002E-04, 3.5000002E-04,5*1.9999999E-04/
     data (cmaq_bc9(13,4,k),k=1,34)&
      /5*1.0000000E-03, 1.0178576E-03,1.0357150E-03,&
      1.1250004E-03,1.2142860E-03,1.2142860E-03,3*1.5000000E-03,&
      2*1.3333330E-03,3*1.1515151E-03,2*9.8148151E-04,&
      9.4629632E-04,9.1111113E-04,8.6481485E-04,2*8.1851851E-04,&
      3*5.0000002E-04, 3.5000002E-04,5*1.9999999E-04/

! for xo2
     data ((cmaq_bc9(14,i,k),k=1,34),i=1,4)&
      /136*0./

! for xo2n
     data ((cmaq_bc9(15,i,k),k=1,34),i=1,4)&
      /136*0./

! for ntr
     data ((cmaq_bc9(16,i,k),k=1,34),i=1,4)&
      /136*0./ 

! for rooh
     data ((cmaq_bc9(17,i,k),k=1,34),i=1,4)&
      /136*0./  

! for form
     data ((cmaq_bc9(18,i,k),k=1,34),i=1,4)&
      /136*0./ 

! for ald2
     data ((cmaq_bc9(19,i,k),k=1,34),i=1,4)&
      /136*0./  

! for aldx
     data ((cmaq_bc9(20,i,k),k=1,34),i=1,4)&
      /136*0./  

! for par
     data ((cmaq_bc9(21,i,k),k=1,34),i=1,4)&
      /136*0./  

! for co
     data (cmaq_bc9(22,1,k),k=1,34)&
      /18*7.0000000E-02,2*6.9537044E-02,6.8657413E-02,&
       6.7777783E-02,6.6620372E-02,2*6.5462962E-02,&
       3*5.9999999E-02,5.7499997E-02,5*5.5000000E-02/
      data (cmaq_bc9(22,2,k),k=1,34)&  
      /13*7.9999998E-02,2*7.8333333E-02,3*7.6515153E-02,&
       2*7.4537039E-02, 7.3657408E-02, 7.2777778E-02,7.1620367E-02,&
       7.0462957E-02,   7.0462957E-02, 6.7500003E-02, 6.7500003E-02,&
       6.7500003E-02,   6.6249996E-02, 5*6.4999998E-02/
      data (cmaq_bc9(22,3,k),k=1,34)&   
      /13*7.9999998E-02,2*7.6666668E-02,3*7.3030300E-02,&
       2*6.9537044E-02,6.8657413E-02,6.7777783E-02,6.6620372E-02,&
       6.5462962E-02, 6.5462962E-02, 5.7500001E-02, 5.7500001E-02,&
       5.7500001E-02, 5.3750001E-02, 5*5.0000001E-02/
      data (cmaq_bc9(22,4,k),k=1,34)&
      /13*7.9999998E-02,2*7.6666668E-02,3*7.3030300E-02,&
       2*6.9537044E-02,6.8657413E-02,6.7777783E-02,6.6620372E-02,&
       6.5462962E-02, 6.5462962E-02, 5.7500001E-02, 5.7500001E-02,&
       5.7500001E-02, 5.3750001E-02, 5*5.0000001E-02/

! for meo2
     data ((cmaq_bc9(23,i,k),k=1,34),i=1,4)&
      /136*0./

! for mepx
      data (cmaq_bc9(24,1,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07, &
       2*1.9537040E-07,1.8657411E-07,1.7777781E-07,1.6620371E-07,&
       2*1.5462960E-07,3*1.2500000E-07,1.1250000E-07,&
       5*1.0000000E-07/  
      data (cmaq_bc9(24,2,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07, &
       2*1.9537040E-07,1.8657411E-07,1.7777781E-07,1.6620371E-07,&
       2*1.5462960E-07,3*1.2500000E-07,1.1250000E-07,&
       5*1.0000000E-07/
      data (cmaq_bc9(24,3,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9074071E-07,1.7314815E-07,1.5555560E-07,1.3240745E-07,&
       1.0925930E-07,1.0925930E-07,3*7.5000003E-08,6.2500000E-08,&
       5*5.0000001E-08/
     data (cmaq_bc9(24,4,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9074071E-07,1.7314815E-07,1.5555560E-07,1.3240745E-07,&
       1.0925930E-07,1.0925930E-07,3*7.5000003E-08,6.2500000E-08,&
       5*5.0000001E-08/

! for meoh
     data ((cmaq_bc9(25,i,k),k=1,34),i=1,4)&
      /136*0./

! for hco3
     data ((cmaq_bc9(26,i,k),k=1,34),i=1,4)&    
      /136*0./ 

! for facd
     data (cmaq_bc9(27,1,k),k=1,34)&
     /2*1.5000001E-06,1.3571430E-06,1.1428572E-06,&
      9.2857141E-07,7.1428570E-07,28*5.0000000E-07/
     data (cmaq_bc9(27,2,k),k=1,34)&
     /2*1.5000001E-06,1.3571430E-06,1.1428572E-06,&
      9.2857141E-07,7.1428570E-07,28*5.0000000E-07/
     data (cmaq_bc9(27,3,k),k=1,34)&
     /5*1.0000000E-06,9.8214264E-07,9.6428539E-07,&
        8.7499978E-07,7.8571418E-07,7.8571418E-07,&
        15*5.0000000E-07,3*2.5000000E-07,1.2500000E-07,&
        5*0./
     data (cmaq_bc9(27,4,k),k=1,34)&
     /5*1.0000000E-06,9.8214264E-07,9.6428539E-07,&
        8.7499978E-07,7.8571418E-07,7.8571418E-07,&
        15*5.0000000E-07,3*2.5000000E-07,1.2500000E-07,&
        5*0./

! for c2o3
     data ((cmaq_bc9(28,i,k),k=1,34),i=1,4)&
     /136*9.9999997E-10/

! for pan
     data ((cmaq_bc9(29,i,k),k=1,34),i=1,2)&
     /68*1.5000000E-05/
     data (cmaq_bc9(29,3,k),k=1,34)&
     /5*9.9999997E-05,9.8214266E-05,9.6428543E-05,&
      8.7499982E-05,7.8571422E-05,7.8571422E-05,8*4.9999999E-05,&
      2*4.6759251E-05, 4.0601852E-05,3.4444449E-05,2.6342594E-05,&
      2*1.8240740E-05, 3*7.4999998E-06,3.7499999E-06,5*0./
     data (cmaq_bc9(29,4,k),k=1,34)&
     /5*9.9999997E-05,9.8214266E-05,9.6428543E-05,&
      8.7499982E-05,7.8571422E-05,7.8571422E-05,8*4.9999999E-05,&
      2*4.6759251E-05, 4.0601852E-05,3.4444449E-05,2.6342594E-05,&
      2*1.8240740E-05, 3*7.4999998E-06,3.7499999E-06,5*0./

! for pacd
     data (cmaq_bc9(30,1,k),k=1,34)&
      /18*9.9999997E-05,2*9.5370357E-05,8.6574073E-05,&
       7.7777782E-05,6.6203706E-05,5.4629629E-05,5.4629629E-05,&
       9*4.9999999E-05/
     data (cmaq_bc9(30,2,k),k=1,34)&
      /34*4.9999999E-05/
     data (cmaq_bc9(30,3,k),k=1,34)&
      /13*2.9999999E-05,2*2.8333330E-05,3*2.6515150E-05,&
        2*2.4537039E-05,2.3657409E-05,2.2777780E-05,2.1620370E-05,&
        2*2.0462960E-05, 3*1.7500000E-05,1.6250000E-05,5*1.5000000E-05/
     data (cmaq_bc9(30,4,k),k=1,34)&
      /13*2.9999999E-05,2*2.8333330E-05,3*2.6515150E-05,&
        2*2.4537039E-05,2.3657409E-05,2.2777780E-05,2.1620370E-05,&
        2*2.0462960E-05, 3*1.7500000E-05,1.6250000E-05,5*1.5000000E-05/

! for aacd
      data (cmaq_bc9(31,1,k),k=1,34)&
        /2*1.5000001E-06, 1.3571430E-06,1.1428572E-06,9.2857141E-07,&
        7.1428570E-07,  28*5.0000000E-07/
      data (cmaq_bc9(31,2,k),k=1,34)&
        /2*1.5000001E-06, 1.3571430E-06,1.1428572E-06,9.2857141E-07,&
        7.1428570E-07,  28*5.0000000E-07/
      data (cmaq_bc9(31,3,k),k=1,34)&
       /5*1.0000000E-06,9.8214264E-07,9.6428539E-07,&
        8.7499978E-07,7.8571418E-07,7.8571418E-07,15*5.0000000E-07,&
        3*2.5000000E-07,1.2500000E-07,5*0/
      data (cmaq_bc9(31,4,k),k=1,34)&
       /5*1.0000000E-06,9.8214264E-07,9.6428539E-07,&
        8.7499978E-07,7.8571418E-07,7.8571418E-07,15*5.0000000E-07,&
        3*2.5000000E-07,1.2500000E-07,5*0/ 

! for cxo3
     data ((cmaq_bc9(32,i,k),k=1,34),i=1,4)&
      /136*9.9999997E-10/

! for panx
     data ((cmaq_bc9(33,i,k),k=1,34),i=1,4)&
      /136*9.9999997E-10/

! for ror
     data ((cmaq_bc9(34,i,k),k=1,34),i=1,4)&
      /136*0./

! for ole
     data ((cmaq_bc9(35,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(35,3,k),k=1,34)&
      /5*2.0000000E-07,1.9642854E-07,1.9285710E-07,&
      1.7499994E-07,2*1.5714279E-07,3*1.0000000E-07,2*6.6666693E-08,&
      3*3.0303010E-08,16*0./
     data (cmaq_bc9(35,4,k),k=1,34)&
      /5*2.0000000E-07,1.9642854E-07,1.9285710E-07,&
      1.7499994E-07,2*1.5714279E-07,3*1.0000000E-07,2*6.6666693E-08,&
      3*3.0303010E-08,16*0./

! for eth
     data ((cmaq_bc9(36,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(36,3,k),k=1,34)&
      /2*4.9999999E-06, 4.7142862E-06,4.2857146E-06,&
      3.8571429E-06,3.3928568E-06,2.9285709E-06,2.7499996E-06,&
      2*2.5714280E-06,3*2.0000000E-06,2*1.6666670E-06,3*1.3030300E-06,&
      7*1.0000000E-06, 3* 5.0000000E-07,2.5000000E-07,5*0./
     data (cmaq_bc9(36,4,k),k=1,34)&
      /2*4.9999999E-06, 4.7142862E-06,4.2857146E-06,&
      3.8571429E-06,3.3928568E-06,2.9285709E-06,2.7499996E-06,&
      2*2.5714280E-06,3*2.0000000E-06,2*1.6666670E-06,3*1.3030300E-06,&
      7*1.0000000E-06, 3* 5.0000000E-07,2.5000000E-07,5*0./

! for iole
     data ((cmaq_bc9(37,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(37,3,k),k=1,34)&
      /5*1.0000000E-07,9.6428536E-08,9.2857071E-08,&
       7.4999960E-08,2*5.7142849E-08,24*0./
     data (cmaq_bc9(37,4,k),k=1,34)&
      /5*1.0000000E-07,9.6428536E-08,9.2857071E-08,&
       7.4999960E-08,2*5.7142849E-08,24*0./

! for tol
     data ((cmaq_bc9(38,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(38,3,k),k=1,34)&
      /18*1.0000000E-06, 2*9.0740730E-07,&
       7.3148146E-07,5.5555557E-07,3.2407408E-07,&
       2*9.2592600E-08,9*0./
     data (cmaq_bc9(38,4,k),k=1,34)&
      /18*1.0000000E-06, 2*9.0740730E-07,&
       7.3148146E-07,5.5555557E-07,3.2407408E-07,&
       2*9.2592600E-08,9*0./

! for cres
     data ((cmaq_bc9(39,i,k),k=1,34),i=1,4)&
      /136*9.9999997E-10/

! for to2
     data ((cmaq_bc9(40,i,k),k=1,34),i=1,4)&
      /136*9.9999997E-10/

! for tolaer1
     data ((cmaq_bc9(41,i,k),k=1,34),i=1,4)&
      /136*0./
! for tolaer2
     data ((cmaq_bc9(42,i,k),k=1,34),i=1,4)&
      /136*0./

! for open
     data (cmaq_bc9(43,1,k),k=1,34)&
      /13*2.5000000E-07, 2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9537040E-07, 1.8657411E-07,1.7777781E-07,1.6620371E-07,&
       2*1.5462960E-07,3*1.2500000E-07,1.1250000E-07,5*1.0000000E-07/
     data (cmaq_bc9(43,2,k),k=1,34)&
      /13*2.5000000E-07, 2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9537040E-07, 1.8657411E-07,1.7777781E-07,1.6620371E-07,&
       2*1.5462960E-07,3*1.2500000E-07,1.1250000E-07,5*1.0000000E-07/
     data (cmaq_bc9(43,3,k),k=1,34)&
      / 13*2.5000000E-07, 2*2.3333330E-07,&
        3*2.1515150E-07, 2*1.9074071E-07,1.7314815E-07,1.5555560E-07,&
        1.3240745E-07,2*1.0925930E-07,&
        3*7.5000003E-08, 6.2500000E-08, 5*5.0000001E-08/
     data (cmaq_bc9(43,4,k),k=1,34)&
      / 13*2.5000000E-07, 2*2.3333330E-07,&
        3*2.1515150E-07, 2*1.9074071E-07,1.7314815E-07,1.5555560E-07,&
        1.3240745E-07,2*1.0925930E-07,&
        3*7.5000003E-08, 6.2500000E-08, 5*5.0000001E-08/

! for cro
      data ((cmaq_bc9(44,i,k),k=1,34),i=1,4)&
      /136*0./

! for cslaer
     data ((cmaq_bc9(45,i,k),k=1,34),i=1,4)&
      /136*0./

! for mgly
     data (cmaq_bc9(46,1,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9537040E-07, 1.8657411E-07, 1.7777781E-07, 1.6620371E-07,&
       2*1.5462960E-07, 3*1.2500000E-07,1.1250000E-07,5*1.0000000E-07/
     data (cmaq_bc9(46,2,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9537040E-07, 1.8657411E-07, 1.7777781E-07, 1.6620371E-07,&
       2*1.5462960E-07, 3*1.2500000E-07,1.1250000E-07,5*1.0000000E-07/
     data (cmaq_bc9(46,3,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9074071E-07,1.7314815E-07,1.5555560E-07,1.3240745E-07,&
       2*1.0925930E-07, 3*7.5000003E-08,6.2500000E-08,&
       5*5.0000001E-08 /
     data (cmaq_bc9(46,4,k),k=1,34)&
      /13*2.5000000E-07,2*2.3333330E-07,3*2.1515150E-07,&
       2*1.9074071E-07,1.7314815E-07,1.5555560E-07,1.3240745E-07,&
       2*1.0925930E-07, 3*7.5000003E-08,6.2500000E-08,&
       5*5.0000001E-08 /

! for xyl
     data ((cmaq_bc9(47,i,k),k=1,34),i=1,2)&
      /68*0./ 
      data (cmaq_bc9(47,3,k),k=1,34)&
      / 5*2.0000000E-07,1.9642854E-07,1.9285710E-07,&
        1.7499994E-07,2*1.5714279E-07,3*1.0000000E-07,&
         2*6.6666693E-08,3*3.0303010E-08, 16*0.0/
      data (cmaq_bc9(47,4,k),k=1,34)&
      / 5*2.0000000E-07,1.9642854E-07,1.9285710E-07,&
        1.7499994E-07,2*1.5714279E-07,3*1.0000000E-07,&
         2*6.6666693E-08,3*3.0303010E-08, 16*0.0/

! for xylaer1
      data ((cmaq_bc9(48,i,k),k=1,34),i=1,4)&
      /136*0./ 

! for xylaer2
      data ((cmaq_bc9(49,i,k),k=1,34),i=1,4)&
      /136*0./

! for isop
     data ((cmaq_bc9(50,i,k),k=1,34),i=1,2)&
      /68*0./  
     data (cmaq_bc9(50,3,k),k=1,34)& 
      /5*1.5000001E-04,1.4464281E-04,1.3928561E-04,&
       1.1249994E-04,2*8.5714273E-05,24*0./
     data (cmaq_bc9(50,4,k),k=1,34)&
      /5*1.5000001E-04,1.4464281E-04,1.3928561E-04,&
       1.1249994E-04,2*8.5714273E-05,24*0./

! for ispd
      data ((cmaq_bc9(51,i,k),k=1,34),i=1,4)&
      /136*0./  

! for isoaer1
      data ((cmaq_bc9(52,i,k),k=1,34),i=1,4)&
      /136*0./

! for isoaer2
      data ((cmaq_bc9(53,i,k),k=1,34),i=1,4)&
      /136*0./

! so2
     data ((cmaq_bc9(54,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(54,3,k),k=1,34)&
      /2*3.0000001E-04,2.8571431E-04, 2.6428571E-04,&
       2.4285710E-04,2.1785710E-04,1.9285710E-04,&
       1.7499995E-04,2*1.5714280E-04,1.5714280E-04,8*9.9999997E-05,&
       2*9.2592578E-05,7.8518511E-05,6.4444437E-05,4.5925924E-05,&
       2*2.7407410E-05,3*1.5000000E-05,5*9.9999997E-06/
     data (cmaq_bc9(54,4,k),k=1,34)&
      /2*3.0000001E-04,2.8571431E-04, 2.6428571E-04,&
       2.4285710E-04,2.1785710E-04,1.9285710E-04,&
       1.7499995E-04,2*1.5714280E-04,1.5714280E-04,8*9.9999997E-05,&
       2*9.2592578E-05,7.8518511E-05,6.4444437E-05,4.5925924E-05,&
       2*2.7407410E-05,3*1.5000000E-05,5*9.9999997E-06/

! sulf
     data ((cmaq_bc9(55,i,k),k=1,34),i=1,4)&
      /136*1.0000000E-30/

! sulfaer
      data ((cmaq_bc9(56,i,k),k=1,34),i=1,4)&
      /136*0./  

! etoh 
      data ((cmaq_bc9(57,i,k),k=1,34),i=1,4)&
      /136*0./

! etha
      data (cmaq_bc9(58,1,k),k=1,34)&
      /13*9.9999997E-05,2*8.9000008E-05,&
       3*7.6999990E-05,16*6.7000001E-05/
      data (cmaq_bc9(58,2,k),k=1,34)&
      /13*9.9999997E-05,2*8.9000008E-05,&
       3*7.6999990E-05,2*6.3851847E-05,&
       5.7870369E-05,5.1888888E-05,4.4018518E-05,&
       2*3.6148151E-05,9*3.3000000E-05/
     data (cmaq_bc9(58,3,k),k=1,34)&
      /5*9.9999997E-05,9.9285702E-05,9.8571407E-05,&
       9.4999989E-05,2*9.1428570E-05,3*7.9999998E-05,&
       2*6.3333340E-05, 3*4.5151501E-05,&
       2*2.8148141E-05, 2.4629626E-05, 2.1111109E-05,1.6481479E-05,&
       2*1.1851850E-05, 3*4.9999999E-06,2.4999999E-06,5*0./  
     data (cmaq_bc9(58,4,k),k=1,34)&
      /5*9.9999997E-05,9.9285702E-05,9.8571407E-05,&
       9.4999989E-05,2*9.1428570E-05,3*7.9999998E-05,&
       2*6.3333340E-05, 3*4.5151501E-05,&
       2*2.8148141E-05, 2.4629626E-05, 2.1111109E-05,1.6481479E-05,&
       2*1.1851850E-05, 3*4.9999999E-06,2.4999999E-06,5*0./

! for terp
     data ((cmaq_bc9(59,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(59,3,k),k=1,34)&
      / 5*3.7500002E-05,3.6160702E-05,3.4821402E-05,&
       2.8124985E-05, 2*2.1428570E-05,24*0./
     data (cmaq_bc9(59,4,k),k=1,34)&
      / 5*3.7500002E-05,3.6160702E-05,3.4821402E-05,&
       2.8124985E-05, 2*2.1428570E-05,24*0./

! terpaer
      data ((cmaq_bc9(60,i,k),k=1,34),i=1,4)&
      /136*0./

! for hum
     data ((cmaq_bc9(61,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(61,3,k),k=1,34)&
      / 5*1.2412500662E-05,1.1969192362E-05,1.1525884062E-05,&
       0.9309370035E-05, 2*0.709285667E-05,24*0./
     data (cmaq_bc9(61,4,k),k=1,34)&
      / 5*1.218262564974E-05,1.174752725874E-05,1.131242886774E-05,&
       0.913696387695E-05, 2*0.69614995359E-05,24*0./

! for humaer
      data ((cmaq_bc9(62,i,k),k=1,34),i=1,4)&
      /136*0./

! for lim
     data ((cmaq_bc9(63,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(63,3,k),k=1,34)&
      / 5*0.506250027E-05,0.488269477E-05,0.470088927E-05,&
       0.3796872975E-05, 2*0.289285695E-05,24*0./
     data (cmaq_bc9(63,4,k),k=1,34)&
      / 5*0.3016087660858E-05,0.2908369101158E-05,0.2800650541458E-05,&
       0.2262064418565E-05, 2*0.172347845653E-05,24*0./

! for limaer1
      data ((cmaq_bc9(64,i,k),k=1,34),i=1,4)&
      /136*0./

! for limaer2
      data ((cmaq_bc9(65,i,k),k=1,34),i=1,4)&
      /136*0./

! for oci
     data ((cmaq_bc9(66,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(66,3,k),k=1,34)&
      / 5*0.8137500434E-05,0.7846872334E-05,0.7556244234E-05,&
       0.6103121745E-05, 2*0.464999969E-05,24*0./
     data (cmaq_bc9(66,4,k),k=1,34)&
      / 5*0.6262500334E-05,0.6038837234E-05,0.5815174134E-05,&
       0.4696872495E-05, 2*0.357857119E-05,24*0./

! for ociaer1
      data ((cmaq_bc9(67,i,k),k=1,34),i=1,4)&
      /136*0./

! for ociaer2
      data ((cmaq_bc9(68,i,k),k=1,34),i=1,4)&
      /136*0./

! for apin
     data ((cmaq_bc9(69,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(69,3,k),k=1,34)&
      / 5*0.331875E-05,0.3200222127E-05,0.3081694077E-05,&
       0.24890611725E-05, 2*0.1896428445E-05,24*0./
     data (cmaq_bc9(69,4,k),k=1,34)&
      / 5*0.69192375E-05,0.6672119608126E-05,0.6425001347226E-05,&
       0.5189425357305E-05, 2*0.395384973641E-05,24*0./

! for apinaer1
      data ((cmaq_bc9(70,i,k),k=1,34),i=1,4)&
      /136*0./

! for apinaer2
      data ((cmaq_bc9(71,i,k),k=1,34),i=1,4)&
      /136*0./

! for apinaer3
      data ((cmaq_bc9(72,i,k),k=1,34),i=1,4)&
      /136*0./

! for apinaer4
      data ((cmaq_bc9(73,i,k),k=1,34),i=1,4)&
      /136*0./

! for bpin
     data ((cmaq_bc9(74,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(74,3,k),k=1,34)&
      / 5*0.2587500138E-05,0.2495088438E-05,0.2402676738E-05,&
       0.1940623965E-05, 2*0.147857133E-05,24*0./
     data (cmaq_bc9(74,4,k),k=1,34)&
      / 5*0.85875E-05,0.8280800758E-05,0.7974101058E-05,&
       0.6440621565E-05, 2*0.490714253E-05,24*0./

! for bpinaer1
      data ((cmaq_bc9(75,i,k),k=1,34),i=1,4)&
      /136*0./

! for bpinaer2
      data ((cmaq_bc9(76,i,k),k=1,34),i=1,4)&
      /136*0./

! for bpinaer3
      data ((cmaq_bc9(77,i,k),k=1,34),i=1,4)&
      /136*0./

! for bpinaer4
      data ((cmaq_bc9(78,i,k),k=1,34),i=1,4)&
      /136*0./

! for bpinaer5
      data ((cmaq_bc9(79,i,k),k=1,34),i=1,4)&
      /136*0./

! for ter

     data ((cmaq_bc9(80,i,k),k=1,34),i=1,2)&
      /68*0./
     data (cmaq_bc9(80,3,k),k=1,34)&
      / 5*0.5925000316E-05,0.5713390916E-05,0.5501781516E-05,&
       0.444374763E-05, 2*0.338571406E-05,24*0./
     data (cmaq_bc9(80,4,k),k=1,34)&
      / 5*0.0532237528386E-05,0.0513228843486E-05,0.0494220158586E-05,&
       0.0399177912105E-05, 2*0.030413569401E-05,24*0./

! for teraer1
      data ((cmaq_bc9(81,i,k),k=1,34),i=1,4)&
      /136*0./

! for teraer2
      data ((cmaq_bc9(82,i,k),k=1,34),i=1,4)&
      /136*0./

! for alkh
      data ((cmaq_bc9(83,i,k),k=1,34),i=1,4)&
      /136*0./

! for alkhaer1
      data ((cmaq_bc9(84,i,k),k=1,34),i=1,4)&
      /136*0./

! for pah
      data ((cmaq_bc9(85,i,k),k=1,34),i=1,4)&
      /136*0./

! for pahaer1
      data ((cmaq_bc9(86,i,k),k=1,34),i=1,4)&
      /136*0./

! for pahaer2
      data ((cmaq_bc9(87,i,k),k=1,34),i=1,4)&
      /136*0./

! for h2
      data ((cmaq_bc9(88,i,k),k=1,34),i=1,4)&
      /136*0.5/

! for ch4
      data ((cmaq_bc9(89,i,k),k=1,34),i=1,4)&
      /136*1.7/

! for cl
      data ((cmaq_bc9(90,i,k),k=1,34),i=1,4)&
      /136*0./

! for hcl
      data ((cmaq_bc9(91,i,k),k=1,34),i=1,4)&
      /136*0./

! for fmcl
      data ((cmaq_bc9(92,i,k),k=1,34),i=1,4)&
      /136*0./
! for hg0
      data ((cmaq_bc9(93,i,k),k=1,34),i=1,4)&
      /136*0./
! for hg2
      data ((cmaq_bc9(94,i,k),k=1,34),i=1,4)&
      /136*0./
! for hocl
      data ((cmaq_bc9(95,i,k),k=1,34),i=1,4)&
      /136*0./
! for clo
      data ((cmaq_bc9(96,i,k),k=1,34),i=1,4)&
      /136*0./
! for cl2
      data ((cmaq_bc9(97,i,k),k=1,34),i=1,4)&
      /136*0./
! for nh3
      data ((cmaq_bc9(98,i,k),k=1,34),i=1,4)&
      /136*0./

       if (id_bdy .eq. 1 ) then
          chem_bv = cmaq_bc9(nch-1,1,kk)
       elseif (id_bdy .eq. 2 ) then
          chem_bv = cmaq_bc9(nch-1,3,kk)
       elseif (id_bdy .eq. 3 ) then
          chem_bv = cmaq_bc9(nch-1,4,kk)
       elseif (id_bdy .eq. 4 ) then
          chem_bv = cmaq_bc9(nch-1,2,kk)
       endif
       return
       end subroutine bdy_chem_value_cb05