/[MITgcm]/MITgcm/pkg/mom_common/mom_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_calc_visc.F

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

revision 1.25 by cnh, Fri Jul 13 16:36:41 2007 UTC revision 1.26 by heimbach, Fri Jul 27 22:18:58 2007 UTC
# Line 259  C-     Initialise to zero gradient of vo Line 259  C-     Initialise to zero gradient of vo
259  C      horizontal gradient of horizontal divergence:  C      horizontal gradient of horizontal divergence:
260    
261  C-       gradient in x direction:  C-       gradient in x direction:
262  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
263           IF (useCubedSphereExchange) THEN           IF (useCubedSphereExchange) THEN
264  C        to compute d/dx(hDiv), fill corners with appropriate values:  C        to compute d/dx(hDiv), fill corners with appropriate values:
265             CALL FILL_CS_CORNER_TR_RL( .TRUE., hDiv, bi,bj, myThid )             CALL FILL_CS_CORNER_TR_RL( .TRUE., hDiv, bi,bj, myThid )
266           ENDIF           ENDIF
267  #endif  cph-exch2#endif
268           DO j=2-Oly,sNy+Oly-1           DO j=2-Oly,sNy+Oly-1
269            DO i=2-Olx,sNx+Olx-1            DO i=2-Olx,sNx+Olx-1
270              divDx(i,j) = (hDiv(i,j)-hDiv(i-1,j))*recip_DXC(i,j,bi,bj)              divDx(i,j) = (hDiv(i,j)-hDiv(i-1,j))*recip_DXC(i,j,bi,bj)
# Line 272  C        to compute d/dx(hDiv), fill cor Line 272  C        to compute d/dx(hDiv), fill cor
272           ENDDO           ENDDO
273    
274  C-       gradient in y direction:  C-       gradient in y direction:
275  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
276           IF (useCubedSphereExchange) THEN           IF (useCubedSphereExchange) THEN
277  C        to compute d/dy(hDiv), fill corners with appropriate values:  C        to compute d/dy(hDiv), fill corners with appropriate values:
278             CALL FILL_CS_CORNER_TR_RL(.FALSE., hDiv, bi,bj, myThid )             CALL FILL_CS_CORNER_TR_RL(.FALSE., hDiv, bi,bj, myThid )
279           ENDIF           ENDIF
280  #endif  cph-exch2#endif
281           DO j=2-Oly,sNy+Oly-1           DO j=2-Oly,sNy+Oly-1
282            DO i=2-Olx,sNx+Olx-1            DO i=2-Olx,sNx+Olx-1
283              divDy(i,j) = (hDiv(i,j)-hDiv(i,j-1))*recip_DYC(i,j,bi,bj)              divDy(i,j) = (hDiv(i,j)-hDiv(i,j-1))*recip_DYC(i,j,bi,bj)
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22