Changeset 9045


Ignore:
Timestamp:
07/18/11 10:51:43 (14 years ago)
Author:
seroussi
Message:

fixed spctemperature

Location:
issm/trunk/test/NightlyRun
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test1301.m

    r8630 r9045  
    1616md.observed_temperature=273.15*ones(md.numberofnodes,1);
    1717md.temperature=md.observed_temperature;
    18 md.spctemperature(:,2)=md.observed_temperature;
     18pos=find(md.nodeonsurface);
     19md.spctemperature(pos)=md.observed_temperature(pos);
    1920md.rheology_B=paterson(md.temperature);
    2021
  • issm/trunk/test/NightlyRun/test1302.m

    r8630 r9045  
    1212
    1313%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;
     14pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3))=10;
     15pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
    1616md.vz=0.1*ones(md.numberofnodes,1);
    1717md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2);
  • issm/trunk/test/NightlyRun/test1303.m

    r8630 r9045  
    1111md=extrude(md,11,2);
    1212md=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;
     13pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3))=10;
     14pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
    1515md.pressure=zeros(md.numberofnodes,1);
    1616
  • issm/trunk/test/NightlyRun/test1304.m

    r8630 r9045  
    1212md=setelementstype(md,'Pattyn','all');
    1313
    14 pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
     14pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
    1515md.pressure=zeros(md.numberofnodes,1);
    1616md.geothermalflux(:)=0.1; %100mW/m^2
Note: See TracChangeset for help on using the changeset viewer.