Changeset 1024


Ignore:
Timestamp:
06/19/09 09:05:06 (16 years ago)
Author:
Mathieu Morlighem
Message:

updategridposition now use S

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

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticHoriz.m

    r121 r1024  
    143143grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    144144grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     145grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    145146grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    146147grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticHutter.m

    r1 r1024  
    102102grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    103103grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     104grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    104105grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    105106grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticStokes.m

    r121 r1024  
    9797grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    9898grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     99grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    99100grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    100101grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticVert.m

    r334 r1024  
    8989grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    9090grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     91grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    9192grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    9293grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorMelting.m

    r51 r1024  
    8686grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    8787grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     88grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    8889grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    8990grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorPrognostic.m

    r32 r1024  
    7979grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    8080grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     81grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    8182grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    8283grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorSlopeCompute.m

    r51 r1024  
    7979grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    8080grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     81grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    8182grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    8283grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m

    r994 r1024  
    103103grids(pos)=SetStructureField(grids(pos),'grid','y',md.y(pos));
    104104grids(pos)=SetStructureField(grids(pos),'grid','z',md.z(pos));
     105grids(pos)=SetStructureField(grids(pos),'grid','s',(md.z(pos)-md.bed(pos))./md.thickness(pos));
    105106grids(pos)=SetStructureField(grids(pos),'grid','onbed',md.gridonbed(pos));
    106107grids(pos)=SetStructureField(grids(pos),'grid','border',bordergrids(pos));

  • issm/trunk/src/m/solutions/ice/UpdateGridPosition.m

    r987 r1024  
    1010count2d=[1:md.numberofgrids2d]';
    1111
    12 %for i=1:md.numlayers,
    13 %       count=count2d+(i-1)*md.numberofgrids2d;
    14 %       for j=1:size(count,1),
    15 %               grids(count(j)).grid.z=bed(count(j))+(i-1)/(md.numlayers-1)*(thickness(count(j)));
    16 %       end
    17 %end
    18 
    19 %Only update the lower and upper layer.
    20 count=md.numberofgrids2d*(md.numlayers-1);
    21 for i=1:md.numberofgrids2d,
    22         grids(i).grid.z=bed(i);
    23         grids(count+i).grid.z=bed(i)+ thickness(i);
     12for i=1:md.numlayers,
     13        count=count2d+(i-1)*md.numberofgrids2d;
     14        for j=1:size(count,1),
     15                grids(count(j)).grid.z=bed(count(j))+grids(count(j)).grid.s*thickness(count(j));
     16        end
    2417end
Note: See TracChangeset for help on using the changeset viewer.