|
From: Aaron D. <aar...@gm...> - 2015-03-05 05:15:33
|
Hi,
You're trying to make an early pre-main sequence star rotate at 100 km/s at
the surface. This is going to be very hard to converge, even for MESA.
Perhaps you should try the options which let you start rotation near the
ZAMS? As Bill suggests:
! You might want to start a run at pre-ms but only turn on
rotation when near zams
! rather than force you to stop the run near zams, change the
inlist, and restart.
! The following will turn on rotation automatically.
! The working definition of "near zams" is
! `L_nuc_burn_total/L_phot >= Lnuc_div_L_upper_limit`
! `Lnuc_div_L_upper_limit` is in the controls part of the inlist.
See mesa/star/defaults/star_job.defaults for more near_zams options.
Aaron
On Thu, Mar 5, 2015 at 4:06 PM, 王江涛 <wan...@16...> wrote:
>
>
> hi!
> I was in mesa/star/work/inlist_project set
> star_job
> ..........................
> new_rotation_flag = .true.
> change_rotation_flag = .false.
> change_initial_rotation_flag = .true.
>
>
> new_surface_rotation_v = 100
> set_surface_rotation_v = .false.
> set_initial_surface_rotation_v = .true
> ................................................
> And
> star_controls
> !# rotation controls
> .............
> but it is show
> read history_columns.list
> read profile_columns.list
> net name basic.net
> rotation_flag = T
> kappa_file_prefix gn93
> kappa_lowT_prefix lowT_fa05_gs98
> eos_file_prefix mesa
> OMP_NUM_THREADS 4
>
>
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 1 1
> -5.602060 0.500000
> 1st model retry log10(dt/yr)
> -5.3010299956639804D+00
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 2 1
> -5.903090 0.500000
> 1st model retry log10(dt/yr)
> -5.6020599913279616D+00
> have reached retry limit so now backup
> 1st model retry log10(dt/yr)
> -7.6020599913279616D+00
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 1 1
> -8.204120 0.500000
> 1st model retry log10(dt/yr)
> -7.9030899869919429D+00
> first model is slow to converge: num tries 70
> first model is slow to converge: num tries 80
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 2 1
> -8.505150 0.500000
> 1st model retry log10(dt/yr)
> -8.2041199826559232D+00
> first model is slow to converge: num tries 90
> first model is slow to converge: num tries 100
> first model is slow to converge: num tries 110
> have reached retry limit so now backup
> 1st model retry log10(dt/yr)
> -1.0806179973983886D+01
> first model is slow to converge: num tries 120
> first model is slow to converge: num tries 130
> first model is slow to converge: num tries 140
> first model is slow to converge: num tries 150
> first model is slow to converge: num tries 160
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 1 1
> -11.408240 0.500000
> 1st model retry log10(dt/yr)
> -1.1107209969647867D+01
> first model is slow to converge: num tries 170
> first model is slow to converge: num tries 180
> first model is slow to converge: num tries 190
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 2 1
> -11.709270 0.500000
> 1st model retry log10(dt/yr)
> -1.1408239965311848D+01
> first model is slow to converge: num tries 200
> first model is slow to converge: num tries 210
> first model is slow to converge: num tries 220
> first model is slow to converge: num tries 230
> first model is slow to converge: num tries 240
> have reached retry limit so now backup
> 1st model retry log10(dt/yr)
> -1.4612359947967773D+01
> first model is slow to converge: num tries 250
> first model is slow to converge: num tries 260
> first model is slow to converge: num tries 270
> first model is slow to converge: num tries 280
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 1 1
> -15.214420 0.500000
> 1st model retry log10(dt/yr)
> -1.4913389943631755D+01
> first model is slow to converge: num tries 290
> first model is slow to converge: num tries 300
> first model is slow to converge: num tries 310
> first model is slow to converge: num tries 320
> first model is slow to converge: num tries 330
> first model is slow to converge: num tries 340
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 2 1
> -15.515450 0.500000
> 1st model retry log10(dt/yr)
> -1.5214419939295736D+01
> first model is slow to converge: num tries 350
> first model is slow to converge: num tries 360
> first model is slow to converge: num tries 370
> first model is slow to converge: num tries 380
> first model is slow to converge: num tries 390
> have reached retry limit so now backup
> 1st model retry log10(dt/yr)
> -1.9020599913279622D+01
> first model is slow to converge: num tries 400
> first model is slow to converge: num tries 410
> first model is slow to converge: num tries 420
> first model is slow to converge: num tries 430
> first model is slow to converge: num tries 440
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 1 1
> -19.622660 0.500000
> 1st model retry log10(dt/yr)
> -1.9321629908943603D+01
> first model is slow to converge: num tries 450
> first model is slow to converge: num tries 460
> first model is slow to converge: num tries 470
> first model is slow to converge: num tries 480
> first model is slow to converge: num tries 490
> hydro_newton_step returned ierr
> retry cnt, step, log10(dt/yr), retry_factor 2 1
> -19.923690 0.500000
> 1st model retry log10(dt/yr)
> -1.9622659904607584D+01
> first model is slow to converge: num tries 500
> first model is slow to converge: num tries 510
> first model is slow to converge: num tries 520
> first model is slow to converge: num tries 530
> first model is slow to converge: num tries 540
> first model is slow to converge: num tries 550
> have reached retry limit so now backup
> dt
> 2.9390817493200311D-17
> min_timestep_limit
> 9.9999999999999998D-17
>
> stopping because of convergence problems dt < min_timestep_limit
>
>
> terminated evolution: convergence problems
>
> termination code: min_timestep_limit
> DATE: 2015-03-03
> TIME: 12:46:54
> [hfsong@localhost wangjt_star_work]$
> I hope you help me!
>
> wang J.T
>
>
>
>
>
>
> ------------------------------------------------------------------------------
> Dive into the World of Parallel Programming The Go Parallel Website,
> sponsored
> by Intel and developed in partnership with Slashdot Media, is your hub for
> all
> things parallel software development, from weekly thought leadership blogs
> to
> news, videos, case studies, tutorials and more. Take a look and join the
> conversation now. http://goparallel.sourceforge.net/
> _______________________________________________
> mesa-users mailing list
> mes...@li...
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
|