[FLASH-USERS] Target expansion, adjustEvolution.F90 for bdry_var

[Aodhan McIlvenny]

Hi Klaus,

many thanks for the help. This makes sense now and the simulation runs 
much more like how I would expect.

All the best,

Aodhan

On 30/04/2020 17:22, Klaus Weide wrote:
> This message is from an external sender. Please take care when responding, clicking links or opening attachments.
>
> On Thu, 30 Apr 2020, Aodhan McIlvenny wrote:
>
>> http://flash.uchicago.edu/pipermail/flash-users/2019-April/002869.html
>>
>> http://flash.uchicago.edu/pipermail/flash-users/2014-March/001427.html
>>
>> So i get a significant uniform expansion of the target over a few ns which is
>> quite unphysical. As mentioned in these replies, a workaround is to set the
>> BDRY_VAR to 1 inside the solid and only set it to -1 when the temperature
>> inside the cell is greater than some threshold. I have initialised the solid
>> as such however I am not sure how to call this information in the
>> adjustEvolution.F90 function to first test the cell temperature and then
>> changing the BDRY to -1 as appropriate.
> Thank you for first checking the archives!
>
> The code for doing what you want could be as simple as
>
>       !  hardwired temperature threshold of 1000 K
>       if (blkPtr(TEMP_VAR,i,j,k) > 1.e3) then
>          blkPtr(BDRY_VAR,i,j,k) = -1.0
>       end if
>
> in Simulation_adjustEvolution.F90 in your simulation directory under
> .../SimulationMain.
>
> If you don't know what to put around those lines, look at other instances
> of Simulation_adjustEvolution.F90 in various simulation directories.
> If you can't find something more suitable,
>    Simulation/SimulationMain/SodSpherical/Simulation_adjustEvolution.F90
> may be a good starting point (still more complicated than you need for
> your purpose, you'll have to delete stuff).
>
> Klaus