1 |
jmc |
1.1 |
C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_main.F,v 1.17 2007/08/30 15:02:51 jmc Exp $ |
2 |
|
|
C $Name: $ |
3 |
|
|
|
4 |
|
|
#ifdef ALLOW_ADDFLUID |
5 |
|
|
C this is a small piece of code to add to S/R thsice_main.F |
6 |
|
|
C to check AddFluid implementation : |
7 |
|
|
C uses addMass instead of EmPmR but solution should not change. |
8 |
|
|
C (tested in global_ocean.cs32x15.thsice set-up) |
9 |
|
|
|
10 |
|
|
IF ( selectAddFluid.NE.0 .AND. hsMax.GT.hiMax ) THEN |
11 |
|
|
IF ( myIter.EQ.nIter0 ) THEN |
12 |
|
|
WRITE(0,*) 'THSICE_MAIN: mv EmPmR to addMass(k=1)' |
13 |
|
|
ENDIF |
14 |
|
|
DO bj=myByLo(myThid),myByHi(myThid) |
15 |
|
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
16 |
|
|
DO j = 1-OLy, sNy+OLy-1 |
17 |
|
|
DO i = 1-OLx, sNx+OLx-1 |
18 |
|
|
addMass(i,j,1,bi,bj) = -EmPmR(i,j,bi,bj) |
19 |
|
|
& *rA(i,j,bi,bj) |
20 |
|
|
& *maskC(i,j,1,bi,bj) |
21 |
|
|
EmPmR(i,j,bi,bj) = 0. _d 0 |
22 |
|
|
ENDDO |
23 |
|
|
ENDDO |
24 |
|
|
ENDDO |
25 |
|
|
ENDDO |
26 |
|
|
ENDIF |
27 |
|
|
#endif /* ALLOW_ADDFLUID */ |