Menu
▾
▴
mesa-users — open list for mesa comments, questions, & answers
This list is closed, nobody may subscribe to it.
| 2007 |
Jan
|
Feb
|
Mar
|
Apr
|
May
(2) |
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
|
Dec
|
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 2008 |
Jan
|
Feb
|
Mar
|
Apr
|
May
(6) |
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
(8) |
Dec
(5) |
| 2009 |
Jan
(5) |
Feb
(1) |
Mar
(3) |
Apr
(4) |
May
(4) |
Jun
|
Jul
|
Aug
|
Sep
(8) |
Oct
|
Nov
(1) |
Dec
(13) |
| 2010 |
Jan
|
Feb
(6) |
Mar
(4) |
Apr
(1) |
May
(10) |
Jun
(43) |
Jul
(37) |
Aug
(3) |
Sep
(6) |
Oct
(26) |
Nov
(17) |
Dec
(29) |
| 2011 |
Jan
(28) |
Feb
(18) |
Mar
(42) |
Apr
(18) |
May
(13) |
Jun
(32) |
Jul
(32) |
Aug
(25) |
Sep
(46) |
Oct
(41) |
Nov
(36) |
Dec
(43) |
| 2012 |
Jan
(92) |
Feb
(120) |
Mar
(40) |
Apr
(75) |
May
(40) |
Jun
(93) |
Jul
(115) |
Aug
(67) |
Sep
(38) |
Oct
(92) |
Nov
(95) |
Dec
(47) |
| 2013 |
Jan
(171) |
Feb
(200) |
Mar
(100) |
Apr
(134) |
May
(112) |
Jun
(142) |
Jul
(123) |
Aug
(66) |
Sep
(175) |
Oct
(236) |
Nov
(141) |
Dec
(98) |
| 2014 |
Jan
(91) |
Feb
(88) |
Mar
(126) |
Apr
(63) |
May
(123) |
Jun
(122) |
Jul
(105) |
Aug
(83) |
Sep
(114) |
Oct
(90) |
Nov
(181) |
Dec
(85) |
| 2015 |
Jan
(111) |
Feb
(120) |
Mar
(161) |
Apr
(95) |
May
(93) |
Jun
(185) |
Jul
(170) |
Aug
(119) |
Sep
(128) |
Oct
(110) |
Nov
(145) |
Dec
(92) |
| 2016 |
Jan
(105) |
Feb
(106) |
Mar
(101) |
Apr
(59) |
May
(96) |
Jun
(168) |
Jul
(110) |
Aug
(183) |
Sep
(85) |
Oct
(79) |
Nov
(87) |
Dec
(86) |
| 2017 |
Jan
(100) |
Feb
(77) |
Mar
(85) |
Apr
(52) |
May
(60) |
Jun
(63) |
Jul
(67) |
Aug
(24) |
Sep
(1) |
Oct
|
Nov
(2) |
Dec
|
|
From: Bill P. <pa...@ki...> - 2011-09-21 16:59:14
|
Hi Sam, On Sep 21, 2011, at 9:36 AM, Samuel Harrold wrote: > I am having difficulty with the short duration of the SourceForge connection timeout. I've set up a mirror repository of MESA for my own use that I would like to keep updated to SourceForge. However, because of a large replacement operation at r2025, the connection to SourceForge always times out even after many tries before the "svnsync synchronize" operation has gathered the revision data for transmission to the mirror repository, even after many tries. > > Is there a workaround or a solution to the short SourceForge connection timeout? (error message and svn log for r2025 below) I've also had problems with SourceForge timeouts, but I've always eventually managed to get things updated. When I start a check out and get a connection timeout, I retry doing multiple svn update's until all of the files finally make. If I had to get a complete checkout to go though it might never make it! I have no experience with setting up mirror repositories, so I'm afraid I have no help to offer on that one. My use of SourceForge has always been limited to simple commits and check-outs and updates. And I've always been able to work around connection timeout problems for those simple operations. You are trying something way beyond my experience!!! Good luck, Bill |
|
From: Samuel H. <ha...@as...> - 2011-09-21 16:36:34
|
Hi, I am having difficulty with the short duration of the SourceForge connection timeout. I've set up a mirror repository of MESA for my own use that I would like to keep updated to SourceForge. However, because of a large replacement operation at r2025, the connection to SourceForge always times out even after many tries before the "svnsync synchronize" operation has gathered the revision data for transmission to the mirror repository, even after many tries. Is there a workaround or a solution to the short SourceForge connection timeout? (error message and svn log for r2025 below) Thank you, Sam __________ Samuel Harrold Graduate Student Dept of Astronomy Univ Texas at Austin _timeout message_ > svnsync initialize --username sam...@gm... https://mesa-mirror.googlecode.com/svn http://mesa.svn.sourceforge.net/svnroot/mesa Copied properties for revision 0. > svnsync synchronize https://mesa-mirror.googlecode.com/svn Committed revision 1. Copied properties for revision 1. Transmitting file data ... ... Committed revision 2024. Copied properties for revision 2024. svnsync: REPORT of 'http://mesa.svn.sourceforge.net/svnroot/mesa': Could not read response body: Connection reset by peer ( http://mesa.svn.sourceforge.net) _svn log for r2025_ > svn log -r 2025 -v http://mesa.svn.sourceforge.net/svnroot/mesa ------------------------------------------------------------------------ r2025 | aaron_dotter | 2010-01-05 19:27:30 -0600 (Tue, 05 Jan 2010) | 1 line Changed paths: A /trunk/README (from /trunk/README:2021) A /trunk/README_mesa_numerics (from /trunk/README_mesa_numerics:2021) A /trunk/alert (from /trunk/alert:2021) ... (2631 files replaced in total) ... R /trunk/website/src/micro_physics.page (from /trunk/website/src/micro_physics.page:2021) R /trunk/website/src/news.page (from /trunk/website/src/news.page:2021) R /trunk/website/src/numerics.page (from /trunk/website/src/numerics.page:2021) fixing the problem caused by paper ------------------------------------------------------------------------ |
|
From: Bill P. <pa...@ki...> - 2011-09-21 14:09:54
|
Hello, I just found and fixed a nasty memory leak in mesa/star. Please get the bug fix by updating to version 3633. The various installation problems are still a mystery. I wish I could provide quick solutions for everyone who runs into trouble, but it is simply beyond my ability. I need your help -- please don't fail to reply to emails on mesa-users if you have relevant experience. Thanks, Bill |
|
From: Samuel H. <ha...@as...> - 2011-09-20 21:28:37
|
Hi, I get the message below when trying to operate on the root of the SourceForge repository. Is there an outdated hook script left in place? _message_ > svn list http://mesa.svn.sourceforge.net/svnroot svn: Repository moved temporarily to ' http://mesa.svn.sourceforge.net/svnroot'; please relocate Thanks, Sam __________ Samuel Harrold Graduate Student Dept of Astronomy Univ Texas at Austin |
|
From: Ken S. <ke...@as...> - 2011-09-19 19:39:24
|
Hi all. I wanted to modify the opacity used in MESA/star in a way
that the current controls can't do, so I tried to edit my
run_star_extras.f as outlined in other_kap.f. Namely, I added these
lines to the extra_controls subroutine in run_star_extras.f:
----------------
s% use_other_kap = .true.
s% other_kap_get_Type1 => my_other_kap_get_Type1
----------------
As a simple test, my_other_kap_get_Type1 looks like this:
----------------
subroutine my_other_kap_get_Type1( &
id, handle, zbar, X, Zbase, log10_rho, log10_T, &
species, chem_id, net_iso, xa, &
lnfree_e, d_lnfree_e_dlnRho, d_lnfree_e_dlnT, &
kap, dln_kap_dlnRho, dln_kap_dlnT, ierr)
use kap_lib, only : kap_get_Type1
! INPUT
integer, intent(in) :: id ! star id if available; 0 otherwise
integer, intent(in) :: handle ! from alloc_other_kap_handle
real*8, intent(in) :: zbar ! average ion charge
real*8, intent(in) :: X ! the hydrogen mass fraction
real*8, intent(in) :: Zbase ! the metallicity
real*8, intent(in) :: log10_rho ! the density
real*8, intent(in) :: log10_T ! the temperature
integer, intent(in) :: species
integer, pointer :: chem_id(:) ! maps species to chem id
! index from 1 to species
! value is between 1 and num_chem_isos
integer, pointer :: net_iso(:) ! maps chem id to species number
! index from 1 to num_chem_isos (defined in chem_def)
! value is 0 if the iso is not in the current net
! else is value between 1 and number of species in current net
real*8, intent(in) :: xa(:) ! mass fractions
double precision, intent(in) :: lnfree_e, d_lnfree_e_dlnRho,
d_lnfree_e_dlnT
! free_e := total combined number per nucleon of free
electrons and positrons
! OUTPUT
real*8, intent(out) :: kap ! opacity
real*8, intent(out) :: dln_kap_dlnRho ! partial derivative at
constant T
real*8, intent(out) :: dln_kap_dlnT ! partial derivative at
constant Rho
integer, intent(out) :: ierr ! 0 means AOK.
call kap_get_Type1( &
handle, zbar, X, Zbase, log10_rho, log10_T, &
lnfree_e, d_lnfree_e_dlnRho, d_lnfree_e_dlnT, &
kap, dln_kap_dlnRho, dln_kap_dlnT, ierr)
end subroutine my_other_kap_get_Type1
----------------
However, these changes cause the run to exit immediately, without
producing any error messages:
----------------
DATE: 2011-09-19
TIME: 12:27:28
read inlist_wd
set_eos_PC_parameters
mass_fraction_limit_for_PC 1.0000000000000000E-02
logRho1_PC_limit 2.0000000000000000E+01
logRho2_PC_limit 2.0000000000000000E+01
log_Gamma_all_HELM 1.6020600000000000E+00
log_Gamma_all_PC 1.9030899999999999E+00
PC_min_Z 9.9900000000000000E-01
load saved model 1_0_x12_0_50.mod
DATE: 2011-09-19
TIME: 12:27:32
----------------
>From inserting debugging messages, I know the lines "s% use_other_kap
= .true." and "s% other_kap_get_Type1 => my_other_kap_get_Type1" are
reached, but my_other_kap_get_Type1 doesn't ever seem to be called.
Any thoughts out there?
Thanks,
Ken
|
|
From: Aaron D. <aar...@gm...> - 2011-09-19 14:06:12
|
Hi, I'll add my tale of woe to Phil's. Everything goes fine until rates. Looks like Phil and I are consistent, at least. /Users/dotter/mesa3611/rates building rates package. make: `librates.a' is up to date. make: `librates.a' is up to date. gfortran -fopenmp -o ../tester test_rates.o -L../../make -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o -L../../make -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas 58c58 < rne20ap_to_mg24 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rne20ap_to_mg24 1.3318310086100496E-08 60c60 < rmg24ap_to_si28 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rmg24ap_to_si28 2.6558494648638122E-06 62c62 < r_c12_to_he4_he4_he4 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_c12_to_he4_he4_he4 1.4771305847143582E-26 63c63 < r_ne20_ga_o16 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_ne20_ga_o16 3.8122575515639085E-16 /Users/dotter/mesa3611/rates/test TEST FAILED -- compare test_output to tmp.txt /Users/dotter/mesa3611/rates ./build_and_test FAILED This is MESA 3611 on Mac OS X 10.6.8 with gfortran 4.5.2. Aaron On Mon, Sep 19, 2011 at 8:30 AM, Phil Arras <pl...@vi...> wrote: > > Greetings, > > I encountered a compilation error for version 3611, for Mac OS X (10.6.8), > with gfortran 4.6.1: > > > -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- > gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o -L../../make > -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum > -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas > > 58c58 > < rne20ap_to_mg24 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 > > rne20ap_to_mg24 1.3318310086100496E-08 > 60c60 > < rmg24ap_to_si28 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 > > rmg24ap_to_si28 2.6558494648638122E-06 > 62c62 > < r_c12_to_he4_he4_he4 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 > > r_c12_to_he4_he4_he4 1.4771305847143582E-26 > 63c63 > < r_ne20_ga_o16 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 > > r_ne20_ga_o16 3.8122575515639085E-16 > > /Users/arras/stars/mesa/rates/test > TEST FAILED -- compare test_output to tmp.txt > > > /Users/arras/stars/mesa/rates > ./build_and_test FAILED > > -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- > > > > Comparing the two files gave: > > > -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- > Phil-Arras-MacBook-Pro:mesa/rates/test] arras% diff test_output tmp.txt > 1c1 > < > --- > > > 5c5 > < > --- > > > 21c21 > < rn15pg_aux 3.0733146900102957E+02 > --- > > rn15pg_aux 3.0733146900102952E+02 > 26c26 > < r_mg24_ag_si28 7.3991076437333121E-04 > --- > > r_mg24_ag_si28 7.3991076437333111E-04 > 31c31 > < > --- > > > 33,34c33,34 > < > < > --- > > > > > 38c38 > < > --- > > > 58c58 > < rne20ap_to_mg24 1.3318310086100496E-08 > --- > > rne20ap_to_mg24 0.0000000000000000E+00 > 60c60 > < rmg24ap_to_si28 2.6558494648638122E-06 > --- > > rmg24ap_to_si28 0.0000000000000000E+00 > 62,64c62,64 > < r_c12_to_he4_he4_he4 1.4771305847143582E-26 > < r_ne20_ga_o16 3.8122575515639085E-16 > < > --- > > r_c12_to_he4_he4_he4 0.0000000000000000E+00 > > r_ne20_ga_o16 0.0000000000000000E+00 > > > 66,69c66,69 > < > < FL_epsnuc_3alf 1.1675035548327446E+21 > < > < > --- > > > > FL_epsnuc_3alf 1.1675035548327443E+21 > > > > > 73c73 > < > --- > > > 75c75 > < > --- > > > 77c77 > < > --- > > > > -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- > > So they only differ by a tiny amount. > > Oddly, the first time I tried the installation it barfed on the mtx > routine. Then just for the fun of it I tried doing ./install again and it > got through mtx, > but then stopped on rates, and is stuck there no matter how many times I > try ./install. > > Help! > > cheers, > > Phil > > > > ------------------------------------------------------------------------------ > BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA > Learn about the latest advances in developing for the > BlackBerry® mobile platform with sessions, labs & more. > See new tools and technologies. Register for BlackBerry® DevCon today! > http://p.sf.net/sfu/rim-devcon-copy1 > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > |
|
From: Phil A. <pl...@vi...> - 2011-09-19 12:30:40
|
Greetings, I encountered a compilation error for version 3611, for Mac OS X (10.6.8), with gfortran 4.6.1: -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o -L../../make -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas 58c58 < rne20ap_to_mg24 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rne20ap_to_mg24 1.3318310086100496E-08 60c60 < rmg24ap_to_si28 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rmg24ap_to_si28 2.6558494648638122E-06 62c62 < r_c12_to_he4_he4_he4 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_c12_to_he4_he4_he4 1.4771305847143582E-26 63c63 < r_ne20_ga_o16 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_ne20_ga_o16 3.8122575515639085E-16 /Users/arras/stars/mesa/rates/test TEST FAILED -- compare test_output to tmp.txt /Users/arras/stars/mesa/rates ./build_and_test FAILED -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Comparing the two files gave: -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Phil-Arras-MacBook-Pro:mesa/rates/test] arras% diff test_output tmp.txt 1c1 < --- > 5c5 < --- > 21c21 < rn15pg_aux 3.0733146900102957E+02 --- > rn15pg_aux 3.0733146900102952E+02 26c26 < r_mg24_ag_si28 7.3991076437333121E-04 --- > r_mg24_ag_si28 7.3991076437333111E-04 31c31 < --- > 33,34c33,34 < < --- > > 38c38 < --- > 58c58 < rne20ap_to_mg24 1.3318310086100496E-08 --- > rne20ap_to_mg24 0.0000000000000000E+00 60c60 < rmg24ap_to_si28 2.6558494648638122E-06 --- > rmg24ap_to_si28 0.0000000000000000E+00 62,64c62,64 < r_c12_to_he4_he4_he4 1.4771305847143582E-26 < r_ne20_ga_o16 3.8122575515639085E-16 < --- > r_c12_to_he4_he4_he4 0.0000000000000000E+00 > r_ne20_ga_o16 0.0000000000000000E+00 > 66,69c66,69 < < FL_epsnuc_3alf 1.1675035548327446E+21 < < --- > > FL_epsnuc_3alf 1.1675035548327443E+21 > > 73c73 < --- > 75c75 < --- > 77c77 < --- > -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- So they only differ by a tiny amount. Oddly, the first time I tried the installation it barfed on the mtx routine. Then just for the fun of it I tried doing ./install again and it got through mtx, but then stopped on rates, and is stuck there no matter how many times I try ./install. Help! cheers, Phil |
|
From: Ming Y. <my...@ma...> - 2011-09-19 03:02:47
|
Dear all, I have updated the MESA to 3611 yesterday. But there is a problem when installing the atm package. I also read the digest of mesa and find some people have the same problems. Are there any solutions? (change the threads number cannot fix this) /home/ymray/works/mesa/atm building atm package. ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../public -I../private -I../../include -warn all -warn nounused -implicitnone -O2 -fp-model precise -prec-div -prec-sqrt -ftz -c -free ../public/atm_def.f ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../public -I../private -I../../include -warn all -warn nounused -implicitnone -O2 -fp-model precise -prec-div -prec-sqrt -ftz -c -free ../private/integrate_atm.f ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../public -I../private -I../../include -warn all -warn nounused -implicitnone -O2 -fp-model precise -prec-div -prec-sqrt -ftz -c -free ../private/table_atm.f ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../public -I../private -I../../include -warn all -warn nounused -implicitnone -O2 -fp-model precise -prec-div -prec-sqrt -ftz -c -free ../private/mod_atm.f ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../public -I../private -I../../include -warn all -warn nounused -implicitnone -O2 -fp-model precise -prec-div -prec-sqrt -ftz -c -free ../public/atm_lib.f ar crs libatm.a integrate_atm.o table_atm.o mod_atm.o atm_lib.o make: “libatm.a”是最新的。 ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../../make -I../../public -I../../../include -check uninit -check pointers -check bounds -check all -g -c -O2 -free ../src/test_atm_support.f ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../../make -I../../public -I../../../include -check uninit -check pointers -check bounds -check all -g -c -O2 -free ../src/test_atm.f ifort -openmp -threads -o ../tester test_atm_support.o test_atm.o -L../../make -latm -L../../../lib -leos -lkap -lnet -lscreen -lrates -lreaclib -lweak -lneu -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../../make -I../../public -I../../../include -check uninit -check pointers -check bounds -check all -g -c -O2 -free ../src/test_atm_quietly.f ifort -openmp -threads -o ../test_quietly test_atm_support.o test_atm_quietly.o -L../../make -latm -L../../../lib -leos -lkap -lnet -lscreen -lrates -lreaclib -lweak -lneu -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas ifort -vec-report0 -traceback -error-limit 6 -openmp -threads -I../../make -I../../public -I../../../include -check uninit -check pointers -check bounds -check all -g -c -O2 -free ../src/plot_atm.f ifort -openmp -threads -o ../plotter test_atm_support.o plot_atm.o -L../../make -latm -L../../../lib -leos -lkap -lnet -lscreen -lrates -lreaclib -lweak -lneu -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas 64c64 < Do_One_Grey_Irradiated: kap 6.0992874184124302E-02 --- field 3 relative error 2.35e+00 > Do_One_Grey_Irradiated: kap 1.8176165663029643E-02 71c71 < tau 5.8037143231695874E+00 --- field 2 relative error 9.12e-01 > tau 1.1100502830879284E+01 75c75 < T 1.4754276464049410E+03 --- field 2 relative error 1.50e-01 > T 1.6974910696457246E+03 76c76 < logT 3.1689179169620103E+00 --- field 2 relative error 1.92e-02 > logT 3.2298074981726947E+00 /home/ymray/works/mesa/atm/test TEST FAILED -- compare test_output to tmp.txt /home/ymray/works/mesa/atm ./build_and_test FAILED Best Regards --------------------------------- Ming Yang Department of Astronomy Beijing Normal University 19 Xinjiekou Wai Street Beijing, 100875 P.R.China E-mail: my...@ma... --------------------------------- |
|
From: chenhl <chl...@gm...> - 2011-09-18 05:49:37
|
Hello all:
I came across a problem when I was installing the newest MESA
3611.Here is the problem:
cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src
../klu_src/btf_order.c -o btf_order.o
cc: error: unrecognized option ‘-no-cpp-precomp’
make: *** [btf_order.o] Error 1
/home/chenhl/mesa/mtx/make
FAILED
BTW: I use the gcc 4.6.0 on my fedora 15.
Any opinions from the specialists?
Chen
|
|
From: Aaron D. <aar...@gm...> - 2011-09-17 17:54:37
|
Hi Ehsan, The ifort traceback information is telling you that the problem occurs in run_star_extras.f on line 175. Is this in a routine that you modified? The preceding error message tells you than the problem is with the array depsnuc_dlnt. I hope that gives you something to work with as you track down the problem. Aaron On Fri, Sep 16, 2011 at 11:29 PM, Ehsan Moravveji <mor...@ia...>wrote: > Dear MESA friends, > Hello. > I get the following error when trying to run a model of massive star. The > inlist used to work quite fine with the last official release. I use ifort > 12.0.4 on a Linux machine. > > > __________________________________________________________________________________________________________________________________________________ > > step lg_Tcntr Teff lg_LH lg_Lnuc Mass > H_rich H_cntr N_cntr Y_surf X_avg eta_cntr pts retry > lg_dt_yr lg_Dcntr lg_R lg_L3a lg_Lneu lg_Mdot > H_poor He_cntr O_cntr Z_surf Y_avg gam_cntr iters bckup > age_yr lg_Pcntr lg_L lg_LZ lg_Psurf lg_Dsurf > He_poor C_cntr Ne_cntr Si_cntr Z_avg v_div_cs dt_limit > > __________________________________________________________________________________________________________________________________________________ > > 1 7.548962 3.343E+04 4.801886 4.865465 18.999998 > 18.999998 0.699120 0.004651 0.279996 0.699596 -6.013394 1009 > 0 > 1.803116 0.692165 0.758225 -28.142923 3.795931 -7.568808 > 0.000000 0.280376 0.009270 0.020000 0.280172 0.012186 3 > 0 > 6.3550E+01 16.435724 4.566671 3.999594 3.921911 -8.841020 > 0.000000 0.000326 0.001853 0.000607 2.023E-02 0.309E-05 > varcontrol > > forrtl: severe (408): fort: (3): Subscript #1 of the array D_EPSNUC_DLNT > has value 1 which is less than the lower bound of 8 > > Image PC Routine Line > Source > star 0000000000ACE6CA Unknown Unknown Unknown > star 0000000000ACD245 Unknown Unknown Unknown > star 0000000000A7CDF6 Unknown Unknown Unknown > star 0000000000A27906 Unknown Unknown Unknown > star 0000000000A27D59 Unknown Unknown Unknown > star 0000000000405B0B run_star_extras_m 175 > run_star_extras.f > star 000000000047C85C profile_mp_do_pro 557 > profile.f > star 0000000000485C58 profile_mp_write_ 494 > profile.f > star 0000000000488FDA do_one_utils_mp_d 152 > do_one_utils.f > star 00000000004A928E evolve_mp_finish_ 827 > evolve.f > star 0000000000407451 star_lib_mp_star_ 370 > star_lib.f > star 0000000000412B63 run_star_support_ 1989 > run_star_support.f > star 0000000000407018 run_star_mp_do_ru 14 > run_star.f > star 000000000040710A MAIN__ 13 run.f > star 00000000004043FC Unknown Unknown Unknown > libc.so.6 00007F68007D5C4D Unknown Unknown Unknown > star 00000000004042F9 Unknown Unknown Unknown > Fri Sep 16 10:07:19 EDT 2011 > > > Any help is highly appreciated. > > Moravveji, Ehsan. > Ph.D student of Astrophysics. > Department of Physics, Institute for Advanced Studies in Basic Sciences > (IASBS), GavaZang Road, > Zanjan 45137-66731, Iran. > > Office: (+98)241-415 2118 > Fax: (+98)241-415 2104 > http://iasbs.ac.ir/students/moravveji > -- > This message has been scanned for viruses and > dangerous content by *MailScanner* <http://www.mailscanner.info/>, and is > believed to be clean. > > ------------------------------------------------------------------------------ > BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA > http://p.sf.net/sfu/rim-devcon-copy2 > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > > |
|
From: Ehsan M. <mor...@ia...> - 2011-09-17 03:41:49
|
<html><div> Dear MESA friends, Hello. I get the following error when trying to run a model of massive star. The inlist used to work quite fine with the last official release. I use ifort 12.0.4 on a Linux machine. __________________________________________________________________________________________________________________________________________________ step lg_Tcntr Teff lg_LH lg_Lnuc Mass H_rich H_cntr N_cntr Y_surf X_avg eta_cntr pts retry lg_dt_yr lg_Dcntr lg_R lg_L3a lg_Lneu lg_Mdot H_poor He_cntr O_cntr Z_surf Y_avg gam_cntr iters bckup age_yr lg_Pcntr lg_L lg_LZ lg_Psurf lg_Dsurf He_poor C_cntr Ne_cntr Si_cntr Z_avg v_div_cs dt_limit __________________________________________________________________________________________________________________________________________________ 1 7.548962 3.343E+04 4.801886 4.865465 18.999998 18.999998 0.699120 0.004651 0.279996 0.699596 -6.013394 1009 0 1.803116 0.692165 0.758225 -28.142923 3.795931 -7.568808 0.000000 0.280376 0.009270 0.020000 0.280172 0.012186 3 0 6.3550E+01 16.435724 4.566671 3.999594 3.921911 -8.841020 0.000000 0.000326 0.001853 0.000607 2.023E-02 0.309E-05 varcontrol forrtl: severe (408): fort: (3): Subscript #1 of the array D_EPSNUC_DLNT has value 1 which is less than the lower bound of 8 Image PC Routine Line Source star 0000000000ACE6CA Unknown Unknown Unknown star 0000000000ACD245 Unknown Unknown Unknown star 0000000000A7CDF6 Unknown Unknown Unknown star 0000000000A27906 Unknown Unknown Unknown star 0000000000A27D59 Unknown Unknown Unknown star 0000000000405B0B run_star_extras_m 175 run_star_extras.f star 000000000047C85C profile_mp_do_pro 557 profile.f star 0000000000485C58 profile_mp_write_ 494 profile.f star 0000000000488FDA do_one_utils_mp_d 152 do_one_utils.f star 00000000004A928E evolve_mp_finish_ 827 evolve.f star 0000000000407451 star_lib_mp_star_ 370 star_lib.f star 0000000000412B63 run_star_support_ 1989 run_star_support.f star 0000000000407018 run_star_mp_do_ru 14 run_star.f star 000000000040710A MAIN__ 13 run.f star 00000000004043FC Unknown Unknown Unknown libc.so.6 00007F68007D5C4D Unknown Unknown Unknown star 00000000004042F9 Unknown Unknown Unknown Fri Sep 16 10:07:19 EDT 2011 Any help is highly appreciated. Moravveji, Ehsan. Ph.D student of Astrophysics. Department of Physics, Institute for Advanced Studies in Basic Sciences (IASBS), GavaZang Road, Zanjan 45137-66731, Iran. Office: (+98)241-415 2118 Fax: (+98)241-415 2104 http://iasbs.ac.ir/students/moravveji -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean. |
|
From: Tuguldur S. <tug...@em...> - 2011-09-15 01:35:26
|
Did you set up the number of threads correctly? I had the same issue with 3611 on ifort 12.1.0 today. Similar failure on tmp.txt but in the mesa/mxt/test/. I was running on macmini with 2 cores so I did: export OMP_NUM_THREADS=2 this fixed my particular problem with tmp.txt but I got bumped in to another failure later. Tuguldur 2011/9/14 Ryan O'Leary <ol...@be...> > Hi, > > I was having the same issue with gfortran 4.6.1 on Ubuntu for mesa > revision 3611: > GNU Fortran (Ubuntu 4.6.0-3~ppa1) 4.6.1 20110409 (prerelease) > > I eventually decided to just use the intel compilers. > > -Ryan. > > > > > > > On Wed, Sep 14, 2011 at 4:01 PM, Tsing Wong > <tsingwong2012@u.northwestern.edu> wrote: > > Hi, > > > > I tried to install the most recent released copy of mesa with gfortran > 4.6.1 > > and stuck at building the rates package. Here is the error message shown > on > > my screen: > > > > /bigdisk/mesa/rates > > building rates package. > > > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../public/rates_def.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../private/ratelib.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../private/pycno.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../private/raw_rates.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../private/rates_initialize.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -ffree-form -fbounds-check -c -O0 > ../private/rates_names.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../private/rates_support.f > > gfortran -fno-range-check -fopenmp -I../public -I../private > > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > > -ffixed-form -ffixed-line-length-132 ../public/rates_lib.f > > ar crs librates.a rates_def.o ratelib.o pycno.o raw_rates.o > > rates_initialize.o rates_names.o rates_support.o rates_lib.o > > make: Warning: File `librates.a' has modification time 6.3 s in the > future > > make: `librates.a' is up to date. > > make: warning: Clock skew detected. Your build may be incomplete. > > gfortran -fno-range-check -fopenmp -I../../make -I../../public > > -I../../../include -fbounds-check -g -c -ffixed-form > > -ffixed-line-length-132 ../src/test_rates.f > > gfortran -fopenmp -o ../tester test_rates.o -L../../make -lrates > > -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum > -lutils > > -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas > > gfortran -fno-range-check -fopenmp -I../../make -I../../public > > -I../../../include -fbounds-check -g -c -ffree-form > > ../src/reaction_rate_from_ > > cache.f > > gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o > -L../../make > > -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d > -lnum > > -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas > > > > 58c58 > > < rne20ap_to_mg24 0.0000000000000000E+00 > > --- field 2 relative error 1.00e+00 > >> rne20ap_to_mg24 1.3318310086100496E-08 > > 60c60 > > < rmg24ap_to_si28 0.0000000000000000E+00 > > --- field 2 relative error 1.00e+00 > >> rmg24ap_to_si28 2.6558494648638122E-06 > > 62c62 > > < r_c12_to_he4_he4_he4 0.0000000000000000E+00 > > --- field 2 relative error 1.00e+00 > >> r_c12_to_he4_he4_he4 1.4771305847143582E-26 > > 63c63 > > < r_ne20_ga_o16 0.0000000000000000E+00 > > --- field 2 relative error 1.00e+00 > >> r_ne20_ga_o16 3.8122575515639085E-16 > > > > /bigdisk/mesa/rates/test > > TEST FAILED -- compare test_output to tmp.txt > > > > > > /bigdisk/mesa/rates > > ./build_and_test FAILED > > > > Does this happen to anyone too? And does anyone have any idea of how to > fix > > that? > > > > Thanks, > > TsingWai > > > ------------------------------------------------------------------------------ > > BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA > > Learn about the latest advances in developing for the > > BlackBerry® mobile platform with sessions, labs & more. > > See new tools and technologies. Register for BlackBerry® DevCon > today! > > http://p.sf.net/sfu/rim-devcon-copy1 > > _______________________________________________ > > mesa-users mailing list > > mes...@li... > > https://lists.sourceforge.net/lists/listinfo/mesa-users > > > > > > > ------------------------------------------------------------------------------ > Doing More with Less: The Next Generation Virtual Desktop > What are the key obstacles that have prevented many mid-market businesses > from deploying virtual desktops? How do next-generation virtual desktops > provide companies an easier-to-deploy, easier-to-manage and more affordable > virtual desktop model.http://www.accelacomm.com/jaw/sfnl/114/51426474/ > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > -- ------------------------------------------------------------------------- "No word from Landau, No thought from Liefshitz" - Landau's wife "Science is interesting, and if you don't agree, you can fuck off" - Editor of New Scientist Magazine Three out of two people have difficulties with fractions. Kenyan proverb: "Take care of the Earth. It was not given to you by your parents; it was loaned to you by your children." "The Bible tells us how to go to heaven, not how the heavens go" Galileo |
|
From: Ryan O'L. <ol...@be...> - 2011-09-15 00:20:30
|
Hi, I was having the same issue with gfortran 4.6.1 on Ubuntu for mesa revision 3611: GNU Fortran (Ubuntu 4.6.0-3~ppa1) 4.6.1 20110409 (prerelease) I eventually decided to just use the intel compilers. -Ryan. On Wed, Sep 14, 2011 at 4:01 PM, Tsing Wong <tsingwong2012@u.northwestern.edu> wrote: > Hi, > > I tried to install the most recent released copy of mesa with gfortran 4.6.1 > and stuck at building the rates package. Here is the error message shown on > my screen: > > /bigdisk/mesa/rates > building rates package. > > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../public/rates_def.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../private/ratelib.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../private/pycno.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../private/raw_rates.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../private/rates_initialize.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -ffree-form -fbounds-check -c -O0 ../private/rates_names.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../private/rates_support.f > gfortran -fno-range-check -fopenmp -I../public -I../private > -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c > -ffixed-form -ffixed-line-length-132 ../public/rates_lib.f > ar crs librates.a rates_def.o ratelib.o pycno.o raw_rates.o > rates_initialize.o rates_names.o rates_support.o rates_lib.o > make: Warning: File `librates.a' has modification time 6.3 s in the future > make: `librates.a' is up to date. > make: warning: Clock skew detected. Your build may be incomplete. > gfortran -fno-range-check -fopenmp -I../../make -I../../public > -I../../../include -fbounds-check -g -c -ffixed-form > -ffixed-line-length-132 ../src/test_rates.f > gfortran -fopenmp -o ../tester test_rates.o -L../../make -lrates > -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils > -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas > gfortran -fno-range-check -fopenmp -I../../make -I../../public > -I../../../include -fbounds-check -g -c -ffree-form > ../src/reaction_rate_from_ > cache.f > gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o -L../../make > -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum > -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas > > 58c58 > < rne20ap_to_mg24 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 >> rne20ap_to_mg24 1.3318310086100496E-08 > 60c60 > < rmg24ap_to_si28 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 >> rmg24ap_to_si28 2.6558494648638122E-06 > 62c62 > < r_c12_to_he4_he4_he4 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 >> r_c12_to_he4_he4_he4 1.4771305847143582E-26 > 63c63 > < r_ne20_ga_o16 0.0000000000000000E+00 > --- field 2 relative error 1.00e+00 >> r_ne20_ga_o16 3.8122575515639085E-16 > > /bigdisk/mesa/rates/test > TEST FAILED -- compare test_output to tmp.txt > > > /bigdisk/mesa/rates > ./build_and_test FAILED > > Does this happen to anyone too? And does anyone have any idea of how to fix > that? > > Thanks, > TsingWai > ------------------------------------------------------------------------------ > BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA > Learn about the latest advances in developing for the > BlackBerry® mobile platform with sessions, labs & more. > See new tools and technologies. Register for BlackBerry® DevCon today! > http://p.sf.net/sfu/rim-devcon-copy1 > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > > |
|
From: Tsing W. <tsingwong2012@u.northwestern.edu> - 2011-09-14 23:58:04
|
Hi, I tried to install the most recent released copy of mesa with gfortran 4.6.1 and stuck at building the rates package. Here is the error message shown on my screen: /bigdisk/mesa/rates building rates package. gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../public/rates_def.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../private/ratelib.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../private/pycno.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../private/raw_rates.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../private/rates_initialize.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -ffree-form -fbounds-check -c -O0 ../private/rates_names.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../private/rates_support.f gfortran -fno-range-check -fopenmp -I../public -I../private -I../../include -Wunused-value -Werror -W -fimplicit-none -O2 -c -ffixed-form -ffixed-line-length-132 ../public/rates_lib.f ar crs librates.a rates_def.o ratelib.o pycno.o raw_rates.o rates_initialize.o rates_names.o rates_support.o rates_lib.o make: Warning: File `librates.a' has modification time 6.3 s in the future make: `librates.a' is up to date. make: warning: Clock skew detected. Your build may be incomplete. gfortran -fno-range-check -fopenmp -I../../make -I../../public -I../../../include -fbounds-check -g -c -ffixed-form -ffixed-line-length-132 ../src/test_rates.f gfortran -fopenmp -o ../tester test_rates.o -L../../make -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas gfortran -fno-range-check -fopenmp -I../../make -I../../public -I../../../include -fbounds-check -g -c -ffree-form ../src/reaction_rate_from_ cache.f gfortran -fopenmp -o ../show_rates reaction_rate_from_cache.o -L../../make -lrates -L../../../lib -lreaclib -lweak -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas 58c58 < rne20ap_to_mg24 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rne20ap_to_mg24 1.3318310086100496E-08 60c60 < rmg24ap_to_si28 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > rmg24ap_to_si28 2.6558494648638122E-06 62c62 < r_c12_to_he4_he4_he4 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_c12_to_he4_he4_he4 1.4771305847143582E-26 63c63 < r_ne20_ga_o16 0.0000000000000000E+00 --- field 2 relative error 1.00e+00 > r_ne20_ga_o16 3.8122575515639085E-16 /bigdisk/mesa/rates/test TEST FAILED -- compare test_output to tmp.txt /bigdisk/mesa/rates ./build_and_test FAILED Does this happen to anyone too? And does anyone have any idea of how to fix that? Thanks, TsingWai |
|
From: Ken S. <ke...@as...> - 2011-09-13 22:30:07
|
Hi all. When trying to install MESA v3611 with ifort v12.1, I run into this problem: ------------------------ dyld: Library not loaded: libiomp5.dylib Referenced from: /Users/kenshen/Documents/programs/mesa/const/test/./tester Reason: image not found ./rn: line 1: 66257 Trace/BPT trap ./tester ------------------------ In the past, I fixed this by setting the appropriate path in my .bash_profile: ------------------------ export DYLD_LIBRARY_PATH="/opt/intel/composerxe/lib" ------------------------ However, this no longer seems to fix it. The tester executable is still being linked without the appropriate library directory, as evidenced by an otool command: ------------------------ Ken-Shens-MacBook-Pro:test kenshen$ otool -L tester tester: libiomp5.dylib (compatibility version 5.0.0, current version 5.0.0) /usr/lib/libgcc_s.1.dylib (compatibility version 1.0.0, current version 830.0.0) /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 125.2.11) ------------------------ What silly error am I committing?? Thanks, Ken |
|
From: Bill P. <pa...@ki...> - 2011-09-12 15:38:13
|
Hi,
On Sep 12, 2011, at 7:33 AM, Ehsan Moravveji wrote:
> I want to know if it is possible to have uneven mesh spacing in some regions of interest in our models during evolution.
> For instance, I like all layers whose interior mass m(r) is less than say 1.2345 Msun have 3 times more dense mesh than the evelope. Or similarly in terms of log T or distance from the core.
> Do we have currently a criteria for this in MESA?
There are lots of controls to adjust the mesh. start by reading the section in star_defaults "gradient controls for mesh adjustment".
The basic idea is that the grid will be adjusted so that none of the mesh functions have any cell-to-cell jumps larger than the value of mesh_delta_coeff.
For example the P_function = P_function_weight*log10(P), with a default P_function_weight = 30.
This means the mesh will be adjusted so that the cell to cell jump of 30*log10(P) <= mesh_delta_coeff.
There are mesh functions for T, R, and M that might be useful.
The ones for grad_ad and omega are more experimental -- I haven't found them to be very useful so far.
There are also mesh functions for abundance mass fractions -- these can be very handy.
And controls to modify the value of mesh_delta_coeff around convection boundaries.
The actual mesh is constrained by more than the mesh functions -- there are also restrictions on the relative sizes of adjacent cells.
Also there are single step limits on how many times a cell can be divided and similarly on how many adjacent cells can be merged.
Try to see if you can get what you want to adjusting the current set of mesh functions.
I can add an option for a user-defined mesh function if necessary.
Cheers,
Bill
p.s. as an example of special mesh adjustments, here's the relevant part of the inlist in test_suite/make_solar
! mesh adjustment
mesh_delta_coeff = 0.9
P_function_weight = 25
T_function1_weight = 75
xtra_coef_czb_full_on = 0
xtra_coef_czb_full_off = 1
xtra_coef_a_l_nb_czb = 0.01 ! above lower nonburn convective boundary
xtra_dist_a_l_nb_czb = 10 ! above lower nonburn convective boundary
xtra_coef_b_l_nb_czb = 0.05 ! below lower nonburn convective boundary
xtra_dist_b_l_nb_czb = 3 ! below lower nonburn convective boundary
xa_function_species(1) = 'he4' ! name of nuclide as defined in chem_def
xa_function_weight(1) = 80
xa_function_param(1) = 1d-2
xa_function_species(2) = 'he3' ! name of nuclide as defined in chem_def
xa_function_weight(2) = 20
xa_function_param(2) = 1d-5
|
|
From: Bill P. <pa...@ki...> - 2011-09-11 15:57:58
|
Hi Zahra, On Sep 10, 2011, at 1:36 AM, zahra altaha motahar wrote: > I want to use the MESA code for modelling the evolution of AGB stars. I wonder if the MESA code supports the input physics and factors that are important for evolution of stars in this epoch. In my modelling I need to include the followings: > > 1) convective core overshooting: overshooting parameter should be set for different initial masses > > 2) adopting parametrization which the third dredge up mixing of the stellar mantle as a result of the H-burning phase, and the enhanced luminosity of massive AGB stars undergoing Hot Bottom Burning (HBB). > > 3) molecular opacity > > 4) predictions for dust production in the winds of LPVs and the associated reddening. > > 5) predictions for the radial pulsation properties (period and mode) as a function of luminosity. Those are big questions! Certainly too big for me to answer directly. Perhaps we can get to specific items that I can answer if you find a case that works for something related to what you want to do, and then ask about details of converting it to do what you really want. For your case, perhaps the "c13_pocket" test suite case would do as a start. Then, for info about features that are currently supported in mesa, you should scan the files "star_defaults" and "run_star_defaults" in mesa/star/public. You'll find sections on overshooting, opacities, and winds. You'll also find controls related to brunt frequencies and ways to output FGONG files for use in getting pulsation modes. Good luck, Bill |
|
From: Frank T. <fx...@ma...> - 2011-09-10 17:32:17
|
dear zahra, i suggest reading the first mesa paper - http://iopscience.iop.org/0067-0049/192/1/3/ fxt On Sep 10, 2011, at 1:36 AM, zahra altaha motahar wrote: > Dear MESA developers , > > I have started a new research project recently. I want to use the MESA code for modelling the evolution of AGB stars. I wonder if the MESA code supports the input physics and factors that are important for evolution of stars in this epoch. In my modelling I need to include the followings: > > 1) convective core overshooting: overshooting parameter should be set for different initial masses > > 2) adopting parametrization which the third dredge up mixing of the stellar mantle as a result of the H-burning phase, and the enhanced luminosity of massive AGB stars undergoing Hot Bottom Burning (HBB). > > 3) molecular opacity > > 4) predictions for dust production in the winds of LPVs and the associated reddening. > > 5) predictions for the radial pulsation properties (period and mode) as a function of luminosity. > > How can I find these aspects in the code? If they are not included in the MESA, Is it possible to add them? > > Best Regards, > > Zahra Motahar |
|
From: zahra a. m. <zam...@ya...> - 2011-09-10 08:36:44
|
Dear MESA developers , I have started a new research project recently. I want to use the MESA code for modelling the evolution of AGB stars. I wonder if the MESA code supports the input physics and factors that are important for evolution of stars in this epoch. In my modelling I need to include the followings: 1) convective coreovershooting: overshooting parameter should be set for different initial masses 2) adopting parametrization which the third dredge up mixing of the stellar mantle as a result of the H-burning phase, and the enhanced luminosity of massive AGB stars undergoing Hot Bottom Burning (HBB). 3) molecular opacity 4) predictions for dust production in the winds of LPVs and the associated reddening. 5) predictions for the radial pulsation properties (period and mode)as a function of luminosity. How can I find these aspects in the code? If they are not included in the MESA, Is it possible to add them? Best Regards, Zahra Motahar |
|
From: Bill P. <pa...@ki...> - 2011-09-08 14:32:42
|
On Sep 8, 2011, at 1:35 AM, 張世昕 wrote:
> Because I combine another simple code with MESA, there are some
> variables I want to add ( of course) that aren't included in log_columns.list .
> Could you teach me how to do ?
Hi,
here's how to do it.
in your working copy of run_star_extras, replace
include 'standard_run_star_extras.dek'
by the contents of the included file from star/public
Don't make any edits to any of the files in star/public.
You do all of this in your private copy of run_star_extras.
before doing anything else, let's make sure your working copy of run_star_extras works.
edit the extras_controls routine
subroutine extras_controls(s, ierr)
type (star_info), pointer :: s
integer, intent(out) :: ierr
ierr = 0
write(*,*) 'hello from extra_controls'
end subroutine extras_controls
then, in your work directory, do ./mk and ./rn to check that it is okay.
assuming that worked, now edit the function how_many_extra_log_columns
to specify the number of columns you'll add.
then edit the subroutine data_for_extra_log_columns to give
the names and values. these 2 routines will be called every time
the system writes a new line in star.log
you can also add columns to profiles by editing the routines
how_many_extra_profile_columns and data_for_extra_profile_columns.
good luck!
-Bill
|
|
From: Max K. <ka...@rp...> - 2011-09-07 14:18:31
|
I too have had some trouble installing 3607 on Ubuntu. I am using version 11.04 but gcc 4.6 is not officially supported until 11.10 is released next month. Nevertheless, I installed the 4.6 binary available at http://gfortran.org/download/x86_64/snapshots/ (it's a different URL for the 32 bit version) and it seems to work relatively fine after pointing the system to the manually downloaded location. Even in 4.6, the 'no-cpp-precomp' warning still occurs but unlike 4.5 does *not *cause the installation to fail. My installation failure seems to be some sort of linking error between gcc and MPFR (I'm guessing because of my choice to manually install it), so as far as I can tell downloading the 'unsupported' version of gcc may be a good place to start. Max Katz Rensselaer Polytechnic Institute On Wed, Sep 7, 2011 at 1:49 AM, Gautschy Alfred < alf...@li...> wrote: > Dear Bill > > I posted the material below on the Help Forum of Mesa that might or might > not take off one day. As it might not attract yet enough attention, I also > send you copy of it for your information. I.e. to let you know that I > continue to have problems with the installation of 3607 despite considerable > changes in the system setup. > > > This posting is related to the one in > > https://sourceforge.net/mailarchive/message.php?msg_id=28040448 > > In the mean time I upgraded my operating system and the compilers > to the current system configuration: Ubuntu 11.04 64-bit with ifort 12.1.0 > 20110811 > > However, as previously on my old 32-bit Ubuntu 10.04 LTS system, > I fail to get the latest Mesa release installed. > The installation process on the 64-bit system dies at another place though: > > ifort -openmp -threads -o ../test_quietly test_ionization_support.o > test_ionization_quitely.o -L../../make -lionization -L../../../lib -leos > -lkap -lnet -lscreen -lrates -lreaclib -lweak -lneu -lchem -linterp_2d > -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack > -lmesablas > > forrtl: severe (174): SIGSEGV, segmentation fault occurred > > Image PC Routine Line Source > > test_quietly 000000000043342A mod_ionizationdo_ 154 mod_ionization.f > test_quietly 00000000004324C2 mod_ionization_mp 86 mod_ionization.f > test_quietly 0000000000431E48 mod_ionization_mp 231 mod_ionization.f > test_quietly 000000000040D10B ionization_lib_mp 93 ionization_lib.f > test_quietly 00000000004044DF test_ionization_s 90 > test_ionization_support.f > test_quietly 0000000000404418 test_ionization_s 77 > test_ionization_support.f > test_quietly 0000000000404230 test_ionization_s 58 > test_ionization_support.f > test_quietly 000000000040C256 MAIN__ 3 > test_ionization_quitely.f > test_quietly 0000000000403DEC Unknown Unknown Unknown > libc.so.6 00007F0298D99EFF Unknown Unknown Unknown > test_quietly 0000000000403CE9 Unknown Unknown Unknown > > The commong dark clouds I observe during both installation attempts (32-bit > and 64-bit) > are via the warnings: > > cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src > ../klu_src/btf_order.c -o btf_order.o cc: unrecognized option > '-no-cpp-precomp' > cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src > ../klu_src/btf_maxtrans.c -o btf_maxtrans.o > cc: unrecognized option '-no-cpp-precomp' cc -O3 -fno-common > -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/btf_strongcomp.c -o > btf_strongcomp.o > cc: unrecognized option '-no-cpp-precomp' cc -O3 -fno-common > -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/colamd_global.c -o > colamd_global.o > > ... and so on through the whole package > > I wonder if this is related to the installation failure. > > It is worthwhile to note that the ubuntu distribution does not provide > access to gcc 4.6.x versions, only gcc 4.5 is supported and offered (could > not find any place where a gcc 4.6 package could be got) > > BTW: I did not change anything in the makefile_header file, as compared > with > older installations/updates (i.e. before 3607) > > Any opinions from the specialists? > Alfred > > > ------------------------------------------------------------------------------ > Using storage to extend the benefits of virtualization and iSCSI > Virtualization increases hardware utilization and delivers a new level of > agility. Learn what those decisions are and how to modernize your storage > and backup environments for virtualization. > http://www.accelacomm.com/jaw/sfnl/114/51434361/ > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > |
|
From: Gautschy A. <alf...@li...> - 2011-09-07 06:02:44
|
Dear Bill I posted the material below on the Help Forum of Mesa that might or might not take off one day. As it might not attract yet enough attention, I also send you copy of it for your information. I.e. to let you know that I continue to have problems with the installation of 3607 despite considerable changes in the system setup. This posting is related to the one in https://sourceforge.net/mailarchive/message.php?msg_id=28040448 In the mean time I upgraded my operating system and the compilers to the current system configuration: Ubuntu 11.04 64-bit with ifort 12.1.0 20110811 However, as previously on my old 32-bit Ubuntu 10.04 LTS system, I fail to get the latest Mesa release installed. The installation process on the 64-bit system dies at another place though: ifort -openmp -threads -o ../test_quietly test_ionization_support.o test_ionization_quitely.o -L../../make -lionization -L../../../lib -leos -lkap -lnet -lscreen -lrates -lreaclib -lweak -lneu -lchem -linterp_2d -linterp_1d -lnum -lutils -lalert -lconst -lmtx -lmesaklu -lmesalapack -lmesablas forrtl: severe (174): SIGSEGV, segmentation fault occurred Image PC Routine Line Source test_quietly 000000000043342A mod_ionizationdo_ 154 mod_ionization.f test_quietly 00000000004324C2 mod_ionization_mp 86 mod_ionization.f test_quietly 0000000000431E48 mod_ionization_mp 231 mod_ionization.f test_quietly 000000000040D10B ionization_lib_mp 93 ionization_lib.f test_quietly 00000000004044DF test_ionization_s 90 test_ionization_support.f test_quietly 0000000000404418 test_ionization_s 77 test_ionization_support.f test_quietly 0000000000404230 test_ionization_s 58 test_ionization_support.f test_quietly 000000000040C256 MAIN__ 3 test_ionization_quitely.f test_quietly 0000000000403DEC Unknown Unknown Unknown libc.so.6 00007F0298D99EFF Unknown Unknown Unknown test_quietly 0000000000403CE9 Unknown Unknown Unknown The commong dark clouds I observe during both installation attempts (32-bit and 64-bit) are via the warnings: cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/btf_order.c -o btf_order.o cc: unrecognized option '-no-cpp-precomp' cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/btf_maxtrans.c -o btf_maxtrans.o cc: unrecognized option '-no-cpp-precomp' cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/btf_strongcomp.c -o btf_strongcomp.o cc: unrecognized option '-no-cpp-precomp' cc -O3 -fno-common -no-cpp-precomp -fexceptions -c -I../klu_src ../klu_src/colamd_global.c -o colamd_global.o ... and so on through the whole package I wonder if this is related to the installation failure. It is worthwhile to note that the ubuntu distribution does not provide access to gcc 4.6.x versions, only gcc 4.5 is supported and offered (could not find any place where a gcc 4.6 package could be got) BTW: I did not change anything in the makefile_header file, as compared with older installations/updates (i.e. before 3607) Any opinions from the specialists? Alfred |
|
From: Giannotti, M. <MGi...@ma...> - 2011-09-05 22:40:23
|
Hi Bill,
Sorry, I still have some problems. I need to use z2bar in neu/private/neu.f. I tried to make the changes that you recommended and they worked (no compilation error). However, when I try to use z2bar in neu.f (in the subroutine "neutrinos") I have a compilation error since, it seems, that subroutine does not know what z2bar is. Any suggestion?
Thanks again,
maurizio
Maurizio Giannotti,
Assistant Professor of Physics,
Physical Sciences, Barry University,
11300 NE 2nd Ave., Miami Shores, FL 33161
mgi...@ma...<mailto:mgi...@ma...>
305-899-4565 (office), 305-899-3479 (fax)
________________________________
From: Bill Paxton [pa...@ki...]
Sent: Monday, September 05, 2011 4:10 PM
To: Giannotti, Maurizio
Cc: mes...@li...
Subject: Re: [mesa-users] use mass fractions in neu.f
Hi Maurizio,
We need to add z2bar to the argument list of neu_get for you.
Then we need to change star to pass that new arg.
I'll change mesa to do that, but if you'd like to try it with
your current version, you can make 2 small changes:
first, in neu/public/neu_lib, add z2bar as another arg for neu_get.
old:
subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, log10_Tlim, flags,
> loss, sources, info)
new:
subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, z2bar, log10_Tlim, flags,
> loss, sources, info)
.....
double precision, intent(in) :: zbar2 ! mean charge squared
then
cd neu/test
./mk
./ck
./export
then in star/private/micro.f
old:
call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), &
log10_Tlim, flags, loss, sources, ierr)
new (with z2bar):
call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), s% z2bar(k), &
log10_Tlim, flags, loss, sources, ierr)
then
cd star/test
./mk
./ck
./export
Let me know how it goes.
Cheers,
Bill
On Sep 5, 2011, at 11:05 AM, Giannotti, Maurizio wrote:
Hello,
I am trying to add a few lines to neu.f (in mesa/neu/private) and I need to know the values of the hydrogen, helium, carbon and oxygen mass fraction. More specifically, I need to calculate
y_e+\sum z_j^2 y_j
where y_e and y_j are the number fractions per baryon of electrons and nuclei. From neu.f I can access directly only y_e.
What should I do (or where should i look)?
Thank you very much,
Maurizio
Maurizio Giannotti,
Assistant Professor of Physics,
Physical Sciences, Barry University,
11300 NE 2nd Ave., Miami Shores, FL 33161
mgi...@ma...<mailto:mgi...@ma...>
305-899-4565 (office), 305-899-3479 (fax)
------------------------------------------------------------------------------
Special Offer -- Download ArcSight Logger for FREE!
Finally, a world-class log management solution at an even better
price-free! And you'll get a free "Love Thy Logs" t-shirt when you
download Logger. Secure your free ArcSight Logger TODAY!
http://p.sf.net/sfu/arcsisghtdev2dev_______________________________________________
mesa-users mailing list
mes...@li...<mailto:mes...@li...>
https://lists.sourceforge.net/lists/listinfo/mesa-users
|
|
From: Bill P. <pa...@ki...> - 2011-09-05 20:10:48
|
Hi Maurizio,
We need to add z2bar to the argument list of neu_get for you.
Then we need to change star to pass that new arg.
I'll change mesa to do that, but if you'd like to try it with
your current version, you can make 2 small changes:
first, in neu/public/neu_lib, add z2bar as another arg for neu_get.
old:
subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, log10_Tlim, flags,
> loss, sources, info)
new:
subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, z2bar, log10_Tlim, flags,
> loss, sources, info)
.....
double precision, intent(in) :: zbar2 ! mean charge squared
then
cd neu/test
./mk
./ck
./export
then in star/private/micro.f
old:
call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), &
log10_Tlim, flags, loss, sources, ierr)
new (with z2bar):
call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), s% z2bar(k), &
log10_Tlim, flags, loss, sources, ierr)
then
cd star/test
./mk
./ck
./export
Let me know how it goes.
Cheers,
Bill
On Sep 5, 2011, at 11:05 AM, Giannotti, Maurizio wrote:
> Hello,
>
> I am trying to add a few lines to neu.f (in mesa/neu/private) and I need to know the values of the hydrogen, helium, carbon and oxygen mass fraction. More specifically, I need to calculate
>
> y_e+\sum z_j^2 y_j
>
> where y_e and y_j are the number fractions per baryon of electrons and nuclei. From neu.f I can access directly only y_e.
>
> What should I do (or where should i look)?
>
> Thank you very much,
> Maurizio
>
> Maurizio Giannotti,
> Assistant Professor of Physics,
> Physical Sciences, Barry University,
> 11300 NE 2nd Ave., Miami Shores, FL 33161
> mgi...@ma...
> 305-899-4565 (office), 305-899-3479 (fax)
>
> ------------------------------------------------------------------------------
> Special Offer -- Download ArcSight Logger for FREE!
> Finally, a world-class log management solution at an even better
> price-free! And you'll get a free "Love Thy Logs" t-shirt when you
> download Logger. Secure your free ArcSight Logger TODAY!
> http://p.sf.net/sfu/arcsisghtdev2dev_______________________________________________
> mesa-users mailing list
> mes...@li...
> https://lists.sourceforge.net/lists/listinfo/mesa-users
|
|
From: Giannotti, M. <MGi...@ma...> - 2011-09-05 18:05:22
|
Hello, I am trying to add a few lines to neu.f (in mesa/neu/private) and I need to know the values of the hydrogen, helium, carbon and oxygen mass fraction. More specifically, I need to calculate y_e+\sum z_j^2 y_j where y_e and y_j are the number fractions per baryon of electrons and nuclei. From neu.f I can access directly only y_e. What should I do (or where should i look)? Thank you very much, Maurizio Maurizio Giannotti, Assistant Professor of Physics, Physical Sciences, Barry University, 11300 NE 2nd Ave., Miami Shores, FL 33161 mgi...@ma...<mailto:mgi...@ma...> 305-899-4565 (office), 305-899-3479 (fax) |
138 messages has been excluded from this view by a project administrator.
