Changeset 630


Ignore:
Timestamp:
05/28/09 14:59:46 (16 years ago)
Author:
Mathieu Morlighem
Message:

cosmetics

Location:
issm/trunk/src/m/solutions/ice
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/ice/StrainRateCompute.m

    r51 r630  
    3232
    3333                        strainratematrix=[strainratevector(1) strainratevector(3)
    34                                                                         strainratevector(3)  strainratevector(2)];
     34                        strainratevector(3)  strainratevector(2)];
    3535
    3636                        %eigen values and vectors
     
    7272
    7373                        strainratematrix=[strainratevector(1) strainratevector(4) strainratevector(5)
    74                                       strainratevector(4)  strainratevector(2)  strainratevector(6)
    75                                       strainratevector(5)  strainratevector(6)  strainratevector(3)];
     74                        strainratevector(4)  strainratevector(2)  strainratevector(6)
     75                        strainratevector(5)  strainratevector(6)  strainratevector(3)];
    7676
    7777                        %eigen values and vectors
    7878                        [directions,value]=eig(strainratematrix);
    7979
    80                         %Plug into global vectors
     80                        %Plug into global vectors
    8181                        strainrate1(n,:)=strainratevector;
    82                         A1(n,1)=value(1,1); A2(n,1)=value(2,2); A3(n,1)=value(3,3);
    83                         Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2); Vx3(n,1)=directions(1,3);
    84                         Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2); Vy3(n,1)=directions(2,3);
    85                         Vz1(n,1)=directions(3,1); Vz2(n,1)=directions(3,2); Vz3(n,1)=directions(3,3);
     82                        A1(n,1)=value(1,1); A2(n,1)=value(2,2); A3(n,1)=value(3,3);
     83                        Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2); Vx3(n,1)=directions(1,3);
     84                        Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2); Vy3(n,1)=directions(2,3);
     85                        Vz1(n,1)=directions(3,1); Vz2(n,1)=directions(3,2); Vz3(n,1)=directions(3,3);
    8686                end
    8787        end
  • issm/trunk/src/m/solutions/ice/StressCompute.m

    r51 r630  
    3333
    3434                        stressmatrix=[stressvector(1) stressvector(3)
    35                                       stressvector(3)  stressvector(2)];
     35                        stressvector(3)  stressvector(2)];
    3636
    3737                        %eigen values and vectors
    3838                        [directions,value]=eig(stressmatrix);
    3939
    40                         %Plug into global vectors
     40                        %Plug into global vectors
    4141                        stress1(n,:)=stressvector;
    42                         A1(n,1)=value(1,1); A2(n,1)=value(2,2);
    43                         Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2);
    44                         Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2);
     42                        A1(n,1)=value(1,1); A2(n,1)=value(2,2);
     43                        Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2);
     44                        Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2);
    4545                end
    4646        end
     
    7474
    7575                        stressmatrix=[stressvector(1) stressvector(4) stressvector(5)
    76                                       stressvector(4)  stressvector(2)  stressvector(6)
    77                                       stressvector(5)  stressvector(6)  stressvector(3)];
     76                        stressvector(4)  stressvector(2)  stressvector(6)
     77                        stressvector(5)  stressvector(6)  stressvector(3)];
    7878
    7979                        %eigen values and vectors
    8080                        [directions,value]=eig(stressmatrix);
    8181
    82                         %Plug into global vectors
     82                        %Plug into global vectors
    8383                        stress1(n,:)=stressvector;
    84                         A1(n,1)=value(1,1); A2(n,1)=value(2,2); A3(n,1)=value(3,3);
    85                         Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2); Vx3(n,1)=directions(1,3);
    86                         Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2); Vy3(n,1)=directions(2,3);
    87                         Vz1(n,1)=directions(3,1); Vz2(n,1)=directions(3,2); Vz3(n,1)=directions(3,3);
     84                        A1(n,1)=value(1,1); A2(n,1)=value(2,2); A3(n,1)=value(3,3);
     85                        Vx1(n,1)=directions(1,1); Vx2(n,1)=directions(1,2); Vx3(n,1)=directions(1,3);
     86                        Vy1(n,1)=directions(2,1); Vy2(n,1)=directions(2,2); Vy3(n,1)=directions(2,3);
     87                        Vz1(n,1)=directions(3,1); Vz2(n,1)=directions(3,2); Vz3(n,1)=directions(3,3);
    8888                end
    8989        end
  • issm/trunk/src/m/solutions/ice/parameters.m

    r34 r630  
    22%PARAMETERS - compute parameters
    33%
    4 %   This routine is caled by the solver to compute the parameters
     4%   This routine is called by the solver to compute the parameters
    55%   that are given in md.outputparameters (stress, strain rate,...)
    66%
Note: See TracChangeset for help on using the changeset viewer.