SUBROUTINE SOLVER58(IJ,IK)

C   ROUTINE TO SOLVE only for transported species for levels 47-58, 
C      using chemical loss at level 46 in the CTLOSS(T$,L$,Z$) array which is in COMMON,
C       only for 46 levels, and is calculated in the regular solver and ncb

#include "com2d.h"

c                             co2sc is the scale factor for CO2 decrease above 100 km (Solomon et al., 1985)
c         REAL co2sc(12)
c
c         DATA co2sc/1, 1, 1, 1, .99, .97, .95, .93, .91, .89, .87, .85/


C  Ox - leave the same
	cn58(1,IJ,IK) = cn58(1,IJ,IK)

C  NOy
       noyp = 0.0
       cn58(2,ij,ik) = (cn58(2,ij,ik)+(noyp-difn(2,ij,ik))*dt)/
     c           (1.+ctloss(2,ij,46)*dt)

C  Cly
       clyp = 0.0
       cn58(3,IJ,IK) = (cn58(3,IJ,IK)+(clyp-difn(3,ij,ik))*dt)/
     c           (1.+ctloss(3,ij,46)*dt)

C  Bry
       bryp = 0.0
       cn58(15,ij,ik) = (cn58(15,ij,ik)+(bryp-difn(15,ij,ik))*dt)/
     c           (1.+ctloss(15,ij,46)*dt)


C  Note: The BIG KAHUNAS of stratospheric chemistry (ie, HNO3, N2O5, ClONO2, solid HNO3, H2O) 
C        are all set to VERY LARGE LOSS RATES of 1.E-4 sec-1  (= lifetime of ~.1 day) above 90 km

C  HNO3
       hno3p = 0.0
       cn58(14,ij,ik) = (cn58(14,ij,ik)+(hno3p-difn(14,ij,ik))*dt)/
     c           (1.+1.E-4*dt)

C  Solid HNO3 
       shno3p = 0.0
       cn58(29,ij,ik) = (cn58(29,ij,ik)+(shno3p-difn(29,ij,ik))*dt)/
     c           (1.+1.E-4*dt)

C  ClONO2
       cntp = 0.0
       cn58(26,ij,ik) = (cn58(26,ij,ik)+(cntp-difn(26,ij,ik))*dt)/
     c           (1.+1.E-4*dt)

C  N2O5
       n2o5p = 0.0
       cn58(27,ij,ik) = (cn58(27,ij,ik)+(n2o5p-difn(27,ij,ik))*dt)/
     c           (1.+1.E-4*dt)

C  Solid H2O 
       sh2op = 0.0
       cn58(30,ij,ik) = (cn58(30,ij,ik)+(sh2op-difn(30,ij,ik))*dt)/
     c           (1.+1.E-4*dt)


C  CO  -- Note: CO production in lower thermosphere only from CO2 photolysis, CO2 is now computed

        COP = (J(41,IJ,46) + J(12,IJ,46))*C58(31,IJ,IK)*TFD(IJ)
   	CN58(11,IJ,IK) = (CN58(11,IJ,IK)+(COP-DIFN(11,IJ,IK))*DT)
     C		/(1.+CTLOSS(11,ij,46)*DT)

C  CO2 -- Now computed, we don't really know what OH is, so we won't include production from CO+OH
C         but we compute CO for levels 47-58, and we'll just keep O as constant MRatio above
C         level 46 (seems reasonable since it slowly peaks around 95-100 km), and reaction K111
c          is constant, so we'll use CO+O+M as CO2 production

	CO2P = (K(111,IJ,46)*C58(11,IJ,IK)*M(IJ,IK)
     c             *C(1,IJ,46)/M(IJ,46)*M(ij,ik))*TFD(IJ)
   	CN58(31,IJ,IK) = (CN58(31,IJ,IK)+(CO2P-DIFN(31,IJ,IK))*DT)
     C		/(1.+CTLOSS(31,ij,46)*DT)

C  CH3OOH CN(26)
        CH3OOHP = 0.0
	CN58(12,IJ,IK) = (CN58(12,IJ,IK)+(CH3OOHP-DIFN(12,IJ,IK))
     C 		*DT)/(1.+CTLOSS(12,ij,46)*DT)

C  N2O
	pn2o=0.0
	cn58(4,IJ,IK) = (cn58(4,IJ,IK)+(pn2o-difn(4,IJ,IK))*dt)/
     c		(1.+CTLOSS(4,ij,46)*dt)

C  CFCL3
	f11p=0.
	cn58(5,IJ,IK) = (cn58(5,IJ,IK)+(f11p-difn(5,IJ,IK))*dt)/
     c		(1.+CTLOSS(5,ij,46)*dt)

