132 |
|
|
133 |
C Calculate density |
C Calculate density |
134 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
135 |
kkey = (ikey-1)*Nr + k |
kkey = (ikey-1)*Nr + k |
136 |
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 |
137 |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
138 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
139 |
CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k, |
CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k, |
217 |
DO j=jMin,jMax |
DO j=jMin,jMax |
218 |
DO i=iMin,iMax |
DO i=iMin,iMax |
219 |
phiHyd(i,j,k)=0. |
phiHyd(i,j,k)=0. |
220 |
|
#ifdef ATMOSPHERIC_LOADING |
221 |
phiHyd(i,j,k)=pload(i,j,bi,bj) |
phiHyd(i,j,k)=pload(i,j,bi,bj) |
222 |
|
#endif |
223 |
ENDDO |
ENDDO |
224 |
ENDDO |
ENDDO |
225 |
ENDIF |
ENDIF |
227 |
C Calculate density |
C Calculate density |
228 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
229 |
kkey = (ikey-1)*Nr + k |
kkey = (ikey-1)*Nr + k |
230 |
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 |
231 |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE sFld (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
232 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
233 |
CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k, |
CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k, |
234 |
& tFld, sFld, |
& tFld, sFld, |
235 |
& alphaRho, myThid) |
& alphaRho, myThid) |
236 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
237 |
|
CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
238 |
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
239 |
|
|
240 |
|
|
241 |
C Hydrostatic pressure at cell centers |
C Hydrostatic pressure at cell centers |
242 |
DO j=jMin,jMax |
DO j=jMin,jMax |