[aspect-devel] Questions on Boundary temperature model/Function

[Timothy Nelson]

Hi,

I am trying to figure out how to assign the function to a specific
boundary. Below is the current function. I would like to have it only
pertain to the bottom boundary. Thank you!


subsection Initial conditions
  set Model name = function

  subsection Function
    set Variable names      = x,z
    set Function constants  = p=0.01, L=3e6, pi=3.1415926536, k=1
    set Function expression = 1673*((1.0-z/3e6) - p*cos(k*pi*x/L)*sin(pi*z))
  end
end



subsection Boundary temperature model
  set Model name          = function

  subsection Function
   set Variable names      = x,y
   set Function constants  = Tcon=1673, dT=100, xstep=0.9, L=3e6, sm=0.02
   set Function expression =
(Tcon-(1-xstep)*dT)+dT*(0.5+(1/pi)*atan((x-xstep*L)/(sm*L)))

  end
end


-- 
Best,

Tim Nelson

*University of California Davis*
Undergraduate
tknelson at ucdavis.edu
(480) 559-4747
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20150220/261bd1ec/attachment.html>