/[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.9 by jmc, Tue Oct 9 00:02:50 2007 UTC revision 1.10 by heimbach, Tue Feb 5 12:49:10 2008 UTC
# Line 132  c--     Read time averages and the month Line 132  c--     Read time averages and the month
132    
133         if ( wmean_sflux .NE. 0. ) then         if ( wmean_sflux .NE. 0. ) then
134            objf_sfluxmm = objf_sfluxmm            objf_sfluxmm = objf_sfluxmm
135       &        + ( (fctiletot/sumtot-1.D0)/wmean_sflux )**2       &        + ( (fctiletot/sumtot)/wmean_sflux )**2
136         else         else
137            objf_sfluxmm = 0. _d 0            objf_sfluxmm = 0. _d 0
138         endif         endif
# 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       &      ((fctiletot/sumtot)/wmean_sflux)**2
145    
146        enddo        enddo
147    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22