/[MITgcm]/MITgcm/model/src/ini_spherical_polar_grid.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_spherical_polar_grid.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by cnh, Mon May 25 18:01:32 1998 UTC revision 1.5 by cnh, Mon Jun 8 21:43:01 1998 UTC
# Line 81  C     Set up my local grid first Line 81  C     Set up my local grid first
81  C     Note: In the spherical polar case delX and delY are given in  C     Note: In the spherical polar case delX and delY are given in
82  C           degrees and are relative to some starting latitude and  C           degrees and are relative to some starting latitude and
83  C           longitude - phiMin and thetaMin.  C           longitude - phiMin and thetaMin.
84          xC0 = thetaMin
85          yC0 = phiMin
86        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
87         jG = myYGlobalLo + (bj-1)*sNy         jG = myYGlobalLo + (bj-1)*sNy
88         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 197  C     Calculate vertical face area Line 199  C     Calculate vertical face area
199          ENDDO          ENDDO
200         ENDDO         ENDDO
201        ENDDO        ENDDO
202        _EXCH_XY_R4(zA, myThid )  
203          DO bj = myByLo(myThid), myByHi(myThid)
204           DO bi = myBxLo(myThid), myBxHi(myThid)
205            DO K=1,Nz
206             DO J=1,sNy
207              DO I=1,sNx
208               IF (HFacC(I,J,K,bi,bj) .NE. 0. D0 ) THEN
209                rHFacC(I,J,K,bi,bj) = 1. D0 / HFacC(I,J,K,bi,bj)
210               ELSE
211                rHFacC(I,J,K,bi,bj) = 0. D0
212               ENDIF
213               IF (HFacW(I,J,K,bi,bj) .NE. 0. D0 ) THEN
214                rHFacW(I,J,K,bi,bj) = 1. D0 / HFacW(I,J,K,bi,bj)
215                maskW(I,J,K,bi,bj) = 1. D0
216               ELSE
217                rHFacW(I,J,K,bi,bj) = 0. D0
218                maskW(I,J,K,bi,bj) = 0.0 D0
219               ENDIF
220               IF (HFacS(I,J,K,bi,bj) .NE. 0. D0 ) THEN
221                rHFacS(I,J,K,bi,bj) = 1. D0 / HFacS(I,J,K,bi,bj)
222                maskS(I,J,K,bi,bj) = 1. D0
223               ELSE
224                rHFacS(I,J,K,bi,bj) = 0. D0
225                maskS(I,J,K,bi,bj) = 0. D0
226               ENDIF
227              ENDDO
228             ENDDO
229            ENDDO
230           ENDDO
231          ENDDO
232    C     Now sync. and get/send edge regions that are shared with
233    C     other threads.
234          _EXCH_XYZ_R4(rHFacC    , myThid )
235          _EXCH_XYZ_R4(rHFacW    , myThid )
236          _EXCH_XYZ_R4(rHFacS    , myThid )
237          _EXCH_XYZ_R4(maskW    , myThid )
238          _EXCH_XYZ_R4(maskS    , myThid )
239          _EXCH_XY_R4 (zA       , myThid )
240    
241    CcnhDebugStarts
242          tanPhiAtU = 0. _d 0
243          tanPhiAtV = 0. _d 0
244          _EXCH_XY_R4 (tanPhiAtU , myThid )
245          _EXCH_XY_R4 (tanPhiAtV , myThid )
246    CcnhDebugEnds
247    
248    
249  C  C
250        RETURN        RETURN
251        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22