/[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.20 by adcroft, Thu Sep 27 18:14:52 2001 UTC revision 1.21 by jmc, Thu Dec 23 23:51:22 2004 UTC
# Line 300  C--     Calculate vertical face area (vo Line 300  C--     Calculate vertical face area (vo
300          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
301           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
302  C         by formula  C         by formula
303            lat=yC(I,J,bi,bj)            lat =0.5 _d 0*(yGloc(I,J)+yGloc(I,J+1))
304            dlon=delX( iGl(I,bi) )            dlon=0.5 _d 0*( delX( iGl(I,bi) ) + delX( iGl(I-1,bi) ) )
305            dlat=0.5*( delY( jGl(J,bj) ) + delY( jGl(J-1,bj) ) )            dlat=0.5 _d 0*( delY( jGl(J,bj) ) + delY( jGl(J-1,bj) ) )
306            rAz(I,J,bi,bj) = rSphere*rSphere*dlon*deg2rad            rAz(I,J,bi,bj) = rSphere*rSphere*dlon*deg2rad
307       &     *abs( sin(lat*deg2rad)-sin((lat-dlat)*deg2rad) )       &     *abs( sin(lat*deg2rad)-sin((lat-dlat)*deg2rad) )
308            IF (abs(lat).GT.90..OR.abs(lat-dlat).GT.90.) rAz(I,J,bi,bj)=0.            IF (abs(lat).GT.90..OR.abs(lat-dlat).GT.90.) rAz(I,J,bi,bj)=0.
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22