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

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

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

revision 1.2 by jmc, Sat Nov 5 01:00:57 2005 UTC revision 1.3 by heimbach, Wed Jun 7 01:55:12 2006 UTC
# Line 66  C     i,j :: Loop counters Line 66  C     i,j :: Loop counters
66        IF ( DIAGNOSTICS_IS_ON('URHOMASS',myThid) ) THEN        IF ( DIAGNOSTICS_IS_ON('URHOMASS',myThid) ) THEN
67          DO j=1,sNy          DO j=1,sNy
68           DO i=1,sNx+1           DO i=1,sNx+1
69             tmpFld(i,j) = uVel(i,j,k,bi,bj)*hFacW(i,j,k,bi,bj)             tmpFld(i,j) = uVel(i,j,k,bi,bj)*_hFacW(i,j,k,bi,bj)
70       &                 *(rhoK(i-1,j)+rhoK(i,j))*0.5 _d 0       &                 *(rhoK(i-1,j)+rhoK(i,j))*0.5 _d 0
71           ENDDO           ENDDO
72          ENDDO          ENDDO
# Line 76  C     i,j :: Loop counters Line 76  C     i,j :: Loop counters
76        IF ( DIAGNOSTICS_IS_ON('VRHOMASS',myThid) ) THEN        IF ( DIAGNOSTICS_IS_ON('VRHOMASS',myThid) ) THEN
77          DO j=1,sNy+1          DO j=1,sNy+1
78           DO i=1,sNx           DO i=1,sNx
79             tmpFld(i,j) = vVel(i,j,k,bi,bj)*hFacS(i,j,k,bi,bj)             tmpFld(i,j) = vVel(i,j,k,bi,bj)*_hFacS(i,j,k,bi,bj)
80       &                 *(rhoK(i,j-1)+rhoK(i,j))*0.5 _d 0       &                 *(rhoK(i,j-1)+rhoK(i,j))*0.5 _d 0
81           ENDDO           ENDDO
82          ENDDO          ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22