C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_init_obcs_variables.F,v 1.10 2014/10/09 00:49:27 gforget Exp $ C $Name: $ #include "CTRL_OPTIONS.h" subroutine CTRL_INIT_OBCS_VARIABLES( mythid ) c ================================================================== c SUBROUTINE ctrl_init_obcs_variables c ================================================================== c c o Set parts of the vector of control variables and initialize the c rest to zero. c c started: heimbach@mit.edu 25-Mar-2002 c c ================================================================== c SUBROUTINE ctrl_init_obcs_variables c ================================================================== implicit none c == global variables == #include "EEPARAMS.h" #include "SIZE.h" #include "PARAMS.h" c#include "GRID.h" #include "ctrl.h" #include "CTRL_OBCS.h" c == routine arguments == integer mythid c == local variables == integer bi,bj integer i,j,k integer itlo,ithi integer jtlo,jthi integer jmin,jmax integer imin,imax integer ntmp integer ivarindex integer iobcs c == end of interface == jtlo = mybylo(mythid) jthi = mybyhi(mythid) itlo = mybxlo(mythid) ithi = mybxhi(mythid) jmin = 1-oly jmax = sny+oly imin = 1-olx imax = snx+olx #ifdef ALLOW_OBCSN_CONTROL do iobcs = 1, nobcs do bj = jtlo,jthi do bi = itlo,ithi do k = 1,nr do i = imin,imax xx_obcsn0(i,k,bi,bj,iobcs) = 0. _d 0 xx_obcsn1(i,k,bi,bj,iobcs) = 0. _d 0 enddo enddo enddo enddo enddo #endif #ifdef ALLOW_OBCSS_CONTROL do iobcs = 1, nobcs do bj = jtlo,jthi do bi = itlo,ithi do k = 1,nr do i = imin,imax xx_obcss0(i,k,bi,bj,iobcs) = 0. _d 0 xx_obcss1(i,k,bi,bj,iobcs) = 0. _d 0 enddo enddo enddo enddo enddo #endif #ifdef ALLOW_OBCSW_CONTROL do iobcs = 1, nobcs do bj = jtlo,jthi do bi = itlo,ithi do k = 1,nr do j = jmin,jmax xx_obcsw0(j,k,bi,bj,iobcs) = 0. _d 0 xx_obcsw1(j,k,bi,bj,iobcs) = 0. _d 0 enddo enddo enddo enddo enddo #endif #ifdef ALLOW_OBCSE_CONTROL do iobcs = 1, nobcs do bj = jtlo,jthi do bi = itlo,ithi do k = 1,nr do j = jmin,jmax xx_obcse0(j,k,bi,bj,iobcs) = 0. _d 0 xx_obcse1(j,k,bi,bj,iobcs) = 0. _d 0 enddo enddo enddo enddo enddo #endif return end