[FLASH-USERS] Initialising the magnetic field, stagard mesh MHD

[Simon Daley-Yates]

Thanks Klaus, that makes things clear.

Cheers,

Simon


On 24/04/2017 17:36, Klaus Weide wrote:
> On Mon, 24 Apr 2017, Simon Daley-Yates wrote:
>
>> There are a number of different examples of setting up a MHD simulation in the
>> source/Simulation/SimulationMain/magnetoMHD folder but they all use the method
>> of setting initial conditions using the solnData(DENS_VAR,i,j,k) array and not
>> the call Grid_putPointData(blockId, CENTER, DENS_VAR, EXTERIOR, axis, rhoZone)
>> function call as recommended in the user guide pdf. 
Which is the correct
>> method? 

> For difference between Grid_putPointData and Grid_getBlkPtr etc., see the
> answer (~ second half or response) here:
>
> http://flash.uchicago.edu/pipermail/flash-users/2016-December/002142.html
>
>> Also, there are no examples for initialising the face centred variables i.e.
>> face centred magnetic fields. Is this also done with call
>> Grid_putPointData(blockId, CENTER, DENS_VAR, EXTERIOR, axis, rhoZone) method
>> and is there an example? 

> One COULD initialize face variables with Grid_putPointData, but all
> existing examples use block pointers for this.
>
> To find examples, grep for "FACEX" in all files
>
>    Simulation/SimulationMain/magnetoHD/*/Simulation_initBlock.F90* .
>
>
> Klaus