[mesa-users] MLT problem

[Michael Zingale]

Fortran is not case-sensitive.  In fact, the original Fortran had
everything in all caps.

On Tue, Jan 19, 2016 at 12:39 PM, Kenny Van <kvan at ualberta.ca> wrote:

> Thanks Bill,
>
> One thing I've also noticed going through the code is that there are
> slightly different rho values used throughout the code. For example near
> the beginning of "set_convective_mixing", there is the following
>
> x = Q*Rho / (2*P)
>
> d_x_dvb = x * (drho_dvb / rho  + dQ_dvb/Q - dP_dvb/P)
>
> From my understanding fortran is a keysensitive language, so are these two
> "rho" values equivalent?
>
> On Tue, Jan 19, 2016 at 10:19 AM, Bill Paxton <paxton at kitp.ucsb.edu>
> wrote:
>
>>
>> On Jan 19, 2016, at 8:47 AM, Kenny Van wrote:
>>
>> I'm not sure I'm properly implementing the ierr = -1. If I understand
>> correctly, I would not only need to include ierr = -1 in the test of
>> (is_bad_num(d_Bcubed_dvb(mlt_dgrada))), but I would also need to define it
>> in the initial parameters of the subroutine correct?
>>
>>
>> ierr /= 0 is used by "lower" level routines to indicate a failure.   it
>> must be an arg to the routine so it can be returned to tell the caller what
>> has happened.
>>
>> "higher" level routine go beyond a simple binary okay or failed to return
>> an integer value from the following set (defined in
>> star/public/star_def.inc)
>>
>>       ! result codes for various evolution routines
>>       integer, parameter :: keep_going = 0
>>       integer, parameter :: redo = 1
>>          ! repeat current step with same timestep
>>          ! for example, after changing a parameter such as mdot
>>       integer, parameter :: retry = 2
>>          ! repeat current step with smaller timestep
>>       integer, parameter :: backup = 3
>>          ! go back to previous model and retry with smaller timestep
>>       integer, parameter :: terminate = 4
>>
>> Routines that are returning such a value may declare ierr as a local
>> variable for use in calls to lower level routines.
>>
>> For an example of this, see 'integer function do_evolve_step_part1' in
>> star/private/evolve
>>
>> -Bill
>>
>>
>>
>>
>
>
> ------------------------------------------------------------------------------
> Site24x7 APM Insight: Get Deep Visibility into Application Performance
> APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
> Monitor end-to-end web transactions and take corrective actions now
> Troubleshoot faster and improve end-user experience. Signup Now!
> http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>


-- 
Michael Zingale
Associate Professor

Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY
11794-3800
*phone*:  631-632-8225
*e-mail*: Michael.Zingale at stonybrook.edu
*web*: http://www.astro.sunysb.edu/mzingale
github: http://github.com/zingale
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160119/fcd0147e/attachment.html>