Changeset 400


Ignore:
Timestamp:
05/13/09 16:13:12 (16 years ago)
Author:
seroussi
Message:

2d determinant needed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/classes/@icefront/CreatePVector.m

    r384 r400  
    273273        complete_list=zeros(3,3); %a third coordinate is needed for the jacobian determinant calculation, here it is zero
    274274        complete_list(1:3,1:2)=xyz_tria(3*i-2:3*i,:);
    275         J(i)=GetJacobianDeterminant3d(tria,complete_list,l1l2l3_tria);
     275        J(i)=GetJacobianDeterminant2d(tria,complete_list,l1l2l3_tria);
    276276         end
    277277
     
    295295                water_pressure_tria=0;
    296296        else
    297                 error('QuadPressureLoad error message: unknow fill type for icefront boundary condition');
     297                error('QuadPressureLoadStokes error message: unknow fill type for icefront boundary condition');
    298298        end
    299299
     
    710710                complete_list=zeros(3,3); %a third coordinate is needed for the jacobian determinant calculation, here it is zero
    711711                complete_list(1:3,1:2)=xyz_tria(3*i-2:3*i,:);
    712                 J(i)=GetJacobianDeterminant3d(tria,complete_list,l1l2l3_tria);
     712                J(i)=GetJacobianDeterminant2d(tria,complete_list,l1l2l3_tria);
    713713                 end
    714714
Note: See TracChangeset for help on using the changeset viewer.