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

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

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

revision 1.6 by cnh, Thu Jun 25 20:43:23 1998 UTC revision 1.7 by adcroft, Thu Jul 2 14:16:24 1998 UTC
# Line 151  C     dxC, dyC is separation between cel Line 151  C     dxC, dyC is separation between cel
151        ENDDO        ENDDO
152        _EXCH_XY_R4(dxC, myThid )        _EXCH_XY_R4(dxC, myThid )
153        _EXCH_XY_R4(dyC, myThid )        _EXCH_XY_R4(dyC, myThid )
 C     Calculate recipricols  
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO J=1,sNy  
          DO I=1,sNx  
           rDxG(I,J,bi,bj)=1.d0/dxG(I,J,bi,bj)  
           rDyG(I,J,bi,bj)=1.d0/dyG(I,J,bi,bj)  
           rDxC(I,J,bi,bj)=1.d0/dxC(I,J,bi,bj)  
           rDyC(I,J,bi,bj)=1.d0/dyC(I,J,bi,bj)  
           rDxF(I,J,bi,bj)=1.d0/dxF(I,J,bi,bj)  
           rDyF(I,J,bi,bj)=1.d0/dyF(I,J,bi,bj)  
           rDxV(I,J,bi,bj)=1.d0/dxV(I,J,bi,bj)  
           rDyU(I,J,bi,bj)=1.d0/dyU(I,J,bi,bj)  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       _EXCH_XY_R4(rDxG, myThid )  
       _EXCH_XY_R4(rDyG, myThid )  
       _EXCH_XY_R4(rDxC, myThid )  
       _EXCH_XY_R4(rDyC, myThid )  
       _EXCH_XY_R4(rDxF, myThid )  
       _EXCH_XY_R4(rDyF, myThid )  
       _EXCH_XY_R4(rDxV, myThid )  
       _EXCH_XY_R4(rDyU, myThid )  
154  C     Calculate vertical face area  C     Calculate vertical face area
155        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
156         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
157          DO J=1,sNy          DO J=1,sNy
158           DO I=1,sNx           DO I=1,sNx
159            zA(I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj)            zA(I,J,bi,bj) = dxF(I,J,bi,bj)*dyF(I,J,bi,bj)
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
   
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO K=1,Nz  
          DO J=1,sNy  
           DO I=1,sNx  
            IF (HFacC(I,J,K,bi,bj) .NE. 0. D0 ) THEN  
             rHFacC(I,J,K,bi,bj) = 1. D0 / HFacC(I,J,K,bi,bj)  
            ELSE  
             rHFacC(I,J,K,bi,bj) = 0. D0  
            ENDIF  
            IF (HFacW(I,J,K,bi,bj) .NE. 0. D0 ) THEN  
             rHFacW(I,J,K,bi,bj) = 1. D0 / HFacW(I,J,K,bi,bj)  
             maskW(I,J,K,bi,bj) = 1. D0  
            ELSE  
             rHFacW(I,J,K,bi,bj) = 0. D0  
             maskW(I,J,K,bi,bj) = 0.0 D0  
            ENDIF  
            IF (HFacS(I,J,K,bi,bj) .NE. 0. D0 ) THEN  
             rHFacS(I,J,K,bi,bj) = 1. D0 / HFacS(I,J,K,bi,bj)  
             maskS(I,J,K,bi,bj) = 1. D0  
            ELSE  
             rHFacS(I,J,K,bi,bj) = 0. D0  
             maskS(I,J,K,bi,bj) = 0. D0  
            ENDIF  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 C     Now sync. and get/send edge regions that are shared with  
 C     other threads.  
       _EXCH_XYZ_R4(rHFacC    , myThid )  
       _EXCH_XYZ_R4(rHFacW    , myThid )  
       _EXCH_XYZ_R4(rHFacS    , myThid )  
       _EXCH_XYZ_R4(maskW    , myThid )  
       _EXCH_XYZ_R4(maskS    , myThid )  
       _EXCH_XY_R4 (zA       , myThid )  
   
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO J=1,sNy  
          DO I=1,sNx  
160            tanPhiAtU(I,J,bi,bj) = 0. _d 0            tanPhiAtU(I,J,bi,bj) = 0. _d 0
161            tanPhiAtV(I,J,bi,bj) = 0. _d 0            tanPhiAtV(I,J,bi,bj) = 0. _d 0
162           ENDDO           ENDDO
163          ENDDO          ENDDO
164         ENDDO         ENDDO
165        ENDDO        ENDDO
166          _EXCH_XY_R4 (zA       , myThid )
167        _EXCH_XY_R4 (tanPhiAtU , myThid )        _EXCH_XY_R4 (tanPhiAtU , myThid )
168        _EXCH_XY_R4 (tanPhiAtV , myThid )        _EXCH_XY_R4 (tanPhiAtV , myThid )
169    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22