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

[Richard Townsend]

On a meta-note: please be sure to post questions first to the mailing list rather than to Bill --- otherwise, he'll spend all his time forwarding emails instead of implementing cool new stuff in MESA!

cheers,

Rich

On Aug 15, 2013, at 1:14 PM, Max Katz <maxpkatz at gmail.com> wrote:

> 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
> 
> 
> ------------------------------------------------------------------------------
> 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/29418cc1/attachment.html>