Changeset 9045
- Timestamp:
- 07/18/11 10:51:43 (14 years ago)
- Location:
- issm/trunk/test/NightlyRun
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1301.m
r8630 r9045 16 16 md.observed_temperature=273.15*ones(md.numberofnodes,1); 17 17 md.temperature=md.observed_temperature; 18 md.spctemperature(:,2)=md.observed_temperature; 18 pos=find(md.nodeonsurface); 19 md.spctemperature(pos)=md.observed_temperature(pos); 19 20 md.rheology_B=paterson(md.temperature); 20 21 -
issm/trunk/test/NightlyRun/test1302.m
r8630 r9045 12 12 13 13 %Thermal boundary conditions 14 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3) ,1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;15 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6) ,1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;14 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3))=10; 15 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0; 16 16 md.vz=0.1*ones(md.numberofnodes,1); 17 17 md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2); -
issm/trunk/test/NightlyRun/test1303.m
r8630 r9045 11 11 md=extrude(md,11,2); 12 12 md=setelementstype(md,'Pattyn','all'); 13 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3) ,1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;14 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6) ,1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;13 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3))=10; 14 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0; 15 15 md.pressure=zeros(md.numberofnodes,1); 16 16 -
issm/trunk/test/NightlyRun/test1304.m
r8630 r9045 12 12 md=setelementstype(md,'Pattyn','all'); 13 13 14 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6) ,1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;14 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0; 15 15 md.pressure=zeros(md.numberofnodes,1); 16 16 md.geothermalflux(:)=0.1; %100mW/m^2
Note:
See TracChangeset
for help on using the changeset viewer.