!
      module ox_module
!
! Calculate matrix coefficients in FS for partitioning of OX into O and O3. 
! This must be called before the composition module (comp_major.F).
!
      implicit none
      contains
!-----------------------------------------------------------------------
      subroutine comp_ox(tn,o2,ox,op,o1,o1d,o3,no,n4s,n2d,o2p,n2p,nplus,
     |  nop,ne,xiop2p,xiop2d,barm,xnmbar,lev0,lev1,lon0,lon1,lat)
!
! Input args:
      integer,intent(in) :: lev0,lev1,lon0,lon1,lat
      real,dimension(lev0:lev1,lon0-2:lon1+2),intent(in) ::
     |  tn,          ! neutral temperature
     |  o2,          ! molecular oxygen
     |  ox,          ! ox from previous step
     |  op,          ! o+ ion
     |  o1,          ! atomic oxygen
     |  o1d,         ! o(1d)
     |  o3,          ! ozone
     |  no,          ! no
     |  n4s,         ! n(4s)
     |  n2d,         ! n(2d)
     |  o2p,         ! o2+ from elden
     |  n2p,         ! n2+ from elden
     |  nplus,       ! n+  from oplus
     |  nop,         ! no+ from elden
     |  ne,          ! ne  from elden
     |  xiop2p,      ! O+(2p) from oplus module
     |  xiop2d,      ! O+(2d) from oplus module
     |  barm,        ! mbar
     |  xnmbar       ! p0e(-z), etc.


      end subroutine comp_ox
      end module ox_module