149 |
& tFld, sFld, |
& tFld, sFld, |
150 |
& alphaRho, myThid) |
& alphaRho, myThid) |
151 |
|
|
152 |
|
#ifdef ALLOW_DIAGNOSTICS |
153 |
|
IF ( useDiagnostics ) |
154 |
|
& CALL DIAGNOSTICS_FILL(alphaRho,'RHOAnoma',k,1,2,bi,bj,myThid) |
155 |
|
#endif |
156 |
|
|
157 |
C Quasi-hydrostatic terms are added in as if they modify the buoyancy |
C Quasi-hydrostatic terms are added in as if they modify the buoyancy |
158 |
IF (quasiHydrostatic) THEN |
IF (quasiHydrostatic) THEN |
159 |
CALL QUASIHYDROSTATICTERMS(bi,bj,k,alphaRho,myThid) |
CALL QUASIHYDROSTATICTERMS(bi,bj,k,alphaRho,myThid) |
238 |
CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
239 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
240 |
|
|
241 |
|
#ifdef ALLOW_DIAGNOSTICS |
242 |
|
IF ( useDiagnostics ) |
243 |
|
& CALL DIAGNOSTICS_FILL(alphaRho,'RHOAnoma',k,1,2,bi,bj,myThid) |
244 |
|
#endif |
245 |
|
|
246 |
C-- Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst |
C-- Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst |
247 |
DO j=jMin,jMax |
DO j=jMin,jMax |
248 |
DO i=iMin,iMax |
DO i=iMin,iMax |