/[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.7 by adcroft, Wed Dec 9 16:11:52 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_CORI( myThid )        SUBROUTINE INI_CORI( myThid )
# Line 8  C     /================================= Line 8  C     /=================================
8  C     | SUBROUTINE INI_CORI                                      |  C     | SUBROUTINE INI_CORI                                      |
9  C     | o Initialise coriolis term.                              |  C     | o Initialise coriolis term.                              |
10  C     \==========================================================/  C     \==========================================================/
11          IMPLICIT NONE
12    
13  C     === Global variables ===  C     === Global variables ===
14  #include "SIZE.h"  #include "SIZE.h"
# Line 36  C--   Initialise coriolis parameter Line 37  C--   Initialise coriolis parameter
37  C      o Constant F case  C      o Constant F case
38         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
39          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
40           DO K=1,Nz           DO K=1,Nr
41            DO J=1,sNy            DO J=1,sNy
42             DO I=1,sNx             DO I=1,sNx
43              fCori(i,j,bi,bj)=f0              fCori(i,j,bi,bj)=f0
# 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,sNy
57             DO I=1,sNx             DO I=1,sNx
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
# Line 65  C      o Spherical case Line 66  C      o Spherical case
66  C        Note in this case we assume yC is in degrees.  C        Note in this case we assume yC is in degrees.
67         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
68          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
69           DO K=1,Nz           DO K=1,Nr
70            DO J=1,sNy            DO J=1,sNy
71             DO I=1,sNx             DO I=1,sNx
72              fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)              fCori(i,j,bi,bj)=
73         &       2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)
74             ENDDO             ENDDO
75            ENDDO            ENDDO
76           ENDDO           ENDDO
# Line 78  C        Note in this case we assume yC Line 80  C        Note in this case we assume yC
80  C      o Special custom form  C      o Special custom form
81         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
82          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
83           DO K=1,Nz           DO K=1,Nr
84            DO J=1,sNy            DO J=1,sNy
85             DO I=1,sNx             DO I=1,sNx
86              fCori(i,j,bi,bj)=0.              fCori(i,j,bi,bj)=0.
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22