--- MITgcm_contrib/submesoscale/code/ini_grid.F	2008/05/31 00:47:06	1.2
+++ 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.2 2008/05/31 00:47:06 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 )