pkg/thsice/thsice_test_addfluid.h

579b39f951 Jean*0001 #ifdef ALLOW_ADDFLUID
                0002 C   this is a small piece of code to add to S/R thsice_main.F
                0003 C   to check AddFluid implementation :
                0004 C   uses addMass instead of EmPmR but solution should not change.
                0005 C   (tested in global_ocean.cs32x15.thsice set-up)
                0006
                0007       IF ( selectAddFluid.NE.0 .AND. hsMax.GT.hiMax ) THEN
                0008        IF ( myIter.EQ.nIter0 ) THEN
                0009         WRITE(0,*) 'THSICE_MAIN: mv EmPmR to addMass(k=1)'
                0010        ENDIF
                0011        DO bj=myByLo(myThid),myByHi(myThid)
                0012         DO bi=myBxLo(myThid),myBxHi(myThid)
                0013          DO j = 1-OLy, sNy+OLy-1
                0014           DO i = 1-OLx, sNx+OLx-1
                0015            addMass(i,j,1,bi,bj) = -EmPmR(i,j,bi,bj)
                0016      &                            *rA(i,j,bi,bj)
                0017      &                            *maskC(i,j,1,bi,bj)
                0018            EmPmR(i,j,bi,bj) = 0. _d 0
                0019           ENDDO
                0020          ENDDO
                0021         ENDDO
                0022        ENDDO
                0023       ENDIF
                0024 #endif /* ALLOW_ADDFLUID */