[mesa-users] Turning off neutrino types

[Bill Paxton]

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