Hi Rob:<br><br>Yes, if you have defined two species with the names FOIL and GAS then there will be two variables called FOIL_SPEC and GAS_SPEC which represent the mass fractions of each species in each cell. In the Simulation_initBlock.F90 file, you just have to set these two variables such that they sum to one. If you create the species using the "species=foil,gas" setup variable, you can set the properties of each species (A, Z, etc...) in the flash.par file. Follow the link below for a description:<br>

<br><a href="http://flash.uchicago.edu/site/flashcode/user_support/flash_ug_devel/node17.html#SECTION05440000000000000000">http://flash.uchicago.edu/site/flashcode/user_support/flash_ug_devel/node17.html#SECTION05440000000000000000</a><br>

<br>This link goes through an example of how to set the species properties in the flash.par file:<br><br><a href="http://flash.uchicago.edu/site/flashcode/user_support/flash_ug_devel/node34.html#SECTION08175000000000000000">http://flash.uchicago.edu/site/flashcode/user_support/flash_ug_devel/node34.html#SECTION08175000000000000000</a><br>

<br>Once you set the values of FOIL_SPEC and GAS_SPEC in the Simulation_initBlock.F90 file and the species properties in the flash.par file, then the EOS unit itself will take care of setting things like YE and SUMY for you - if you use an appropriate EOS, such as Multigamma or MTMMMT.<br>

<br>Is that helpful?<br><br>Thanks<br>Milad<br><br><div class="gmail_quote">On Wed, Aug 29, 2012 at 12:27 PM, Robert Crowston <span dir="ltr"><<a href="mailto:rhc504@york.ac.uk" target="_blank">rhc504@york.ac.uk</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi, when using the unstaggered mesh MHD solver, is it possible to specify that different species start in different places in space? I am trying to model the explosion of an aluminium (say) foil into a magnetised hydrogen gas. I've been looking at the LaserSlab example problem to understand how to configure different species. As I understand it, the calls to Grid_putPointData( ..., ..., YE_MSCALAR, ..., zbar/abar) and similarly with SUMY_MSCALAR are setting up the mass fractions in these blocks. Is that right? But I don't know how to get these constants defined with USM.<br>



<br>Thanks for your help! <br>Rob Crowston<br>
</blockquote></div><br>