/[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.32 by jmc, Mon Apr 19 15:10:05 2010 UTC revision 1.33 by jmc, Fri Nov 12 03:16:16 2010 UTC
# Line 30  C     !USES: Line 30  C     !USES:
30  #endif  #endif
31    
32  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
33    C     myThid  :: my Thread Id number
34        INTEGER myThid        INTEGER myThid
35  CEOP  CEOP
36    
# Line 54  C     facGrid :: Factor for grid to mete Line 55  C     facGrid :: Factor for grid to mete
55    
56    
57  C     Initialise coriolis parameter  C     Initialise coriolis parameter
58        IF     ( useConstantF ) THEN        IF     ( selectCoriMap.EQ.0 ) THEN
59  C       Constant F case  C       Constant F case
60          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
61            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 62  C       Constant F case Line 63  C       Constant F case
63                DO i=1-Olx,sNx+Olx                DO i=1-Olx,sNx+Olx
64                  fCori(i,j,bi,bj)  = f0                  fCori(i,j,bi,bj)  = f0
65                  fCoriG(i,j,bi,bj) = f0                  fCoriG(i,j,bi,bj) = f0
66                  fCoriCos(i,j,bi,bj)=0. _d 0                  fCoriCos(i,j,bi,bj)=fPrime
67                ENDDO                ENDDO
68              ENDDO              ENDDO
69            ENDDO            ENDDO
70          ENDDO          ENDDO
71        ELSEIF ( useBetaPlaneF ) THEN        ELSEIF ( selectCoriMap.EQ.1 ) THEN
72  C       Beta plane case  C       Beta plane case
73          facGrid = 1. _d 0          facGrid = 1. _d 0
74          IF ( usingSphericalPolarGrid          IF ( usingSphericalPolarGrid
# Line 78  C       Beta plane case Line 79  C       Beta plane case
79                DO i=1-Olx,sNx+Olx                DO i=1-Olx,sNx+Olx
80                  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
81                  fCoriG(i,j,bi,bj) = f0+beta* yG(i,j,bi,bj)*facGrid                  fCoriG(i,j,bi,bj) = f0+beta* yG(i,j,bi,bj)*facGrid
82                  fCoriCos(i,j,bi,bj)=0. _d 0                  fCoriCos(i,j,bi,bj)=fPrime
83                ENDDO                ENDDO
84              ENDDO              ENDDO
85            ENDDO            ENDDO
86          ENDDO          ENDDO
87        ELSEIF ( useSphereF ) THEN        ELSEIF ( selectCoriMap.EQ.2 ) THEN
88  C       Spherical case  C       Spherical case
89  C       Note in this case we assume yC is in degrees.  C       Note in this case we assume yC is in degrees.
90          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
# Line 104  c       CALL WRITE_FLD_XY_RL('fCoriC',' Line 105  c       CALL WRITE_FLD_XY_RL('fCoriC','
105  c       CALL WRITE_FLD_XY_RL('fCoriG',' ',fCoriG , 0,myThid)  c       CALL WRITE_FLD_XY_RL('fCoriG',' ',fCoriG , 0,myThid)
106  c       CALL WRITE_FLD_XY_RL('fCorCs',' ',fCoriCos,0,myThid)  c       CALL WRITE_FLD_XY_RL('fCorCs',' ',fCoriCos,0,myThid)
107        ELSE        ELSE
108  C       Special custom form  C       Initialise to zero
109          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
110            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
111              DO j=1-Oly,sNy+Oly              DO j=1-Oly,sNy+Oly
# Line 116  C       Special custom form Line 117  C       Special custom form
117              ENDDO              ENDDO
118            ENDDO            ENDDO
119          ENDDO          ENDDO
120          ENDIF
121    
122          IF ( selectCoriMap.EQ.3 ) THEN
123    C     Special custom form: read from files
124          CALL READ_REC_XY_RS( 'fCoriC.bin', fCori,   1, 0, myThid )          CALL READ_REC_XY_RS( 'fCoriC.bin', fCori,   1, 0, myThid )
125          CALL READ_REC_XY_RS( 'fCorCs.bin', fCoriCos,1, 0, myThid )          CALL READ_REC_XY_RS( 'fCorCs.bin', fCoriCos,1, 0, myThid )
126          IF ( .NOT.useCubedSphereExchange ) THEN          IF ( .NOT.useCubedSphereExchange ) THEN
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22