Changeset 4953


Ignore:
Timestamp:
08/03/10 17:23:56 (15 years ago)
Author:
Mathieu Morlighem
Message:

updated par files

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

Legend:

Unmodified
Added
Removed

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

    r4949 r4953  
    11%Start defining model parameters here
    22
    3 %Create geometry
    4 load('79Northdata.mat');
    5 md.surface=surface;
    6 md.bed=bed;
    7 md.thickness=thickness;
    8 md.vx=vx;
    9 md.vy=vy;
     3%Geometry and observation
     4load('./../Data/79North.data','-mat');
     5md.vx       =InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0);
     6md.vy       =InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0);
     7md.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y,0);
     8md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y,0);
     9md.bed=md.surface-md.thickness;
     10clear surface thickness vx vy x y index;
    1011
    11 %Put melting on ice shelves
     12%Materials
     13md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     14md.rheology_B=paterson(md.observed_temperature);
     15md.rheology_n=3*ones(md.numberofelements,1);
     16
     17%Friction
     18pos=find(md.elementoniceshelf);
     19md.drag_type=2;
     20md.drag_coefficient=50*ones(md.numberofgrids,1);
     21md.drag_coefficient(md.elements(pos,:))=0;
     22md.drag_p=ones(md.numberofelements,1);
     23md.drag_q=ones(md.numberofelements,1);
     24
     25%Ice shelf melting
    1226md.melting_rate=zeros(md.numberofgrids,1);
    1327pos=find(md.gridoniceshelf);
    1428md.melting_rate(pos)=10;
    1529
    16 %Create drag
    17 md.drag_type=2;
    18 md.drag_coefficient=50*ones(md.numberofgrids,1);
    19 %Take care of iceshelves: no basal drag
    20 pos=find(md.elementoniceshelf);
    21 md.drag_coefficient(md.elements(pos,:))=0;
    22 md.drag_p=ones(md.numberofelements,1);
    23 md.drag_q=ones(md.numberofelements,1);
     30%Numerical parameters
    2431md.viscosity_overshoot=0.3;
     32md.artificial_diffusivity=1;
     33md.verbose=0;
     34md.waitonlock=30;
     35md.np=3;
     36md.eps_res=0.05;
     37md.eps_rel=0.005;
     38md.eps_abs=NaN;
    2539
    26 %Create temperature and rheology
    27 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
    28 md.rheology_B=paterson(md.observed_temperature);
    29 md.rheology_n=3*ones(md.numberofelements,1);
    30 
    31 %Deal with boundary conditions:
     40%Boundary conditions:
    3241md=SetMarineIceSheetBC(md);
    3342pos=find(md.gridonboundary);
    3443md.spcthickness(pos,1)=1;
    3544md.spcthickness(pos,2)=md.thickness(pos);
    36 
    37 %Parallel options
    38 md.np=7;
    39 md.time=50;
    40 md.waitonlock=90;
    41 
    42 %Misc
    43 md.verbose=0;
    44 md.eps_abs=NaN;
    45 md.eps_rel=0.05;

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

    r4949 r4953  
    11%Start defining model parameters here
    22
    3 load('Pigdata.mat');
    4 md.surface=surface;
    5 md.bed=bed;
    6 md.thickness=thickness;
    7 md.vx_obs=vx_obs;
    8 md.vy_obs=vy_obs;
     3%Geometry and observation
     4load('./../Data/Pig.data','-mat');
     5md.vx_obs   =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y,0);
     6md.vy_obs   =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y,0);
     7md.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y,0);
     8md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y,0);
     9md.bed=md.surface-md.thickness;
     10clear surface thickness vx_obs vy_obs x y index;
    911
    10 %Create drag
    11 md.drag_type=2;
    12 md.drag_coefficient=50*ones(md.numberofgrids,1);
    13 %Take care of iceshelves: no basal drag
    14 pos=find(md.elementoniceshelf);
    15 md.drag_coefficient(md.elements(pos,:))=0;
    16 md.drag_p=ones(md.numberofelements,1);
    17 md.drag_q=ones(md.numberofelements,1);
    18 md.viscosity_overshoot=0.3;
    19 
    20 %Create temperature and rheology
     12%Materials
    2113md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
    2214md.rheology_B=paterson(md.observed_temperature);
    2315md.rheology_n=3*ones(md.numberofelements,1);
    2416
    25 %Deal with boundary conditions:
     17%Friction
     18pos=find(md.elementoniceshelf);
     19md.drag_type=2;
     20md.drag_coefficient=50*ones(md.numberofgrids,1);
     21md.drag_coefficient(md.elements(pos,:))=0;
     22md.drag_p=ones(md.numberofelements,1);
     23md.drag_q=ones(md.numberofelements,1);
     24
     25%Numerical parameters
     26md.viscosity_overshoot=0.3;
     27md.artificial_diffusivity=1;
     28md.verbose=0;
     29md.waitonlock=30;
     30md.np=3;
     31md.eps_res=0.05;
     32md.eps_rel=0.005;
     33md.eps_abs=NaN;
     34
     35%Boundary conditions:
    2636md=SetMarineIceSheetBC(md);
    27 
    28 %Parallel options
    29 md.np=3;
    30 md.time=50;
    31 md.waitonlock=60;
    32 
    33 %Misc
    34 md.verbose=0;
    35 md.eps_abs=NaN;
    36 md.eps_rel=0.005;

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

    r4949 r4953  
    11%Start defining model parameters here
    22
    3 %dynamics
    4 md.verbose=0;
    5 md.artificial_diffusivity=1;
    6 
     3%Geometry
    74hmin=300;
    85hmax=1000;
    … …  
    1310md.surface=md.bed+md.thickness;
    1411
    15 %Create drag
    16 md.drag_type=2; %0 none 1 plastic 2 viscous
    17 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
    18 %Take care of iceshelves: no basal drag
    19 pos=find(md.elementoniceshelf);
    20 md.drag_coefficient(md.elements(pos,:))=0;
    21 md.drag_p=ones(md.numberofelements,1);
    22 md.drag_q=ones(md.numberofelements,1);
    23 md.viscosity_overshoot=0.3;
     12%Initial velocity
     13load('./../Data/SquareShelf.data','-mat');
     14md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0);
     15md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0);
     16clear vx vy x y index;
    2417
    25 %Create temperature and rheology
     18%Materials
    2619md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
    2720md.rheology_B=paterson(md.observed_temperature);
    2821md.rheology_n=3*ones(md.numberofelements,1);
    2922
    30 %Deal with boundary conditions:
     23%Friction
     24pos=find(md.elementoniceshelf);
     25md.drag_type=2;
     26md.drag_coefficient=20*ones(md.numberofgrids,1);
     27md.drag_coefficient(md.elements(pos,:))=0;
     28md.drag_p=ones(md.numberofelements,1);
     29md.drag_q=ones(md.numberofelements,1);
     30
     31%Numerical parameters
     32md.viscosity_overshoot=0.3;
     33md.artificial_diffusivity=1;
     34md.verbose=0;
     35md.waitonlock=30;
     36md.np=3;
     37md.eps_res=0.05;
     38md.eps_rel=0.005;
     39md.eps_abs=NaN;
     40
     41%Boundary conditions:
    3142md=SetIceSheetBC(md);
    32 
    33 %Parallel options
    34 md.np=3;
    35 md.time=50;
    36 md.waitonlock=30;

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

    r4952 r4953  
    2323%Friction
    2424pos=find(md.elementoniceshelf);
     25md.drag_type=2;
    2526md.drag_coefficient=20*ones(md.numberofgrids,1);
    2627md.drag_coefficient(md.elements(pos,:))=0;
    … …  
    3031%Numerical parameters
    3132md.viscosity_overshoot=0.3;
     33md.artificial_diffusivity=1;
     34md.waitonlock=30;
    3235md.verbose=0;
    33 md.artificial_diffusivity=1;
    3436md.np=3;
    35 md.waitonlock=30;
     37md.eps_res=0.10;
     38md.eps_rel=0.02;
    3639md.eps_abs=NaN;
    37 md.eps_rel=0.02;
    38 md.eps_res=1.10;
    3940
    4041%Boundary conditions:

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

    r4949 r4953  
    11%Start defining model parameters here
    22
    3 %dynamics
    4 md.verbose=0;
    5 md.artificial_diffusivity=1;
    6 
     3%Geometry
    74hmin=300;
    85hmax=1000;
    … …  
    1310md.surface=md.bed+md.thickness;
    1411
    15 md.drag_type=2; %0 none 1 plastic 2 viscous
    16 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
    17 %Take care of iceshelves: no basal drag
    18 pos=find(md.elementoniceshelf);
    19 md.drag_coefficient(md.elements(pos,:))=0;
    20 md.drag_p=ones(md.numberofelements,1);
    21 md.drag_q=ones(md.numberofelements,1);
    22 md.viscosity_overshoot=0.3;
     12%Initial velocity
     13load('./../Data/SquareShelfConstrained.data','-mat');
     14md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0);
     15md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0);
     16clear vx vy x y index;
    2317
    24 %Create temperature and rheology
     18%Materials
    2519md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
    2620md.rheology_B=paterson(md.observed_temperature);
    2721md.rheology_n=3*ones(md.numberofelements,1);
    2822
     23%Friction
     24pos=find(md.elementoniceshelf);
     25md.drag_type=2;
     26md.drag_coefficient=20*ones(md.numberofgrids,1);
     27md.drag_coefficient(md.elements(pos,:))=0;
     28md.drag_p=ones(md.numberofelements,1);
     29md.drag_q=ones(md.numberofelements,1);
     30
     31%Numerical parameters
     32md.viscosity_overshoot=0.3;
     33md.artificial_diffusivity=1;
     34md.verbose=0;
     35md.waitonlock=30;
     36md.np=3;
     37md.eps_res=0.05;
     38md.eps_rel=0.005;
     39md.eps_abs=NaN;
     40
    2941%Deal with boundary conditions:
    30 md=SetIceSheetBC(md);
    31 
    32 %Parallel options
    33 md.np=3;
    34 md.time=50;
    35 md.waitonlock=30;
     42md=SetIceShelfBC(md);
Note: See TracChangeset for help on using the changeset viewer.