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

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

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

revision 1.6 by jmc, Tue Oct 9 00:02:50 2007 UTC revision 1.7 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_hflux .NE. 0. ) then         if ( wmean_hflux .NE. 0. ) then
134            objf_hfluxmm = objf_hfluxmm            objf_hfluxmm = objf_hfluxmm
135       &        + (fctiletot/sumtot/wmean_hflux)**2       &        + ((fctiletot/sumtot+1.D0)/wmean_hflux)**2
136         else         else
137            objf_hfluxmm = 0. _d 0            objf_hfluxmm = 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_hfluxmm    =', irec,       &      ' --> bal_hfluxmm    =', irec,
143       &      fctiletot/sumtot,       &      fctiletot/sumtot,
144       &      (fctiletot/sumtot/wmean_hflux)**2       &      ((fctiletot/sumtot+1.D0)/wmean_hflux)**2
145    
146        enddo        enddo
147    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22