Changeset 17769


Ignore:
Timestamp:
04/18/14 10:00:36 (11 years ago)
Author:
cborstad
Message:

CHG: reverting previous change, isdamage=0 is specified by default in the damage class, problem was a typo elsewhere...

Location:
issm/trunk-jpl/test
Files:
41 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test271.m

    r17764 r17769  
    1515md.damage.stabilization=1;
    1616
    17 md.damage.requested_outputs={'default','DamageF'}
     17md.damage.requested_outputs={'default','DamageF'};
    1818
    1919md=setflowequation(md,'SSA','all');
    … …  
    2222
    2323%Fields and tolerances to track changes
    24 field_names={'D','F'}
     24field_names={'D','F'};
    2525field_tolerances={1e-13,1e-13};
    2626field_values={...

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    r17767 r17769  
    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
    24 md.damage.isdamage=0;
    2522
    2623disp('      creating surface mass balance');

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

    r17767 r17769  
    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
    27 md.damage.isdamage=0
    2825
    2926print "      creating surface mass balance"

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

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

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

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

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

    r17767 r17769  
    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
    28 md.damage.isdamage=0;
    2926
    3027disp('      creating surface mass balance');

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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