/[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.13 by jmc, Tue May 24 20:48:28 2011 UTC revision 1.16 by jmc, Tue Sep 18 20:21:23 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_OPTIONS.h"
5  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
6  # include "OBCS_OPTIONS.h"  # include "OBCS_OPTIONS.h"
7  #endif  #endif
# Line 37  c     == global variables == Line 37  c     == global variables ==
37  c#include "OBCS_PARAMS.h"  c#include "OBCS_PARAMS.h"
38  #include "OBCS_GRID.h"  #include "OBCS_GRID.h"
39  #include "OBCS_FIELDS.h"  #include "OBCS_FIELDS.h"
40    #include "CTRL_SIZE.h"
41  #include "ctrl.h"  #include "ctrl.h"
42  #include "ctrl_dummy.h"  #include "ctrl_dummy.h"
43  #include "optim.h"  #include "optim.h"
# Line 130  c--   Get the counters, flags, and the i Line 130  c--   Get the counters, flags, and the i
130            if (iobcs .gt. 2) then            if (iobcs .gt. 2) then
131             do j = jmin,jmax             do j = jmin,jmax
132              i = OB_Iw(j,bi,bj)              i = OB_Iw(j,bi,bj)
133                IF ( i.EQ.OB_indexNone ) i = 1
134  cih    Determine number of open vertical layers.  cih    Determine number of open vertical layers.
135              nz = 0              nz = 0
136              do k = 1,Nr              do k = 1,Nr
# Line 198  cgg   & Line 199  cgg   &
199            if (iobcs .gt. 2) then            if (iobcs .gt. 2) then
200             do j = jmin,jmax             do j = jmin,jmax
201              i = OB_Iw(j,bi,bj)              i = OB_Iw(j,bi,bj)
202                IF ( i.EQ.OB_indexNone ) i = 1
203  cih    Determine number of open vertical layers.  cih    Determine number of open vertical layers.
204              nz = 0              nz = 0
205              do k = 1,Nr              do k = 1,Nr
# Line 249  c--   Calculate mask for tracer cells (0 Line 251  c--   Calculate mask for tracer cells (0
251           do k = 1,nr           do k = 1,nr
252            do j = 1,sny            do j = 1,sny
253             i = OB_Iw(j,bi,bj)             i = OB_Iw(j,bi,bj)
254               IF ( i.EQ.OB_indexNone ) i = 1
255             if (iobcs .EQ. 1) then             if (iobcs .EQ. 1) then
256              OBWt(j,k,bi,bj) = OBWt (j,k,bi,bj)              OBWt(j,k,bi,bj) = OBWt (j,k,bi,bj)
257       &           + obcswfac            *xx_obcsw0(j,k,bi,bj,iobcs)       &           + obcswfac            *xx_obcsw0(j,k,bi,bj,iobcs)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22