Changeset 630
- Timestamp:
- 05/28/09 14:59:46 (16 years ago)
- Location:
- issm/trunk/src/m/solutions/ice
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/ice/StrainRateCompute.m
r51 r630 32 32 33 33 strainratematrix=[strainratevector(1) strainratevector(3) 34 34 strainratevector(3) strainratevector(2)]; 35 35 36 36 %eigen values and vectors … … 72 72 73 73 strainratematrix=[strainratevector(1) strainratevector(4) strainratevector(5) 74 75 74 strainratevector(4) strainratevector(2) strainratevector(6) 75 strainratevector(5) strainratevector(6) strainratevector(3)]; 76 76 77 77 %eigen values and vectors 78 78 [directions,value]=eig(strainratematrix); 79 79 80 80 %Plug into global vectors 81 81 strainrate1(n,:)=strainratevector; 82 83 84 85 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); 86 86 end 87 87 end -
issm/trunk/src/m/solutions/ice/StressCompute.m
r51 r630 33 33 34 34 stressmatrix=[stressvector(1) stressvector(3) 35 35 stressvector(3) stressvector(2)]; 36 36 37 37 %eigen values and vectors 38 38 [directions,value]=eig(stressmatrix); 39 39 40 40 %Plug into global vectors 41 41 stress1(n,:)=stressvector; 42 43 44 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); 45 45 end 46 46 end … … 74 74 75 75 stressmatrix=[stressvector(1) stressvector(4) stressvector(5) 76 77 76 stressvector(4) stressvector(2) stressvector(6) 77 stressvector(5) stressvector(6) stressvector(3)]; 78 78 79 79 %eigen values and vectors 80 80 [directions,value]=eig(stressmatrix); 81 81 82 82 %Plug into global vectors 83 83 stress1(n,:)=stressvector; 84 85 86 87 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); 88 88 end 89 89 end -
issm/trunk/src/m/solutions/ice/parameters.m
r34 r630 2 2 %PARAMETERS - compute parameters 3 3 % 4 % This routine is cal ed by the solver to compute the parameters4 % This routine is called by the solver to compute the parameters 5 5 % that are given in md.outputparameters (stress, strain rate,...) 6 6 %
Note:
See TracChangeset
for help on using the changeset viewer.