/[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.2 by cnh, Fri Apr 24 02:05:41 1998 UTC revision 1.3 by cnh, Mon Jun 8 21:43:01 1998 UTC
# Line 31  C     facGrid - Factor for grid to meter Line 31  C     facGrid - Factor for grid to meter
31        INTEGER  I,  J, K        INTEGER  I,  J, K
32        _RL facGrid        _RL facGrid
33    
34  C--   Initialise beta plane f.  C--   Initialise coriolis parameter
35        facGrid = 1. _d 0        IF     ( useConstantF ) THEN
36        IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere  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,Nz
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+beta*yC(i,j,bi,bj)*facGrid              fCori(i,j,bi,bj)=f0
43               ENDDO
44            ENDDO            ENDDO
45           ENDDO           ENDDO
46          ENDDO          ENDDO
47         ENDDO         ENDDO
48        ENDDO        ELSEIF ( useBetaPlaneF ) THEN
49    C      o Beta plane case
50           facGrid = 1. _d 0
51           IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere
52           DO bj = myByLo(myThid), myByHi(myThid)
53            DO bi = myBxLo(myThid), myBxHi(myThid)
54             DO K=1,Nz
55              DO J=1,sNy
56               DO I=1,sNx
57                fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid
58               ENDDO
59              ENDDO
60             ENDDO
61            ENDDO
62           ENDDO
63          ELSEIF ( useSphereF ) THEN
64    C      o Spherical case
65    C        Note in this case we assume yC is in degrees.
66           DO bj = myByLo(myThid), myByHi(myThid)
67            DO bi = myBxLo(myThid), myBxHi(myThid)
68             DO K=1,Nz
69              DO J=1,sNy
70               DO I=1,sNx
71                fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)
72               ENDDO
73              ENDDO
74             ENDDO
75            ENDDO
76           ENDDO
77          ELSE
78    C      o Special custom form
79           DO bj = myByLo(myThid), myByHi(myThid)
80            DO bi = myBxLo(myThid), myBxHi(myThid)
81             DO K=1,Nz
82              DO J=1,sNy
83               DO I=1,sNx
84                fCori(i,j,bi,bj)=0.
85               ENDDO
86              ENDDO
87             ENDDO
88            ENDDO
89           ENDDO
90          ENDIF
91  C  C
92        _EXCH_XY_R4(fCori , myThid )        _EXCH_XY_R4(fCori , myThid )
93    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22