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

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

      it=int(Iftm/Ifprd)+1

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

            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)+

       INTEGER I,J, kLev, it

        do it=1,12

         call TRACER_MEANAREA(myThid,po4obsl1,it,

     &                    po4avm(it,1))

         call TRACER_MEANAREA(myThid,po4obsl2,it,

     &                    po4avm(it,2))

         call TRACER_MEANAREA(myThid,po4obsl3,it,

     &                    po4avm(it,3))

         call TRACER_MEANAREA(myThid,o2obsl1,it,

     &                    o2avm(it,1))

         call TRACER_MEANAREA(myThid,o2obsl2,it,

     &                    o2avm(it,2))

         call TRACER_MEANAREA(myThid,o2obsl3,it,

     &                    o2avm(it,3))

          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(bi,bj,it)=0.d0

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

           do it=1,12

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

            po4avm(it,k)=0.d0

            o2avm(it,k)=0.d0

            po4varm(it,k)=0.d0

            o2varm(it,k)=0.d0

          do it=1,10

       INTEGER I,J, kLev, it

     &          /o2varm(it,1)

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

     &          /o2varm(it,2)

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

     &          /o2varm(it,3)

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

           do it=1,12

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

     &          /po4varm(it,1)

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

     &          /po4varm(it,2)

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

     &          /po4varm(it,3)

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

      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)

      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)

      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  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, it, ntot, nloop, ni

      DO it = 1,nloop

      INTEGER I, j, it, ntot, nloop, ni

      DO it = 1,nloop

      INTEGER I, j, it, ntot, nloop, ni

      DO it = 1,nloop

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

      INTEGER it