model/src/taueddy_init_varia.F

C $Header: /u/gcmpack/MITgcm/model/src/taueddy_init_varia.F,v 1.2 2008/05/30 02:45:43 gforget Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: TAUEDDY_INIT_VARIA
C     !INTERFACE:
      SUBROUTINE TAUEDDY_INIT_VARIA( myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE TAUEDDY_INIT_VARIA                                    
C     | o Initialise to zero some active arrays  
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "FFIELDS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance of INI_UVEL 
      INTEGER myThid

#ifdef ALLOW_EDDYPSI

C     !LOCAL VARIABLES:
C     == Local variables ==
C     bi,bj  - Loop counters
C     I,J,K
      INTEGER bi, bj
      INTEGER  I, J, K
CEOP

C--   Over all tiles
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO K=1,Nr
         DO J=1-Oly,sNy+Oly
          DO I=1-Olx,sNx+Olx
           eddyPsiX(i,j,k,bi,bj)=0.
           eddyPsiY(i,j,k,bi,bj)=0.
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      IF ( eddyPsiXFile .NE. ' ' ) THEN
       CALL MDSREADFIELD(eddyPsiXFile,readBinaryPrec,
     &                  'RS',Nr,eddyPsiX,1,myThid)
      ENDIF
      IF ( eddyPsiYFile .NE. ' ' ) THEN
       CALL MDSREADFIELD(eddyPsiYFile,readBinaryPrec,
     &                  'RS',Nr,eddyPsiY,1,myThid)
      ENDIF

      CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)

#endif /* ALLOW_AUTODIFF_TAMC */

      END
1	dfer	1.3	C $Header: /u/gcmpack/MITgcm/model/src/taueddy_init_varia.F,v 1.2 2008/05/30 02:45:43 gforget Exp $
2	heimbach	1.1	C $Name: $
3
4			#include "PACKAGES_CONFIG.h"
5			#include "CPP_OPTIONS.h"
6
7			CBOP
8			C !ROUTINE: TAUEDDY_INIT_VARIA
9			C !INTERFACE:
10			SUBROUTINE TAUEDDY_INIT_VARIA( myThid )
11			C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| SUBROUTINE TAUEDDY_INIT_VARIA
14			C \| o Initialise to zero some active arrays
15			C ==========================================================
16			C \ev
17
18			C !USES:
19			IMPLICIT NONE
20			C === Global variables ===
21			#include "SIZE.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "GRID.h"
25			#include "FFIELDS.h"
26
27			C !INPUT/OUTPUT PARAMETERS:
28			C == Routine arguments ==
29			C myThid - Number of this instance of INI_UVEL
30			INTEGER myThid
31
32	gforget	1.2	#ifdef ALLOW_EDDYPSI
33	heimbach	1.1
34			C !LOCAL VARIABLES:
35			C == Local variables ==
36			C bi,bj - Loop counters
37			C I,J,K
38			INTEGER bi, bj
39			INTEGER I, J, K
40			CEOP
41
42			C-- Over all tiles
43			DO bj = myByLo(myThid), myByHi(myThid)
44			DO bi = myBxLo(myThid), myBxHi(myThid)
45			DO K=1,Nr
46			DO J=1-Oly,sNy+Oly
47			DO I=1-Olx,sNx+Olx
48	gforget	1.2	eddyPsiX(i,j,k,bi,bj)=0.
49			eddyPsiY(i,j,k,bi,bj)=0.
50	heimbach	1.1	ENDDO
51			ENDDO
52			ENDDO
53			ENDDO
54			ENDDO
55
56	dfer	1.3	IF ( eddyPsiXFile .NE. ' ' ) THEN
57			CALL MDSREADFIELD(eddyPsiXFile,readBinaryPrec,
58	gforget	1.2	& 'RS',Nr,eddyPsiX,1,myThid)
59	heimbach	1.1	ENDIF
60	dfer	1.3	IF ( eddyPsiYFile .NE. ' ' ) THEN
61			CALL MDSREADFIELD(eddyPsiYFile,readBinaryPrec,
62	gforget	1.2	& 'RS',Nr,eddyPsiY,1,myThid)
63	heimbach	1.1	ENDIF
64
65	gforget	1.2	CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
66	heimbach	1.1
67			#endif /* ALLOW_AUTODIFF_TAMC */
68
69			END