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

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

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

revision 1.3 by cnh, Mon Jun 8 21:43:01 1998 UTC revision 1.9.2.2 by adcroft, Mon Apr 9 19:13:42 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE INI_CORI( myThid )        SUBROUTINE INI_CORI( myThid )
# Line 8  C     /================================= Line 9  C     /=================================
9  C     | SUBROUTINE INI_CORI                                      |  C     | SUBROUTINE INI_CORI                                      |
10  C     | o Initialise coriolis term.                              |  C     | o Initialise coriolis term.                              |
11  C     \==========================================================/  C     \==========================================================/
12          IMPLICIT NONE
13    
14  C     === Global variables ===  C     === Global variables ===
15  #include "SIZE.h"  #include "SIZE.h"
# Line 22  C     myThid -  Number of this instance Line 24  C     myThid -  Number of this instance
24  CEndOfInterface  CEndOfInterface
25    
26  C     == Local variables ==  C     == Local variables ==
 C     iG, jG - Global coordinate index  
27  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
28  C     I,J,K  C     I,J,K
29  C     facGrid - Factor for grid to meter conversion  C     facGrid - Factor for grid to meter conversion
       INTEGER iG, jG  
30        INTEGER bi, bj        INTEGER bi, bj
31        INTEGER  I,  J, K        INTEGER  I,  J, K
32        _RL facGrid        _RL facGrid
# Line 36  C--   Initialise coriolis parameter Line 36  C--   Initialise coriolis parameter
36  C      o Constant F case  C      o Constant F case
37         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
38          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
39           DO K=1,Nz           DO K=1,Nr
40            DO J=1,sNy            DO J=1-Oly,sNy+Oly
41             DO I=1,sNx             DO I=1-Olx,sNx+Olx
42              fCori(i,j,bi,bj)=f0              fCori(i,j,bi,bj)=f0
43                fCoriG(i,j,bi,bj)=f0
44             ENDDO             ENDDO
45            ENDDO            ENDDO
46           ENDDO           ENDDO
# Line 51  C      o Beta plane case Line 52  C      o Beta plane case
52         IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere         IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere
53         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
54          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
55           DO K=1,Nz           DO K=1,Nr
56            DO J=1,sNy            DO J=1-Oly,sNy+Oly
57             DO I=1,sNx             DO I=1-Olx,sNx+Olx
58              fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid              fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid
59                fCoriG(i,j,bi,bj)=f0+beta*yG(i,j,bi,bj)*facGrid
60             ENDDO             ENDDO
61            ENDDO            ENDDO
62           ENDDO           ENDDO
# Line 65  C      o Spherical case Line 67  C      o Spherical case
67  C        Note in this case we assume yC is in degrees.  C        Note in this case we assume yC is in degrees.
68         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
69          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
70           DO K=1,Nz           DO K=1,Nr
71            DO J=1,sNy            DO J=1-Oly,sNy+Oly
72             DO I=1,sNx             DO I=1-Olx,sNx+Olx
73              fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)              fCori(i,j,bi,bj)=
74         &       2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)
75                fCoriG(i,j,bi,bj)=
76         &       2. _d 0*omega*sin(yG(i,j,bi,bj)*deg2rad)
77             ENDDO             ENDDO
78            ENDDO            ENDDO
79           ENDDO           ENDDO
# Line 78  C        Note in this case we assume yC Line 83  C        Note in this case we assume yC
83  C      o Special custom form  C      o Special custom form
84         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
85          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
86           DO K=1,Nz           DO K=1,Nr
87            DO J=1,sNy            DO J=1-Oly,sNy+Oly
88             DO I=1,sNx             DO I=1-Olx,sNx+Olx
89              fCori(i,j,bi,bj)=0.              fCori(i,j,bi,bj)=0.
90                fCoriG(i,j,bi,bj)=0.
91             ENDDO             ENDDO
92            ENDDO            ENDDO
93           ENDDO           ENDDO
# Line 89  C      o Special custom form Line 95  C      o Special custom form
95         ENDDO         ENDDO
96        ENDIF        ENDIF
97  C  C
98        _EXCH_XY_R4(fCori , myThid )  c     _EXCH_XY_R4(fCori,myThid)
99    c     CALL EXCH_Z_XY_RS(fCoriG,myThid)
100    
101        RETURN        RETURN
102        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.9.2.2
  ViewVC Help
Powered by ViewVC 1.1.22