call COLUMN(m,n,np,pa,dt,dco2,sco3,scopre,scotem)

        call COLUMN(m,n,np,pa,dt,do3,sco3,scopre,scotem)

        call H2OEXPS(ib,m,n,np,dh2o,pa,dt,xkw,aw,bw,pm,mw,h2oexp)

        call CO2EXPS(m,n,np,dco2,pa,dt,co2exp)

        call N2OEXPS(ib,m,n,np,dn2o,pa,dt,n2oexp)

        call CH4EXPS(ib,m,n,np,dch4,pa,dt,ch4exp)

        call B10EXPS(m,n,np,dh2o,dcont,dco2,dn2o,pa,dt

      subroutine COLUMN (m,n,np,pa,dt,sabs0,sabs,spre,stem)

      _RL pa(m,n,np),dt(m,n,np),sabs0(m,n,np)

           spre(i,j,k+1)=spre(i,j,k)+pa(i,j,k)*sabs0(i,j,k)

      subroutine H2OEXPS(ib,m,n,np,dh2o,pa,dt,xkw,aw,bw,pm,mw,h2oexp)

      _RL dh2o(m,n,np),pa(m,n,np),dt(m,n,np)

           xh = dh2o(i,j,k)*(pa(i,j,k)/500.)**pm(ib)

      subroutine CO2EXPS(m,n,np,dco2,pa,dt,co2exp)

      _RL dco2(m,n,np),pa(m,n,np),dt(m,n,np)

           xc = dco2(i,j,k)*(pa(i,j,k)/300.0)**0.5

           xc = dco2(i,j,k)*(pa(i,j,k)/30.0)**0.85

      subroutine N2OEXPS(ib,m,n,np,dn2o,pa,dt,n2oexp)

      _RL dn2o(m,n,np),pa(m,n,np),dt(m,n,np)

           xc=dn2o(i,j,k)*(pa(i,j,k)/500.0)**0.48

      subroutine CH4EXPS(ib,m,n,np,dch4,pa,dt,ch4exp)

      _RL dch4(m,n,np),pa(m,n,np),dt(m,n,np)

           xc=dch4(i,j,k)*(pa(i,j,k)/500.0)**0.65

      subroutine B10EXPS(m,n,np,dh2o,dcont,dco2,dn2o,pa,dt

      _RL dco2(m,n,np),pa(m,n,np),dt(m,n,np)

           xx=dh2o(i,j,k)*(pa(i,j,k)/500.0)

           xx=dco2(i,j,k)*(pa(i,j,k)/300.0)**0.5

      _RL x1,x2,x3,we,pe,fw,fp,pa,pb,pc,ax,ba,bb,t1,ca,cb,t2

        pa = coef1(ip,iw-1,nn)*(1.-fp)+coef1(ip+1,iw-1,nn)*fp

        ax = (-pa*(1.-fw)+pc*(1.+fw)) *fw*0.5 + pb*(1.-fw*fw)

      _RL pa(m,n,np),dt(m,n,np)

         pa(i,j,k)=0.5*(pl(i,j,k)+pl(i,j,k+1))

         xx=pa(i,j,k)*0.001618*wa(i,j,k)*wa(i,j,k)*dp

        call COLUMN(m,n,np,pa,dt,dh2o,sh2o,swpre,swtem)