/[MITgcm]/MITgcm/model/src/taueddy_init_varia.F
ViewVC logotype

Diff of /MITgcm/model/src/taueddy_init_varia.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by dfer, Fri Jun 5 14:27:43 2009 UTC revision 1.4 by jmc, Wed Sep 9 15:31:42 2009 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
 #include "PACKAGES_CONFIG.h"  
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CBOP  CBOP
# Line 10  C     !INTERFACE: Line 9  C     !INTERFACE:
9        SUBROUTINE TAUEDDY_INIT_VARIA( myThid )        SUBROUTINE TAUEDDY_INIT_VARIA( myThid )
10  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
11  C     *==========================================================*  C     *==========================================================*
12  C     | SUBROUTINE TAUEDDY_INIT_VARIA                                      C     | SUBROUTINE TAUEDDY_INIT_VARIA
13  C     | o Initialise to zero some active arrays    C     | o Initialise to zero some active arrays
14  C     *==========================================================*  C     *==========================================================*
15  C     \ev  C     \ev
16    
# Line 21  C     === Global variables === Line 20  C     === Global variables ===
20  #include "SIZE.h"  #include "SIZE.h"
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
 #include "GRID.h"  
23  #include "FFIELDS.h"  #include "FFIELDS.h"
24    
25  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
26  C     == Routine arguments ==  C     == Routine arguments ==
27  C     myThid -  Number of this instance of INI_UVEL  C     myThid :: my Thread Id number
28        INTEGER myThid        INTEGER myThid
29    
30  #ifdef ALLOW_EDDYPSI  #ifdef ALLOW_EDDYPSI
# Line 34  C     myThid -  Number of this instance Line 32  C     myThid -  Number of this instance
32  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
33  C     == Local variables ==  C     == Local variables ==
34  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
35  C     I,J,K  C     i,j,k
36        INTEGER bi, bj        INTEGER bi, bj
37        INTEGER  I, J, K        INTEGER i, j, k
38  CEOP  CEOP
39    
40  C--   Over all tiles  C--   Over all tiles
41        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
42         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
43          DO K=1,Nr          DO k=1,Nr
44           DO J=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
45            DO I=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
46             eddyPsiX(i,j,k,bi,bj)=0.             eddyPsiX(i,j,k,bi,bj)=0.
47             eddyPsiY(i,j,k,bi,bj)=0.             eddyPsiY(i,j,k,bi,bj)=0.
48            ENDDO            ENDDO
# Line 54  C--   Over all tiles Line 52  C--   Over all tiles
52        ENDDO        ENDDO
53    
54        IF ( eddyPsiXFile .NE. ' ' ) THEN        IF ( eddyPsiXFile .NE. ' ' ) THEN
55         CALL MDSREADFIELD(eddyPsiXFile,readBinaryPrec,         CALL READ_FLD_XYZ_RS( eddyPsiXFile, ' ',
56       &                  'RS',Nr,eddyPsiX,1,myThid)       &                       eddyPsiX, nIter0, myThid )
57        ENDIF        ENDIF
58        IF ( eddyPsiYFile .NE. ' ' ) THEN        IF ( eddyPsiYFile .NE. ' ' ) THEN
59         CALL MDSREADFIELD(eddyPsiYFile,readBinaryPrec,         CALL READ_FLD_XYZ_RS( eddyPsiYFile, ' ',
60       &                  'RS',Nr,eddyPsiY,1,myThid)       &                       eddyPsiY, nIter0, myThid )
61        ENDIF        ENDIF
62    
63        CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)        CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
64    
65  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_EDDYPSI */
66    
67          RETURN
68        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22