pkg/rbcs/rbcs_init_varia.F

c3afacc940 Jean*0001 #include "RBCS_OPTIONS.h"
429c85036b Patr*0002
                0003 C !INTERFACE: ==========================================================
c3afacc940 Jean*0004       SUBROUTINE RBCS_INIT_VARIA( myThid )
429c85036b Patr*0005
                0006 C !DESCRIPTION:
                0007 C calls subroutines that initialized variables for relaxed
                0008 c boundary conditions
                0009
                0010 C !USES: ===============================================================
                0011       IMPLICIT NONE
                0012 #include "SIZE.h"
                0013 #include "EEPARAMS.h"
                0014 #include "PARAMS.h"
                0015 #ifdef ALLOW_PTRACERS
                0016 #include "PTRACERS_SIZE.h"
                0017 #endif
a16c4403c6 Jean*0018 #include "RBCS_SIZE.h"
                0019 c#include "RBCS_PARAMS.h"
                0020 #include "RBCS_FIELDS.h"
429c85036b Patr*0021
                0022 C !INPUT PARAMETERS: ===================================================
                0023 C  myThid               :: thread number
                0024       INTEGER myThid
                0025 CEOP
                0026
                0027 #ifdef ALLOW_RBCS
                0028
                0029 C     !LOCAL VARIABLES:
                0030 C     i,j,k,bi,bj,iTracer  :: loop indices
                0031       INTEGER i,j,k,bi,bj
                0032 #ifdef ALLOW_PTRACERS
                0033       INTEGER iTracer
                0034 #endif
                0035
869c678fd8 Patr*0036       DO bj = myByLo(myThid), myByHi(myThid)
                0037         DO bi = myBxLo(myThid), myBxHi(myThid)
72a8e44ea5 Jean*0038           rbcsLdRec(bi,bj) = 0
869c678fd8 Patr*0039           DO k=1,Nr
                0040             DO j=1-Oly,sNy+OLy
                0041               DO i=1-Olx,sNx+Olx
a16c4403c6 Jean*0042 #ifndef DISABLE_RBCS_MOM
                0043                 rbcu0(i,j,k,bi,bj) = 0. _d 0
                0044                 rbcv0(i,j,k,bi,bj) = 0. _d 0
                0045                 rbcu1(i,j,k,bi,bj) = 0. _d 0
                0046                 rbcv1(i,j,k,bi,bj) = 0. _d 0
                0047                 RBCuVel(i,j,k,bi,bj) = 0. _d 0
                0048                 RBCvVel(i,j,k,bi,bj) = 0. _d 0
                0049 #endif
c3afacc940 Jean*0050                 rbct0(i,j,k,bi,bj) = 0. _d 0
                0051                 rbcs0(i,j,k,bi,bj) = 0. _d 0
                0052                 rbct1(i,j,k,bi,bj) = 0. _d 0
                0053                 rbcs1(i,j,k,bi,bj) = 0. _d 0
                0054                 RBCtemp(i,j,k,bi,bj) = 0. _d 0
                0055                 RBCsalt(i,j,k,bi,bj) = 0. _d 0
429c85036b Patr*0056               ENDDO
869c678fd8 Patr*0057             ENDDO
429c85036b Patr*0058           ENDDO
869c678fd8 Patr*0059         ENDDO
                0060       ENDDO
429c85036b Patr*0061
                0062 #ifdef ALLOW_PTRACERS
                0063 C     Loop over tracers
                0064       DO iTracer = 1, PTRACERS_num
                0065         DO bj = myByLo(myThid), myByHi(myThid)
c3afacc940 Jean*0066          DO bi = myBxLo(myThid), myBxHi(myThid)
                0067            DO k=1,Nr
                0068              DO j=1-Oly,sNy+OLy
                0069                DO i=1-Olx,sNx+Olx
                0070                 rbcptr0(i,j,k,bi,bj,iTracer) = 0. _d 0
                0071                 rbcptr1(i,j,k,bi,bj,iTracer) = 0. _d 0
                0072                 RBC_ptracers(i,j,k,bi,bj,iTracer) = 0. _d 0
429c85036b Patr*0073                ENDDO
                0074               ENDDO
                0075             ENDDO
                0076           ENDDO
                0077         ENDDO
                0078 C       end of Tracer loop
                0079       ENDDO
c3afacc940 Jean*0080 #endif /* ALLOW_PTRACERS */
429c85036b Patr*0081
                0082 #endif /* ALLOW_RBCS */
                0083
                0084       RETURN
                0085       END