C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_test_addfluid.h,v 1.1 2008/08/24 21:50:12 jmc Exp $
C $Name:  $

#ifdef ALLOW_ADDFLUID
C   this is a small piece of code to add to S/R thsice_main.F
C   to check AddFluid implementation : 
C   uses addMass instead of EmPmR but solution should not change.
C   (tested in global_ocean.cs32x15.thsice set-up)

      IF ( selectAddFluid.NE.0 .AND. hsMax.GT.hiMax ) THEN
       IF ( myIter.EQ.nIter0 ) THEN
        WRITE(0,*) 'THSICE_MAIN: mv EmPmR to addMass(k=1)'
       ENDIF
       DO bj=myByLo(myThid),myByHi(myThid)
        DO bi=myBxLo(myThid),myBxHi(myThid)
         DO j = 1-OLy, sNy+OLy-1
          DO i = 1-OLx, sNx+OLx-1
           addMass(i,j,1,bi,bj) = -EmPmR(i,j,bi,bj)
     &                            *rA(i,j,bi,bj)
     &                            *maskC(i,j,1,bi,bj)
           EmPmR(i,j,bi,bj) = 0. _d 0
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDIF
#endif /* ALLOW_ADDFLUID */