/[MITgcm]/MITgcm/pkg/ctrl/ctrl_obcsbal.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_obcsbal.F

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

revision 1.4 by mlosch, Fri Dec 3 00:48:57 2004 UTC revision 1.5 by jmc, Tue Oct 9 00:00:00 2007 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
5  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
# Line 14  c     ================================== Line 16  c     ==================================
16  c     SUBROUTINE ctrl_obcsbal  c     SUBROUTINE ctrl_obcsbal
17  c     ==================================================================  c     ==================================================================
18  c  c
19  c     o volumetrically balance the control vector contribution.  c     o volumetrically balance the control vector contribution.
20  c     o Assume the calendar is identical  c     o Assume the calendar is identical
21  c       for all open boundaries. Need to save the barotropic adjustment  c       for all open boundaries. Need to save the barotropic adjustment
22  c       velocity so it can be used in all ctrl_getobcs files.  c       velocity so it can be used in all ctrl_getobcs files.
23  c     o WARNING: eastern boundary (not defined) filenames have been a  c     o WARNING: eastern boundary (not defined) filenames have been a
24  c                problem in the past.  c                problem in the past.
25  c  c
26  c     - started G. Gebbie, MIT-WHOI, 15-June-2002  c     - started G. Gebbie, MIT-WHOI, 15-June-2002
# Line 159  c--   Do the global summation. Line 161  c--   Do the global summation.
161    
162          shiftvel(2) = volflux / area          shiftvel(2) = volflux / area
163          print*,'volflux,area',volflux,area          print*,'volflux,area',volflux,area
164        endif          endif
165  cgg    End of the obcsnfirst loop.  cgg    End of the obcsnfirst loop.
166    
167        if ( ( obcsnfirst) .or. (obcsnchanged)) then        if ( ( obcsnfirst) .or. (obcsnchanged)) then
168    
169  cgg     Swap the value.  cgg     Swap the value.
170          shiftvel(1) = shiftvel(2)          shiftvel(1) = shiftvel(2)
171    
172          volflux = 0. d0          volflux = 0. d0
173          area= 0. d0          area= 0. d0
174    
# Line 179  c--   Do the global summation. Line 181  c--   Do the global summation.
181          shiftvel(2) = volflux /area          shiftvel(2) = volflux /area
182          print*,'volflux,area',volflux,area          print*,'volflux,area',volflux,area
183    
184        endif          endif
185  cgg    End of the obcsnfirst, obcsnchanged loop.  cgg    End of the obcsnfirst, obcsnchanged loop.
186    
187  #endif  #endif
188    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22