152 |
ENDDO |
ENDDO |
153 |
|
|
154 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
155 |
|
CADJ STORE SlopeX(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
156 |
|
CADJ STORE SlopeY(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
157 |
CADJ STORE SlopeSqr(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE SlopeSqr(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
158 |
|
CADJ STORE dsigmadrreal(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
159 |
|
CADJ STORE taperfct(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
160 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
161 |
|
|
162 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
220 |
ENDDO |
ENDDO |
221 |
ENDDO |
ENDDO |
222 |
ENDIF |
ENDIF |
223 |
|
cph( NEW |
224 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
225 |
|
CADJ STORE VisbeckK(:,:,bi,bj) = comlev1_bibj, key=igmkey, byte=isbyte |
226 |
|
#endif |
227 |
|
cph) |
228 |
#endif /* GM_VISBECK_VARIABLE_K */ |
#endif /* GM_VISBECK_VARIABLE_K */ |
229 |
|
|
230 |
|
|
238 |
|
|
239 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
240 |
kkey = (igmkey-1)*Nr + k |
kkey = (igmkey-1)*Nr + k |
241 |
#ifdef GM_NON_UNITY_DIAGONAL |
#if (defined (GM_NON_UNITY_DIAGONAL) || \ |
242 |
|
defined (GM_VISBECK_VARIABLE_K)) |
243 |
CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
244 |
CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
245 |
CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
251 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
252 |
Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K |
Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K |
253 |
#ifdef GM_VISBECK_VARIABLE_K |
#ifdef GM_VISBECK_VARIABLE_K |
254 |
& + VisbeckK(i,j,bi,bj)*(1.+GM_skewflx) |
& + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx) |
255 |
#endif |
#endif |
256 |
Kwx(i,j,k,bi,bj)= Kgm_tmp*Kwx(i,j,k,bi,bj) |
Kwx(i,j,k,bi,bj)= Kgm_tmp*Kwx(i,j,k,bi,bj) |
257 |
Kwy(i,j,k,bi,bj)= Kgm_tmp*Kwy(i,j,k,bi,bj) |
Kwy(i,j,k,bi,bj)= Kgm_tmp*Kwy(i,j,k,bi,bj) |
294 |
O SlopeSqr, taperFct, |
O SlopeSqr, taperFct, |
295 |
I bi, bj, myThid ) |
I bi, bj, myThid ) |
296 |
|
|
297 |
|
cph( NEW |
298 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
299 |
|
cph( |
300 |
|
CADJ STORE taperfct(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
301 |
|
cph) |
302 |
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
303 |
|
cph) |
304 |
|
|
305 |
#ifdef GM_NON_UNITY_DIAGONAL |
#ifdef GM_NON_UNITY_DIAGONAL |
306 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
307 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
315 |
ENDDO |
ENDDO |
316 |
ENDDO |
ENDDO |
317 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
318 |
# ifndef GM_TAPER_ORIG_CLIPPING |
# ifdef GM_EXCLUDE_CLIPPING |
319 |
CADJ STORE Kux(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE Kux(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
320 |
# endif |
# endif |
321 |
#endif |
#endif |
375 |
O SlopeSqr, taperFct, |
O SlopeSqr, taperFct, |
376 |
I bi, bj, myThid ) |
I bi, bj, myThid ) |
377 |
|
|
378 |
|
cph( |
379 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
380 |
|
cph( |
381 |
|
CADJ STORE taperfct(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
382 |
|
cph) |
383 |
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
384 |
|
cph) |
385 |
|
|
386 |
#ifdef GM_NON_UNITY_DIAGONAL |
#ifdef GM_NON_UNITY_DIAGONAL |
387 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
388 |
DO i=1-Olx+1,sNx+Olx-1 |
DO i=1-Olx+1,sNx+Olx-1 |
396 |
ENDDO |
ENDDO |
397 |
ENDDO |
ENDDO |
398 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
399 |
# ifndef GM_TAPER_ORIG_CLIPPING |
# ifdef GM_EXCLUDE_CLIPPING |
400 |
CADJ STORE Kvy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE Kvy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
401 |
# endif |
# endif |
402 |
#endif |
#endif |