C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/nesting_sannino/nest_child/nest_child_init_varia.F,v 1.3 2010/11/28 02:01:49 jmc Exp $
C $Name:  $

#include "NEST_CHILD_OPTIONS.h"

CBOP 0
C     !ROUTINE: NEST_CHILD_INIT_VARIA
C     !INTERFACE:

      SUBROUTINE NEST_CHILD_INIT_VARIA( myThid )

C     !DESCRIPTION:
C     Routine to initialize NEST_CHILD variables.

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "NEST_CHILD.h"
#include "NEST_CHILD_PARAMS.h"

C     !INPUT PARAMETERS:
C     myThid   :: my Thread Id number
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
C     bi,bj   :: tile indices
C     i,j,k   :: loop index
      INTEGER bi, bj
      INTEGER i, j, k

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

C--   Initialise passage counter
      _BEGIN_MASTER(myThid)
       PASSI = 0
      _END_MASTER(myThid)

C--   Initialise all arrays in common blocks (from NEST_CHILD.h)

C-    Vertical Section arrays (i=1 : West section , i=2 : Eastern section)
C     note: tile indices are missing
      DO i=1,2
        DO k=1,Nr
         DO j=1-Oly,sNy+Oly
           U_F1      (j,k,i) = 0. _d 0
           V_F1      (j,k,i) = 0. _d 0
           T_F1      (j,k,i) = 0. _d 0
           S_F1      (j,k,i) = 0. _d 0
           ETA_F1    (j,k,i) = 0. _d 0
           dETAhdt_F1(j,k,i) = 0. _d 0
           gU_F1     (j,k,i) = 0. _d 0
           gV_F1     (j,k,i) = 0. _d 0
           gT_F1     (j,k,i) = 0. _d 0
           gS_F1     (j,k,i) = 0. _d 0
           guNm1_F1  (j,k,i) = 0. _d 0
           gvNm1_F1  (j,k,i) = 0. _d 0
           gtNm1_F1  (j,k,i) = 0. _d 0
           gsNm1_F1  (j,k,i) = 0. _d 0
         ENDDO
        ENDDO
      ENDDO

C--   Over all tiles
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)

C-      3D arrays
        DO k=1,Nr
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
           guNm1_MEMO(i,j,k,bi,bj) = 0. _d 0
           gvNm1_MEMO(i,j,k,bi,bj) = 0. _d 0
           gtNm1_MEMO(i,j,k,bi,bj) = 0. _d 0
           gsNm1_MEMO(i,j,k,bi,bj) = 0. _d 0
          ENDDO
         ENDDO
        ENDDO

C-      2D arrays
        DO j=1-Oly,sNy+Oly
         DO i=1-Olx,sNx+Olx
           dEtaHdt_MEMO(i,j,bi,bj) = 0. _d 0
         ENDDO
        ENDDO

C-- end bi,bj loops
       ENDDO
      ENDDO

      RETURN
      END