Index: /issm/trunk-jpl/src/m/classes/solidearthsettings.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/solidearthsettings.py	(revision 25765)
+++ /issm/trunk-jpl/src/m/classes/solidearthsettings.py	(revision 25766)
@@ -23,4 +23,5 @@
         self.runfrequency           = 1 # How many time steps we skip before we run grd_core
         self.computesealevelchange  = 0 # Will grd_core compute sea level?
+        self.isgrd                  = 1 # Will GRD patterns be computed?
         self.degacc                 = 0 # Degree increment for resolution of Green tables
         self.horiz                  = 0 # Compute horizontal displacement?
@@ -42,4 +43,5 @@
         s += '{}\n'.format(fielddisplay(self, 'ocean_area_scaling', 'correction for model representation of ocean area [default: No correction]'))
         s += '{}\n'.format(fielddisplay(self, 'computesealevelchange', 'compute sealevel change from GRD in addition to steric?) default 0'))
+        s += '{}\n'.format(fielddisplay(self, 'isgrd', 'compute GRD patterns (default: 1'))
         s += '{}\n'.format(fielddisplay(self, 'runfrequency', 'how many time steps we skip before we run solidearthsettings solver during transient (default: 1)'))
         s += '{}\n'.format(fielddisplay(self, 'rigid', 'rigid earth graviational potential perturbation'))
@@ -63,4 +65,5 @@
         self.rotation = 1
         self.ocean_area_scaling = 0
+        self.isgrd = 1
         self.computesealevelchange = 0
 
@@ -90,5 +93,5 @@
 
         # A coupler to planet model is provided
-        if self.computesealevelchange:
+        if self.isgrd:
             if md.transient.iscoupler:
                 # We are good
@@ -115,4 +118,5 @@
         WriteData(fid, prefix, 'object', self, 'fieldname', 'horiz', 'name', 'md.solidearth.settings.horiz', 'format', 'Integer')
         WriteData(fid, prefix, 'object', self, 'fieldname', 'computesealevelchange', 'name', 'md.solidearth.settings.computesealevelchange', 'format', 'Integer')
+        WriteData(fid, prefix, 'object', self, 'fieldname','isgrd', 'name', 'md.solidearth.settings.isgrd', 'format', 'Integer')
         WriteData(fid, prefix, 'object', self, 'fieldname','glfraction', 'name', 'md.solidearth.settings.glfraction', 'format', 'Integer')
     #}}}
