/[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.5 by jmc, Mon Nov 22 16:27:05 2004 UTC revision 1.6 by jmc, Wed Jun 22 00:25:32 2005 UTC
# Line 83  C-    Internal exchange for calculations Line 83  C-    Internal exchange for calculations
83         ENDDO         ENDDO
84        ENDDO        ENDDO
85    
 C- jmc: Patrick, will this work for TAF ?  
86        IF (K.EQ.1) THEN        IF (K.EQ.1) THEN
87         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
88          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
# Line 94  C- jmc: Patrick, will this work for TAF Line 93  C- jmc: Patrick, will this work for TAF
93         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
94          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
95           sigmaR(i,j,k)= maskC(i,j,k,bi,bj)           sigmaR(i,j,k)= maskC(i,j,k,bi,bj)
96       &                *recip_drC(k)*rkFac       &                *recip_drC(k)*rkSign
97       &                *(sigKm1(i,j)-sigKp1(i,j))       &                *(sigKp1(i,j)-sigKm1(i,j))
98          ENDDO          ENDDO
99         ENDDO         ENDDO
100        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  
101    
102        RETURN        RETURN
103        END        END
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22