Changeset 9425


Ignore:
Timestamp:
08/25/11 15:18:21 (14 years ago)
Author:
Mathieu Morlighem
Message:

no more md.observed_temperature

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

Legend:

Unmodified
Added
Removed

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

    r9080 r9425  
    1111
    1212%Materials
    13 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    14 md.rheology_B=paterson(md.observed_temperature);
     13md.temperature=(273-20)*ones(md.numberofnodes,1);
     14md.rheology_B=paterson(md.temperature);
    1515md.rheology_n=3*ones(md.numberofelements,1);
    1616

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

    r8824 r9425  
    1515
    1616%Materials
    17 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    18 md.rheology_B=paterson(md.observed_temperature);
     17md.temperature=(273-20)*ones(md.numberofnodes,1);
     18md.rheology_B=paterson(md.temperature);
    1919md.rheology_n=3*ones(md.numberofelements,1);
    20 md.temperature=md.observed_temperature;
     20md.temperature=md.temperature;
    2121
    2222%Friction

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

    r8999 r9425  
    1515st=1.67*10^-2/1000; %k/m;
    1616radius=sqrt((md.x).^2+(md.y).^2);
    17 md.observed_temperature=(tmin+st*radius);
    18 md.temperature=md.observed_temperature;
     17md.temperature=(tmin+st*radius);
    1918md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
    2019

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

    r9389 r9425  
    4040
    4141%Materials
    42 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    43 md.rheology_B=paterson(md.observed_temperature);
     42md.temperature=(273-20)*ones(md.numberofnodes,1);
     43md.rheology_B=paterson(md.temperature);
    4444md.rheology_n=3*ones(md.numberofelements,1);
    45 md.temperature=md.observed_temperature;
    4645
    4746%Surface mass balance and basal melting

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

    r8999 r9425  
    2121tmin=238.15; %K
    2222st=1.67*10^-2/1000; %k/m;
    23 md.observed_temperature=(tmin+st*radius);
    24 md.temperature=md.observed_temperature;
     23md.temperature=(tmin+st*radius);
    2524md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
    2625

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

    r9196 r9425  
    1919
    2020disp('      creating temperature');
    21 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
     21md.temperature=(273-20)*ones(md.numberofnodes,1);
    2222
    2323disp('      creating flow law paramter');
    24 %md.B=paterson(md.observed_temperature);
    2524md.rheology_B=1.7687*10^8*ones(md.numberofnodes,1);
    2625md.rheology_n=3*ones(md.numberofelements,1);

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

    r8305 r9425  
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    22 md.rheology_B=paterson(md.observed_temperature);
     21md.temperature=(273-20)*ones(md.numberofnodes,1);
     22md.rheology_B=paterson(md.temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    24 md.temperature=md.observed_temperature;
    2524
    2625%Friction

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

    r8999 r9425  
    2222
    2323%Materials
    24 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    25 md.rheology_B=paterson(md.observed_temperature);
     24md.temperature=(273-20)*ones(md.numberofnodes,1);
     25md.rheology_B=paterson(md.temperature);
    2626md.rheology_n=3*ones(md.numberofelements,1);
    27 md.temperature=md.observed_temperature;
    2827
    2928%Accumulation and melting

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

    r8305 r9425  
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    22 md.rheology_B=paterson(md.observed_temperature);
     21md.temperature=(273-20)*ones(md.numberofnodes,1);
     22md.rheology_B=paterson(md.temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    24 md.temperature=md.observed_temperature;
    2524
    2625%Friction

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

    r8999 r9425  
    1919
    2020%Materials
    21 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    22 md.rheology_B=paterson(md.observed_temperature);
     21md.temperature=(273-20)*ones(md.numberofnodes,1);
     22md.rheology_B=paterson(md.temperature);
    2323md.rheology_n=3*ones(md.numberofelements,1);
    24 md.temperature=md.observed_temperature;
    2524
    2625%Surface mass balance and basal melting

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

    r9046 r9425  
    2525
    2626disp('      creating temperatures');
    27 md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
     27md.temperature=(273-20)*ones(md.numberofnodes,1);
    2828
    2929disp('      creating flow law paramter');
    30 md.rheology_B=paterson(md.observed_temperature);
     30md.rheology_B=paterson(md.temperature);
    3131md.rheology_n=3*ones(md.numberofelements,1);
    3232
    … …  
    4141
    4242disp('      boundary conditions for thermal model');
    43 md.temperature=md.observed_temperature;
    44 md.spctemperature(:)=md.observed_temperature;
     43md.spctemperature(:)=md.temperature;
    4544md.geothermalflux=zeros(md.numberofnodes,1);
    4645pos=find(md.elementonicesheet);md.geothermalflux(md.elements(pos,:))=1*10^-3; %1 mW/m^2
Note: See TracChangeset for help on using the changeset viewer.