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

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

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

revision 1.1 by heimbach, Thu Nov 6 22:10:07 2003 UTC revision 1.2 by jmc, Tue Oct 9 00:02:50 2007 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
5  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
# Line 116  c     Number of records to be used. Line 118  c     Number of records to be used.
118        sumvol = 0. _d 0        sumvol = 0. _d 0
119        wobcsvol = .01        wobcsvol = .01
120  cgg   Acceptable volume flux is 10^-3. Corresponds to 5 mm change over a year.  cgg   Acceptable volume flux is 10^-3. Corresponds to 5 mm change over a year.
121  cgg   Added a factor of 1000 because its very important to me.  cgg   Added a factor of 1000 because its very important to me.
122        wobcsvol = 1./(wobcsvol * wobcsvol)        wobcsvol = 1./(wobcsvol * wobcsvol)
123    
124  #ifdef ECCO_VERBOSE  #ifdef ECCO_VERBOSE
# Line 140  cgg   Added a factor of 1000 because its Line 142  cgg   Added a factor of 1000 because its
142        if (optimcycle .ge. 0) then        if (optimcycle .ge. 0) then
143  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
144          ilfld=ilnblnk( xx_obcsn_file )          ilfld=ilnblnk( xx_obcsn_file )
145          write(fnamefldn(1:80),'(2a,i10.10)')          write(fnamefldn(1:80),'(2a,i10.10)')
146       &        xx_obcsn_file(1:ilfld),'.', optimcycle       &        xx_obcsn_file(1:ilfld),'.', optimcycle
147  #endif  #endif
148  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
149          ilfld=ilnblnk( xx_obcss_file )          ilfld=ilnblnk( xx_obcss_file )
150          write(fnameflds(1:80),'(2a,i10.10)')          write(fnameflds(1:80),'(2a,i10.10)')
151       &        xx_obcss_file(1:ilfld),'.',optimcycle       &        xx_obcss_file(1:ilfld),'.',optimcycle
152  #endif  #endif
153  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
154          ilfld=ilnblnk( xx_obcsw_file )          ilfld=ilnblnk( xx_obcsw_file )
155          write(fnamefldw(1:80),'(2a,i10.10)')          write(fnamefldw(1:80),'(2a,i10.10)')
156       &        xx_obcsw_file(1:ilfld),'.',optimcycle       &        xx_obcsw_file(1:ilfld),'.',optimcycle
157  #endif  #endif
158  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
159          ilfld=ilnblnk( xx_obcse_file )          ilfld=ilnblnk( xx_obcse_file )
160          write(fnameflde(1:80),'(2a,i10.10)')          write(fnameflde(1:80),'(2a,i10.10)')
161       &        xx_obcse_file(1:ilfld),'.',optimcycle       &        xx_obcse_file(1:ilfld),'.',optimcycle
162  #endif  #endif
163        else        else
# Line 185  cgg   think of a more general way to do Line 187  cgg   think of a more general way to do
187       &                       , xx_obcsn_dummy )       &                       , xx_obcsn_dummy )
188    
189  cgg  At this point, do not be concerned with the overlap halos.  cgg  At this point, do not be concerned with the overlap halos.
190  cgg  From experience, there is no control contribution in the  cgg  From experience, there is no control contribution in the
191  cgg  velocity points outside the boundaries. This has something  cgg  velocity points outside the boundaries. This has something
192  cgg  to do with the computational stencil, and the fact that we  cgg  to do with the computational stencil, and the fact that we
193  cgg  are diagonally offset. Could check later by employing both  cgg  are diagonally offset. Could check later by employing both
# Line 228  cgg   think of a more general way to do Line 230  cgg   think of a more general way to do
230       &                       , xx_obcss_dummy )       &                       , xx_obcss_dummy )
231    
232  cgg  At this point, do not be concerned with the overlap halos.  cgg  At this point, do not be concerned with the overlap halos.
233  cgg  From experience, there is no control contribution in the  cgg  From experience, there is no control contribution in the
234  cgg  velocity points outside the boundaries. This has something  cgg  velocity points outside the boundaries. This has something
235  cgg  to do with the computational stencil, and the fact that we  cgg  to do with the computational stencil, and the fact that we
236  cgg  are diagonally offset. Could check later by employing both  cgg  are diagonally offset. Could check later by employing both
# Line 271  cgg   think of a more general way to do Line 273  cgg   think of a more general way to do
273       &                       , xx_obcsw_dummy )       &                       , xx_obcsw_dummy )
274    
275  cgg  At this point, do not be concerned with the overlap halos.  cgg  At this point, do not be concerned with the overlap halos.
276  cgg  From experience, there is no control contribution in the  cgg  From experience, there is no control contribution in the
277  cgg  velocity points outside the boundaries. This has something  cgg  velocity points outside the boundaries. This has something
278  cgg  to do with the computational stencil, and the fact that we  cgg  to do with the computational stencil, and the fact that we
279  cgg  are diagonally offset. Could check later by employing both  cgg  are diagonally offset. Could check later by employing both
# Line 314  cgg   think of a more general way to do Line 316  cgg   think of a more general way to do
316       &                       , xx_obcse_dummy )       &                       , xx_obcse_dummy )
317    
318  cgg  At this point, do not be concerned with the overlap halos.  cgg  At this point, do not be concerned with the overlap halos.
319  cgg  From experience, there is no control contribution in the  cgg  From experience, there is no control contribution in the
320  cgg  velocity points outside the boundaries. This has something  cgg  velocity points outside the boundaries. This has something
321  cgg  to do with the computational stencil, and the fact that we  cgg  to do with the computational stencil, and the fact that we
322  cgg  are diagonally offset. Could check later by employing both  cgg  are diagonally offset. Could check later by employing both
# Line 351  c--   End of loop over records. Line 353  c--   End of loop over records.
353    
354  c--   Do the global summation.  c--   Do the global summation.
355        _GLOBAL_SUM_R8( sumvol, mythid )        _GLOBAL_SUM_R8( sumvol, mythid )
356        objf_obcsvol =  wobcsvol * sumvol* sumvol        objf_obcsvol =  wobcsvol * sumvol* sumvol
357    
358  #endif  #endif
359  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22