/[MITgcm]/MITgcm/model/src/calc_phi_hyd.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_phi_hyd.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.22 by adcroft, Thu Nov 7 21:51:15 2002 UTC revision 1.23 by heimbach, Fri Nov 15 03:01:21 2002 UTC
# Line 132  C       P(z=eta) = P(atmospheric_loading Line 132  C       P(z=eta) = P(atmospheric_loading
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,
# Line 217  CADJ GENERAL Line 217  CADJ GENERAL
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
# Line 225  CADJ GENERAL Line 227  CADJ GENERAL
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
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22