Changeset 14106


Ignore:
Timestamp:
12/07/12 11:09:49 (12 years ago)
Author:
jschierm
Message:

NEW: Parameter files for python tests 1200.

Location:
issm/trunk-jpl/test/Par
Files:
3 added
3 edited

Legend:

Unmodified
Added
Removed

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

    r13137 r14106  
    11%Ok, start defining model parameters here
    22disp('      creating thickness');
    3 md.geometry.thickness=10*ones(md.mesh.numberofvertices,1);
     3md.geometry.thickness=10.*ones(md.mesh.numberofvertices,1);
    44md.geometry.bed=zeros(md.mesh.numberofvertices,1);
    55md.geometry.surface=md.geometry.bed+md.geometry.thickness;
    66
    77disp('      creating drag');
    8 md.friction.coefficient=20*ones(md.mesh.numberofvertices,1); %q=1. %no drag is specified in the analytical solution
     8md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1); %q=1. no drag is specified in the analytical solution
    99md.friction.p=ones(md.mesh.numberofelements,1);
    1010md.friction.q=ones(md.mesh.numberofelements,1);
    … …  
    1212disp('      creating temperatures');
    1313tmin=238.15; %K
    14 st=1.67*10^-2/1000; %k/m;
     14st=1.67*10^-2/1000.; %k/m
    1515radius=sqrt((md.mesh.x).^2+(md.mesh.y).^2);
    1616md.initialization.temperature=(tmin+st*radius);
    1717md.basalforcings.geothermalflux=4.2*10^-2*ones(md.mesh.numberofvertices,1);
    1818
    19 disp('      creating flow law paramter');
    20 md.materials.rheology_B=6.81*10^(7)*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
    21 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     19disp('      creating flow law parameter');
     20md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
     21md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2222
    2323disp('      creating surface mass balance');
    2424smb_max=0.5; %m/yr
    25 sb=10^-2/1000; %m/yr/m
    26 rel=450*1000; %m
     25sb=10^-2/1000.; %m/yr/m
     26rel=450.*1000.; %m
    2727md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
    2828
    2929disp('      creating velocities');
    3030constant=0.3;
    31 md.inversion.vx_obs=constant/2*md.mesh.x.*(md.geometry.thickness).^-1;
    32 md.inversion.vy_obs=constant/2*md.mesh.y.*(md.geometry.thickness).^-1;
     31md.inversion.vx_obs=constant/2.*md.mesh.x.*(md.geometry.thickness).^-1;
     32md.inversion.vy_obs=constant/2.*md.mesh.y.*(md.geometry.thickness).^-1;
    3333md.inversion.vel_obs=(sqrt((md.inversion.vx_obs).^2+(md.inversion.vy_obs).^2));
    3434md.initialization.vx=zeros(md.mesh.numberofvertices,1);
    … …  
    3838
    3939%Deal with boundary conditions:
    40 disp('      boundary conditions for diagnostic model: ');
     40disp('      boundary conditions for diagnostic model:');
    4141md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp');
    4242
    43 radius=sqrt((md.mesh.x).*md.mesh.x+(md.mesh.y).*md.mesh.y);
     43radius=sqrt((md.mesh.x).^2+(md.mesh.y).^2);
    4444pos=find(radius==min(radius));
    45 md.mesh.x(pos)=0; md.mesh.y(pos)=0; %the closest node to the center is changed to be exactly at the center
     45md.mesh.x(pos)=0.; md.mesh.y(pos)=0.; %the closest node to the center is changed to be exactly at the center
    4646
    47 md.diagnostic.spcvx(pos)=0;
    48 md.diagnostic.spcvy(pos)=0;
    49 md.diagnostic.spcvz(pos)=0;
     47md.diagnostic.spcvx(pos)=0.;
     48md.diagnostic.spcvy(pos)=0.;
     49md.diagnostic.spcvz(pos)=0.;
    5050
    5151%parallel options
    52 md.timestepping.final_time=50000;
     52md.timestepping.final_time=50000.;
    5353
    5454%Constants
    55 md.materials.rho_ice=910;
     55md.materials.rho_ice=910.;
    5656md.materials.thermalconductivity=2.1;
    5757md.materials.latentheat=3.35*10^5;
    5858md.materials.beta=8.66*10^-4/(md.materials.rho_ice*md.constants.g); %conversion from K/m to K/Pa
    59 md.constants.yts=31556926;
     59md.constants.yts=31556926.;

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

    r13137 r14106  
    22hmin=0.01;
    33hmax=2756.7;
    4 radius=(sqrt((md.mesh.x).^2+(md.mesh.y).^2));
     4radius=sqrt((md.mesh.x).^2+(md.mesh.y).^2);
    55radiusmax=max(radius);
    6 md.geometry.thickness=hmin*ones(size(md.mesh.x,1),1)+hmax*(4*((1/2)^(4/3)*ones(size(md.mesh.x,1),1)-((radius)./(2*radiusmax)).^(4/3))).^(3/8);
    7 md.geometry.bed=0*md.geometry.thickness;
     6md.geometry.thickness=hmin*ones(size(md.mesh.x,1),1)+hmax*(4.*((1./2.)^(4./3.)*ones(size(md.mesh.x,1),1)-((radius)./(2.*radiusmax)).^(4./3.))).^(3./8.);
     7md.geometry.bed=0.*md.geometry.thickness;
    88md.geometry.surface=md.geometry.bed+md.geometry.thickness;
    99
    1010disp('      creating drag');
    11 md.friction.coefficient=20*ones(md.mesh.numberofvertices,1); %q=1. %no drag is specified in the analytical solution
     11md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1); %q=1. no drag is specified in the analytical solution
    1212%Take care of iceshelves: no basal drag
    1313pos=find(md.mask.elementonfloatingice);
    14 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     14md.friction.coefficient(md.mesh.elements(pos,:))=0.;
    1515md.friction.p=ones(md.mesh.numberofelements,1);
    1616md.friction.q=ones(md.mesh.numberofelements,1);
    … …  
    1818disp('      creating temperatures');
    1919tmin=238.15; %K
    20 st=1.67*10^-2/1000; %k/m;
    21 md.initialization.temperature=(tmin+st*radius);
     20st=1.67*10^-2/1000.; %k/m
     21md.initialization.temperature=tmin+st*radius;
    2222md.basalforcings.geothermalflux=4.2*10^-2*ones(md.mesh.numberofvertices,1);
    2323
    24 disp('      creating flow law paramter');
    25 md.materials.rheology_B=6.81*10^(7)*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
    26 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     24disp('      creating flow law parameter');
     25md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
     26md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2727
    2828disp('      creating surface mass balance');
    2929smb_max=0.5; %m/yr
    30 sb=10^-2/1000; %m/yr/m
    31 rel=450*1000; %m
     30sb=10^-2/1000.; %m/yr/m
     31rel=450.*1000.; %m
    3232md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
    3333
    3434disp('      creating velocities');
    3535constant=0.3;
    36 md.inversion.vx_obs=constant/2*md.mesh.x.*(md.geometry.thickness).^-1;
    37 md.inversion.vy_obs=constant/2*md.mesh.y.*(md.geometry.thickness).^-1;
    38 md.inversion.vel_obs=(sqrt((md.inversion.vx_obs).^2+(md.inversion.vy_obs).^2));
     36md.inversion.vx_obs=constant/2.*md.mesh.x.*(md.geometry.thickness).^-1;
     37md.inversion.vy_obs=constant/2.*md.mesh.y.*(md.geometry.thickness).^-1;
     38md.inversion.vel_obs=sqrt((md.inversion.vx_obs).^2+(md.inversion.vy_obs).^2);
    3939md.initialization.vx=zeros(md.mesh.numberofvertices,1);
    4040md.initialization.vy=zeros(md.mesh.numberofvertices,1);
    … …  
    4343
    4444%Deal with boundary conditions:
    45 disp('      boundary conditions for diagnostic model: ');
     45disp('      boundary conditions for diagnostic model:');
    4646md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp');
    4747
    48 radius=sqrt((md.mesh.x).*md.mesh.x+(md.mesh.y).*md.mesh.y);
     48radius=sqrt((md.mesh.x).^2+(md.mesh.y).^2);
    4949pos=find(radius==min(radius));
    50 md.mesh.x(pos)=0; md.mesh.y(pos)=0; %the closest node to the center is changed to be exactly at the center
     50md.mesh.x(pos)=0.; md.mesh.y(pos)=0.; %the closest node to the center is changed to be exactly at the center
    5151
    52 md.diagnostic.spcvx(pos)=0;
    53 md.diagnostic.spcvy(pos)=0;
    54 md.diagnostic.spcvz(pos)=0;
     52md.diagnostic.spcvx(pos)=0.;
     53md.diagnostic.spcvy(pos)=0.;
     54md.diagnostic.spcvz(pos)=0.;

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

    r13137 r14106  
    44ymin=min(md.mesh.y);
    55ymax=max(md.mesh.y);
    6 md.geometry.thickness=500*ones(md.mesh.numberofvertices,1);
     6md.geometry.thickness=500.*ones(md.mesh.numberofvertices,1);
    77md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness;
    88md.geometry.surface=md.geometry.bed+md.geometry.thickness;
    99
    1010disp('      creating drag');
    11 md.friction.coefficient=200*ones(md.mesh.numberofvertices,1); %q=1.
     11md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1); %q=1.
    1212%Take care of iceshelves: no basal drag
    1313pos=find(md.mask.elementonfloatingice);
    14 md.friction.coefficient(md.mesh.elements(pos,:))=0;
     14md.friction.coefficient(md.mesh.elements(pos,:))=0.;
    1515md.friction.p=ones(md.mesh.numberofelements,1);
    1616md.friction.q=ones(md.mesh.numberofelements,1);
    1717
    1818disp('      creating initial values');
    19 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
     19md.initialization.temperature=(273.-20.)*ones(md.mesh.numberofvertices,1);
    2020md.initialization.vx=zeros(md.mesh.numberofvertices,1);
    2121md.initialization.vy=zeros(md.mesh.numberofvertices,1);
    … …  
    2424md.initialization.pressure=zeros(md.mesh.numberofvertices,1);
    2525
    26 disp('      creating flow law paramter');
     26disp('      creating flow law parameter');
    2727md.materials.rheology_B=1.7687*10^8*ones(md.mesh.numberofvertices,1);
    28 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     28md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2929
    3030disp('      creating surface mass balance');
    3131md.surfaceforcings.mass_balance=0.2*ones(md.mesh.numberofvertices,1); %0m/a
    32 md.basalforcings.melting_rate=0*ones(md.mesh.numberofvertices,1); %0m/a
     32md.basalforcings.melting_rate=0.*ones(md.mesh.numberofvertices,1); %0m/a
    3333
    34 disp('      boundary conditions ');
     34disp('      boundary conditions');
    3535md=SetMarineIceSheetBC(md,'../Exp/SquareFrontEISMINT.exp');
    3636
    3737%Evolution of the ice shelf
    38 pos=find(md.mesh.y==200000); %nodes on the upper boundary condition
     38pos=find(md.mesh.y==200000.); %nodes on the upper boundary condition
    3939md.balancethickness.spcthickness=NaN*ones(md.mesh.numberofvertices,1);
    40 md.balancethickness.spcthickness(pos)=500;
     40md.balancethickness.spcthickness(pos)=500.;
    4141md.prognostic.spcthickness=NaN*ones(md.mesh.numberofvertices,1);
    42 md.prognostic.spcthickness(pos)=500;
     42md.prognostic.spcthickness(pos)=500.;
    4343md.prognostic.stabilization=0; %Better result with no artificial diffusivity
    44 md.thermal.stabilization=0; 
    45 md.timestepping.final_time=500;
     44md.thermal.stabilization=0;
     45md.timestepping.final_time=500.;
    4646md.timestepping.time_step=1;
Note: See TracChangeset for help on using the changeset viewer.