Changeset 17767


Ignore:
Timestamp:
04/18/14 09:46:20 (11 years ago)
Author:
cborstad
Message:

BUG: forgot to set isdamage=0 in par files

Location:
issm/trunk-jpl/test/Par
Files:
40 edited

Legend:

Unmodified
Added
Removed

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

    r17762 r17767  
    2020md.materials.rheology_B=paterson(md.initialization.temperature);
    2121md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     22
     23%Damage
     24md.damage.isdamage=0;
    2225
    2326%Friction

  • issm/trunk-jpl/test/Par/79North.py

    r17762 r17767  
    3333md.initialization.temperature=md.initialization.temperature
    3434
     35#Damage
     36md.damage.isdamage=0
     37
    3538#Friction
    3639md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    4444md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    4545
     46%Damage
     47md.damage.isdamage=0;
     48
    4649%Friction
    4750md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);

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

    r17762 r17767  
    4343md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    4444
     45%Damage
     46md.damage.isdamage=0;
     47
    4548%Friction
    4649md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);

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

    r17762 r17767  
    1616md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    1717
     18%Damage
     19md.damage.isdamage=0;
     20
    1821disp('      boundary conditions for stressbalance model');
    1922%Create node on boundary first (because we cannot use mesh)

  • issm/trunk-jpl/test/Par/ISMIPA.py

    r17762 r17767  
    1919md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    2020
     21#Damage
     22md.damage.isdamage=0
     23
    2124print "      boundary conditions for stressbalance model"
    2225#Create node on boundary first (because we cannot use mesh)

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

    r17762 r17767  
    1616md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    1717
     18%Damage
     19md.damage.isdamage=0;
     20
    1821disp('      boundary conditions for stressbalance model');
    1922%Create node on boundary first (because we cannot use mesh)

  • issm/trunk-jpl/test/Par/ISMIPB.py

    r17762 r17767  
    1919md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    2020
     21#Damage
     22md.damage.isdamage=0
     23
    2124print "      boundary conditions for stressbalance model"
    2225#Create node on boundary first (because we cannot use mesh)

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

    r17762 r17767  
    1717md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    1818
     19%Damage
     20md.damage.isdamage=0;
     21
    1922disp('      boundary conditions for stressbalance model:');
    2023%Create node on boundary first (because we can not use mesh)

  • issm/trunk-jpl/test/Par/ISMIPC.py

    r17762 r17767  
    2020md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    2121
     22#Damage
     23md.damage.isdamage=0
     24
    2225print "      boundary conditions for stressbalance model:"
    2326#Create node on boundary first (because we can not use mesh)

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

    r17762 r17767  
    1616md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    1717
     18%Damage
     19md.damage.isdamage=0;
     20
    1821disp('      boundary conditions for stressbalance model:');
    1922%Create node on boundary first (because we can not use mesh)

  • issm/trunk-jpl/test/Par/ISMIPD.py

    r17762 r17767  
    1919md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    2020
     21#Damage
     22md.damage.isdamage=0
     23
    2124print "      boundary conditions for stressbalance model:"
    2225#Create node on boundary first (because we can not use mesh)

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

    r17762 r17767  
    2626md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2727
     28%Damage
     29md.damage.isdamage=0;
     30
    2831disp('      boundary conditions for stressbalance model:');
    2932%Create node on boundary first (because we can not use mesh)

  • issm/trunk-jpl/test/Par/ISMIPE.py

    r17762 r17767  
    3131md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    3232
     33#Damage
     34md.damage.isdamage=0
     35
    3336print "      boundary conditions for stressbalance model:"
    3437#Create node on boundary first (because we can not use mesh)

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

    r17762 r17767  
    1717md.materials.rheology_n=1.*ones(md.mesh.numberofelements,1);
    1818md.materials.rheology_law='None';
     19
     20%Damage
     21md.damage.isdamage=0;
    1922
    2023disp('      boundary conditions for stressbalance model');

  • issm/trunk-jpl/test/Par/ISMIPF.py

    r17762 r17767  
    2020md.materials.rheology_n=1.*numpy.ones((md.mesh.numberofelements,1))
    2121md.materials.rheology_law='None'
     22
     23#Damage
     24md.damage.isdamage=0
    2225
    2326print "      boundary conditions for stressbalance model"

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

    r17762 r17767  
    2626md.initialization.temperature=md.initialization.temperature;
    2727
     28%Damage
     29md.damage.isdamage=0;
     30
    2831%Friction
    2932md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/Pig.py

    r17762 r17767  
    3737md.initialization.temperature=md.initialization.temperature
    3838
     39#Damage
     40md.damage.isdamage=0
     41
    3942#Friction
    4043md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2020md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
    2121md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
     22
     23%Damage
     24md.damage.isdamage=0;
    2225
    2326disp('      creating surface mass balance');

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

    r17762 r17767  
    2323md.materials.rheology_B=6.81*10**7*numpy.ones((md.mesh.numberofvertices,1))    #to have the same B as the analytical solution
    2424md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
     25
     26#Damage
     27md.damage.isdamage=0
    2528
    2629print "      creating surface mass balance"

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

    r17762 r17767  
    6161md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    6262
     63%Damage
     64md.damage.isdamage=0;
     65
    6366%Surface mass balance and basal melting
    6467md.surfaceforcings.mass_balance=-10.*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/RoundSheetShelf.py

    r17762 r17767  
    6868md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    6969
     70#Damage
     71md.damage.isdamage=0
     72
    7073#Surface mass balance and basal melting
    7174md.surfaceforcings.mass_balance=-10.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2424md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution
    2525md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
     26
     27%Damage
     28md.damage.isdamage=0;
    2629
    2730disp('      creating surface mass balance');

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

    r17762 r17767  
    2828md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    2929
     30#Damage
     31md.damage.isdamage=0
     32
    3033print "      creating surface mass balance"
    3134smb_max=0.5    #m/yr

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

    r17762 r17767  
    2626md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2727
     28%Damage
     29md.damage.isdamage=0;
     30
    2831disp('      creating surface mass balance');
    2932md.surfaceforcings.mass_balance=0.2*ones(md.mesh.numberofvertices,1); %0m/a

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

    r17762 r17767  
    2929md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    3030
     31#Damage
     32md.damage.isdamage=0
     33
    3134print "      creating surface mass balance"
    3235md.surfaceforcings.mass_balance=0.2*numpy.ones((md.mesh.numberofvertices,1))    #0m/a

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

    r17762 r17767  
    2929md.materials.rheology_B=paterson(md.initialization.temperature);
    3030md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
     31
     32%Damage
     33md.damage.isdamage=0;
    3134
    3235md.masstransport.calvingrate=zeros(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/SquareSheetConstrained.py

    r17762 r17767  
    4040md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    4141
     42#Damage
     43md.damage.isdamage=0
     44
    4245#Masstransport
    4346md.masstransport.calvingrate=0.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    3232md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    3333
     34%Damage
     35md.damage.isdamage=0;
     36
    3437%Accumulation and melting
    3538md.surfaceforcings.mass_balance=10.*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/SquareSheetShelf.py

    r17762 r17767  
    4343md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    4444
     45#Damage
     46md.damage.isdamage=0
     47
    4548#Accumulation and melting
    4649md.surfaceforcings.mass_balance=10.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2929md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
    3030
     31%Damage
     32md.damage.isdamage=0;
     33
    3134%Friction
    3235md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/SquareShelf.py

    r17762 r17767  
    5959md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1))
    6060
     61#Damage
     62md.damage.isdamage=0
     63
    6164#Friction
    6265md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2929md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
    3030
     31%Damage
     32md.damage.isdamage=0;
     33
    3134%Friction
    3235md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/SquareShelf2.py

    r17762 r17767  
    5959md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1))
    6060
     61#Damage
     62md.damage.isdamage=0
     63
    6164#Friction
    6265md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2929md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
    3030
     31%Damage
     32md.damage.isdamage=0;
     33
    3134%Surface mass balance and basal melting
    3235md.surfaceforcings.mass_balance=10*ones(md.mesh.numberofvertices,1);

  • issm/trunk-jpl/test/Par/SquareShelfConstrained.py

    r17762 r17767  
    4040md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    4141
     42#Damage
     43md.damage.isdamage=0
     44
    4245#Surface mass balance and basal melting
    4346md.surfaceforcings.mass_balance=10.*numpy.ones((md.mesh.numberofvertices,1))

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

    r17762 r17767  
    2727md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
    2828
     29%Damage
     30md.damage.isdamage=0;
     31
    2932disp('      creating surface mass balance');
    3033md.surfaceforcings.mass_balance=ones(md.mesh.numberofvertices,1)/md.constants.yts; %1m/a

  • issm/trunk-jpl/test/Par/SquareThermal.py

    r17762 r17767  
    3131md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
    3232
     33#Damage
     34md.damage.isdamage=0
     35
    3336print "      creating surface mass balance"
    3437md.surfaceforcings.mass_balance=numpy.ones((md.mesh.numberofvertices,1))/md.constants.yts    #1m/a

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

    r17762 r17767  
    4040md.materials.rheology_B = paterson(md.initialization.temperature);
    4141md.materials.rheology_n = 3.*ones(md.mesh.numberofelements,1);
     42
     43%Damage
     44md.damage.isdamage=0;
    4245
    4346%Thermal

  • issm/trunk-jpl/test/Par/ValleyGlacierShelf.py

    r17762 r17767  
    4747md.materials.rheology_B=paterson(md.initialization.temperature)
    4848md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
     49
     50#Damage
     51md.damage.isdamage=0
    4952
    5053#Thermal
Note: See TracChangeset for help on using the changeset viewer.