pkg/ctrl/ctrl_init_obcs_variables.F

7bfe6112e8 Jean*0001 #include "CTRL_OPTIONS.h"
7109a141b2 Patr*0002
9f5240b52a Jean*0003       subroutine ctrl_init_obcs_variables( myThid )
7109a141b2 Patr*0004
                0005 c     ==================================================================
                0006 c     SUBROUTINE ctrl_init_obcs_variables
                0007 c     ==================================================================
                0008 c
                0009 c     o Set parts of the vector of control variables and initialize the
                0010 c       rest to zero.
                0011 c
                0012 c     started: heimbach@mit.edu 25-Mar-2002
                0013 c
                0014 c     ==================================================================
                0015 c     SUBROUTINE ctrl_init_obcs_variables
                0016 c     ==================================================================
                0017
                0018       implicit none
                0019
                0020 c     == global variables ==
                0021
                0022 #include "EEPARAMS.h"
                0023 #include "SIZE.h"
                0024 #include "PARAMS.h"
7ef594fa1c Jean*0025 c#include "GRID.h"
5cf4364659 Mart*0026 #include "CTRL_SIZE.h"
4d72283393 Mart*0027 #include "CTRL.h"
e612621177 Gael*0028 #include "CTRL_OBCS.h"
7109a141b2 Patr*0029
                0030 c     == routine arguments ==
9f5240b52a Jean*0031       integer myThid
7109a141b2 Patr*0032
7c50f07931 Mart*0033 #ifdef ALLOW_OBCS
7109a141b2 Patr*0034 c     == local variables ==
                0035       integer bi,bj
                0036       integer i,j,k
                0037       integer iobcs
                0038
                0039 c     == end of interface ==
                0040
                0041 #ifdef ALLOW_OBCSN_CONTROL
                0042       do iobcs = 1, nobcs
9f5240b52a Jean*0043         do bj = myByLo(myThid), myByHi(myThid)
                0044           do bi = myBxLo(myThid), myBxHi(myThid)
                0045             do k = 1,Nr
4d72283393 Mart*0046               do i = 1-OLx, sNx+OLx
7109a141b2 Patr*0047                 xx_obcsn0(i,k,bi,bj,iobcs) = 0. _d 0
                0048                 xx_obcsn1(i,k,bi,bj,iobcs) = 0. _d 0
                0049               enddo
                0050             enddo
                0051           enddo
                0052         enddo
                0053       enddo
                0054 #endif
                0055
                0056 #ifdef ALLOW_OBCSS_CONTROL
                0057       do iobcs = 1, nobcs
9f5240b52a Jean*0058         do bj = myByLo(myThid), myByHi(myThid)
                0059           do bi = myBxLo(myThid), myBxHi(myThid)
                0060             do k = 1,Nr
4d72283393 Mart*0061               do i = 1-OLx, sNx+OLx
7109a141b2 Patr*0062                 xx_obcss0(i,k,bi,bj,iobcs) = 0. _d 0
                0063                 xx_obcss1(i,k,bi,bj,iobcs) = 0. _d 0
                0064               enddo
                0065             enddo
                0066           enddo
                0067         enddo
                0068       enddo
                0069 #endif
                0070
                0071 #ifdef ALLOW_OBCSW_CONTROL
                0072       do iobcs = 1, nobcs
9f5240b52a Jean*0073         do bj = myByLo(myThid), myByHi(myThid)
                0074           do bi = myBxLo(myThid), myBxHi(myThid)
                0075             do k = 1,Nr
4d72283393 Mart*0076               do j = 1-OLy, sNy+OLy
7109a141b2 Patr*0077                 xx_obcsw0(j,k,bi,bj,iobcs) = 0. _d 0
                0078                 xx_obcsw1(j,k,bi,bj,iobcs) = 0. _d 0
                0079               enddo
                0080             enddo
                0081           enddo
                0082         enddo
                0083       enddo
                0084 #endif
                0085
                0086 #ifdef ALLOW_OBCSE_CONTROL
                0087       do iobcs = 1, nobcs
9f5240b52a Jean*0088         do bj = myByLo(myThid), myByHi(myThid)
                0089           do bi = myBxLo(myThid), myBxHi(myThid)
                0090             do k = 1,Nr
4d72283393 Mart*0091               do j = 1-OLy, sNy+OLy
7109a141b2 Patr*0092                 xx_obcse0(j,k,bi,bj,iobcs) = 0. _d 0
                0093                 xx_obcse1(j,k,bi,bj,iobcs) = 0. _d 0
                0094               enddo
                0095             enddo
                0096           enddo
                0097         enddo
                0098       enddo
                0099 #endif
7c50f07931 Mart*0100 #endif /* ALLOW_OBCS */
7109a141b2 Patr*0101
                0102       return
                0103       end