[mesa-users] Fwd: Suggesting run-time parameter update
Bill Paxton
paxton at kitp.ucsb.edu
Mon Sep 17 13:17:47 EDT 2012
Hi,
Again, can anyone else help Ehsan debug this.
Thanks,
Bill
Begin forwarded message:
> From: "Ehsan Moravveji" <moravveji at iasbs.ac.ir>
> Date: September 16, 2012 11:11:25 PM PDT
> To: "Bill Paxton" <paxton at kitp.ucsb.edu>
> Subject: Re: Suggesting run-time parameter update
>
> Thanks Bill for your guideline. I followed it exactly. I compile fine.
> However, during the run-time, I get this error:
> __________________________________________________________________________________________________________________________________________________
>
> 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 lg_Pcntr lg_L lg_LZ lg_Psurf lg_Dsurf
> He_poor C_cntr Ne_cntr Z_cntr Z_avg v_div_cs
> dt_limit
> __________________________________________________________________________________________________________________________________________________
>
>
> save profile LOGS/log1.data 1
>
> pgstar_out/main1.png
>
> tried to increase number of mesh points beyond max allowed nz 80000
>
> mesh_plan problem
> doing mesh_call_number 206
> s% model_number 7
>
> failed in do_mesh
>
>
>
> terminated evolution: adjust_mesh_failed
>
> 7 5.476649 3561.274 -21.697809 -21.697809 10.000000
> 10.000000 0.719989 0.000725 0.266011 0.719989 -12.291545
> 44530 0
> -5.405369 -5.148974 2.370404 -99.000000 -23.095677 -99.000000
> 0.000000 0.266011 0.006003 0.014000 0.266011 0.016070 2
> 0
> 5.4271E-05 8.490208 3.900712 -99.000000 3.371681 -7.991517
> 0.000000 0.002476 0.001316 1.400E-02 1.400E-02 0.157E+00
> varcontrol
>
> I presume more has to be done when tweaking mesh adjustment controls. True?
> Is there any conflict?
>
> I have attached all files that you need to reproduce the error. The three
> modules that start with user_*.f have to be put at any directory.
> You can also see how I have implemented the runtime control update in my
> extras_finish_step.
>
> Thanks again.
> Ehsan.
>
>
>> Hi,
>>
>> I've pondered the option of run-time changes to &controls myself,
>> but I've decided against doing it wholesale, and I don't want to
>> get into having controls that say which controls are run-time
>> and which are start-time only.
>>
>> So here's my suggestion for you -- do it yourself in your run_star_extras
>> in extras_finish_step with your own namelist file for things you
>> want to be able to modify at runtime. E.g., your "runtime" namelist
>> could have values for mesh_delta_coeff and varcontrol_target.
>> Then at each step you read in that namelist and change
>> the values for s% mesh_delta_coeff and s% varcontrol_target.
>>
>> Share it with us once you have it working!
>>
>> Cheers,
>> Bill
>>
>> p.s.
>>
>> here's how you define your namelist (at the top of your run_star_extras)
>>
>> real(dp) :: mesh_delta_coeff, varcontrol_target
>> namelist /runtime_controls/ mesh_delta_coeff, varcontrol_target
>>
>> here's how you should open a namelist file (assuming the file name is
>> inlist_runtime_controls)
>> unit = ... << get an iounit in your favorite way
>> open(unit=unit, file='inlist_runtime_controls', action='read',
>> delim='quote', iostat=ierr)
>>
>> here's how you should read it (assuming your namelist is called
>> 'runtime_controls')
>> read(unit, nml= runtime_controls, iostat=ierr)
>>
>> after doing the read, check ierr to make sure it worked, then do
>> s% mesh_delta_coeff = mesh_delta_coeff
>> s% varcontrol_target = varcontrol_target
>>
>> don't forget to close the file
>> close(unit)
>>
>> the inlist_runtime_controls file should look something like this
>>
>> &runtime_controls ! read these at each timestep
>> mesh_delta_coeff = 0.5
>> varcontrol_target = 0.5d-4
>> / ! end of runtime_controls namelist
>>
>>
>>
>>
>>
>> On Sep 15, 2012, at 9:09 AM, Ehsan Moravveji wrote:
>>
>>> Dear Bill,
>>> I came up with a minor suggestion, yet, it is yours to decide upon.
>>>
>>> The &pgstar namelist is modifiable during the run-time. It means, once
>>> the
>>> user edits and saves the &pgstar namelist during the run-time, the
>>> effects
>>> are simultaneous.
>>>
>>> The same idea could (probably and/or possibly) be pursued for some other
>>> parameters in &controls namelist. As we know, this is
>>> an-already-implemented capability for restarts photos.
>>>
>>> Given the fact that the speed and precision of the evolution is partly
>>> (read critically) governed by
>>> mesh_delta_coeff = 0.4
>>> varcontrol_target = 1d-4
>>> what if these two parameters be also given the capability of run-time
>>> modification?
>>>
>>> It is up to you and the MESA board to decide upon implementing this.
>>>
>>> Best wishes.
>>> Ehsan.
>>>
>>>
>>> --
>>> This message has been scanned for viruses and
>>> dangerous content by MailScanner, and is
>>> believed to be clean.
>>>
>>
>>
>> --
>> This message has been scanned for viruses and
>> dangerous content by MailScanner, and is
>> believed to be clean.
>>
>>
>
> --
> This message has been scanned for viruses and
> dangerous content by MailScanner, and is
> believed to be clean.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: makefile
Type: application/octet-stream
Size: 1701 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist
Type: application/octet-stream
Size: 395 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_wr
Type: application/octet-stream
Size: 25984 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0002.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0003.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_runtime_controls
Type: application/octet-stream
Size: 314 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0003.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0004.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: application/octet-stream
Size: 16907 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0004.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0005.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: user_evol_state.f
Type: application/octet-stream
Size: 10481 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0005.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0006.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: user_rotation.f
Type: application/octet-stream
Size: 2841 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0006.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0007.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: user_wolf_rayet.f
Type: application/octet-stream
Size: 5645 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0007.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/49c143e1/attachment-0008.html>
More information about the Mesa-users
mailing list