[FLASH-USERS] Issue with Abundances in particle file, using FLASH4.5

[Ji, Suoqing]

Hi Vishal,

Hope you’re doing well!

The default number of particle attributes is 10, so I guess that’s why you cannot get proper values beyond 10 fluid-related quantities. You might try to increase the default number of particle attributes, by adding a line to your Config file:

PPDEFINE PT_MAX_ATTRIBUTES 29

Since you have a total of 29 particle attributes.

Best regards,
—
Suoqing Ji
Sherman Fairchild Fellow
TAPIR & Walter Burke Institute for Theoretical Physics
California Institute of Technology
http://www.tapir.caltech.edu/~suoqing

On Aug 18, 2020, at 3:30 PM, Vishal Tiwari <vtiwari at umassd.edu<mailto:vtiwari at umassd.edu>> wrote:

Hi all,

We are facing some issues with the abundances of the tracer particles in our setup (using FLASH4.5). All the abundances for the aprox19 nuclides in all "part" files are coming out to "-1.0". But, the hydro quantities like density, temperature, velocities are correct as expected. (We are using hd5dump to look at the tracer particle abundances from the "part" file generated by the setup/run.)

The "plt" files have the right abundances of the nuclides.

Please find the Config file relevant to the particle setup below. I kindly request you if you could help in figuring out the issue. Thank you!


USESETUPVARS withParticles

IF withParticles

  REQUIRES Grid/GridParticles/
  REQUESTS Grid/GridParticles/GridParticlesMove/Sieve/BlockMatch

  PARTICLETYPE passive_m INITMETHOD CellMass MAPMETHOD WEIGHTED ADVMETHOD rungekutta
  REQUIRES Particles/ParticlesMain
  REQUESTS Particles/ParticlesMapping/meshWeighting/CIC  REQUESTS Particles/ParticlesInitialization/WithDensity/CellMassBins   # by density
  REQUESTS IO/IOParticles/hdf5/parallel

PARTICLEPROP    dens            REAL
PARTICLEPROP    temp            REAL
PARTICLEPROP    eint            REAL
PARTICLEPROP    ener            REAL
PARTICLEPROP    enuc            REAL
PARTICLEPROP    gamc            REAL
PARTICLEPROP    game            REAL
PARTICLEPROP    gpol            REAL
PARTICLEPROP    gpot            REAL
PARTICLEPROP    pres            REAL
#PARTICLEPROP MOMX REAL
PARTICLEPROP    h1              REAL
PARTICLEPROP    he3             REAL
PARTICLEPROP    n14             REAL
PARTICLEPROP    fe54            REAL
PARTICLEPROP    neut            REAL
PARTICLEPROP    prot            REAL
PARTICLEPROP    ar36            REAL
PARTICLEPROP    c12             REAL
PARTICLEPROP    ca40            REAL
PARTICLEPROP    cr48            REAL
PARTICLEPROP    fe52            REAL
PARTICLEPROP    he4             REAL
PARTICLEPROP    mg24            REAL
PARTICLEPROP    ne20            REAL
PARTICLEPROP    ni56            REAL
PARTICLEPROP    o16             REAL
PARTICLEPROP    s32             REAL
PARTICLEPROP    si28            REAL
PARTICLEPROP    ti44            REAL

PARTICLEMAP  TO dens  FROM VARIABLE   dens
PARTICLEMAP  TO temp  FROM VARIABLE   temp
PARTICLEMAP  TO eint  FROM VARIABLE   eint
PARTICLEMAP  TO ener  FROM VARIABLE   ener
PARTICLEMAP  TO enuc  FROM VARIABLE   enuc
PARTICLEMAP  TO gamc  FROM VARIABLE   gamc
PARTICLEMAP  TO game  FROM VARIABLE   game
PARTICLEMAP  TO gpol  FROM VARIABLE   gpol
PARTICLEMAP  TO gpot  FROM VARIABLE   gpot
PARTICLEMAP  TO pres  FROM VARIABLE   pres

# particle map of species
PARTICLEMAP  TO h1    FROM SPECIES   H1
PARTICLEMAP  TO he3   FROM SPECIES   HE3
PARTICLEMAP  TO n14   FROM SPECIES   N14
PARTICLEMAP  TO fe54  FROM SPECIES   FE54
PARTICLEMAP  TO neut  FROM SPECIES   NEUT
PARTICLEMAP  TO prot  FROM SPECIES   PROT
PARTICLEMAP  TO ar36  FROM SPECIES   AR36
PARTICLEMAP  TO c12   FROM SPECIES   C12
PARTICLEMAP  TO ca40  FROM SPECIES   CA40
PARTICLEMAP  TO cr48  FROM SPECIES   CR48
PARTICLEMAP  TO fe52  FROM SPECIES   FE52
PARTICLEMAP  TO he4   FROM SPECIES   HE4
PARTICLEMAP  TO mg24  FROM SPECIES   MG24
PARTICLEMAP  TO ne20  FROM SPECIES   NE20
PARTICLEMAP  TO ni56  FROM SPECIES   NI56
PARTICLEMAP  TO o16   FROM SPECIES   O16
PARTICLEMAP  TO s32   FROM SPECIES   S32
PARTICLEMAP  TO si28  FROM SPECIES   SI28
PARTICLEMAP  TO ti44  FROM SPECIES   TI44

PARAMETER particle_attribute_1 STRING "dens"
PARAMETER particle_attribute_2 STRING "temp"
PARAMETER particle_attribute_3 STRING "eint"
PARAMETER particle_attribute_4 STRING "ener"
PARAMETER particle_attribute_5 STRING "enuc"
PARAMETER particle_attribute_6 STRING "gamc"
PARAMETER particle_attribute_7 STRING "game"
PARAMETER particle_attribute_8 STRING "gpol"
PARAMETER particle_attribute_9 STRING "gpot"
PARAMETER particle_attribute_10 STRING "pres"

PARAMETER particle_attribute_11 STRING "H1"
PARAMETER particle_attribute_12 STRING "HE4"
PARAMETER particle_attribute_13 STRING "N14"
PARAMETER particle_attribute_14 STRING "FE54"
PARAMETER particle_attribute_15 STRING "NEUT"
PARAMETER particle_attribute_16 STRING "PROT"
PARAMETER particle_attribute_17 STRING "AR36"
PARAMETER particle_attribute_18 STRING "C12"
PARAMETER particle_attribute_19 STRING "CA40"
PARAMETER particle_attribute_20 STRING "CR48"
PARAMETER particle_attribute_21 STRING "FE52"
PARAMETER particle_attribute_22 STRING "HE4"
PARAMETER particle_attribute_23 STRING "MG24"
PARAMETER particle_attribute_24 STRING "NE20"
PARAMETER particle_attribute_25 STRING "NI56"
PARAMETER particle_attribute_26 STRING "O16"
PARAMETER particle_attribute_27 STRING "S32"
PARAMETER particle_attribute_28 STRING "SI28"
PARAMETER particle_attribute_29 STRING "TI44"


ENDIF

Best,
Vishal Tiwari
Physics Department
University of Massachusetts, Dartmouth
https://vishaltiwari.github.io/

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20200819/3aeea574/attachment-0001.htm>