/[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.9 by cnh, Wed Oct 28 03:11:37 1998 UTC revision 1.10 by cnh, Fri Nov 6 22:44:47 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_SPHERICAL_POLAR_GRID( myThid )        SUBROUTINE INI_SPHERICAL_POLAR_GRID( myThid )
# Line 102  C           longitude - phiMin and theta Line 102  C           longitude - phiMin and theta
102            xc(I,J,bi,bj)  = xG + delX(iG+i-1)*0.5 _d 0            xc(I,J,bi,bj)  = xG + delX(iG+i-1)*0.5 _d 0
103            yc(I,J,bi,bj)  = yG + delY(jG+j-1)*0.5 _d 0            yc(I,J,bi,bj)  = yG + delY(jG+j-1)*0.5 _d 0
104            xG = xG + delX(iG+I-1)            xG = xG + delX(iG+I-1)
105            dxF(I,J,bi,bj) = delX(iG+i-1)*deg2rad*            dxF(I,J,bi,bj) = delX(iG+i-1)*deg2rad
106       &     rSphere*COS(yc(I,J,bi,bj)*deg2rad)       &    *rSphere*COS(yc(I,J,bi,bj)*deg2rad)
107            dyF(I,J,bi,bj) = delY(jG+j-1)*deg2rad*rSphere            dyF(I,J,bi,bj) = delY(jG+j-1)*deg2rad*rSphere
108           ENDDO           ENDDO
109           yG = yG + delY(jG+J-1)           yG = yG + delY(jG+J-1)
# Line 171  C     Calculate vertical face area and t Line 171  C     Calculate vertical face area and t
171            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0            latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0
172            rA(I,J,bi,bj) = dyF(I,J,bi,bj)            rA(I,J,bi,bj) = dyF(I,J,bi,bj)
173       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))       &    *rSphere*(SIN(latN*deg2rad)-SIN(latS*deg2rad))
174    C         Area cannot be zero but sin can be if lat if < -90.
175              IF ( rA(I,J,bi,bj) .LT. 0. ) rA(I,J,bi,bj) = -rA(I,J,bi,bj)
176            tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad)            tanPhiAtU(i,j,bi,bj)=tan(_yC(i,j,bi,bj)*deg2rad)
177            tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad)            tanPhiAtV(i,j,bi,bj)=tan(latS*deg2rad)
178           ENDDO           ENDDO
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22