Using collpase

in-woo

Dear developers,
I have a question with collapse function.

I want to solve the temperature field in 3D for ice rigidity, and changed 3D ice rigidity values into 2D.

   md = extrude(md,10,2);

   disp('   -- Solve 3D thermal model');
   md = solve(md,'th');

   disp('   -- Set ice viscosity as 3D thermal model ');
   md.initialization.temperature =md.results.ThermalSolution.Temperature;
   md.materials.rheology_B = cuffey(md.results.ThermalSolution.Temperature);

   disp('   -- Collapse 3D viscosity field into 2D');
   md = collapse(md);

Does collapse function provide the depth-average temperature field from 3D? I'm not sure that collapse function really supports the average for depth values. Or additional process is required to average the 3D values into depth average?

Thank you,
Inwoo

soroush

Hi Inwoo,
It returns the depth-averaged value of initial temperature,

if ~isnan(md.initialization.temperature),md.initialization.temperature=DepthAverage(md,md.initialization.temperature);
end;

But the temperature boundary conditions (spctemperature) are just projected to a 2d mesh by
md.thermal.spctemperature=project2d(md,md.thermal.spctemperature,md.mesh.numberoflayers);

You can also calculate the depth-averaged value of any variable of interest yourself using the same DepthAverage() function.

Check out the model.m file in your /trunk/bin/ directory, where function md = collapse(md) is declared (in case you already haven't). You can see the details of what collapse function does.

hope this somewhat helps,
Soroush