48 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
49 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
50 |
C- state variables : |
C- state variables : |
51 |
DO j=1-Oly,sNy+Oly |
DO j=1-OLy,sNy+OLy |
52 |
DO i=1-Olx,sNx+Olx |
DO i=1-OLx,sNx+OLx |
53 |
iceMask(i,j,bi,bj) = 0. _d 0 |
iceMask(i,j,bi,bj) = 0. _d 0 |
54 |
iceHeight(i,j,bi,bj)= 0. _d 0 |
iceHeight(i,j,bi,bj)= 0. _d 0 |
55 |
snowHeight(i,j,bi,bj)=0. _d 0 |
snowHeight(i,j,bi,bj)=0. _d 0 |
62 |
ENDDO |
ENDDO |
63 |
ENDDO |
ENDDO |
64 |
C- fluxes : |
C- fluxes : |
65 |
DO j=1-Oly,sNy+Oly |
DO j=1-OLy,sNy+OLy |
66 |
DO i=1-Olx,sNx+Olx |
DO i=1-OLx,sNx+OLx |
67 |
sHeating(i,j,bi,bj) = 0. _d 0 |
sHeating(i,j,bi,bj) = 0. _d 0 |
68 |
flxCndBt(i,j,bi,bj) = 0. _d 0 |
flxCndBt(i,j,bi,bj) = 0. _d 0 |
69 |
snowPrc(i,j,bi,bj) = 0. _d 0 |
snowPrc(i,j,bi,bj) = 0. _d 0 |
71 |
#ifdef COUPLE_MODEL |
#ifdef COUPLE_MODEL |
72 |
dFdT(i,j,bi,bj) = 0. _d 0 |
dFdT(i,j,bi,bj) = 0. _d 0 |
73 |
#endif |
#endif |
74 |
|
icFlxSW (i,j,bi,bj) = 0. _d 0 |
75 |
|
icFlxAtm(i,j,bi,bj) = 0. _d 0 |
76 |
|
icFrwAtm(i,j,bi,bj) = 0. _d 0 |
77 |
ENDDO |
ENDDO |
78 |
ENDDO |
ENDDO |
79 |
C- oceanic mixed layer state : |
C- oceanic mixed layer state : |
80 |
v2Loc = vMxL_default*vMxL_default |
v2Loc = vMxL_default*vMxL_default |
81 |
DO j=1-Oly,sNy+Oly |
DO j=1-OLy,sNy+OLy |
82 |
DO i=1-Olx,sNx+Olx |
DO i=1-OLx,sNx+OLx |
83 |
hOceMxL(i,j,bi,bj) = hMxL_default |
hOceMxL(i,j,bi,bj) = hMxL_default |
84 |
tOceMxL(i,j,bi,bj) = 0. _d 0 |
tOceMxL(i,j,bi,bj) = 0. _d 0 |
85 |
sOceMxL(i,j,bi,bj) = sMxL_default |
sOceMxL(i,j,bi,bj) = sMxL_default |