Changeset 2890

Timestamp:

01/22/10 07:24:45 (15 years ago)

Author:

Eric.Larour

Message:

Added aliases, similar to bash aliases.
meshread: can not read quads.
BamgCall: cosmetic.
pargenerate: generate idl style parameter file for processed images (Rignot)
plotbasins: replaced basinzoom by basin.

Location:

issm/trunk/src/m/utils

Files:

• Alias (added)
• Alias/cda.m (added)
• Alias/cdg.m (added)
• Basins (added)
• Basins/plotbasins.m (added)
• LatLong/pargenerate.m (added)
• Mesh/BamgCall.m (modified) (3 diffs)
• Mesh/GetAreas.m (modified) (1 diff)
• Mesh/meshread.m (modified) (1 diff)
• Plot/basinzoom.m (modified) (1 diff)

issm/trunk/src/m/utils/Mesh/BamgCall.m

@@ -1 @@
-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 +10 @@
 %
 %   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 +63 @@
 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']);

issm/trunk/src/m/utils/Mesh/GetAreas.m

@@ -44 +44 @@
 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

issm/trunk/src/m/utils/Mesh/meshread.m

@@ -29 +29 @@
                 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

issm/trunk/src/m/utils/Plot/basinzoom.m

@@ -128 +128 @@
         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);