Changeset 9425
- Timestamp:
- 08/25/11 15:18:21 (14 years ago)
- Location:
- issm/trunk/test/Par
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Par/79North.par
r9080 r9425 11 11 12 12 %Materials 13 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);14 md.rheology_B=paterson(md. observed_temperature);13 md.temperature=(273-20)*ones(md.numberofnodes,1); 14 md.rheology_B=paterson(md.temperature); 15 15 md.rheology_n=3*ones(md.numberofelements,1); 16 16 -
issm/trunk/test/Par/Pig.par
r8824 r9425 15 15 16 16 %Materials 17 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);18 md.rheology_B=paterson(md. observed_temperature);17 md.temperature=(273-20)*ones(md.numberofnodes,1); 18 md.rheology_B=paterson(md.temperature); 19 19 md.rheology_n=3*ones(md.numberofelements,1); 20 md.temperature=md. observed_temperature;20 md.temperature=md.temperature; 21 21 22 22 %Friction -
issm/trunk/test/Par/RoundSheetEISMINT.par
r8999 r9425 15 15 st=1.67*10^-2/1000; %k/m; 16 16 radius=sqrt((md.x).^2+(md.y).^2); 17 md.observed_temperature=(tmin+st*radius); 18 md.temperature=md.observed_temperature; 17 md.temperature=(tmin+st*radius); 19 18 md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1); 20 19 -
issm/trunk/test/Par/RoundSheetShelf.par
r9389 r9425 40 40 41 41 %Materials 42 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);43 md.rheology_B=paterson(md. observed_temperature);42 md.temperature=(273-20)*ones(md.numberofnodes,1); 43 md.rheology_B=paterson(md.temperature); 44 44 md.rheology_n=3*ones(md.numberofelements,1); 45 md.temperature=md.observed_temperature;46 45 47 46 %Surface mass balance and basal melting -
issm/trunk/test/Par/RoundSheetStaticEISMINT.par
r8999 r9425 21 21 tmin=238.15; %K 22 22 st=1.67*10^-2/1000; %k/m; 23 md.observed_temperature=(tmin+st*radius); 24 md.temperature=md.observed_temperature; 23 md.temperature=(tmin+st*radius); 25 24 md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1); 26 25 -
issm/trunk/test/Par/SquareEISMINT.par
r9196 r9425 19 19 20 20 disp(' creating temperature'); 21 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);21 md.temperature=(273-20)*ones(md.numberofnodes,1); 22 22 23 23 disp(' creating flow law paramter'); 24 %md.B=paterson(md.observed_temperature);25 24 md.rheology_B=1.7687*10^8*ones(md.numberofnodes,1); 26 25 md.rheology_n=3*ones(md.numberofelements,1); -
issm/trunk/test/Par/SquareSheetConstrained.par
r8305 r9425 19 19 20 20 %Materials 21 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);22 md.rheology_B=paterson(md. observed_temperature);21 md.temperature=(273-20)*ones(md.numberofnodes,1); 22 md.rheology_B=paterson(md.temperature); 23 23 md.rheology_n=3*ones(md.numberofelements,1); 24 md.temperature=md.observed_temperature;25 24 26 25 %Friction -
issm/trunk/test/Par/SquareSheetShelf.par
r8999 r9425 22 22 23 23 %Materials 24 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);25 md.rheology_B=paterson(md. observed_temperature);24 md.temperature=(273-20)*ones(md.numberofnodes,1); 25 md.rheology_B=paterson(md.temperature); 26 26 md.rheology_n=3*ones(md.numberofelements,1); 27 md.temperature=md.observed_temperature;28 27 29 28 %Accumulation and melting -
issm/trunk/test/Par/SquareShelf.par
r8305 r9425 19 19 20 20 %Materials 21 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);22 md.rheology_B=paterson(md. observed_temperature);21 md.temperature=(273-20)*ones(md.numberofnodes,1); 22 md.rheology_B=paterson(md.temperature); 23 23 md.rheology_n=3*ones(md.numberofelements,1); 24 md.temperature=md.observed_temperature;25 24 26 25 %Friction -
issm/trunk/test/Par/SquareShelfConstrained.par
r8999 r9425 19 19 20 20 %Materials 21 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);22 md.rheology_B=paterson(md. observed_temperature);21 md.temperature=(273-20)*ones(md.numberofnodes,1); 22 md.rheology_B=paterson(md.temperature); 23 23 md.rheology_n=3*ones(md.numberofelements,1); 24 md.temperature=md.observed_temperature;25 24 26 25 %Surface mass balance and basal melting -
issm/trunk/test/Par/SquareThermal.par
r9046 r9425 25 25 26 26 disp(' creating temperatures'); 27 md. observed_temperature=(273-20)*ones(md.numberofnodes,1);27 md.temperature=(273-20)*ones(md.numberofnodes,1); 28 28 29 29 disp(' creating flow law paramter'); 30 md.rheology_B=paterson(md. observed_temperature);30 md.rheology_B=paterson(md.temperature); 31 31 md.rheology_n=3*ones(md.numberofelements,1); 32 32 … … 41 41 42 42 disp(' boundary conditions for thermal model'); 43 md.temperature=md.observed_temperature; 44 md.spctemperature(:)=md.observed_temperature; 43 md.spctemperature(:)=md.temperature; 45 44 md.geothermalflux=zeros(md.numberofnodes,1); 46 45 pos=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.