Re: Problem with source.f

From: Francesco Cerutti (
Date: Wed Nov 01 2006 - 10:18:07 CET

  • Next message: Joseph Comfort: "Re: Problem with source.f"

    Hi Joseph

    I think that in fact your problem is not concerning your source, but your
    geometry instead. If you try to debug it just around the point where you
    place your primary particle, using for instance

    GEOEND 0.0 0.0 10.0 -0.5 -0.1 5.0DEBUG

    (instead of GEOEND alone -> see the manual)

    you can immediately see that there is a huge amount of points not
    belonging to any region.
    The reason is that the operator "+" means intersection and not union. So
    the region R2air (where I guess that your primary particle was supposed to
    start) is in fact empty, due to the intersection of bodies not overlapping
    at all.

    (Moreover, following the recent suggestion of Alberto, in the
    source.f statement


    replace 1.0 by 1.0d0 or oneone (double precision required)).



    Francesco Cerutti
    CH-1211 Geneva 23
    tel. ++41 22 7678962

  • Next message: Joseph Comfort: "Re: Problem with source.f"

    This archive was generated by hypermail 2.1.6 : Thu Nov 02 2006 - 10:10:56 CET