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

        do it=2,nt

           c1 (ip,iw,it)= c1(ip,iw,1)

           c2 (ip,iw,it)= c2(ip,iw,1)

           c3 (ip,iw,it)= c3(ip,iw,1)

           h11(ip,iw,it)=h11(ip,iw,1)

           h12(ip,iw,it)=h12(ip,iw,1)

           h13(ip,iw,it)=h13(ip,iw,1)

           h21(ip,iw,it)=h21(ip,iw,1)

           h22(ip,iw,it)=h22(ip,iw,1)

           h23(ip,iw,it)=h23(ip,iw,1)

           h81(ip,iw,it)=h81(ip,iw,1)

           h82(ip,iw,it)=h82(ip,iw,1)

           h83(ip,iw,it)=h83(ip,iw,1)

        do it=2,nt

           o1 (ip,iw,it)= o1(ip,iw,1)

           o2 (ip,iw,it)= o2(ip,iw,1)

           o3 (ip,iw,it)= o3(ip,iw,1)

      integer i,j,k,im,it,ia,kk

           it=int(ft+1.5)

           it=max(it,2)

           it=min(it,nt-1)

           ft=ft-float(it-1)

           xai=    (-caib(im-1,it,ia)*(1.-fm)+

     *      caib(im+1,it,ia)*(1.+fm))*fm*.5+caib(im,it,ia)*(1.-fm*fm)

           xai=xai+(-caib(im,it-1,ia)*(1.-ft)+

     *      caib(im,it+1,ia)*(1.+ft))*ft*.5+caib(im,it,ia)*(1.-ft*ft)

           xai=xai+(-caib(im,it,ia-1)*(1.-fa)+

     *     caib(im,it,ia+1)*(1.+fa))*fa*.5+caib(im,it,ia)*(1.-fa*fa)

           xai= xai-2.*caib(im,it,ia)

           xai=    (-caif(it-1,ia)*(1.-ft)+

     *      caif(it+1,ia)*(1.+ft))*ft*.5+caif(it,ia)*(1.-ft*ft)

           xai=xai+(-caif(it,ia-1)*(1.-fa)+

     *      caif(it,ia+1)*(1.+fa))*fa*.5+caif(it,ia)*(1.-fa*fa)

           xai= xai-caif(it,ia)

  as an error then it means you have not placed your configuration

      INTEGER it

      integer  nTimePnts, Iunit, it

      DO it = 1,nTimePnts

         READ(Iunit,*) year(it),ACFC11(it,1),ACFC12(it,1),

     &        ACFC11(it,2),ACFC12(it,2)

     &        ACFC11(it,1),ACFC12(it,1),

     &        ACFC11(it,2),ACFC12(it,2)

    as an error then it means you have not placed your coupler 

    as an error then it means you have not placed your coupler

    as an error then it means you have not placed your coupler 

      integer  nTimePnts, iUnit, it

      DO it = 1,nTimePnts

         READ(iUnit,*) year(it),ACFC11(it,1),ACFC12(it,1),

     &        ACFC11(it,2),ACFC12(it,2)

     &        ACFC11(it,1),ACFC12(it,1),

     &        ACFC11(it,2),ACFC12(it,2)

      INTEGER i, j, bi, bj, k, it

      it=int(Iftm/Ifprd)+1

         OBSM_timetave(it,bi,bj) = OBSM_timetave(it,bi,bj)+deltaTclock

            po4lev1(i,j,it,bi,bj)=po4lev1(i,j,it,bi,bj)+

            po4lev2(i,j,it,bi,bj)=po4lev2(i,j,it,bi,bj)+

            po4lev3(i,j,it,bi,bj)=po4lev3(i,j,it,bi,bj)+

            o2lev1(i,j,it,bi,bj)=o2lev1(i,j,it,bi,bj)+

            o2lev2(i,j,it,bi,bj)=o2lev2(i,j,it,bi,bj)+

            o2lev3(i,j,it,bi,bj)=o2lev3(i,j,it,bi,bj)+

        do it=1,MIN(12,Nr)

         call TRACER_MEANAREA(po4obsl1,it,

     &                    po4avm(it,1),myThid)

         call TRACER_MEANAREA(po4obsl2,it,

     &                    po4avm(it,2),myThid)

         call TRACER_MEANAREA(po4obsl3,it,

     &                    po4avm(it,3),myThid)

         call TRACER_MEANAREA(o2obsl1,it,

     &                    o2avm(it,1),myThid)

         call TRACER_MEANAREA(o2obsl2,it,

     &                    o2avm(it,2),myThid)

         call TRACER_MEANAREA(o2obsl3,it,

     &                    o2avm(it,3),myThid)

          DO it=1,12

           po4varm(it,1)=po4varm(it,1)+

     &                (po4obsl1(i,j,it,bi,bj)-po4avm(it,1))**2

           po4varm(it,2)=po4varm(it,2)+

     &                (po4obsl2(i,j,it,bi,bj)-po4avm(it,2))**2

           po4varm(it,3)=po4varm(it,3)+

     &                (po4obsl3(i,j,it,bi,bj)-po4avm(it,3))**2

           o2varm(it,1)=o2varm(it,1)+

     &                (o2obsl1(i,j,it,bi,bj)-o2avm(it,1))**2

           o2varm(it,2)=o2varm(it,2)+

     &                (o2obsl2(i,j,it,bi,bj)-o2avm(it,2))**2

           o2varm(it,3)=o2varm(it,3)+

     &                (o2obsl3(i,j,it,bi,bj)-o2avm(it,3))**2

         Do it=1,12

           po4varm(it,k)=po4varm(it,k)/volvar(k)

           o2varm(it,k)=o2varm(it,k)/volvar(k)

         do it=1,12

           OBSM_Timetave(it,bi,bj) = 0. _d 0

      INTEGER i, j, bi, bj, k, it

           do it=1,12

            OBSM_Timetave(it,bi,bj) = 0. _d 0

            po4avm(it,k) = 0. _d 0

            o2avm(it,k)  = 0. _d 0

            po4varm(it,k)= 0. _d 0

            o2varm(it,k) = 0. _d 0

      INTEGER i, j, bi, bj, k, it

             DO it=1,10

     I        i,j,kLev,bi,bj, it ,  myThid )

       INTEGER i,j, kLev, it

            it = exch2_neighbourId (np(i),ip(i))

            jt = exch2_myFace(it)

            IF ( exch2_jLo(nt,it).EQ.exch2_jHi(nt,it) )

     &        ii = 2 - MIN(1,exch2_jHi(nt,it))

            IF ( exch2_iLo(nt,it).EQ.exch2_iHi(nt,it) )

     &        ii = 4 - MIN(1,exch2_iHi(nt,it))

      it = 1

        IF ( exch2_nNeighbours(is).GT.exch2_nNeighbours(it) ) it = is

      IF ( it.NE.0 ) THEN

     &  'current Max.Nb.Neighbours (e.g., on tile',it,

     &  ' ) =', exch2_nNeighbours(it)

      INTEGER i, j, js, jp, jt, ii, is, it, ns, nt, k, tx, ty

            it = exch2_neighbourId (np(i),ip(i))

            jt = exch2_myFace(it)

            IF ( exch2_jLo(nt,it).EQ.exch2_jHi(nt,it) )

     &        ii = 2 - MIN(1,exch2_jHi(nt,it))

            IF ( exch2_iLo(nt,it).EQ.exch2_iHi(nt,it) )

     &        ii = 4 - MIN(1,exch2_iHi(nt,it))

            it = 2*ii - 3

               it = is + it*nbTx

               it = is + it

            IF ( exch2_myFace(it).NE.0 ) THEN

               exch2_neighbourId(ns,is) = it

             it = facet_owns(1,jt) + tx + ty*nbTx

             IF ( exch2_myFace(it).NE.0 ) THEN

               exch2_neighbourId(ns,is)  = it

               itbd1 = MIN( MAX( ibnd1, exch2_tBasex(it)+1 ),

     &                                  exch2_tBasex(it)+tNx )

               itbd2 = MIN( MAX( ibnd2, exch2_tBasex(it)+1 ),

     &                                  exch2_tBasex(it)+tNx )

               jtbd1 = MIN( MAX( jbnd1, exch2_tBasey(it)+1 ),

     &                                  exch2_tBasey(it)+tNy )

               jtbd2 = MIN( MAX( jbnd2, exch2_tBasey(it)+1 ),

     &                                  exch2_tBasey(it)+tNy )

      it = 0

          it = is

     &  'current Max.Nb.Neighbours (e.g., on tile',it,' ) =',nbNeighb

        it = exch2_neighbourId(ns,is)

        DO nt=1,exch2_nNeighbours(it)

          ii = MOD(tile_edge2edge(nt,it),10)

          IF ( exch2_neighbourId(nt,it).EQ.is .AND. ii.EQ.i ) THEN

            WRITE(msgBuf,'(A,I5,5(A,I3))') ' with tile', it, ' :',

     &      nn,' (',tile_edge2edge(nn,it),' ) and',

     &      nt,' (',tile_edge2edge(nt,it),' )'

     &      ns,' (',tile_edge2edge(ns,is),' ) no connection from',it

        it = exch2_neighbourId(ns,is)

        ii = exch2_neighbourId(nt,it)

        nn = exch2_opposingSend(nt,it)

          WRITE(msgBuf,'(A,I5,5(A,I3))') ' with tile', it, ' :',

     &      nt,' (',tile_edge2edge(nt,it),' )'

          WRITE(msgBuf,'(A,I5,2(A,I3),A)') ' but',it,' neighb:',

     &      nt,' (',tile_edge2edge(nt,it),' ) connected'

      INTEGER is, js, ns, it, jt, nt, tx, ty