Index: /issm/trunk-jpl/src/m/classes/slr.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/slr.py	(revision 23018)
+++ /issm/trunk-jpl/src/m/classes/slr.py	(revision 23019)
@@ -15,26 +15,28 @@
 	
 	def __init__(self): # {{{
-		self.deltathickness    = float('NaN')
-		self.sealevel          = float('NaN')
-		self.maxiter           = 0
-		self.reltol            = 0
-		self.abstol            = 0
-		self.love_h            = 0 #provided by PREM model()
-		self.love_k            = 0 #ideam
-		self.love_l            = 0 #ideam
-		self.tide_love_h       = 0
-		self.tide_love_k       = 0
-		self.fluid_love        = 0; 
-		self.equatorial_moi    = 0; 
-		self.polar_moi	       = 0; 
-		self.angular_velocity  = 0;
-		self.rigid             = 0
-		self.elastic           = 0
-		self.rotation          = 0
-		self.ocean_area_scaling = 0;
-		steric_rate					= 0; #rate of ocean expansion from steric effects.
-		self.degacc            = 0
-		self.requested_outputs = []
-		self.transitions       = []
+		self.deltathickness         = float('NaN')
+		self.sealevel               = float('NaN')
+		self.maxiter                = 0
+		self.reltol                 = 0
+		self.abstol                 = 0
+		self.love_h                 = 0 #provided by PREM model()
+		self.love_k                 = 0 #ideam
+		self.love_l                 = 0 #ideam
+		self.tide_love_h            = 0
+		self.tide_love_k            = 0
+		self.fluid_love             = 0; 
+		self.equatorial_moi         = 0; 
+		self.polar_moi	             = 0; 
+		self.angular_velocity       = 0;
+		self.rigid                  = 0
+		self.elastic                = 0
+		self.rotation               = 0
+		self.ocean_area_scaling     = 0;
+		steric_rate					    = 0; #rate of ocean expansion from steric effects.
+		self.geodetic_run_frequency = 1; #how many time steps we skip before we run the geodetic part of the solver during transient
+		self.geodetic               = 0; #compute geodetic SLR? (in addition to steric?)
+		self.degacc                 = 0
+		self.requested_outputs      = []
+		self.transitions            = []
 		
 		#set defaults
@@ -78,4 +80,5 @@
 
 		#computational flags: 
+		self.geodetic=0
 		self.rigid=1
 		self.elastic=1
@@ -87,10 +90,10 @@
 		self.tide_love_k=0.3055; #degree 2
 		
-                #secular fluid love number: 
+      #secular fluid love number: 
 		self.fluid_love=0.942; 
 		
 		#moment of inertia: 
 		self.equatorial_moi=8.0077*10**37; # [kg m^2] 
-		self.polar_moi	   =8.0345*10**37; # [kg m^2] 
+		self.polar_moi	    =8.0345*10**37; # [kg m^2] 
 		
 		#mean rotational velocity of earth 
@@ -102,4 +105,7 @@
 		#steric:
 		self.steric_rate=0;
+
+		#how many time steps we skip before we run SLR solver during transient
+		self.geodetic_run_frequency=1;
 		
 		#output default:
@@ -133,4 +139,5 @@
 		md = checkfield(md,'fieldname','slr.abstol','size',[1,1])
 		md = checkfield(md,'fieldname','slr.maxiter','size',[1,1],'>=',1)
+		md = checkfield(md,'fieldname','slr.geodetic_run_frequency','size',[1,1],'>=',1)
 		md = checkfield(md,'fieldname','slr.steric_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices]);
 		md = checkfield(md,'fieldname','slr.degacc','size',[1,1],'>=',1e-10)
@@ -164,7 +171,9 @@
 		WriteData(fid,prefix,'object',self,'fieldname','rotation','format','Boolean')
 		WriteData(fid,prefix,'object',self,'fieldname','ocean_area_scaling','format','Boolean')
+		WriteData(fid,prefix,'object',self,'fieldname','geodetic_run_frequency','format','Integer')
 		WriteData(fid,prefix,'object',self,'fieldname','steric_rate','format','DoubleMat','mattype',1,'scale',1e-3/md.constants.yts);
 		WriteData(fid,prefix,'object',self,'fieldname','degacc','format','Double')
 		WriteData(fid,prefix,'object',self,'fieldname','transitions','format','MatArray')
+		WriteData(fid,prefix,'object',self,'fieldname','geodetic','format','Integer')
 	
 		#process requested outputs