/[MITgcm]/MITgcm/pkg/cost/cost_atlantic_heat.F
ViewVC logotype

Diff of /MITgcm/pkg/cost/cost_atlantic_heat.F

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

revision 1.4 by heimbach, Mon Sep 27 18:00:19 2004 UTC revision 1.5 by heimbach, Fri Dec 10 20:15:10 2004 UTC
# Line 100  c Line 100  c
100       &                      *maskS(i,j,k,bi,bj)       &                      *maskS(i,j,k,bi,bj)
101       &                      *maskC(i,j,k,bi,bj)*maskC(i,j-1,k,bi,bj)       &                      *maskC(i,j,k,bi,bj)*maskC(i,j-1,k,bi,bj)
102               theta_bar(k) = theta_bar(k) +               theta_bar(k) = theta_bar(k) +
 cph     &                      theta(i,j,k,bi,bj)*dxG(i,j,bi,bj)  
 cph     &                      *maskC(i,j,k,bi,bj)  
103       &            0.5*( cMeanTheta(i,j,k,bi,bj)       &            0.5*( cMeanTheta(i,j,k,bi,bj)
104       &                 +cMeanTheta(i,j-1,k,bi,bj) )       &                 +cMeanTheta(i,j-1,k,bi,bj) )
105       &                 *maskS(i,j,k,bi,bj)*dxG(i,j,bi,bj)       &                 *maskS(i,j,k,bi,bj)*dxG(i,j,bi,bj)
# Line 146  c Line 144  c
144       &                      *maskW(i,j,k,bi,bj)       &                      *maskW(i,j,k,bi,bj)
145       &                      *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)       &                      *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
146               theta_bar(k) = theta_bar(k) +               theta_bar(k) = theta_bar(k) +
 cph     &                      + theta(i,j,k,bi,bj)*dyG(i,j,bi,bj)  
 cph     &                      *maskC(i,j,k,bi,bj)  
147       &            0.5*( cMeanTheta(i,j,k,bi,bj)       &            0.5*( cMeanTheta(i,j,k,bi,bj)
148       &                 +cMeanTheta(i-,j,k,bi,bj) )       &                 +cMeanTheta(i-,j,k,bi,bj) )
149       &                 *maskW(i,j,k,bi,bj)*dyG(i,j,bi,bj)       &                 *maskW(i,j,k,bi,bj)*dyG(i,j,bi,bj)

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22