Changeset 18216

Timestamp:

07/02/14 09:26:48 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: improved help

Location:

issm/trunk-jpl/src/m

Files:

• classes/inversion.m (modified) (1 diff)
• classes/model.m (modified) (1 diff)
• classes/model.py (modified) (1 diff)
• extrusion/project2d.m (modified) (2 diffs)
• extrusion/project2d.py (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/inversion.m

@@ -206 +206 @@
                         WriteData(fid,'object',obj,'fieldname','vy_obs','format','DoubleMat','mattype',1,'scale',1./yts);
                         WriteData(fid,'object',obj,'fieldname','vz_obs','format','DoubleMat','mattype',1,'scale',1./yts);
-                        WriteData(fid,'object',obj,'fieldname','thickness_obs','format','DoubleMat','mattype',1);
+                        if(numel(obj.thickness_obs)==md.mesh.numberofelements),
+                                mattype=2;
+                        else
+                                mattype=1;
+                        end
+                        WriteData(fid,'object',obj,'class','inversion','fieldname','thickness_obs','format','DoubleMat','mattype',mattype);
+
 
                         %process control parameters

issm/trunk-jpl/src/m/classes/model.m

@@ -552 +552 @@
                         %
                         %   Example:
-                        %      md=extrude(md,8,3);
-                        %      md=extrude(md,8,3,2);
+                        %      md=extrude(md,15,1.3);
+                        %      md=extrude(md,15,1.3,1.2);
                         %      md=extrude(md,[0 0.2 0.5 0.7 0.9 0.95 1]);
                         %

issm/trunk-jpl/src/m/classes/model.py

@@ -462 +462 @@
 
                    Example:
-                      md=extrude(md,8,3)
-                      md=extrude(md,8,3,2)
-                      md=extrude(md,[0 0.2 0.5 0.7 0.9 0.95 1])
+                                md=extrude(md,15,1.3);
+                                md=extrude(md,15,1.3,1.2);
+                                md=extrude(md,[0 0.2 0.5 0.7 0.9 0.95 1])
 
                    See also: MODELEXTRACT, COLLAPSE

issm/trunk-jpl/src/m/extrusion/project2d.m

@@ -2 +2 @@
 %PROJECT2D - returns the value of a field for a given layer of the mesh
 %
-%   project 'value' vector taken at layer 'layer' from extruded 2d-3d mesh onto 2d mesh 
-%   used to do the extrusion. This routine is used to compare values between a 2d-3d mesh
-%   at a certain layer, and the equivalent value (if it exists), on the original 2d mesh. 
-%   This routine relies heavily on projections (contained in 3d model md) recored during 
-%   the extrude operation.
+%   returns the value of a vector for a given layer from extruded mesh onto the 2d mesh 
+%   used to do the extrusion. This function is used to compare values between different
+%   layers of a 3d mesh.
 %
 %   Usage:
@@ -12 +10 @@
 %
 %   Example:
-%      vel2=project2d(md3d,md3d.vel,2);
+%      vel2=project2d(md3d,md3d.initialization.vel,2);
+%      returns the velocity of the second layer (1 is the base)
 
 %some checks on list of arguments

issm/trunk-jpl/src/m/extrusion/project2d.py

@@ -5 +5 @@
         returns the value of a field for a given layer of the mesh
         
-        project 'value' vector taken at layer 'layer' from extruded 2d-3d
-        mesh onto 2d mesh used to do the extrusion. This routine is used to
-        compare values between a 2d-3d mesh at a certain layer, and the
-        equivalent value (if it exists), on the original 2d mesh.       
 
-        Note: for consistency with usage in ISSM, layers are indexed starting
-        from one.  This routine accounts for the fact that python indexes from
-        zero.  In other words, when using this function consider the bottom
-        layer as layer=1.  
+   returns the value of a vector for a given layer from extruded mesh onto the 2d mesh 
+   used to do the extrusion. This function is used to compare values between different
+   layers of a 3d mesh.
 
    Usage:
@@ -19 +14 @@
 
    Example:
-      vel2=project2d(md3d,md3d.vel,2)
+      vel2=project2d(md3d,md3d.initialization.vel,2);
+      returns the velocity of the second layer (1 is the base)
         '''