Ignore:
Timestamp:
07/06/11 14:30:08 (14 years ago)
Author:
Eric.Larour
Message:

Split spcvelocity into spcvx, spcvy and spcvz

File:
1 edited

Legend:

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

    r8630 r8824  
    1111%Impose a non zero velocity on the upper boundary condition (y=max(y))
    1212pos=find(md.y==max(md.y));
    13 md.spcvelocity(pos,5)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
     13md.spcvy(pos)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
    1414
    1515%Compute solution for MacAyeal's model
Note: See TracChangeset for help on using the changeset viewer.