42 |
INTEGER i, j, k, bi, bj |
INTEGER i, j, k, bi, bj |
43 |
_RL PSTAR |
_RL PSTAR |
44 |
_RS mask_uice |
_RS mask_uice |
45 |
INTEGER myIter, myTile |
INTEGER myIter, myTile, kSurface |
|
|
|
46 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
47 |
INTEGER I_obc, J_obc, kSurface |
INTEGER I_obc, J_obc |
48 |
|
#endif /* ALLOW_OBCS */ |
49 |
|
|
50 |
IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN |
IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN |
51 |
kSurface = Nr |
kSurface = Nr |
52 |
ELSE |
ELSE |
53 |
kSurface = 1 |
kSurface = 1 |
54 |
ENDIF |
ENDIF |
|
#endif /* ALLOW_OBCS */ |
|
55 |
|
|
56 |
C-- Initialise all variables in common blocks: |
C-- Initialise all variables in common blocks: |
57 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
367 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
368 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
369 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
370 |
HSALT(i,j,bi,bj)=HEFF(i,j,bi,bj)*salt(i,j,1,bi,bj)* |
HSALT(i,j,bi,bj)=HEFF(i,j,bi,bj)*salt(i,j,kSurface,bi,bj)* |
371 |
& ICE2WATR*rhoConstFresh*SEAICE_salinity |
& ICE2WATR*rhoConstFresh*SEAICE_salinity |
372 |
ENDDO |
ENDDO |
373 |
ENDDO |
ENDDO |