/[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.9 by cnh, Sun Feb 4 14:38:47 2001 UTC revision 1.10 by adcroft, Tue May 29 14:01:37 2001 UTC
# Line 37  C      o Constant F case Line 37  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,Nr           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 52  C      o Beta plane case Line 53  C      o Beta plane case
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,Nr           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 66  C        Note in this case we assume yC Line 68  C        Note in this case we assume yC
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,Nr           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)=              fCori(i,j,bi,bj)=
74       &       2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)       &       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 80  C      o Special custom form Line 84  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,Nr           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 90  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.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22