Index: /issm/trunk-jpl/src/m/classes/esa.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/esa.js	(revision 22351)
+++ /issm/trunk-jpl/src/m/classes/esa.js	(revision 22352)
@@ -10,5 +10,8 @@
 		//numerical discretization accuracy
 		this.degacc=.01;
-		
+	
+		//computational flags:
+		this.hemisphere=0;
+
 		//output default:
 		this.requested_outputs=['default'];
@@ -26,4 +29,5 @@
 			md = checkfield(md,'fieldname','esa.love_h','NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','esa.love_l','NaN',1,'Inf',1);
+			md = checkfield(md,'fieldname','esa.hemisphere','NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','esa.degacc','size',[1, 1],'>=',1e-10);
 			md = checkfield(md,'fieldname','esa.requested_outputs','stringrow',1);
@@ -48,4 +52,5 @@
 		fielddisplay(this,'love_h','load Love number for radial displacement');
 		fielddisplay(this,'love_l','load Love number for horizontal displacements'); 
+		fielddisplay(this,'hemisphere','North-south, East-west components of 2-D horiz displacement vector: -1 south, 1 north'); 
 		fielddisplay(this,'degacc',"accuracy (default .01 deg) for numerical discretization of the Green's functions");
 		fielddisplay(this,'transitions','indices into parts of the mesh that will be icecaps');
@@ -57,4 +62,5 @@
 			WriteData(fid,prefix,'object',this,'fieldname','love_h','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',this,'fieldname','love_l','format','DoubleMat','mattype',1);
+			WriteData(fid,prefix,'object',this,'fieldname','hemisphere','format','Integer');
 			WriteData(fid,prefix,'object',this,'fieldname','degacc','format','Double');
 			WriteData(fid,prefix,'object',this,'fieldname','transitions','format','MatArray');
@@ -75,4 +81,5 @@
 			this.love_h=NullFix(this.love_h,NaN);
 			this.love_l=NullFix(this.love_l,NaN);
+			this.hemisphere=NullFix(this.hemisphere,NaN);
 			this.degacc=NullFix(this.degacc,NaN);
 		}//}}}
@@ -82,4 +89,5 @@
 	this.love_h         = 0; //provided by PREM model
 	this.love_l         = 0; //idam
+	this.hemisphere     = 0;
 	this.degacc         = 0;
 	this.requested_outputs = [];
Index: /issm/trunk-jpl/src/m/classes/esa.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/esa.m	(revision 22351)
+++ /issm/trunk-jpl/src/m/classes/esa.m	(revision 22352)
@@ -44,4 +44,5 @@
 			md = checkfield(md,'fieldname','esa.love_h','NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','esa.love_l','NaN',1,'Inf',1);
+			md = checkfield(md,'fieldname','esa.hemisphere','NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','esa.degacc','size',[1 1],'>=',1e-10);
 			md = checkfield(md,'fieldname','esa.requested_outputs','stringrow',1);
Index: /issm/trunk-jpl/src/m/classes/esa.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/esa.py	(revision 22351)
+++ /issm/trunk-jpl/src/m/classes/esa.py	(revision 22352)
@@ -18,4 +18,5 @@
 		self.love_h            = 0 #provided by PREM model()
 		self.love_l            = 0 #ideam
+		self.hemisphere        = 0
 		self.degacc            = 0
 		self.requested_outputs = []
@@ -30,4 +31,5 @@
 			string="%s\n%s"%(string,fielddisplay(self,'love_h','load Love number for radial displacement'))
 			string="%s\n%s"%(string,fielddisplay(self,'love_l','load Love number for horizontal displaements'))
+                        string="%s\n%s"%(string,fielddisplay(self,'hemisphere','North-south, East-west components of 2-D horiz displacement vector: -1 south, 1 north'))
 			string="%s\n%s"%(string,fielddisplay(self,'degacc','accuracy (default .01 deg) for numerical discretization of the Green''s functions'))
 			string="%s\n%s"%(string,fielddisplay(self,'transitions','indices into parts of the mesh that will be icecaps'))
@@ -40,5 +42,8 @@
 		#numerical discretization accuracy
 		self.degacc=.01
-		
+	
+                #computational flags:
+                self.hemisphere=0;
+
 		#output default:
 		self.requested_outputs=['default']
@@ -60,4 +65,5 @@
 		md = checkfield(md,'fieldname','esa.love_h','NaN',1,'Inf',1)
 		md = checkfield(md,'fieldname','esa.love_l','NaN',1,'Inf',1)
+		md = checkfield(md,'fieldname','esa.hemisphere','NaN',1,'Inf',1)
 		md = checkfield(md,'fieldname','esa.degacc','size',[1,1],'>=',1e-10)
 		md = checkfield(md,'fieldname','esa.requested_outputs','stringrow',1)
@@ -75,4 +81,5 @@
 		WriteData(fid,prefix,'object',self,'fieldname','love_h','format','DoubleMat','mattype',1)
 		WriteData(fid,prefix,'object',self,'fieldname','love_l','format','DoubleMat','mattype',1)
+		WriteData(fid,prefix,'object',self,'fieldname','hemisphere','format','Integer')
 		WriteData(fid,prefix,'object',self,'fieldname','degacc','format','Double')
 		WriteData(fid,prefix,'object',self,'fieldname','transitions','format','MatArray')