h2otbl=high.and.(ib.eq.1.or.ib.eq.2.or.ib.eq.8)

       conbnd=ib.ge.3.and.ib.le.6

       co2bnd=ib.eq.3

       oznbnd=ib.eq.5

       n2obnd=ib.eq.6.or.ib.eq.7

       ch4bnd=ib.eq.6.or.ib.eq.7

       combnd=ib.eq.4.or.ib.eq.5

       f11bnd=ib.eq.4.or.ib.eq.5

       f12bnd=ib.eq.4.or.ib.eq.6

       f22bnd=ib.eq.4.or.ib.eq.6

       b10bnd=ib.eq.10

     *                 *(ta(i,j,k)*cb(5,ib)+cb(4,ib))+cb(3,ib))

     *                 +cb(2,ib))+cb(1,ib)

     *                    *(ts(i,j)*cb(5,ib)+cb(4,ib))+cb(3,ib))

     *                    +cb(2,ib))+cb(1,ib) )*emiss(i,j,ib)

         dbs(i,j)=(ts(i,j)*(ts(i,j)*(ts(i,j)*4.*cb(5,ib)

     *           +3.*cb(4,ib))+2.*cb(3,ib))+cb(2,ib))*emiss(i,j,ib)

          taucl(i,j,k,1)=cwp(i,j,k,1)*(aib(1,ib)+aib(2,ib)/

     *      reff(i,j,k,1)**aib(3,ib))

          taucl(i,j,k,2)=cwp(i,j,k,2)*(awb(1,ib)+(awb(2,ib)+(awb(3,ib)

     *      +awb(4,ib)*reff(i,j,k,2))*reff(i,j,k,2))*reff(i,j,k,2))

           w1=taucl(i,j,k,1)*(aiw(1,ib)+(aiw(2,ib)+(aiw(3,ib)

     *       +aiw(4,ib)*reff1)*reff1)*reff1)

           w2=taucl(i,j,k,2)*(aww(1,ib)+(aww(2,ib)+(aww(3,ib)

     *       +aww(4,ib)*reff2)*reff2)*reff2)

           g1=w1*(aig(1,ib)+(aig(2,ib)+(aig(3,ib)

     *      +aig(4,ib)*reff1)*reff1)*reff1)

           g2=w2*(awg(1,ib)+(awg(2,ib)+(awg(3,ib)

     *      +awg(4,ib)*reff2)*reff2)*reff2)

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

        call CONEXPS(ib,m,n,np,dcont,xke,conexp)

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

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

        call COMEXPS(ib,m,n,np,dco2,dt,comexp)

        call CFCEXPS(ib,m,n,np,a1,b1,fk1,a2,b2,fk2,df11,dt,f11exp)

        call CFCEXPS(ib,m,n,np,a1,b1,fk1,a2,b2,fk2,df12,dt,f12exp)

        call CFCEXPS(ib,m,n,np,a1,b1,fk1,a2,b2,fk2,df22,dt,f22exp)

          if (ib.eq.1) then

          if (ib.eq.2) then

          if (ib.eq.8) then

        call H2OKDIS(ib,m,n,np,k2-1,fkw,gkw,ne,h2oexp,conexp,

          call N2OKDIS(ib,m,n,np,k2-1,n2oexp,tn2o,trant)

          call CH4KDIS(ib,m,n,np,k2-1,ch4exp,tch4,trant)

          call COMKDIS(ib,m,n,np,k2-1,comexp,tcom,trant)

         ff=0.5+(0.3739+(0.0076+0.1185*asyal(i,j,k2-1,ib))

     *      *asyal(i,j,k2-1,ib))*asyal(i,j,k2-1,ib)

         taux=taual(i,j,k2-1,ib)*(1.-ssaal(i,j,k2-1,ib)*ff)

     *           flcd(i,j,np+1)*trantcr(i,j,k)*(1.-emiss(i,j,ib))

     *           flxd(i,j,np+1)*transfc(i,j,k)*(1.-emiss(i,j,ib))

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

      integer ib,m,n,np,i,j,k,ik

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

     1        * ( 1.+(aw(ib)+bw(ib)* dt(i,j,k))*dt(i,j,k) )

           h2oexp(i,j,k,1) = exp(-xh*xkw(ib))

         if (mw(ib).eq.6) then

        elseif (mw(ib).eq.8) then

        elseif (mw(ib).eq.9) then

      subroutine CONEXPS(ib,m,n,np,dcont,xke,conexp)

      integer ib,m,n,np,i,j,k,iq

           conexp(i,j,k,1) = exp(-dcont(i,j,k)*xke(ib))

       if (ib .eq. 3) then

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

      integer ib,m,n,np,i,j,k

          if (ib.eq.6) then

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

      integer ib,m,n,np,i,j,k

          if (ib.eq.6) then

      subroutine COMEXPS(ib,m,n,np,dcom,dt,comexp)

      integer ib,m,n,np,i,j,k

          if (ib.eq.4) then

      subroutine CFCEXPS(ib,m,n,np,a1,b1,fk1,a2,b2,fk2,dcfc,dt,cfcexp)

      integer ib,m,n,np,i,j,k

          if (ib.eq.4) then

      subroutine H2OKDIS(ib,m,n,np,k,fkw,gkw,ne,h2oexp,conexp,

      integer ib,m,n,np,k,i,j

      if (ne(ib).eq.0) then

           trnth2o      =(fkw(1,ib)*th2o(i,j,1)

     *                  + fkw(2,ib)*th2o(i,j,2)

     *                  + fkw(3,ib)*th2o(i,j,3)

     *                  + fkw(4,ib)*th2o(i,j,4)

     *                  + fkw(5,ib)*th2o(i,j,5)

     *                  + fkw(6,ib)*th2o(i,j,6))

      elseif (ne(ib).eq.1) then

           trnth2o      =(fkw(1,ib)*th2o(i,j,1)

     *                  + fkw(2,ib)*th2o(i,j,2)

     *                  + fkw(3,ib)*th2o(i,j,3)

     *                  + fkw(4,ib)*th2o(i,j,4)

     *                  + fkw(5,ib)*th2o(i,j,5)

     *                  + fkw(6,ib)*th2o(i,j,6))*tcon(i,j,1)

      subroutine N2OKDIS(ib,m,n,np,k,n2oexp,tn2o,tran)

      integer ib,m,n,np,k,i,j

          if (ib.eq.6) then

      subroutine CH4KDIS(ib,m,n,np,k,ch4exp,tch4,tran)

      integer ib,m,n,np,k,i,j

          if (ib.eq.6) then

      subroutine COMKDIS(ib,m,n,np,k,comexp,tcom,tran)

      integer ib,m,n,np,k,i,j

           if (ib.eq.4) then

      integer i,j,k,ip,iw,it,ib,ik,iq,isb,k1,k2

      do 1000 ib=1,10

                  ib(i) = 0

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLy + 1

                      ib(ndim)  = OLy + 1

                    ib(ndim)  = 1

                    ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

     &                 dim, dn, ib, ie, myThid)

      ib(1)       = OLx + 1

      ib(2)       = OLy + 1

      ib(3)       = 1

      CALL MNC_CW_ADD_GNAME('fizhi_veg', ndim, dim, dn,ib,ie, myThid)

      ib(1)       = 1

      CALL MNC_CW_ADD_GNAME('fizhi_chip', ndim, dim, dn,ib,ie, myThid)

      ib(1)       = 1

      ib(2)       = 1

      CALL MNC_CW_ADD_GNAME('fizhi_chiplev',ndim,dim,dn,ib,ie,myThid)

      integer dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)

      integer i,L,nc,ib,nt

       tx1(i) = cmb2pa *  (prs(i,ib+1) - prs(i,ib))/(grav*dt)

       tx1(i) = rhfrac(i)*rknob * cmass(i,ib) / tx1(i)

       if(cmass(i,K).gt.0.) clf(i,ib) = clf(i,ib) + tx1(i)

       if( clf(i,ib).gt.1.) clf(i,ib) = 1.

      do L=ib,k

      do L=ib,k

      detrain(i,ib) = detrain(i,ib) + cmass(i,ib)

       rains(I,ib) = rains(I,ib) + PCU(I)*dt * rhfrac(i)

      integer ib,ik,i,j,k

      do 100 ib=1,nband

           w1=(1.-(aia(ib,1)+(aia(ib,2)+

     *         aia(ib,3)*reff1)*reff1))*taucld(i,j,k,1)

           w2=(1.-(awa(ib,1)+(awa(ib,2)+

     *         awa(ib,3)*reff2)*reff2))*taucld(i,j,k,2)

           g1=(aig(ib,1)+(aig(ib,2)+aig(ib,3)*reff1)*reff1)*w1

           g2=(awg(ib,1)+(awg(ib,2)+awg(ib,3)*reff2)*reff2)*w2

             ssatau=ssaal(ib)*taual(i,j,k)

             asysto=asyal(ib)*ssaal(ib)*taual(i,j,k)

          flx(i,j,k) = flx(i,j,k)+fall(i,j,k)*hk(ib,ik)

          flc(i,j,k) = flc(i,j,k)+fclr(i,j,k)*hk(ib,ik)

          fdirir(i,j) = fdirir(i,j)+fsdir(i,j)*hk(ib,ik)

          fdifir(i,j) = fdifir(i,j)+fsdif(i,j)*hk(ib,ik)

      integer i,j,k,ib

      do 100 ib=1,nband

          taurs=ry(ib)*dp(i,j,k)

          tauoz=xk(ib)*oh(i,j,k)

          ssatau=ssaal(ib)*taual(i,j,k)+taurs

          asysto=asyal(ib)*ssaal(ib)*taual(i,j,k)

          flx(i,j,k)=flx(i,j,k)+fall(i,j,k)*hk(ib)

          flc(i,j,k)=flc(i,j,k)+fclr(i,j,k)*hk(ib)

       if(ib.eq.nband) then

           fdirpar(i,j)=fsdir(i,j)*hk(ib)

           fdifpar(i,j)=fsdif(i,j)*hk(ib)

           fdiruv(i,j)=fdiruv(i,j)+fsdir(i,j)*hk(ib)

           fdifuv(i,j)=fdifuv(i,j)+fsdif(i,j)*hk(ib)

      integer ia,ib,L,K

      integer ia,ib,L,K

      integer ia,ib,L,K

      subroutine STRIPIT(a,b,irun,ia,ib,l,k)

      integer ia,ib,irun,l,k

      _RL a(ia,l), b(ib,l)

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      if(Lena.eq.ib) then

        do 400 i=1,ib-Lena

      subroutine STRIPITINT(a,b,irun,ia,ib,l,k)

      integer ia,ib,irun,l,k,a(ia,l),b(ib,l)

      offset = ib*(k-1)

      Lena = min(ib,irun-offset)

      if(Lena.eq.ib) then

        do 400 i=1,ib-Lena

      subroutine PASTIT(b,a,ib,ia,irun,L,k)

      integer ib,ia,L,k,irun,Lena,offset

      _RL a(ia,l), b(ib,l)

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      subroutine PSTBITINT(b,a,ib,ia,irun,l,k)

      integer ib,ia,L,k,irun,Lena,offset

      integer b(ib,l)

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      subroutine PSTBMPIT(b,a,ib,ia,irun,l,k)

      integer ib,ia,L,k,irun,Lena,offset

      _RL a(ia,l), b(ib,l)

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      subroutine STRIP2TILE(a,indx,b,irun,ia,ib,levs,npeice)

      integer ia,ib,irun,levs,npeice

      _RL a(ia,levs), b(ib,levs)

      offset = ib*(npeice-1)

      Lena    = min(ib,irun-offset)

      if(Lena.eq.ib) then

        do 400 i=1,ib-Lena

      subroutine PASTE2GRD_OLD(b,indx,chfr,ib,numpts,a,ia,levs,npeice)

      integer ia,ib,levs,numpts,npeice

      _RL a(ia,levs), b(ib,levs), chfr(ib)

      offset = ib*(npeice-1)

      Lena    = min(ib,numpts-offset)

      subroutine PASTE2GRD (b,indx,chfr,ib,numpts,a,ia,levs,npeice,

      integer ia,ib,levs,numpts,npeice

      _RL a(ia,levs), b(ib,levs), chfr(ib)

      offset = ib*(npeice-1)

      Lena    = min(ib,numpts-offset)

      integer ia,ib,L,K

      integer ia,ib,L,K

                  ib(i) = 0

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLy + 1

                      ib(ndim)  = OLy + 1

                    ib(ndim)  = 1

                    ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

     &                 dim, dn, ib, ie, myThid)

      integer dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)

                  ib(i) = 0

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLx + 1

                      ib(ndim)  = OLy + 1

                      ib(ndim)  = OLy + 1

                    ib(ndim)  = 1

                    ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

                  ib(ndim)  = 1

     &                 '  | ', (ib(i), i=1,5),

     &                 dim, dn, ib, ie, myThid)

      integer dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)

      INTEGER dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)

        ib(i) = 0

      ib(ndim)  = 1

      ib(ndim)  = 1

      ib(ndim)  = 1

      CALL MNC_CW_ADD_GNAME('AIM_xyt',  ndim, dim,dn,ib,ie, myThid)

      ib(ndim)  = 1

      ib(ndim)  = 1

      CALL MNC_CW_ADD_GNAME('AIM_xyzt', ndim, dim,dn,ib,ie, myThid)

      INTEGER dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)

      INTEGER ib,ibj1,ibj2  ! runoff band variables

      DO ib=1,numBands

        IF (ib.GT.1) ibj1= rband(ib-1)+1

        IF (ib.LT.numBands) ibj2= rband(ib)

        run_b(ib)=0. _d 0

          run_b(ib)=run_b(ib) +

      DO ib=1,numBands

        if (ib.GT.1) ibj1=  rband(ib-1)+1

        if (ib.LT.numBands) ibj2= rband(ib)

        totrun_b(ib)=0.d0

            totrun_b(ib)=totrun_b(ib)+runoffVal(i,j)*maskC(i,j,1,1,1)

          runIndex(j)= ib     ! for lookup of rband as fn. of latitude

            runoffVal(i,j)=runoffVal(i,j)*maskC(i,j,1,1,1)/totrun_b(ib)

      INTEGER ib, ibj1, ibj2  ! runoff band loop counters

          offset = ib*(npeice-1)

          Lena   = MIN(ib,numpts-offset)

     &                field,indx,chfr,ib,numpts,npeice,

           offset = ib*(npeice-1)

           Lena    = MIN(ib,numpts-offset)

      INTEGER ib,numpts,npeice

      _RL field(ib,nlevs), chfr(ib)

        ib(1)  = 1

     &       dim, dn, ib, ie, myThid)

          CALL MNC_CW_ADD_GNAME(dn(1), 1, dim, dn, ib, ie, myThid)

      INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)

         ib(1)        = OLx + 1

         ib(2)        = OLy + 1

         ib(3)  = 1

         ib(4)  = 1

     &                          dim, dn, ib, ie, myThid )

      INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)

          ib(1)      = 1

          ib(2)      = 1

     &         dim, dn, ib, ie, myThid)

      INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)

          ib(1)      = OLx + 1

          ib(2)      = OLy + 1

          ib(3)      = 1

          ib(4)      = 1

     &         dim, dn, ib, ie, myThid)

          ib(1)      = OLx + 1

          ib(2)      = OLy + 1

          ib(3)      = 1

          ib(4)      = 1

     &         dim, dn, ib, ie, myThid)

          ib(1)      = 1

          ib(2)      = 1

     &         dim, dn, ib, ie, myThid)

      INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)

        ib(1)  = 1

     &       dim, dn, ib, ie, myThid)

          CALL MNC_CW_ADD_GNAME(dn(1), 1, dim, dn, ib, ie, myThid)

        ib(1)  = 1

        ib(2)      = 1

        ib(3)      = 1

     &       dim, dn, ib, ie, myThid)

     &       dim, dn, ib, ie, myThid)

      INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)