Index: /issm/trunk-jpl/test/MITgcm/coupling_rules.txt
===================================================================
--- /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 22705)
+++ /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 22706)
@@ -14,45 +14,31 @@
 
 >>>> INITIALIZATION FROM SCRATCH
-
-Initialization:
 IceModelTime   = 0
 OceanModelTime = 0
 
-  I1. ISSM sends CouplingTime
-      (the interval at which we couple)
-
+Initialization:
+  I1. ISSM sends CouplingTime, the interval at which we couple
   I2. MITgcm sends grid size (NX and NY)
-
   I3. MITgcm sends grid coordinates of center of cells
       (longitude -180 <= XC < 180 and latitude YC)
-
   I4. ISSM sends ice model time IceTimeTag
-
   I5. MITgcm sends ocean model time OceanTimeTag
-
   I6. ISSM sends IceModelTime=0 base to MITgcm
       (base should be 9999 in areas not covered by ice grid)
+  I7o. MITgcm integrates from OceanModelTime=0
+                           to OceanModelTime=CouplingTime
 
-  I7. MITgcm integrates from OceanModelTime=0
-                          to OceanModelTime=CouplingTime
-
-First Coupling Exchange:
+First Recurrent Coupling Exchange:
 IceModelTime   = 0
 OceanModelTime = CouplingTime
-
   C1. ISSM sends ice model time IceTimeTag
-
   C2. MITgcm sends ocean model time OceanTimeTag
-
   C3. MITgcm sends
       0 <= OceanModelTime < CouplingTime
       time-mean melt rate to ISSM
-  
   C4. ISSM sends IceModelTime=0 base to MITgcm
       (this is duplicate of Initialization)
-
   C5i. ISSM integrates from IceModelTime=0
                          to IceModelTime=CouplingTime
-
   C5o. MITgcm integrates from OceanModelTime=CouplingTime
                            to OceanModelTime=2*CouplingTime
@@ -63,18 +49,12 @@
 IceModelTime   = (N-1) * CouplingTime
 OceanModelTime =   N   * CouplingTime
-
   C1. ISSM sends ice model time IceTimeTag
-
   C2. MITgcm sends ocean model time OceanTimeTag
-
   C3. MITgcm sends
       (N-1)*CouplingTime <= OceanModelTime < N*CouplingTime
       time-mean melt rate to ISSM
-  
   C4. ISSM sends IceModelTime=(N-1)*CouplingTime base to MITgcm
-
   C5i. ISSM integrates from IceModelTime=(N-1)*CouplingTime
                          to IceModelTime=N*CouplingTime
-
   C5o. MITgcm integrates from OceanModelTime=N*CouplingTime
                            to OceanModelTime=(N+1)*CouplingTime
@@ -83,16 +63,11 @@
 IceModelTime   =   N   * CouplingTime
 OceanModelTime = (N+1) * CouplingTime
-
   C1. ISSM sends ice model time IceTimeTag
-
   C2. MITgcm sends ocean model time OceanTimeTag
-
   C3. MITgcm sends
       N*CouplingTime <= OceanModelTime < (N+1)*CouplingTime
       time-mean melt rate to ISSM
-  
   C4i. ISSM saves to restart files the IceModelTime=N state and
        the N*CouplingTime to (N+1)*CouplingTime mean melt rate
-
   C4o. MITgcm saves to restart files
        the OceanModelTime=(N+1)*CouplingTime state
@@ -100,50 +75,36 @@
 
 >>>> RESTART FROM SAVED STATE
-
-Initialization:
 IceModelTime   =   N   * CouplingTime
 OceanModelTime = (N+1) * CouplingTime
 
+  - ISSM reads from restart files the IceModelTime=N state and
+      the N*CouplingTime to (N+1)*CouplingTime mean melt rate
+      
+  - MITgcm reads from restart files
+      the OceanModelTime=(N+1)*CouplingTime state
+
+Initialization:
   I1. ISSM sends CouplingTime
-
   I2. MITgcm sends grid size (NX and NY)
-
   I3. MITgcm sends grid coordinates (XC and YC)
-
-  I4i. ISSM reads from restart files the IceModelTime=N state and
-       the N*CouplingTime to (N+1)*CouplingTime mean melt rate
-
-  I4o. MITgcm reads from restart files
-       the OceanModelTime=(N+1)*CouplingTime state
-
-  I5. ISSM sends ice model time IceTimeTag
-
-  I6. MITgcm sends ocean model time OceanTimeTag
-
-  I7. ISSM sends IceModelTime=N base to MITgcm
-
-  I8i. ISSM integrates from IceModelTime=N*CouplingTime
+  I4. ISSM sends ice model time IceTimeTag
+  I5. MITgcm sends ocean model time OceanTimeTag
+  I6. ISSM sends IceModelTime=N base to MITgcm
+  I7i. ISSM integrates from IceModelTime=N*CouplingTime
                          to IceModelTime=(N+1)*CouplingTime
-
-  I8o. MITgcm integrates from OceanModelTime=(N+1)*CouplingTime
+  I7o. MITgcm integrates from OceanModelTime=(N+1)*CouplingTime
                           to OceanModelTime=(N+2)*CouplingTime
 
-First Coupling Exchange:
+First Recurrent Coupling Exchange:
 IceModelTime   = (N+1)   * CouplingTime
 OceanModelTime = (N+2) * CouplingTime
-
   C1. ISSM sends ice model time IceTimeTag
-
   C2. MITgcm sends ocean model time OceanTimeTag
-
   C3. MITgcm sends
       (N+1)*CouplingTime <= OceanModelTime < (N+2)*CouplingTime
       time-mean melt rate to ISSM
-  
   C4. ISSM sends IceModelTime=(N+1)*CouplingTime base to MITgcm
-
   C5i. ISSM integrates from IceModelTime=(N+1)*CouplingTime
                          to IceModelTime=(N+2)*CouplingTime
-
   C5o. MITgcm integrates from OceanModelTime=(N+2)*CouplingTime
                            to OceanModelTime=(N+3)*CouplingTime