Revision dd8daa98413988b11c594e5a3bff8f5659d7ecf7 authored by Florent Renaud on 26 June 2014, 09:41:19 UTC, committed by Florent Renaud on 26 June 2014, 09:41:19 UTC
Originally released on 18 March 2013 Based on Nbody6 downloaded on 29 January 2013
1 parent edab1b7
trdot2.f
***
SUBROUTINE TRDOT2(KW,AGE,TM,TN,TSCLS,DTM,DTR)
*
*
* Time-scale for evolution changes.
* ---------------------------------
*
REAL*8 AGE,TM,TN,TSCLS(20),DTM,DTR
REAL*8 PTS1,PTS2
PARAMETER(PTS1=0.05D0,PTS2=0.02D0)
*
* Base new time scale for changes in radius & mass on stellar type.
if(kw.le.1)then
dtm = pts1*tm
dtr = tm - age
elseif(kw.ge.10)then
dtm = 1.0d+02
dtr = dtm
elseif(kw.eq.2)then
dtm = pts1*(tscls(1) - tm)
dtr = tscls(1) - age
elseif(kw.eq.3)then
if(age.lt.tscls(6))then
dtm = pts2*(tscls(4) - age)
else
dtm = pts2*(tscls(5) - age)
endif
dtr = MIN(tscls(2),tn) - age
elseif(kw.eq.4)then
dtm = pts1*tscls(3)
dtr = MIN(tn,tscls(2) + tscls(3)) - age
elseif(kw.eq.5)then
if(age.lt.tscls(9))then
dtm = pts2*(tscls(7) - age)
else
dtm = pts2*(tscls(8) - age)
endif
dtr = MIN(tn,tscls(13)) - age
elseif(kw.eq.6)then
if(age.lt.tscls(12))then
dtm = pts2*(tscls(10) - age)
else
dtm = pts2*(tscls(11) - age)
endif
dtm = MIN(dtm,0.005d0)
dtr = tn - age
elseif(kw.eq.7)then
dtm = pts1*tm
dtr = tm - age
elseif(kw.eq.8.or.kw.eq.9)then
if(age.lt.tscls(6))then
dtm = pts2*(tscls(4) - age)
else
dtm = pts2*(tscls(5) - age)
endif
dtr = tn - age
endif
dtm = MIN(dtm,dtr)
dtm = MAX(dtm,1.0d-07)
*
RETURN
END
***
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...