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

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

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

revision 1.1 by gforget, Wed Oct 20 22:06:54 2010 UTC revision 1.5 by gforget, Fri Aug 10 14:34:41 2012 UTC
# Line 29  c     == routine arguments == Line 29  c     == routine arguments ==
29        _RL     mytime        _RL     mytime
30    
31  #ifdef ALLOW_ECCO_EVOLUTION  #ifdef ALLOW_ECCO_EVOLUTION
32          
33  #ifdef ALLOW_COST_ATLANTIC  #ifdef ALLOW_AUTODIFF_TAMC
34  c--   Compute meridional heat transport  # ifdef ALLOW_SSH_COST_CONTRIBUTION
35  #ifndef DISABLE_DEBUGMODE  #  ifdef ALLOW_PROFILES
36        IF ( debugLevel .GE. debLevB )  #   include "profiles.h"
37       &    CALL DEBUG_CALL('cost_atlantic',myThid)  #  endif
38  #endif  # endif
39        call timer_start('cost_atlantic      [ECCO SPIN-DOWN]', mythid)  # ifdef ALLOW_MEAN_HFLUX_COST_CONTRIBUTION
40        call cost_atlantic(  mytime, myiter,mythid )  #  include "ecco_cost.h"
41        call timer_stop ('cost_atlantic      [ECCO SPIN-DOWN]', mythid)  # endif
42  #endif  #endif
43    
44  c--   Compute the cost function contribution of the boundary forcing,  c--   Compute the cost function contribution of the boundary forcing,
45  c--   i.e. heat flux, salt flux, zonal and meridional wind stress.  c--   i.e. heat flux, salt flux, zonal and meridional wind stress.
46  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
47        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_forcing',myThid)
      &    CALL DEBUG_CALL('cost_forcing',myThid)  
48  #endif  #endif
49        call timer_start('cost_forcing       [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_forcing       [ECCO SPIN-DOWN]', mythid)
50        call cost_forcing( myiter, mytime, mythid )        call cost_forcing( myiter, mytime, mythid )
# Line 70  cph) Line 69  cph)
69    
70  c--   Compute cost function contribution of SSH.  c--   Compute cost function contribution of SSH.
71  #ifdef ALLOW_SSH_COST_CONTRIBUTION  #ifdef ALLOW_SSH_COST_CONTRIBUTION
72  # ifndef DISABLE_DEBUGMODE  # ifdef ALLOW_DEBUG
73        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_ssh',myThid)
      &    CALL DEBUG_CALL('cost_ssh',myThid)  
74  # endif  # endif
75  # if (defined(ALLOW_SSHV4_COST))  # if (defined(ALLOW_SSHV4_COST))
76        call timer_start('cost_sshv4         [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_sshv4         [ECCO SPIN-DOWN]', mythid)
# Line 82  c--   Compute cost function contribution Line 80  c--   Compute cost function contribution
80        call timer_start('cost_ssh_new       [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_ssh_new       [ECCO SPIN-DOWN]', mythid)
81        call cost_ssh_new( myiter, mytime, mythid )        call cost_ssh_new( myiter, mytime, mythid )
82        call timer_stop ('cost_ssh_new       [ECCO SPIN-DOWN]', mythid)        call timer_stop ('cost_ssh_new       [ECCO SPIN-DOWN]', mythid)
83  # else  # elif (!defined (ALLOW_GENCOST_SSHV4))
84        call timer_start('cost_ssh           [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_ssh           [ECCO SPIN-DOWN]', mythid)
85        call cost_ssh( myiter, mytime, mythid )        call cost_ssh( myiter, mytime, mythid )
86        call timer_stop ('cost_ssh           [ECCO SPIN-DOWN]', mythid)        call timer_stop ('cost_ssh           [ECCO SPIN-DOWN]', mythid)
# Line 102  c--   Compute bottom pressure cost Line 100  c--   Compute bottom pressure cost
100  #endif  #endif
101    
102  c--   Compute cost function contribution of Temperature and Salinity.  c--   Compute cost function contribution of Temperature and Salinity.
103  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
104        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_hyd',myThid)
      &    CALL DEBUG_CALL('cost_hyd',myThid)  
105  #endif  #endif
106        call timer_start('cost_hyd           [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_hyd           [ECCO SPIN-DOWN]', mythid)
107        call cost_hyd( myiter, mytime, mythid )        call cost_hyd( myiter, mytime, mythid )
108        call timer_stop ('cost_hyd           [ECCO SPIN-DOWN]', mythid)        call timer_stop ('cost_hyd           [ECCO SPIN-DOWN]', mythid)
109    
110  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
111  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
112        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('seaice_cost_driver',myThid)
      &    CALL DEBUG_CALL('seaice_cost_driver',myThid)  
