Re: [FLASH-USERS] selfgravity 2d spherical?

The 2d spherical geometry is supported by the multipole solver. Please
send us your setup, init_block, Makefile.h and logfile so we can take a
look at the problem.

On Thu, 5 May 2005, Matthias Rempel wrote:

> Hi,
>
> I tried to setup a selfgravity problem in 2d spherical geometry
> and it didn't work. The first question is, if this is supported
> in Flash2.5 (according to the manual the multipole solver should
> be able to handle this geometry). I put a few write statements
> into GetGravAccelOneRow to write out the potential and gravitational
> acceleration. For most blocks the values are OK , but there are a few
> blocks that look like the following:
>
> gpot: 6.0139153975266924E-154 6.0134700169990651E-154
> 6.0134700169990651E-154 6.0134700169990651E-154
> 6.0134700169990651E-154
> 6.0134700169990651E-154 6.0134700169990651E-154
> 6.0134700169990651E-154
> 6.0134700169990651E-154 6.0134700169990651E-154
> 7.4570280168308933E-090
> 5.0397108423046450E+223 6.0134700169991797E-154
> 6.0134700169990651E-154
> 6.0134700169990651E-154 6.0134700169990651E-154
>
> grav: 3.5630442209916460E-174 3.5630442209916460E-174
> 0.0000000000000000 0.0000000000000000
> 0.0000000000000000
> 0.0000000000000000 0.0000000000000000
> 0.0000000000000000
> 0.0000000000000000 -5.9656224134647147E-106
> -4.0317686738437160E+207
> 5.9656224134647147E-106 4.0317686738437160E+207
> 8.2671729307096520E-184
> 0.0000000000000000 0.0000000000000000
>
> Matthias
>

-- 
Anshu Dubey
Astronomy & Astrophysics
University of Chicago
5640 S. Ellis Ave.  Chicago IL 60637
Tel : (773) 834 2999
Fax : (773) 834 3230