143 |
C--- Calculate density |
C--- Calculate density |
144 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
145 |
kkey = (ikey-1)*Nr + k |
kkey = (ikey-1)*Nr + k |
146 |
CADJ STORE tFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE tFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte, |
147 |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & kind = isbyte |
148 |
|
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte, |
149 |
|
CADJ & kind = isbyte |
150 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
151 |
CALL FIND_RHO_2D( |
CALL FIND_RHO_2D( |
152 |
I iMin, iMax, jMin, jMax, k, |
I iMin, iMax, jMin, jMax, k, |
260 |
C-- Calculate density |
C-- Calculate density |
261 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
262 |
kkey = (ikey-1)*Nr + k |
kkey = (ikey-1)*Nr + k |
263 |
CADJ STORE tFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE tFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte, |
264 |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & kind = isbyte |
265 |
|
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte, |
266 |
|
CADJ & kind = isbyte |
267 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
268 |
CALL FIND_RHO_2D( |
CALL FIND_RHO_2D( |
269 |
I iMin, iMax, jMin, jMax, k, |
I iMin, iMax, jMin, jMax, k, |
272 |
O alphaRho, |
O alphaRho, |
273 |
I k, bi, bj, myThid ) |
I k, bi, bj, myThid ) |
274 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
275 |
CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte, |
276 |
|
CADJ & kind = isbyte |
277 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
278 |
ELSE |
ELSE |
279 |
DO j=jMin,jMax |
DO j=jMin,jMax |