[mesa-users] Fwd: quick question on definitions in run_star_extras.f

[Max Katz]

Manos,

If you look at the top of run_star_extras you'll notice that it includes
the const_def module, which can be found in mesa/const/public. You can see
all of the constants that MESA includes there, which does include the ones
you want. Since we're just including the module directory, you don't need
any special s% notation, you can just write 'pi' and it will work.


On Thu, Aug 15, 2013 at 2:01 PM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:

> Can someone help Manos?
>
> thanks,
> bill
>
>
>
> Begin forwarded message:
>
> *From: *Manos Chatzopoulos <manolis at astro.as.utexas.edu>
> *Date: *August 15, 2013 10:54:04 AM PDT
> *To: *Bill Paxton <paxton at kitp.ucsb.edu>
> *Subject: **quick question on definitions in run_star_extras.f*
>
> Hi Bill,
>
> I apologize to bother you, I am adding my own column/variable in
> run_star_extras.f and it looks like
> this:
>
>      subroutine data_for_extra_profile_columns(s, id, id_extra, n, nz,
> names, vals, ierr)
>         type (star_info), pointer :: s
>         integer, intent(in) :: id, id_extra, n, nz
>         character (len=maxlen_profile_column_name) :: names(n)
>         real(dp) :: vals(nz,n),lnLambda
>         integer, intent(out) :: ierr
>         integer :: k
>         ierr = 0
>
>     names(1) = 'mri_criterion'
>
>     do k = 2,nz
>        lnLambda = log(2/(3*qe**3)*
>     +      sqrt(mp*(s% abar(k))*(kerg*(s% T(k)))**3/(pi*(s% rho(k))*
>     +      (max(1d-10,s% zbar(k)))**5)))
>
>        vals(k,1) = 7.0d11*lnLambda/(s% T(k))**1.5d0
>
>     enddo
>
>     vals(1,1) = 0.
>
> Do I have to re-define the constants (qe, pi, kerg, mp) or the code will
> recognize them? Or there
> is a special way I have to code them in the equation in the do-loop above?
> Do I also have to use
> the s% idiom like I do with the variables? (for example s% kerg) or things
> will work the way they are?
>
> Thanks!
>
> Manos
>
> --
> Dr. Manos Chatzopoulos
> Department of Astronomy
> University of Texas at Austin
>
>
>
>
> ------------------------------------------------------------------------------
> Get 100% visibility into Java/.NET code with AppDynamics Lite!
> It's a free troubleshooting tool designed for production.
> Get down to code-level detail for bottlenecks, with <2% overhead.
> Download for free and get started troubleshooting in minutes.
> http://pubads.g.doubleclick.net/gampad/clk?id=48897031&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/20130815/5dafed0d/attachment.html>