[mesa-users] Proper use of "other_winds"

[Chris Mankovich]

Hi Mathieu,

Your star_info pointer has been initialized but hasn't yet been mapped to
your current star_info structure.  Your subroutine takes the star id as
input, so you just need to use get_star_ptr from star_def:

call get_star_ptr(id,s,ierr)

then you can access anything you need.

I can't speak to whether there's a better way to call the wind routine you
want; I don't see a public interface for it.

Cheers,
Chris



On Mon, May 26, 2014 at 11:42 AM, Mathieu <mathren90 at gmail.com> wrote:

> Hi everyone,
>
> I am trying to implement a wind scheme in MESA 6208 using the "other_wind"
> hook (see the attached run_star_extras).
> To compile it, I have copied  $MESA_DIR/star/make/mod_kuma.mod into
> $MESA_DIR/include, to be able to call the routine "eval_Kudritzki_wind"
> (defined in /$MESA_DIR/star/private/winds.f )
> I would like to use the surface hydrogen mass fraction, so I tried to
> access this variable with:
>
>     type (star_info), pointer :: s
>     real(dp) :: xsurf
>     xsurf = s% xa(2,1)
>
> but when I run an initial model I get:
>
> Fortran runtime error: Index '1' of dimension 2 of array 's%xa' below
> lower bound of 47946330539336
>
> which looks like s% xa is not initialized. Can someone shed some light on
> how should I access to the surface H mass fraction?
> And maybe there is a better way to call the routine "eval_Kudritzki_wind"
> than copying a module in the $MESA_DIR/include directory?
>
> Eventually I would also like to use the x_cntrl variable to set the
> efficiency of my wind from the inlist_project.
>
> Thanks in advance!
>
> Mathieu
>
>
>
>
>
>
>
> ------------------------------------------------------------------------------
> The best possible search technologies are now affordable for all companies.
> Download your FREE open source Enterprise Search Engine today!
> Our experts will assist you in its installation for $59/mo, no commitment.
> Test it for FREE on our Cloud platform anytime!
>
> http://pubads.g.doubleclick.net/gampad/clk?id=145328191&iu=/4140/ostg.clktrk
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140526/d4388d99/attachment.html>