Changeset 384


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

now have 2d and 3d determinant for tria in ice package

Location:
issm/trunk/src/m/classes
Files:
2 added
1 deleted
9 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/m/classes/@icefront/CreatePVector.m

    r1 r384  
    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)=GetJacobianDeterminant(tria,complete_list,l1l2l3_tria);
     275        J(i)=GetJacobianDeterminant3d(tria,complete_list,l1l2l3_tria);
    276276         end
    277277
    … …  
    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)=GetJacobianDeterminant(tria,complete_list,l1l2l3_tria);
     712                J(i)=GetJacobianDeterminant3d(tria,complete_list,l1l2l3_tria);
    713713                 end
    714714

  • issm/trunk/src/m/classes/@pentaelem/CreateKMatrix.m

    r328 r384  
    9999               
    100100                %Get the Jacobian determinant
    101                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     101                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list,gauss_coord);
    102102
    103103                %Get L matrix if viscous basal drag present:
    … …  
    303303                %Get the Jacobian determinant
    304304                xyz_list_tria=xyz_list(1:3,:);
    305                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord_tria);
     305                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list_tria,gauss_coord_tria);
    306306
    307307                %Get L and Lprime matrix if viscous basal drag present:
    … …  
    437437                       
    438438                        %Get the Jacobian determinant
    439                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     439                        Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list_tria,gauss_coord);
    440440
    441441                        %Get L matrix if viscous basal drag present:
    … …  
    597597                        %Get the Jacobian determinant
    598598                        xyz_list_tria=xyz_list(1:3,:);
    599                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     599                        Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list_tria,gauss_coord);
    600600                       
    601601                        %Get nodal functions value
    … …  
    859859                        %Get the Jacobian determinant
    860860                        xyz_list_tria=xyz_list(1:3,:);
    861                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     861                        Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list_tria,gauss_coord);
    862862
    863863                        %Get L if viscous basal drag present:
    … …  
    931931        if pentaelem.artificial_diffusivity,
    932932                %Get the Jacobian determinant
    933                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,[1/3 1/3 1/3]);
     933                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,[1/3 1/3 1/3]);
    934934
    935935                %Build K matrix (artificial diffusivity matrix)
    … …  
    948948
    949949                %Get the Jacobian determinant
    950                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     950                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    951951
    952952                %Get L matrix if viscous basal drag present:
    … …  
    10441044               
    10451045                %Get the Jacobian determinant
    1046                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     1046                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list,gauss_coord);
    10471047
    10481048                %Get L matrix if viscous basal drag present:

  • issm/trunk/src/m/classes/@pentaelem/CreatePVector.m

    r328 r384  
    258258       
    259259                %Get the Jacobian determinant
    260                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     260                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list_tria,gauss_coord);
    261261
    262262                %Get bed at gaussian point
    … …  
    366366
    367367                %Get the Jacobian determinant
    368                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     368                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list,gauss_coord);
    369369
    370370                %Get L:
    … …  
    515515
    516516                        %Get the Jacobian determinant
    517                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     517                        Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list,gauss_coord);
    518518
    519519                        %Get L:
    … …  
    561561
    562562                                %Get the Jacobian determinant
    563                                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     563                                Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list,gauss_coord);
    564564
    565565                                %Build gaussian vector
    … …  
    741741                        %Get the Jacobian determinant
    742742                        xyz_list_tria=xyz_list(1:3,:);
    743                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     743                        Jdettria=GetJacobianDeterminant3d(triaelem,xyz_list_tria,gauss_coord);
    744744                       
    745745                        %Get nodal functions value
    … …  
    818818                        %Get the Jacobian determinant
    819819                        xyz_list_tria=xyz_list(1:3,:);
    820                         Jdettria=GetJacobianDeterminant(triaelem,xyz_list_tria,gauss_coord);
     820                        Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list_tria,gauss_coord);
    821821                       
    822822                        %Get geothermal flux and basal friction
    … …  
    11601160
    11611161                %Get the Jacobian determinant
    1162                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     1162                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    11631163
    11641164                %Get L:

  • issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m

    r229 r384  
    7474        if triaelem.artificial_diffusivity,
    7575                %Get the Jacobian determinant
    76                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,[1/3 1/3 1/3]);
     76                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,[1/3 1/3 1/3]);
    7777
    7878                %Build K matrix (artificial diffusivity matrix)
    … …  
    9191
    9292                %Get the Jacobian determinant
    93                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     93                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    9494
    9595                %Get L matrix if viscous basal drag present:
    … …  
    178178               
    179179                %Get the Jacobian determinant
    180                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     180                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    181181
    182182                %Get L matrix if viscous basal drag present:
    … …  
    400400
    401401        %Get Jacobian determinant:
    402         Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_l1l2l3);
     402        Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_l1l2l3);
    403403
    404404        % Build the D matrix: we plug the gaussian weight, the thickness, the viscosity, and the jacobian determinant

  • issm/trunk/src/m/classes/@triaelem/CreatePVector.m

    r45 r384  
    5454        xyz_list=getgriddata(triaelem,grids);
    5555
    56         %Just keep the first 3 grids
    57         xyz_list=xyz_list(1:3,:);
    58 
    5956        %recover extra inputs
    6057        if strcmpi(analysis_type,'bed_slope_compute_x') | strcmpi(analysis_type,'bed_slope_compute_y') ,
    … …  
    10097
    10198                %Get the Jacobian determinant
    102                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     99                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    103100
    104101                %Get L:
    … …  
    146143        %Get all element grid data:
    147144        xyz_list=getgriddata(triaelem,grids);
     145
     146        %It is a 2d model so get rid of z coordinates even in 3d collapse
     147        xyz_list(:,3)=0;
    148148
    149149        %initialize extra inputs
    … …  
    219219
    220220                %Get Jacobian determinant:
    221                 Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_l1l2l3);
     221                Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_l1l2l3);
    222222
    223223                %Get nodal functions
    … …  
    331331
    332332                %Get the Jacobian determinant
    333                 Jdettria=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     333                Jdettria=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    334334
    335335                %Get L:

  • issm/trunk/src/m/classes/@triaelem/DeviatoricStress.m

    r1 r384  
    6767
    6868                %Get Jacobian determinant:
    69                 Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     69                Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    7070
    7171                deviatoricstress=deviatoricstress+deviatoricstress_g*Jdet*gauss_weight;

  • issm/trunk/src/m/classes/@triaelem/StrainRate.m

    r1 r384  
    5454
    5555                %Get Jacobian determinant:
    56                 Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     56                Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    5757
    5858                strainrate=strainrate+epsilon*Jdet*gauss_weight;

  • issm/trunk/src/m/classes/@triaelem/Stress.m

    r1 r384  
    7474
    7575                %Get Jacobian determinant:
    76                 Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     76                Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    7777
    7878                stress=stress+stress_g*Jdet*gauss_weight;

  • issm/trunk/src/m/classes/@triaelem/ViscousHeating.m

    r1 r384  
    7979
    8080                %Get Jacobian determinant:
    81                 Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord);
     81                Jdet=GetJacobianDeterminant2d(triaelem,xyz_list,gauss_coord);
    8282
    8383                viscousheating=viscousheating+phi*Jdet*gauss_weight;
Note: See TracChangeset for help on using the changeset viewer.