C  CF2CL2
	f12p=0.
	cn58(6,IJ,IK) = (cn58(6,IJ,IK)+(f12p-difn(6,IJ,IK))*dt)/
     c		(1.+CTLOSS(6,IJ,46)*dt)

C  CCL4
	ccl4p=0.
	cn58(7,IJ,IK) = (cn58(7,IJ,IK)+(ccl4p-difn(7,IJ,IK))*dt)
     c 		/(1.+CTLOSS(7,IJ,46)*dt)

C CN(37)  IS CH3CL
	ch3clp=0.
	cn58(8,IJ,IK) = (cn58(8,IJ,IK)+(ch3clP-difn(8,IJ,IK))*dt)
     c		/(1.+ctloss(8,ij,46)*dt)

C CN(18)  IS CH4
	ch4prod=0.0
	cn58(9,IJ,IK) = (cn58(9,IJ,IK)+(ch4prod-difn(9,IJ,IK))*dt)
     c		/(1.+ctloss(9,ij,46)*dt)

C CN(17) IS H2, Production above 90 km due to H2O photolysis, use c58 array (previous days' H2O), j46

        h2prod = j(25,IJ,46)*c58(28,IJ,IK)*tfd(IJ)
       	cn58(10,IJ,IK) = (cn58(10,IJ,IK)+(h2prod-difn(10,IJ,IK))*dt)
     c		/(1.+ctloss(10,ij,46)*dt)


C CN(40)  IS CH3CCL3
	CMPROD = 0.0E0
	CN58(13,IJ,IK) = (CN58(13,IJ,IK)+(CMPROD-DIFN(13,IJ,IK))*DT)
     c          /(1.+ctloss(13,ij,46)*DT)

c ch3br
	ch3brp=0.
	cn58(16,IJ,IK) = (cn58(16,IJ,IK)+(ch3brp-difn(16,IJ,IK))*dt)
     C		/(1.+ctloss(16,ij,46)*dt)

c CHClF2
	chclf2p=0.
	cn58(17,IJ,IK) = (cn58(17,IJ,IK)+(chclf2p-difn(17,IJ,IK))*dt)
     C		/(1.+ctloss(17,ij,46)*dt)

c C2Cl3F3
	c2cl3f3p=0.
	cn58(18,IJ,IK) = (cn58(18,IJ,IK)+(c2cl3f3p-difn(18,IJ,IK))*dt)
     C		/(1.+ctloss(18,ij,46)*dt)

c C2Cl2F4
	c2cl2f4p=0.
	cn58(19,IJ,IK) = (cn58(19,IJ,IK)+(c2cl2f4p-difn(19,IJ,IK))*dt)
     C		/(1.+ctloss(19,ij,46)*dt)

c C2ClF5
	c2clf5p=0.
	cn58(20,IJ,IK) = (cn58(20,IJ,IK)+(c2clf5p-difn(20,IJ,IK))*dt)
     C		/(1.+ctloss(20,ij,46)*dt)

c CBrClF2
	cbrclf2p=0.
	cn58(21,IJ,IK) = (cn58(21,IJ,IK)+(cbrclf2p-difn(21,IJ,IK))*dt)
     C		/(1.+ctloss(21,ij,46)*dt)

c CBrF3
	cbrf3p=0.
	cn58(22,IJ,IK) = (cn58(22,IJ,IK)+(cbrf3p-difn(22,IJ,IK))*dt)
     C		/(1.+ctloss(22,ij,46)*dt)

c HF
	hfp = 0.0
	cn58(23,IJ,IK) = (cn58(23,IJ,IK)+(hfp-difn(23,IJ,IK))*dt)
     c		/(1.+ctloss(23,ij,46)*dt)

c CClFO
	CClFOp = 0.0
 	cn58(24,IJ,IK) = (cn58(24,IJ,IK)+(cclfop-difn(24,IJ,IK))*dt)
     c		/(1.+ctloss(24,ij,46)*dt)

c COF2
	COF2p = 0.0
	cn58(25,IJ,IK) = (cn58(25,IJ,IK)+(cof2p-difn(25,IJ,IK))*dt)
     c		/(1.+ctloss(25,ij,46)*dt)

C H2O
        h2op = 0.0
        cn58(28,ij,ik) = (cn58(28,ij,ik)+(h2op-difn(28,ij,ik))*dt)
     c         /(1.+ctloss(28,ij,46)*dt)

c
c  Set limit to 1.E-12 number density
C    (m at 115 km ~1.e12, and ~3.e19 at the ground, so the minimum mixing ratio will be ~3.E-31)

       DO 7382 III=1,ITRANS
       IF (CN58(III,IJ,IK) .LT. 1.E-12) CN58(III,IJ,IK) = 1.E-12
7382       CONTINUE


      SAVE
      RETURN
      END