/[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.21 by jmc, Wed Jan 26 00:45:53 2005 UTC revision 1.22 by jmc, Fri Nov 4 01:19:24 2005 UTC
# Line 35  C     bi,bj  - Loop counters Line 35  C     bi,bj  - Loop counters
35  C     I,J,K  C     I,J,K
36  C     facGrid - Factor for grid to meter conversion  C     facGrid - Factor for grid to meter conversion
37        INTEGER bi, bj        INTEGER bi, bj
38        INTEGER  I,  J, K        INTEGER  I, J
39        _RL facGrid        _RL facGrid
40    
41  C     Initialise coriolis parameter  C     Initialise coriolis parameter
# Line 43  C     Initialise coriolis parameter Line 43  C     Initialise coriolis parameter
43  C       Constant F case  C       Constant F case
44          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
45            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
 c           DO K=1,Nr  
46                DO J=1-Oly,sNy+Oly                DO J=1-Oly,sNy+Oly
47                  DO I=1-Olx,sNx+Olx                  DO I=1-Olx,sNx+Olx
48                    fCori(i,j,bi,bj)=f0                    fCori(i,j,bi,bj)=f0
# Line 51  c           DO K=1,Nr Line 50  c           DO K=1,Nr
50                    fCoriCos(i,j,bi,bj)=0.                    fCoriCos(i,j,bi,bj)=0.
51                  ENDDO                  ENDDO
52                ENDDO                ENDDO
 c           ENDDO  
53            ENDDO            ENDDO
54          ENDDO          ENDDO
55        ELSEIF ( useBetaPlaneF ) THEN        ELSEIF ( useBetaPlaneF ) THEN
# Line 60  C       Beta plane case Line 58  C       Beta plane case
58          IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere          IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere
59          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
60            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
 c           DO K=1,Nr  
61                DO J=1-Oly,sNy+Oly                DO J=1-Oly,sNy+Oly
62                  DO I=1-Olx,sNx+Olx                  DO I=1-Olx,sNx+Olx
63                    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 68  c           DO K=1,Nr Line 65  c           DO K=1,Nr
65                    fCoriCos(i,j,bi,bj)=0.                    fCoriCos(i,j,bi,bj)=0.
66                  ENDDO                  ENDDO
67                ENDDO                ENDDO
 c           ENDDO  
68            ENDDO            ENDDO
69          ENDDO          ENDDO
70        ELSEIF ( useSphereF ) THEN        ELSEIF ( useSphereF ) THEN
# Line 76  C       Spherical case Line 72  C       Spherical case
72  C       Note in this case we assume yC is in degrees.  C       Note in this case we assume yC is in degrees.
73          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
74            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
 c           DO K=1,Nr  
75                DO J=1-Oly,sNy+Oly                DO J=1-Oly,sNy+Oly
76                  DO I=1-Olx,sNx+Olx                  DO I=1-Olx,sNx+Olx
77                    fCori(i,j,bi,bj)=                    fCori(i,j,bi,bj)=
# Line 87  c           DO K=1,Nr Line 82  c           DO K=1,Nr
82       &                 2. _d 0*omega*cos(_yC(i,j,bi,bj)*deg2rad)       &                 2. _d 0*omega*cos(_yC(i,j,bi,bj)*deg2rad)
83                  ENDDO                  ENDDO
84                ENDDO                ENDDO
 c           ENDDO  
85            ENDDO            ENDDO
86          ENDDO          ENDDO
87  c       CALL WRITE_FLD_XY_RL('fCoriC',' ',fCori ,  0,myThid)  c       CALL WRITE_FLD_XY_RL('fCoriC',' ',fCori ,  0,myThid)
# Line 97  c       CALL WRITE_FLD_XY_RL('fCorCs',' Line 91  c       CALL WRITE_FLD_XY_RL('fCorCs','
91  C       Special custom form  C       Special custom form
92          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
93            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
 c           DO K=1,Nr  
94                DO J=1-Oly,sNy+Oly                DO J=1-Oly,sNy+Oly
95                  DO I=1-Olx,sNx+Olx                  DO I=1-Olx,sNx+Olx
96                    fCori(i,j,bi,bj)=0.                    fCori(i,j,bi,bj)=0.
# Line 105  c           DO K=1,Nr Line 98  c           DO K=1,Nr
98                    fCoriCos(i,j,bi,bj)=0.                    fCoriCos(i,j,bi,bj)=0.
99                  ENDDO                  ENDDO
100                ENDDO                ENDDO
 c           ENDDO  
101            ENDDO            ENDDO
102          ENDDO          ENDDO
103          CALL READ_REC_XY_RS( 'fCoriC.bin', fCori,   1, 0, myThid )          CALL READ_REC_XY_RS( 'fCoriC.bin', fCori,   1, 0, myThid )
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22