Index: /issm/trunk/src/m/utils/Alias/cda.m
===================================================================
--- /issm/trunk/src/m/utils/Alias/cda.m	(revision 2890)
+++ /issm/trunk/src/m/utils/Alias/cda.m	(revision 2890)
@@ -0,0 +1,1 @@
+cd /u/astrid-r1b/larour/Testing/Models.devel/Antarctica/Agu09
Index: /issm/trunk/src/m/utils/Alias/cdg.m
===================================================================
--- /issm/trunk/src/m/utils/Alias/cdg.m	(revision 2890)
+++ /issm/trunk/src/m/utils/Alias/cdg.m	(revision 2890)
@@ -0,0 +1,1 @@
+cd /u/astrid-r1b/larour/Testing/Models.devel/Greenland/Icebridge
Index: /issm/trunk/src/m/utils/Basins/plotbasins.m
===================================================================
--- /issm/trunk/src/m/utils/Basins/plotbasins.m	(revision 2890)
+++ /issm/trunk/src/m/utils/Basins/plotbasins.m	(revision 2890)
@@ -0,0 +1,18 @@
+%display all the domain outlines in a directory
+
+basins=listfiles;
+
+hold on
+for i=1:length(basins), 
+	%check whether this is a .exp file
+	basin=basins{i};
+	if strcmpi(basin(end-3:end),'.exp'),
+
+		contour=expread(basin,0);
+		x=contour(1).x;
+		y=contour(1).y;
+		x0=mean(x); y0=mean(y);
+		text(x0,y0,basin(1:end-4),'Fontsize',14);
+		expdisp(basin);
+	end
+end
Index: /issm/trunk/src/m/utils/LatLong/pargenerate.m
===================================================================
--- /issm/trunk/src/m/utils/LatLong/pargenerate.m	(revision 2890)
+++ /issm/trunk/src/m/utils/LatLong/pargenerate.m	(revision 2890)
@@ -0,0 +1,24 @@
+function pargenerate(filename,title,data_format,nlines,ncols,secant_lat,meridian_long,corner_north,corner_east,posting_north,posting_east)
+%PARGENERATE generate parameter file for images ala Rignot
+%
+% Usage: 
+%              pargenerate(filename,title,data_format,nlines,ncols,secant_lat,meridian_long,corner_north,corner_east,posting_north,posting_east)
+%
+%
+
+fid=fopen(filename,'wt');
+
+fprintf(fid,'%s\n','ISSM gridded dataset parameter file');
+fprintf(fid,'%s%s\n','title : ',title);
+fprintf(fid,'%s\n','DEM_projection :    PS');
+fprintf(fid,'%s%s\n','data_format :   ',data_format);
+fprintf(fid,'%s%d\n','width:      ',ncols);
+fprintf(fid,'%s%d\n','nlines:     ',nlines);
+fprintf(fid,'%s%g%s\n','PS_secant_lat:    ',secant_lat,'   decimal degrees');
+fprintf(fid,'%s%g%s\n','PS_meridian_long:    ',meridian_long,'   decimal degrees');
+fprintf(fid,'%s%g%s\n','PS_corner_north:    ',corner_north,'  m');
+fprintf(fid,'%s%g%s\n','PS_corner_east:    ',corner_east,'   m');
+fprintf(fid,'%s%g%s\n','PS_post_north:    ',posting_north,'   m');
+fprintf(fid,'%s%g%s\n','PS_post_east::    ',posting_east,'   m');
+
+fclose(fid);
Index: /issm/trunk/src/m/utils/Mesh/BamgCall.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/BamgCall.m	(revision 2889)
+++ /issm/trunk/src/m/utils/Mesh/BamgCall.m	(revision 2890)
@@ -1,7 +1,7 @@
-function md=YamsCall(md,field,hmin,hmax,gradation,epsilon),
-%YAMSCALL - call yams
+function md=BamgCall(md,field,hmin,hmax,gradation,epsilon),
+%BAMGCALL - call bam
 %
 %   build a metric using the Hessian of the given field
-%   call Yams and the output mesh is plugged onto the model
+%   call Bamg and the output mesh is plugged onto the model
 %   -hmin = minimum edge length (m)
 %   -hmax = maximum edge length (m)
@@ -10,8 +10,8 @@
 %
 %   Usage:
-%      md=YamsCall(md,field,hmin,hmax,gradation,epsilon);
+%      md=BamgCall(md,field,hmin,hmax,gradation,epsilon);
 %
 %   Example:
-%      md=YamsCall(md,md.vel_obs,1500,10^8,1.3,0.9);
+%      md=BamgCall(md,md.vel_obs,1500,10^8,1.3,0.9);
 
 global ISSM_DIR
@@ -63,5 +63,5 @@
 t2=clock;fprintf('%s\n',[' done (' num2str(etime(t2,t1)) ' seconds)']);
 
-%call yams
+%call bamg
 fprintf('%s\n','      call Bamg...');
 system(['bamg -ratio ' num2str(gradation) ' -splitpbedge -nbv 1000000 -M carre0.met -b carre0.mesh -o carre1.mesh']);
Index: /issm/trunk/src/m/utils/Mesh/GetAreas.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/GetAreas.m	(revision 2889)
+++ /issm/trunk/src/m/utils/Mesh/GetAreas.m	(revision 2890)
@@ -44,8 +44,8 @@
 if nargin==3,
 	%compute the surface of the triangle
-	areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)));
+	areas=(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)));
 else
 	%V=area(triangle)*1/3(z1+z2+z3)
 	thickness=mean(z(index(:,4:6)),2)-mean(z(index(:,1:3)),2);
-	areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1))).*thickness;
+	areas=(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1))).*thickness;
 end
Index: /issm/trunk/src/m/utils/Mesh/meshread.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/meshread.m	(revision 2889)
+++ /issm/trunk/src/m/utils/Mesh/meshread.m	(revision 2890)
@@ -29,4 +29,8 @@
 		Struct.index=A(1:3,:)';
 
+	elseif strcmp(A,'Quadrilaterals'),
+		Struct.nels=fscanf(fid,'%i',1);
+		A=fscanf(fid,'%i %i %i %i',[5 Struct.nels]);
+		Struct.index=A(1:4,:)';
 	else
 		%do nothing
Index: /issm/trunk/src/m/utils/Plot/basinzoom.m
===================================================================
--- /issm/trunk/src/m/utils/Plot/basinzoom.m	(revision 2889)
+++ /issm/trunk/src/m/utils/Plot/basinzoom.m	(revision 2890)
@@ -128,4 +128,5 @@
 	regions=AddAvailableRegion(regions,'pig',-1.72*10^6,-1.45*10^6,-3.5*10^5,.5*10^5);
 	regions=AddAvailableRegion(regions,'ronneiceshelf',-1.5*10^6,-.7*10^6,1*10^5,9.3*10^5);
+	regions=AddAvailableRegion(regions,'ris',-1.5*10^6,-.7*10^6,1*10^5,9.3*10^5);
 	regions=AddAvailableRegion(regions,'roosevelt',-6.2797*10^5,-1.5414*10^5,-1.2734*10^6,-.9609*10^6);
 	regions=AddAvailableRegion(regions,'rossiceshelf',-8*10^5,4.5*10^5,-13.8*10^5,-5*10^5);