106 |
ENDDO |
ENDDO |
107 |
ENDDO |
ENDDO |
108 |
ENDDO |
ENDDO |
109 |
|
#ifdef ALLOW_SHIFWFLX_CONTROL |
110 |
|
C maskSHI is a hack to play along with the general ctrl-package |
111 |
|
C infrastructure, where only the k=1 layer of a 3D mask is used |
112 |
|
C for 2D fields. We cannot use maskInC instead, because routines |
113 |
|
C like ctrl_get_gen and ctrl_set_unpack_xy require 3D masks. |
114 |
|
DO K=1,Nr |
115 |
|
DO J=1-Oly,sNy+Oly |
116 |
|
DO I=1-Olx,sNx+Olx |
117 |
|
maskSHI(I,J,K,bi,bj) = 0. _d 0 |
118 |
|
ENDDO |
119 |
|
ENDDO |
120 |
|
ENDDO |
121 |
|
DO K=1,Nr |
122 |
|
DO J=1-Oly,sNy+Oly |
123 |
|
DO I=1-Olx,sNx+Olx |
124 |
|
IF ( ABS(R_shelfice(I,J,bi,bj)) .GT. 0. _d 0 |
125 |
|
& .AND. hFacC(I,J,K,bi,bj) .NE. 0. _d 0 ) THEN |
126 |
|
maskSHI(I,J,K,bi,bj) = 1. _d 0 |
127 |
|
maskSHI(I,J,1,bi,bj) = 1. _d 0 |
128 |
|
ENDIF |
129 |
|
ENDDO |
130 |
|
ENDDO |
131 |
|
ENDDO |
132 |
|
#endif |
133 |
C - end bi,bj loops. |
C - end bi,bj loops. |
134 |
ENDDO |
ENDDO |
135 |
ENDDO |
ENDDO |