/[MITgcm]/MITgcm/pkg/ecco/cost_mean_saltflux.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_mean_saltflux.F

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

revision 1.12 by jmc, Tue Apr 28 18:13:28 2009 UTC revision 1.13 by mlosch, Wed Sep 30 15:48:31 2009 UTC
# Line 115  c--     Read time averages and the month Line 115  c--     Read time averages and the month
115                  tmpx=tmpfld2d(i,j,bi,bj)                  tmpx=tmpfld2d(i,j,bi,bj)
116                  if (maskC(i,j,kk,bi,bj) .ne. 0.) then                  if (maskC(i,j,kk,bi,bj) .ne. 0.) then
117                     fctiletot = fctiletot                     fctiletot = fctiletot
118       &                + tmpx*cos(yc(i,j,bi,bj)*deg2rad)/rhoConstFresh       &                + tmpx* _rA(i,j,bi,bj)/rhoConstFresh
119                     sumtot = sumtot                     sumtot = sumtot
120       &                + cos(yc(i,j,bi,bj)*deg2rad)       &                + _rA(i,j,bi,bj)
121                     num_sfluxmm(bi,bj) = num_sfluxmm(bi,bj) + 1                     num_sfluxmm(bi,bj) = num_sfluxmm(bi,bj) + 1
122                  endif                  endif
123                enddo                enddo
# Line 141  c-- diagnostic: imbalance per year: Line 141  c-- diagnostic: imbalance per year:
141         write(standardmessageunit,'(A,I5,2(X,D22.14))')         write(standardmessageunit,'(A,I5,2(X,D22.14))')
142       &      ' --> bal_sfluxmm    =', irec,       &      ' --> bal_sfluxmm    =', irec,
143       &      fctiletot/sumtot,       &      fctiletot/sumtot,
144       &      ((fctiletot/sumtot)/wmean_sflux)**2       &      objf_sfluxmm
145    
146        enddo        enddo
147    
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22