--- MITgcm/model/src/ini_spherical_polar_grid.F	1998/11/30 23:45:25	1.11
+++ MITgcm/model/src/ini_spherical_polar_grid.F	1999/03/16 16:37:42	1.13
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_spherical_polar_grid.F,v 1.11 1998/11/30 23:45:25 adcroft Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_spherical_polar_grid.F,v 1.13 1999/03/16 16:37:42 adcroft Exp $
 
 #include "CPP_OPTIONS.h"
 
@@ -33,6 +33,7 @@
 C     | in X and Y are in degrees. Distance in Z are in m or Pa  |
 C     | depending on the vertical gridding mode.                 |
 C     \==========================================================/
+      IMPLICIT NONE
 
 C     === Global variables ===
 #include "SIZE.h"
@@ -171,6 +172,7 @@
         DO J=1,sNy
          DO I=1,sNx
           jG = myYGlobalLo + (bj-1)*sNy + J-1
+          iG = myXGlobalLo + (bi-1)*sNx + I-1
           latS = yc(i,j,bi,bj)-delY(jG)*0.5 _d 0
           latN = yc(i,j,bi,bj)+delY(jG)*0.5 _d 0
 #ifdef OLD_UV_GEOMETRY