--- MITgcm_contrib/submesoscale/code/ini_grid.F	2008/05/31 00:13:54	1.1
+++ MITgcm_contrib/submesoscale/code/ini_grid.F	2010/03/19 19:24:16	1.3
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/submesoscale/code/ini_grid.F,v 1.1 2008/05/31 00:13:54 dimitri Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/submesoscale/code/ini_grid.F,v 1.3 2010/03/19 19:24:16 jmc Exp $
 C $Name:  $
 
 #include "PACKAGES_CONFIG.h"
@@ -15,6 +15,15 @@
 C     by each thread and initializes only the region of the domain it is
 C     "responsible" for.
 
+C     !CALLING SEQUENCE:
+C     INI_GRID
+C      |   -- LOAD_GRID_SPACING
+C      |   -- INI_VERTICAL_GRID
+C      |    / INI_CARTESIAN_GRID
+C      |   /  INI_SPHERICAL_POLAR_GRID
+C      |   \  INI_CURVILINEAR_GRID
+C      |    \ INI_CYLINDER_GRID
+
 C     !USES:
       IMPLICIT NONE
 #include "SIZE.h"
@@ -47,11 +56,6 @@
 C     load grid spacing (vector) from files
       CALL LOAD_GRID_SPACING( myThid )
 
-C     Set up reference vertical profile (vector) for state variables
-C jmc: this call will not stay here but will move to S/R INITIALISE_FIXED
-C      (a better place) once anelatic initialisation is cleaned-up.
-      CALL INI_REFERENCE_STATE( myThid )
-
 C     Set up vertical grid and coordinate system
       CALL INI_VERTICAL_GRID( myThid )
 
@@ -89,6 +93,13 @@
         _END_MASTER(myThid)
       ENDIF
 
+#ifdef ALLOW_OBCS
+      IF ( useOBCS ) THEN
+C     Fix overlap regions to avoid discontinuities in dxc, dyc, etc.
+       CALL OBCS_CHECK_GRID( myThid )
+      ENDIF
+#endif /* ALLOW_OBCS */
+
 #ifdef ALLOW_MONITOR
       IF ( MASTER_CPU_IO(myThid) ) THEN
 C--   only the master thread is allowed to switch On/Off mon_write_stdout