Changeset 9481


Ignore:
Timestamp:
08/28/11 17:45:50 (14 years ago)
Author:
Mathieu Morlighem
Message:

no more firn property in model

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

Legend:

Unmodified
Added
Removed

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

    r9455 r9481  
    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.numberofnodes,1);
    87
    98disp('      creating drag');

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

    r9455 r9481  
    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.numberofnodes,1);
    87
    98disp('      creating drag');

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

    r9455 r9481  
    55md.bed=md.surface-1000;
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofnodes,1);
    87
    98disp('      creating drag');

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

    r9455 r9481  
    55md.bed=md.surface-1000;
    66md.thickness=md.surface-md.bed;
    7 md.firn_layer=0*ones(md.numberofnodes,1);
    87
    98disp('      creating drag');

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

    r9455 r9481  
    1515end
    1616md.thickness=md.surface-md.bed;
    17 md.firn_layer=0*ones(md.numberofnodes,1);
    1817md.thickness(find(~md.thickness))=0.01;
    1918md.bed=md.surface-md.thickness;

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

    r9465 r9481  
    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.numberofnodes,1);
    109
    1110disp('      creating drag');

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

    r9455 r9481  
    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.numberofnodes,1);
    87md.bed=0*md.thickness;
    98md.surface=md.bed+md.thickness;

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

    r9455 r9481  
    55ymax=max(md.y);
    66md.thickness=500*ones(md.numberofnodes,1);
    7 md.firn_layer=0*ones(md.numberofnodes,1);
    87md.bed=-md.rho_ice/md.rho_water*md.thickness;
    98md.surface=md.bed+md.thickness;

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

    r9455 r9481  
    66h=1000;
    77md.thickness=h*ones(md.numberofnodes,1);
    8 md.firn_layer=10*ones(md.numberofnodes,1);
    98md.bed=-1000*ones(md.numberofnodes,1);
    109md.surface=md.bed+md.thickness;
Note: See TracChangeset for help on using the changeset viewer.