/[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.5 by cnh, Mon Jun 8 21:43:01 1998 UTC revision 1.6 by adcroft, Mon Jun 22 15:26:25 1998 UTC
# Line 185  C     Calculate recipricols Line 185  C     Calculate recipricols
185        _EXCH_XY_R4(rDyF, myThid )        _EXCH_XY_R4(rDyF, myThid )
186        _EXCH_XY_R4(rDxV, myThid )        _EXCH_XY_R4(rDxV, myThid )
187        _EXCH_XY_R4(rDyU, myThid )        _EXCH_XY_R4(rDyU, myThid )
188  C     Calculate vertical face area  C     Calculate vertical face area and trigonometric terms
189        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
190         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
191          DO J=1,sNy          DO J=1,sNy
# Line 195  C     Calculate vertical face area Line 195  C     Calculate vertical face area
195            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0
196            zA(I,J,bi,bj) = dyF(I,J,bi,bj)            zA(I,J,bi,bj) = dyF(I,J,bi,bj)
197       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))
198              tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad)
199              tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad)
200           ENDDO           ENDDO
201          ENDDO          ENDDO
202         ENDDO         ENDDO
# Line 237  C     other threads. Line 239  C     other threads.
239        _EXCH_XYZ_R4(maskW    , myThid )        _EXCH_XYZ_R4(maskW    , myThid )
240        _EXCH_XYZ_R4(maskS    , myThid )        _EXCH_XYZ_R4(maskS    , myThid )
241        _EXCH_XY_R4 (zA       , myThid )        _EXCH_XY_R4 (zA       , myThid )
   
 CcnhDebugStarts  
       tanPhiAtU = 0. _d 0  
       tanPhiAtV = 0. _d 0  
242        _EXCH_XY_R4 (tanPhiAtU , myThid )        _EXCH_XY_R4 (tanPhiAtU , myThid )
243        _EXCH_XY_R4 (tanPhiAtV , myThid )        _EXCH_XY_R4 (tanPhiAtV , myThid )
 CcnhDebugEnds  
   
244    
245  C  C
246        RETURN        RETURN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22