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

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

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

revision 1.15 by jmc, Fri Aug 10 19:38:57 2012 UTC revision 1.17 by gforget, Thu Oct 9 00:49:26 2014 UTC
# Line 40  c#include "OBCS_PARAMS.h" Line 40  c#include "OBCS_PARAMS.h"
40  #include "CTRL_SIZE.h"  #include "CTRL_SIZE.h"
41  #include "ctrl.h"  #include "ctrl.h"
42  #include "ctrl_dummy.h"  #include "ctrl_dummy.h"
43    #include "CTRL_OBCS.h"
44  #include "optim.h"  #include "optim.h"
45  #endif /* ALLOW_OBCSW_CONTROL */  #endif /* ALLOW_OBCSW_CONTROL */
46    
# Line 130  c--   Get the counters, flags, and the i Line 131  c--   Get the counters, flags, and the i
131            if (iobcs .gt. 2) then            if (iobcs .gt. 2) then
132             do j = jmin,jmax             do j = jmin,jmax
133              i = OB_Iw(j,bi,bj)              i = OB_Iw(j,bi,bj)
134                IF ( i.EQ.OB_indexNone ) i = 1
135  cih    Determine number of open vertical layers.  cih    Determine number of open vertical layers.
136              nz = 0              nz = 0
137              do k = 1,Nr              do k = 1,Nr
# Line 198  cgg   & Line 200  cgg   &
200            if (iobcs .gt. 2) then            if (iobcs .gt. 2) then
201             do j = jmin,jmax             do j = jmin,jmax
202              i = OB_Iw(j,bi,bj)              i = OB_Iw(j,bi,bj)
203                IF ( i.EQ.OB_indexNone ) i = 1
204  cih    Determine number of open vertical layers.  cih    Determine number of open vertical layers.
205              nz = 0              nz = 0
206              do k = 1,Nr              do k = 1,Nr
# Line 249  c--   Calculate mask for tracer cells (0 Line 252  c--   Calculate mask for tracer cells (0
252           do k = 1,nr           do k = 1,nr
253            do j = 1,sny            do j = 1,sny
254             i = OB_Iw(j,bi,bj)             i = OB_Iw(j,bi,bj)
255               IF ( i.EQ.OB_indexNone ) i = 1
256             if (iobcs .EQ. 1) then             if (iobcs .EQ. 1) then
257              OBWt(j,k,bi,bj) = OBWt (j,k,bi,bj)              OBWt(j,k,bi,bj) = OBWt (j,k,bi,bj)
258       &           + obcswfac            *xx_obcsw0(j,k,bi,bj,iobcs)       &           + obcswfac            *xx_obcsw0(j,k,bi,bj,iobcs)
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22