[mesa-users] Turning off neutrino types

[Shane Frewen]

Perfect! The routine runs fine when I define my own flags. Thank you both
for your help.

- Shane


On Wed, Mar 20, 2013 at 6:56 PM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:

>
> On Mar 20, 2013, at 6:50 PM, Josiah Schwab wrote:
>
> >
> > Ah, yes.  I did not realize that flags is intent(in) in other_neu!
> >
> > Bill, is there a technical reason that the other_neu routine shouldn't
> > modify flags?
>
> no -- that's my mistake.  I'll remove that restriction for the next
> release.
>
> In the meantime, declare logical :: my_flags(num_neu_types) in the
> other_neu routine,
> then
>         my_flags = flags
>         my_flags(pair_neu_type) = .false.
>         call neu_get( .....   my_flags ..... )
>
> no need to touch the mesa sources.
>
>
> -Bill
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130320/512a1d17/attachment.html>