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

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

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

revision 1.4 by jmc, Sun Nov 21 16:06:05 2004 UTC revision 1.8 by jmc, Thu Aug 16 02:13:40 2007 UTC
# Line 13  C     !INTERFACE: Line 13  C     !INTERFACE:
13       I             myThid )       I             myThid )
14  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
15  C     *==========================================================*  C     *==========================================================*
16  C     | SUBROUTINE CALC_ISOSLOPES                                  C     | SUBROUTINE GRAD_SIGMA
17  C     | o Calculate isoneutral gradients                            C     | o Calculate isoneutral gradients
18  C     *==========================================================*  C     *==========================================================*
19  C     \ev  C     \ev
20    
# Line 56  C-    safer to work on a local copy of r Line 56  C-    safer to work on a local copy of r
56        ENDDO        ENDDO
57    
58  C-    Internal exchange for calculations in X  C-    Internal exchange for calculations in X
59    cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
60        IF ( useCubedSphereExchange ) THEN        IF ( useCubedSphereExchange ) THEN
61          CALL FILL_CS_CORNER_TR_RL( .TRUE., rhoLoc, bi,bj, myThid )          CALL FILL_CS_CORNER_TR_RL( .TRUE., .FALSE.,
62         &                             rhoLoc, bi,bj, myThid )
63        ENDIF        ENDIF
64    cph-exch2#endif
65        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
66         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
67          sigmaX(i,j,k)=_maskW(i,j,k,bi,bj)          sigmaX(i,j,k)=_maskW(i,j,k,bi,bj)
# Line 68  C-    Internal exchange for calculations Line 71  C-    Internal exchange for calculations
71        ENDDO        ENDDO
72    
73  C-    Internal exchange for calculations in Y  C-    Internal exchange for calculations in Y
74    cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
75        IF ( useCubedSphereExchange ) THEN        IF ( useCubedSphereExchange ) THEN
76          CALL FILL_CS_CORNER_TR_RL( .FALSE., rhoLoc, bi,bj, myThid )          CALL FILL_CS_CORNER_TR_RL(.FALSE., .FALSE.,
77         &                             rhoLoc, bi,bj, myThid )
78        ENDIF        ENDIF
79    cph-exch2#endif
80        DO j=1-Oly+1,sNy+Oly        DO j=1-Oly+1,sNy+Oly
81         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
82          sigmaY(i,j,k)=_maskS(i,j,k,bi,bj)          sigmaY(i,j,k)=_maskS(i,j,k,bi,bj)
# Line 79  C-    Internal exchange for calculations Line 85  C-    Internal exchange for calculations
85         ENDDO         ENDDO
86        ENDDO        ENDDO
87    
 C- jmc: Patrick, will this work for TAF ?  
88        IF (K.EQ.1) THEN        IF (K.EQ.1) THEN
89         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
90          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
# Line 90  C- jmc: Patrick, will this work for TAF Line 95  C- jmc: Patrick, will this work for TAF
95         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
96          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
97           sigmaR(i,j,k)= maskC(i,j,k,bi,bj)           sigmaR(i,j,k)= maskC(i,j,k,bi,bj)
98       &                *recip_drC(k)*rkFac       &                *recip_drC(k)*rkSign
99       &                *(sigKm1(i,j)-sigKp1(i,j))       &                *(sigKp1(i,j)-sigKm1(i,j))
100          ENDDO          ENDDO
101         ENDDO         ENDDO
102        ENDIF        ENDIF
 C- jmc: leave the old code commented (in case ...)  
 c     DO j=1-Oly,sNy+Oly  
 c      DO i=1-Olx,sNx+Olx  
 c       IF (K.NE.1 .AND. hfacC(i,j,k,bi,bj).NE.0.) THEN  
 c        sigmaR(i,j,k)=recip_drC(k)*rkFac*(sigKm1(i,j)-sigKp1(i,j))  
 c       ELSE  
 c        sigmaR(i,j,k)=0.  
 c       ENDIF  
 c      ENDDO  
 c     ENDDO  
103    
104        RETURN        RETURN
105        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22