[aspect-devel] update on memory pbs

[Wolfgang Bangerth]

On 02/15/2013 01:57 PM, Wolfgang Bangerth wrote:
> I've run your vkk.prm with global refinement 9 on my machine

In the interest of making things comparable, the input file is attached.

Best
  W.

-- 
------------------------------------------------------------------------
Wolfgang Bangerth               email:            bangerth at math.tamu.edu
                                 www: http://www.math.tamu.edu/~bangerth/

-------------- next part --------------
############### Global parameters
# We use a 3d setup. Since we are only interested
# in a steady state solution, we set the end time
# equal to the start time to force a single time
# step before the program terminates.

set Dimension                              = 2
set Start time                             = 0
set End time                               = 2000
set Use years in output instead of seconds = false
set CFL number                             = 0.5
set Output directory                       = output


subsection Geometry model
  set Model name = box
  subsection Box
    set X extent  = 0.9142
    set Y extent  = 1.0000
  end
end

subsection Model settings
  set Tangential velocity boundary indicators = 0,1
  set Zero velocity boundary indicators       = 2,3
end

subsection Material model
  set Model name = simple
  subsection Simple model
    set Reference density             = 1010
    set Viscosity                     = 1e2
    set Thermal expansion coefficient = 0
  end
end

subsection Gravity model
  set Model name = vertical
  subsection Vertical
    set Magnitude = 10
  end
end


############### Parameters describing the temperature field
# As above, there is no need to set anything for the
# temperature boundary conditions.

subsection Boundary temperature model
  set Model name = box
end

subsection Initial conditions
  set Model name = function
  subsection Function
    set Variable names      = x,z
    set Function constants  = pi=3.14159
    set Function expression = if( (z>0.2+0.02*cos(pi*x/0.9142)) , 0 , 1 )
  end
end


############### Parameters describing the compositional field

subsection Compositional fields
  set Number of fields = 1
end

subsection Compositional initial conditions
  set Model name = function
  subsection Function
    set Variable names      = x,z
    set Function constants  = pi=3.14159
    set Function expression = if( (z>0.2+0.02*cos(pi*x/0.9142)) , 0 , 1 )
  end
end

subsection Material model
  subsection Simple model
    set Density differential for compositional field 1 = -10
  end
end

############### Parameters describing the discretization

subsection Mesh refinement
  set Initial adaptive refinement        = 0
  set Strategy                           = composition
  set Initial global refinement          = 9
  set Time steps between mesh refinement = 0
  set Coarsening fraction                = 0.9
  set Refinement fraction                = 0.9
end

############### Parameters describing the what to do with the solution

subsection Postprocess
  set List of postprocessors = visualization, velocity statistics, composition statistics
  subsection Visualization
    set List of output variables = all
    set Output format                 = vtu
    set Time between graphical output = 25
  end
end