/[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.4 by cnh, Sat Aug 22 17:51:08 1998 UTC
# 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,sNy
41             DO I=1,sNx             DO I=1,sNx
42              fCori(i,j,bi,bj)=f0              fCori(i,j,bi,bj)=f0
# Line 51  C      o Beta plane case Line 51  C      o Beta plane case
51         IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere         IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere
52         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
53          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
54           DO K=1,Nz           DO K=1,Nr
55            DO J=1,sNy            DO J=1,sNy
56             DO I=1,sNx             DO I=1,sNx
57              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 65  C      o Spherical case
65  C        Note in this case we assume yC is in degrees.  C        Note in this case we assume yC is in degrees.
66         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
67          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
68           DO K=1,Nz           DO K=1,Nr
69            DO J=1,sNy            DO J=1,sNy
70             DO I=1,sNx             DO I=1,sNx
71              fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)              fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)
# Line 78  C        Note in this case we assume yC Line 78  C        Note in this case we assume yC
78  C      o Special custom form  C      o Special custom form
79         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
80          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
81           DO K=1,Nz           DO K=1,Nr
82            DO J=1,sNy            DO J=1,sNy
83             DO I=1,sNx             DO I=1,sNx
84              fCori(i,j,bi,bj)=0.              fCori(i,j,bi,bj)=0.
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22