IF ( iceMask(i,j,bi,bj).GT.0. _d 0

     &             *(iceMask(i-1,j,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0

     &             *(iceMask(i,j-1,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0

        CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,iceMask,   1,

        CALL MNC_CW_RL_W('D',fn,0,0,'iceMask',   iceMask,   myThid)

        IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN 

     &                   iceMask(i,j,bi,bj), hIce, hSnow, Tsf

       CALL TIMEAVE_CUMULATE(ice_fract_Ave, iceMask, 

       CALL TIMEAVE_CUMUL_2V( ice_iceH_Ave,  iceheight, iceMask,

       CALL TIMEAVE_CUMUL_2V( ice_snowH_Ave, snowheight, iceMask,

       CALL TIMEAVE_CUMUL_2V( ice_Tsrf_Ave,  Tsrf,  iceMask,

       CALL TIMEAVE_CUMUL_2V( ice_Tice1_Ave, Tice1, iceMask,

       CALL TIMEAVE_CUMUL_2V( ice_Tice2_Ave, Tice2, iceMask,

       CALL TIMEAVE_CUMUL_2V( ice_snowPr_Ave,snowPrc, iceMask,

     I                iceheight, iceMask, maskH, rA, yC, locDr,

     I                snowheight, iceMask, maskH, rA, yC, locDr,

     I                1, Qice1, iceMask, iceheight, rA, locDr,

     I                1, Qice2, iceMask, iceheight, rA, locDr,

     I                Tsrf, iceMask, maskH, rA, yC, locDr,

     I                Tice1, iceMask, iceheight, rA, yC, locDr,

     I                Tice2, iceMask, iceheight, rA, yC, locDr,

         IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN 

          icFrac  = iceMask(i,j,bi,bj)

     I                   snowPrc,   iceMask,tmpFac,1,'SIsnwPrc',

        compact= iceMask(i,j,bi,bj)

        IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN 

     &  iceMask(i,j,bi,bj),flxAtm(i,j),evpAtm(i,j),-Lfresh*snowPr

          icFrac = iceMask(i,j,bi,bj)

          icFrac = iceMask(i,j,bi,bj)

        esurp  = frzmltMxL - Fbot*iceMask(i,j,bi,bj)

          iceMask(i,j,bi,bj)=compact

          iceMask(i,j,bi,bj)  = 0. _d 0

     &                        )*iceMask(i,j,bi,bj)

     &       iceMask, iceHeight, snowHeight,

      _RL iceMask(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

        CALL DIAGNOSTICS_FILL(iceMask,'SI_Fract',0,1,1,bi,bj,myThid)

     I                   iceheight, iceMask,tmpFac,1,'SI_Thick',

     I                   snowheight,iceMask,tmpFac,1,'SI_SnowH',

     I                   Tsrf,      iceMask,tmpFac,1,'SI_Tsrf ',

     I                   Tice1,     iceMask,tmpFac,1,'SI_Tice1',

     I                   Tice2,     iceMask,tmpFac,1,'SI_Tice2',

          tmpFld(i,j) = iceMask(i,j,bi,bj)*iceheight(i,j,bi,bj)

         CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)

        _EXCH_XY_R8(iceMask,myThid)

            iceMask(i,j,bi,bj)  = 0. _d 0

        CALL MDSREADFIELD(fn,prec,'RL',1,iceMask,   1,myThid)

        CALL MNC_CW_RL_R('D',fn,0,0,'iceMask',   iceMask,   myThid)

          CALL WRITE_FLD_XY_RL('ice_fract.',  suff,iceMask,

          CALL MNC_CW_RL_W('D',mncf,0,0,'iceMask',   iceMask,   myThid)

          IF ( iceMask(i,j,bi,bj) .GT. 0. _d 0 ) THEN

          siceAlb(i,j,bi,bj) = iceMask(i,j,bi,bj)*aimAlb(i,j)

           aimSIfrc(i,j) = iceMask(i,j,bi,bj)

           iceMask(i,j,bi,bj) = aimSIfrc(i,j)

         icFrac = iceMask(i,j,bi,bj)

         IF ( iceMask(i,j,bi,bj).GT.0. _d 0 ) THEN

      IF (useThSIce)  _EXCH_XY_R8( iceMask, myThid )

     &                   )*iceMask(i,j,bi,bj)*cplTimeFraction

     &       *(1. _d 0 - iceMask(i,j-1,bi,bj)*stressReduction)

     &       *(1. _d 0 - iceMask(i, j ,bi,bj)*stressReduction)

     &       *(1. _d 0 - iceMask(i-1,j,bi,bj)*stressReduction)

     &       *(1. _d 0 - iceMask( i ,j,bi,bj)*stressReduction)