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

[Klaus Weide]

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