Changeset 8305


Ignore:
Timestamp:
05/16/11 15:50:32 (14 years ago)
Author:
seroussi
Message:

grid in par files

Location:
issm/trunk/test/Par
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/Par/79North.par

    r7332 r8305  
    1111
    1212%Materials
    13 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     13md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    1414md.rheology_B=paterson(md.observed_temperature);
    1515md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    1818pos=find(md.elementoniceshelf);
    1919md.drag_type=2;
    20 md.drag_coefficient=50*ones(md.numberofgrids,1);
     20md.drag_coefficient=50*ones(md.numberofnodes,1);
    2121md.drag_coefficient(md.elements(pos,:))=0;
    2222md.drag_p=ones(md.numberofelements,1);
    … …  
    2424
    2525%Ice shelf melting and accumulation
    26 md.melting_rate=zeros(md.numberofgrids,1);
    27 pos=zeros(md.numberofgrids,1);
     26md.melting_rate=zeros(md.numberofnodes,1);
     27pos=zeros(md.numberofnodes,1);
    2828pos(md.elements(find(md.elementoniceshelf),:))=1;
    2929md.melting_rate(find(pos))=10;
    30 md.accumulation_rate=15*ones(md.numberofgrids,1);
     30md.accumulation_rate=15*ones(md.numberofnodes,1);
    3131
    3232%Numerical parameters

  • issm/trunk/test/Par/ISMIPA.par

    r5966 r8305  
    55md.bed=md.surface-1000+500*sin(md.x*2*pi/max(md.x)).*sin(md.y*2*pi/max(md.x));
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofgrids,1);
     7md.firn_layer=0*ones(md.numberofnodes,1);
    88
    99disp('      creating drag');
    1010md.drag_type=2; %0 none 1 plastic 2 viscous
    11 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
     11md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
    1212%Take care of iceshelves: no basal drag
    1313pos=find(md.elementoniceshelf);
    … …  
    1717
    1818disp('      creating flow law paramter');
    19 md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
     19md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
    2020md.rheology_n=3*ones(md.numberofelements,1);
    2121

  • issm/trunk/test/Par/ISMIPB.par

    r5966 r8305  
    55md.bed=md.surface-1000+500*sin(md.x*2*pi/max(md.x));
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofgrids,1);
     7md.firn_layer=0*ones(md.numberofnodes,1);
    88
    99disp('      creating drag');
    1010md.drag_type=2; %0 none 1 plastic 2 viscous
    11 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
     11md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
    1212%Take care of iceshelves: no basal drag
    1313pos=find(md.elementoniceshelf);
    … …  
    1717
    1818disp('      creating flow law paramter');
    19 md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
     19md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
    2020md.rheology_n=3*ones(md.numberofelements,1);
    2121

  • issm/trunk/test/Par/ISMIPC.par

    r7877 r8305  
    55md.bed=md.surface-1000;
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofgrids,1);
     7md.firn_layer=0*ones(md.numberofnodes,1);
    88
    99disp('      creating drag');
    … …  
    1818
    1919disp('      creating flow law paramter');
    20 md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
     20md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
    2121md.rheology_n=3*ones(md.numberofelements,1);
    2222

  • issm/trunk/test/Par/ISMIPD.par

    r7963 r8305  
    55md.bed=md.surface-1000;
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofgrids,1);
     7md.firn_layer=0*ones(md.numberofnodes,1);
    88
    99disp('      creating drag');
    … …  
    1717
    1818disp('      creating flow law paramter');
    19 md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
     19md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
    2020md.rheology_n=3*ones(md.numberofelements,1);
    2121

  • issm/trunk/test/Par/ISMIPE.par

    r5966 r8305  
    44data=load('../Data/ISMIPE.data','-mat');
    55data=data.data;
    6 md.surface=zeros(md.numberofgrids,1);
    7 md.bed=zeros(md.numberofgrids,1);
    8 for i=1:md.numberofgrids
     6md.surface=zeros(md.numberofnodes,1);
     7md.bed=zeros(md.numberofnodes,1);
     8for i=1:md.numberofnodes
    99        y=md.y(i);
    1010        point1=floor(y/100)+1;
    … …  
    1515end
    1616md.thickness=md.surface-md.bed;
    17 md.firn_layer=0*ones(md.numberofgrids,1);
     17md.firn_layer=0*ones(md.numberofnodes,1);
    1818md.thickness(find(~md.thickness))=0.01;
    1919md.bed=md.surface-md.thickness;
    … …  
    2121disp('      creating drag');
    2222md.drag_type=2; %0 none 1 plastic 2 viscous
    23 md.drag_coefficient=zeros(md.numberofgrids,1);
     23md.drag_coefficient=zeros(md.numberofnodes,1);
    2424md.drag_p=ones(md.numberofelements,1);
    2525md.drag_q=ones(md.numberofelements,1);
    2626
    2727disp('      creating flow law paramter');
    28 md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
     28md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
    2929md.rheology_n=3*ones(md.numberofelements,1);
    3030

  • issm/trunk/test/Par/ISMIPF.par

    r8244 r8305  
    77md.bed=md.surface-1000+100*exp(-((md.x-max(md.x)/2).^2+(md.y-max(md.y)/2).^2)/(10000^2));
    88md.thickness=md.surface-md.bed;
    9 md.firn_layer=0*ones(md.numberofgrids,1);
     9md.firn_layer=0*ones(md.numberofnodes,1);
    1010
    1111disp('      creating drag');
    1212md.drag_type=2; %0 none 1 plastic 2 viscous
    13 md.drag_coefficient=sqrt(md.yts/(2.140373*10^-7*1000))*ones(md.numberofgrids,1);
     13md.drag_coefficient=sqrt(md.yts/(2.140373*10^-7*1000))*ones(md.numberofnodes,1);
    1414md.drag_p=ones(md.numberofelements,1);
    1515md.drag_q=zeros(md.numberofelements,1);
    1616
    1717disp('      creating flow law paramter');
    18 md.rheology_B=1.4734*10^14*ones(md.numberofgrids,1);
     18md.rheology_B=1.4734*10^14*ones(md.numberofnodes,1);
    1919md.rheology_n=1*ones(md.numberofelements,1);
    2020md.rheology_law=NoneEnum;
    … …  
    2323%Create grid on boundary fist (because we cannot use mesh)
    2424md=SetIceSheetBC(md);
    25 md.spcvelocity(:,4)=100*ones(md.numberofgrids,1);
    26 md.vx=zeros(md.numberofgrids,1);
    27 md.vy=zeros(md.numberofgrids,1);
    28 md.vz=zeros(md.numberofgrids,1);
    29 md.vel=zeros(md.numberofgrids,1);
    30 md.pressure=zeros(md.numberofgrids,1);
    31 md.temperature=255*ones(md.numberofgrids,1);
     25md.spcvelocity(:,4)=100*ones(md.numberofnodes,1);
     26md.vx=zeros(md.numberofnodes,1);
     27md.vy=zeros(md.numberofnodes,1);
     28md.vz=zeros(md.numberofnodes,1);
     29md.vel=zeros(md.numberofnodes,1);
     30md.pressure=zeros(md.numberofnodes,1);
     31md.temperature=255*ones(md.numberofnodes,1);
    3232pos=find(md.x==min(md.x) | md.x==max(md.x) | md.y==min(md.y) | md.y==max(md.y));
    3333md.spcthickness(pos,1)=1;
    3434md.spcthickness(pos,2)=md.thickness(pos);
    35 md.spctemperature=[ones(md.numberofgrids,1) 255*ones(md.numberofgrids,1)];
    36 md.geothermalflux=0.4*ones(md.numberofgrids,1);
     35md.spctemperature=[ones(md.numberofnodes,1) 255*ones(md.numberofnodes,1)];
     36md.geothermalflux=0.4*ones(md.numberofnodes,1);
    3737
    3838%Parallel options

  • issm/trunk/test/Par/Pig.par

    r6307 r8305  
    1111md.vx=md.vx_obs;
    1212md.vy=md.vy_obs;
    13 md.vz=zeros(md.numberofgrids,1);
    14 md.pressure=zeros(md.numberofgrids,1);
     13md.vz=zeros(md.numberofnodes,1);
     14md.pressure=zeros(md.numberofnodes,1);
    1515
    1616%Materials
    17 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     17md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    1818md.rheology_B=paterson(md.observed_temperature);
    1919md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    2323pos=find(md.elementoniceshelf);
    2424md.drag_type=2;
    25 md.drag_coefficient=50*ones(md.numberofgrids,1);
     25md.drag_coefficient=50*ones(md.numberofnodes,1);
    2626md.drag_coefficient(md.elements(pos,:))=0;
    2727md.drag_p=ones(md.numberofelements,1);

  • issm/trunk/test/Par/RoundSheetEISMINT.par

    r5966 r8305  
    11%Ok, start defining model parameters here
    22disp('      creating thickness');
    3 md.thickness=10*ones(md.numberofgrids,1);
    4 md.bed=zeros(md.numberofgrids,1);
     3md.thickness=10*ones(md.numberofnodes,1);
     4md.bed=zeros(md.numberofnodes,1);
    55md.surface=md.bed+md.thickness;
    66
    77disp('      creating drag');
    88md.drag_type=2; %0 none 1 plastic 2 viscous
    9 md.drag_coefficient=20*ones(md.numberofgrids,1); %q=1. %no drag is specified in the analytical solution
     9md.drag_coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
    1010md.drag_p=ones(md.numberofelements,1);
    1111md.drag_q=ones(md.numberofelements,1);
    … …  
    1717md.observed_temperature=(tmin+st*radius);
    1818md.temperature=md.observed_temperature;
    19 md.geothermalflux=4.2*10^-2*ones(md.numberofgrids,1);
     19md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
    2020
    2121disp('      creating flow law paramter');
    22 md.rheology_B=6.81*10^(7)*ones(md.numberofgrids,1); %to have the same B as the analytical solution
     22md.rheology_B=6.81*10^(7)*ones(md.numberofnodes,1); %to have the same B as the analytical solution
    2323md.rheology_n=3*ones(md.numberofelements,1);
    2424
    … …  
    3434md.vy_obs=constant/2*md.y.*(md.thickness).^-1;
    3535md.vel_obs=(sqrt((md.vx_obs).^2+(md.vy_obs).^2));
    36 md.vx=zeros(md.numberofgrids,1);
    37 md.vy=zeros(md.numberofgrids,1);
    38 md.vz=zeros(md.numberofgrids,1);
    39 md.pressure=zeros(md.numberofgrids,1);
     36md.vx=zeros(md.numberofnodes,1);
     37md.vy=zeros(md.numberofnodes,1);
     38md.vz=zeros(md.numberofnodes,1);
     39md.pressure=zeros(md.numberofnodes,1);
    4040
    4141%Deal with boundary conditions:

  • issm/trunk/test/Par/RoundSheetShelf.par

    r8244 r8305  
    3434
    3535%Initial velocity
    36 md.vx=zeros(md.numberofgrids,1);
    37 md.vy=zeros(md.numberofgrids,1);
    38 md.vz=zeros(md.numberofgrids,1);
    39 md.pressure=zeros(md.numberofgrids,1);
     36md.vx=zeros(md.numberofnodes,1);
     37md.vy=zeros(md.numberofnodes,1);
     38md.vz=zeros(md.numberofnodes,1);
     39md.pressure=zeros(md.numberofnodes,1);
    4040
    4141%Materials
    42 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     42md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    4343md.rheology_B=paterson(md.observed_temperature);
    4444md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    4646
    4747%Accumulation and melting
    48 md.accumulation_rate=-10*ones(md.numberofgrids,1);
    49 md.melting_rate=zeros(md.numberofgrids,1);
     48md.accumulation_rate=-10*ones(md.numberofnodes,1);
     49md.melting_rate=zeros(md.numberofnodes,1);
    5050pos=find(md.gridoniceshelf);md.melting_rate(pos)=10;
    5151
    … …  
    5353pos=find(md.elementoniceshelf);
    5454md.drag_type=2;
    55 md.drag_coefficient=20*ones(md.numberofgrids,1);
     55md.drag_coefficient=20*ones(md.numberofnodes,1);
    5656md.drag_coefficient(md.elements(pos,:))=0;
    5757md.drag_p=ones(md.numberofelements,1);
    … …  
    7575
    7676%Deal with boundary conditions:
    77 md.spcvelocity=zeros(md.numberofgrids,6);
     77md.spcvelocity=zeros(md.numberofnodes,6);
    7878pos=find(md.x==0 & md.y==0);
    7979md.spcvelocity(pos,1:2)=1;
    … …  
    8181pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))+AirEnum*md.elementonicesheet(pressureload(:,end))];
    8282md.pressureload=pressureload;
    83 md.spcthickness=zeros(md.numberofgrids,2);
    84 md.diagnostic_ref=NaN*ones(md.numberofgrids,6);
     83md.spcthickness=zeros(md.numberofnodes,2);
     84md.diagnostic_ref=NaN*ones(md.numberofnodes,6);
    8585
    8686%Change name so that no test have the same name

  • issm/trunk/test/Par/RoundSheetStaticEISMINT.par

    r5966 r8305  
    55radiusmax=max(radius);
    66md.thickness=hmin*ones(size(md.x,1),1)+hmax*(4*((1/2)^(4/3)*ones(size(md.x,1),1)-((radius)./(2*radiusmax)).^(4/3))).^(3/8);
    7 md.firn_layer=10*ones(md.numberofgrids,1);
     7md.firn_layer=10*ones(md.numberofnodes,1);
    88md.bed=0*md.thickness;
    99md.surface=md.bed+md.thickness;
    … …  
    1111disp('      creating drag');
    1212md.drag_type=2; %0 none 1 plastic 2 viscous
    13 md.drag_coefficient=20*ones(md.numberofgrids,1); %q=1. %no drag is specified in the analytical solution
     13md.drag_coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
    1414%Take care of iceshelves: no basal drag
    1515pos=find(md.elementoniceshelf);
    … …  
    2323md.observed_temperature=(tmin+st*radius);
    2424md.temperature=md.observed_temperature;
    25 md.geothermalflux=4.2*10^-2*ones(md.numberofgrids,1);
     25md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
    2626
    2727disp('      creating flow law paramter');
    28 md.rheology_B=6.81*10^(7)*ones(md.numberofgrids,1); %to have the same B as the analytical solution
     28md.rheology_B=6.81*10^(7)*ones(md.numberofnodes,1); %to have the same B as the analytical solution
    2929md.rheology_n=3*ones(md.numberofelements,1);
    3030
    … …  
    4040md.vy_obs=constant/2*md.y.*(md.thickness).^-1;
    4141md.vel_obs=(sqrt((md.vx_obs).^2+(md.vy_obs).^2));
    42 md.vx=zeros(md.numberofgrids,1);
    43 md.vy=zeros(md.numberofgrids,1);
    44 md.vz=zeros(md.numberofgrids,1);
    45 md.pressure=zeros(md.numberofgrids,1);
     42md.vx=zeros(md.numberofnodes,1);
     43md.vy=zeros(md.numberofnodes,1);
     44md.vz=zeros(md.numberofnodes,1);
     45md.pressure=zeros(md.numberofnodes,1);
    4646
    4747%Deal with boundary conditions:

  • issm/trunk/test/Par/SquareEISMINT.par

    r6245 r8305  
    44ymin=min(md.y);
    55ymax=max(md.y);
    6 md.thickness=500*ones(md.numberofgrids,1);
    7 md.firn_layer=0*ones(md.numberofgrids,1);
     6md.thickness=500*ones(md.numberofnodes,1);
     7md.firn_layer=0*ones(md.numberofnodes,1);
    88md.bed=-md.rho_ice/md.rho_water*md.thickness;
    99md.surface=md.bed+md.thickness;
    … …  
    1111disp('      creating drag');
    1212md.drag_type=2; %0 none 1 plastic 2 viscous
    13 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
     13md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
    1414%Take care of iceshelves: no basal drag
    1515pos=find(md.elementoniceshelf);
    … …  
    1919
    2020disp('      creating temperature');
    21 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     21md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    2222
    2323disp('      creating flow law paramter');
    2424%md.B=paterson(md.observed_temperature);
    25 md.rheology_B=1.7687*10^8*ones(md.numberofgrids,1);
     25md.rheology_B=1.7687*10^8*ones(md.numberofnodes,1);
    2626md.rheology_n=3*ones(md.numberofelements,1);
    2727
    2828disp('      creating accumulation rates');
    29 md.accumulation_rate=0.2*ones(md.numberofgrids,1); %0m/a
    30 md.melting_rate=0*ones(md.numberofgrids,1); %0m/a
     29md.accumulation_rate=0.2*ones(md.numberofnodes,1); %0m/a
     30md.melting_rate=0*ones(md.numberofnodes,1); %0m/a
    3131
    3232disp('      boundary conditions ');
    … …  
    3737
    3838%Evolution of the ice shelf
    39 pos=find(md.y==200000); %grids on the upper boundary condition
     39pos=find(md.y==200000); %nodes on the upper boundary condition
    4040md.spcthickness(pos,1)=1;
    4141md.spcthickness(pos,2)=500;

  • issm/trunk/test/Par/SquareSheetConstrained.par

    r7542 r8305  
    1515md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
    1616clear vx vy x y index;
    17 md.vz=zeros(md.numberofgrids,1);
    18 md.pressure=zeros(md.numberofgrids,1);
     17md.vz=zeros(md.numberofnodes,1);
     18md.pressure=zeros(md.numberofnodes,1);
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     21md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    2222md.rheology_B=paterson(md.observed_temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    2727pos=find(md.elementoniceshelf);
    2828md.drag_type=2;
    29 md.drag_coefficient=20*ones(md.numberofgrids,1);
     29md.drag_coefficient=20*ones(md.numberofnodes,1);
    3030md.drag_coefficient(md.elements(pos,:))=0;
    3131md.drag_p=ones(md.numberofelements,1);

  • issm/trunk/test/Par/SquareSheetShelf.par

    r6307 r8305  
    1818md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
    1919clear vx vy x y index;
    20 md.vz=zeros(md.numberofgrids,1);
    21 md.pressure=zeros(md.numberofgrids,1);
     20md.vz=zeros(md.numberofnodes,1);
     21md.pressure=zeros(md.numberofnodes,1);
    2222
    2323%Materials
    24 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     24md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    2525md.rheology_B=paterson(md.observed_temperature);
    2626md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    2828
    2929%Accumulation and melting
    30 md.accumulation_rate=10*ones(md.numberofgrids,1);
    31 md.melting_rate=5*ones(md.numberofgrids,1);
     30md.accumulation_rate=10*ones(md.numberofnodes,1);
     31md.melting_rate=5*ones(md.numberofnodes,1);
    3232
    3333%Friction
    3434pos=find(md.elementoniceshelf);
    3535md.drag_type=2;
    36 md.drag_coefficient=20*ones(md.numberofgrids,1);
     36md.drag_coefficient=20*ones(md.numberofnodes,1);
    3737md.drag_coefficient(md.elements(pos,:))=0;
    3838md.drag_p=ones(md.numberofelements,1);

  • issm/trunk/test/Par/SquareShelf.par

    r7542 r8305  
    1515md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
    1616clear vx vy x y index;
    17 md.vz=zeros(md.numberofgrids,1);
    18 md.pressure=zeros(md.numberofgrids,1);
     17md.vz=zeros(md.numberofnodes,1);
     18md.pressure=zeros(md.numberofnodes,1);
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     21md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    2222md.rheology_B=paterson(md.observed_temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    2727pos=find(md.elementoniceshelf);
    2828md.drag_type=2;
    29 md.drag_coefficient=20*ones(md.numberofgrids,1);
     29md.drag_coefficient=20*ones(md.numberofnodes,1);
    3030md.drag_coefficient(md.elements(pos,:))=0;
    3131md.drag_p=ones(md.numberofelements,1);

  • issm/trunk/test/Par/SquareShelfConstrained.par

    r7542 r8305  
    1515md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
    1616clear vx vy x y index;
    17 md.vz=zeros(md.numberofgrids,1);
    18 md.pressure=zeros(md.numberofgrids,1);
     17md.vz=zeros(md.numberofnodes,1);
     18md.pressure=zeros(md.numberofnodes,1);
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     21md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    2222md.rheology_B=paterson(md.observed_temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    2525
    2626%Accumulation and melting
    27 md.accumulation_rate=10*ones(md.numberofgrids,1);
    28 md.melting_rate=5*ones(md.numberofgrids,1);
     27md.accumulation_rate=10*ones(md.numberofnodes,1);
     28md.melting_rate=5*ones(md.numberofnodes,1);
    2929
    3030%Friction
    3131pos=find(md.elementoniceshelf);
    3232md.drag_type=2;
    33 md.drag_coefficient=20*ones(md.numberofgrids,1);
     33md.drag_coefficient=20*ones(md.numberofnodes,1);
    3434md.drag_coefficient(md.elements(pos,:))=0;
    3535md.drag_p=ones(md.numberofelements,1);

  • issm/trunk/test/Par/SquareThermal.par

    r5966 r8305  
    77disp('      creating thickness');
    88h=1000;
    9 md.thickness=h*ones(md.numberofgrids,1);
    10 md.firn_layer=10*ones(md.numberofgrids,1);
    11 md.bed=-1000*ones(md.numberofgrids,1);
     9md.thickness=h*ones(md.numberofnodes,1);
     10md.firn_layer=10*ones(md.numberofnodes,1);
     11md.bed=-1000*ones(md.numberofnodes,1);
    1212md.surface=md.bed+md.thickness;
    1313
    1414disp('      creating velocities');
    15 md.vx=zeros(md.numberofgrids,1);
    16 md.vy=zeros(md.numberofgrids,1);
    17 md.vz=zeros(md.numberofgrids,1);
     15md.vx=zeros(md.numberofnodes,1);
     16md.vy=zeros(md.numberofnodes,1);
     17md.vz=zeros(md.numberofnodes,1);
    1818
    1919disp('      creating drag');
    2020md.drag_type=2; %0 none 1 plastic 2 viscous
    21 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
     21md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
    2222%Take care of iceshelves: no basal drag
    2323pos=find(md.elementoniceshelf);
    … …  
    2727
    2828disp('      creating temperatures');
    29 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     29md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    3030
    3131disp('      creating flow law paramter');
    … …  
    3434
    3535disp('      creating accumulation rates');
    36 md.accumulation_rate=ones(md.numberofgrids,1)/md.yts; %1m/a
    37 md.melting_rate=0*ones(md.numberofgrids,1)/md.yts; %1m/a
     36md.accumulation_rate=ones(md.numberofnodes,1)/md.yts; %1m/a
     37md.melting_rate=0*ones(md.numberofnodes,1)/md.yts; %1m/a
    3838
    3939%Deal with boundary conditions:
    … …  
    4444disp('      boundary conditions for thermal model');
    4545md.temperature=md.observed_temperature;
    46 md.spctemperature(:,1)=ones(md.numberofgrids,1); %surface temperature
     46md.spctemperature(:,1)=ones(md.numberofnodes,1); %surface temperature
    4747md.spctemperature(:,2)=md.observed_temperature;
    48 md.geothermalflux=zeros(md.numberofgrids,1);
     48md.geothermalflux=zeros(md.numberofnodes,1);
    4949pos=find(md.elementonicesheet);md.geothermalflux(md.elements(pos,:))=1*10^-3; %1 mW/m^2
Note: See TracChangeset for help on using the changeset viewer.