113  #endif  #endif
114        IF ( useSeaice) THEN        IF ( useSeaice) THEN
115        call timer_start('seaice_cost_driver [ECCO SPIN-DOWN]', mythid)        call timer_start('seaice_cost_driver [ECCO SPIN-DOWN]', mythid)
# Line 123  c--   Compute cost function contribution Line 119  c--   Compute cost function contribution
119  #endif  #endif
120    
121  #ifdef ALLOW_OBCS_COST_CONTRIBUTION  #ifdef ALLOW_OBCS_COST_CONTRIBUTION
122  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
123        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_obcs',myThid)
      &    CALL DEBUG_CALL('cost_obcs',myThid)  
124  #endif  #endif
125        call timer_start('cost_obcs          [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_obcs          [ECCO SPIN-DOWN]', mythid)
126        call cost_obcs( myiter, mytime, mythid )        call cost_obcs( myiter, mytime, mythid )
# Line 133  c--   Compute cost function contribution Line 128  c--   Compute cost function contribution
128  #endif  #endif
129    
130  #ifdef ALLOW_CURMTR_COST_CONTRIBUTION  #ifdef ALLOW_CURMTR_COST_CONTRIBUTION
131  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
132        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_curmtr',myThid)
      &    CALL DEBUG_CALL('cost_curmtr',myThid)  
133  #endif  #endif
134        call timer_start('cost_curmtr        [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_curmtr        [ECCO SPIN-DOWN]', mythid)
135        call cost_curmtr( myiter, mytime, mythid )        call cost_curmtr( myiter, mytime, mythid )
# Line 145  c--   Compute cost function contribution Line 139  c--   Compute cost function contribution
139  c--   Compute cost function contribution of drifter velocities.  c--   Compute cost function contribution of drifter velocities.
140  #ifdef ALLOW_DRIFTER_COST_CONTRIBUTION  #ifdef ALLOW_DRIFTER_COST_CONTRIBUTION
141  #ifndef ALLOW_AUTODIFF_WHTAPEIO  #ifndef ALLOW_AUTODIFF_WHTAPEIO
142  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
143        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_drifter',myThid)
      &    CALL DEBUG_CALL('cost_drifter',myThid)  
144  #endif  #endif
145        call timer_start('cost_drifter         [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_drifter         [ECCO SPIN-DOWN]', mythid)
146        call cost_drifter( myiter, mytime, mythid )        call cost_drifter( myiter, mytime, mythid )
# Line 158  c--   Compute cost function contribution Line 151  c--   Compute cost function contribution
151  c--   Compute cost function contribution of wind stress observations.  c--   Compute cost function contribution of wind stress observations.
152  #if (defined (ALLOW_SCAT_COST_CONTRIBUTION) || \  #if (defined (ALLOW_SCAT_COST_CONTRIBUTION) || \
153       defined (ALLOW_DAILYSCAT_COST_CONTRIBUTION) )       defined (ALLOW_DAILYSCAT_COST_CONTRIBUTION) )
154  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
155        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_scat',myThid)
      &    CALL DEBUG_CALL('cost_scat',myThid)  
156  #endif  #endif
157        call timer_start('cost_scat            [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_scat            [ECCO SPIN-DOWN]', mythid)
158        call cost_scat( myiter, mytime, mythid )        call cost_scat( myiter, mytime, mythid )
# Line 170  c--   Compute cost function contribution Line 162  c--   Compute cost function contribution
162  c--   Compute cost function contribution of drift between the first  c--   Compute cost function contribution of drift between the first
163  c     and the last year.  c     and the last year.
164  #ifdef ALLOW_DRIFT_COST_CONTRIBUTION  #ifdef ALLOW_DRIFT_COST_CONTRIBUTION
165  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
166        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_drift',myThid)
      &    CALL DEBUG_CALL('cost_drift',myThid)  
167  #endif  #endif
168        call timer_start('cost_drift         [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_drift         [ECCO SPIN-DOWN]', mythid)
169        call cost_drift( myiter, mytime, mythid )        call cost_drift( myiter, mytime, mythid )
170        call timer_stop ('cost_drift         [ECCO SPIN-DOWN]', mythid)        call timer_stop ('cost_drift         [ECCO SPIN-DOWN]', mythid)
171  #endif  #endif
172  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION  #ifdef ALLOW_DRIFTW_COST_CONTRIBUTION
173  #ifndef DISABLE_DEBUGMODE  #ifdef ALLOW_DEBUG
174        IF ( debugLevel .GE. debLevB )        IF (debugMode) CALL DEBUG_CALL('cost_driftw',myThid)
      &    CALL DEBUG_CALL('cost_driftw',myThid)  
175  #endif  #endif
176        call timer_start('cost_driftw        [ECCO SPIN-DOWN]', mythid)        call timer_start('cost_driftw        [ECCO SPIN-DOWN]', mythid)
177        call cost_driftw( myiter, mytime, mythid )        call cost_driftw( myiter, mytime, mythid )
# Line 218  c--   Compute user defined cost function Line 208  c--   Compute user defined cost function
208    
209        return        return
210        end        end
   
   
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22