156 |
offsets(i) = 0 |
offsets(i) = 0 |
157 |
ENDDO |
ENDDO |
158 |
offsets(3) = k |
offsets(3) = k |
159 |
C write(*,*) 'offsets = ',(offsets(i),i=1,9) |
c write(*,*) 'offsets = ',(offsets(i),i=1,9) |
160 |
ENDIF |
ENDIF |
161 |
#endif /* ALLOW_MNC */ |
#endif /* ALLOW_MNC */ |
162 |
|
|
163 |
C Initialise intermediate terms |
C-- Initialise intermediate terms |
164 |
DO J=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
165 |
DO I=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
166 |
vF(i,j) = 0. |
vF(i,j) = 0. |
167 |
vrF(i,j) = 0. |
vrF(i,j) = 0. |
168 |
uCf(i,j) = 0. |
uCf(i,j) = 0. |
176 |
vort3(i,j) = 0. |
vort3(i,j) = 0. |
177 |
omega3(i,j)= 0. |
omega3(i,j)= 0. |
178 |
KE(i,j) = 0. |
KE(i,j) = 0. |
179 |
c hDiv(i,j) = 0. |
C- need to initialise hDiv for MOM_VI_DEL2UV(call FILL_CS_CORNER_TR_RL) |
180 |
|
hDiv(i,j) = 0. |
181 |
viscAh_Z(i,j) = 0. |
viscAh_Z(i,j) = 0. |
182 |
viscAh_D(i,j) = 0. |
viscAh_D(i,j) = 0. |
183 |
viscA4_Z(i,j) = 0. |
viscA4_Z(i,j) = 0. |
237 |
C and mask relative vorticity (free-slip case): |
C and mask relative vorticity (free-slip case): |
238 |
|
|
239 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
240 |
CADJ STORE vort3(:,:) = |
CADJ STORE vort3(:,:) = |
241 |
CADJ & comlev1_bibj_k, key = imomkey, byte = isbyte |
CADJ & comlev1_bibj_k, key = imomkey, byte = isbyte |
242 |
#endif |
#endif |
243 |
|
|