From andrew.szymkowiak at yale.edu  Thu Oct  1 09:14:34 2009
From: andrew.szymkowiak at yale.edu (Andrew Szymkowiak)
Date: Thu, 01 Oct 2009 10:14:34 -0400
Subject: [FLASH-USERS] Problems with refinement - nevermind...
Message-ID: <4AC4B94A.4040602@yale.edu>

Sorry, I found my error (for those that might see something similar, it 
was not due to me copying the outflow into wrong regions, but from BC 
routine returning before setting the guardcells for all the variables in 
the regions that did not include the outflow - evidently the "old" 
values in the reused blocks would sometimes contain the values from 
under the outflow).


From mateuszr at umich.edu  Mon Oct  5 17:03:21 2009
From: mateuszr at umich.edu (Mateusz Ruszkowski)
Date: Mon, 5 Oct 2009 18:03:21 -0400 (EDT)
Subject: [FLASH-USERS] H5Screate_simple: zero sized dimension
Message-ID: <Pine.LNX.4.64.0910051758030.19545@lem.astro.lsa.umich.edu>



  Hi all,

[apologies if you are getting this message more then once; we had some 
trouble sending the e-mail]

We are getting a runtime error related to hdf5. The crash occurs the 
moment the code is attempting to write the initial conditions. The 
standard output includes the following error messages:

  .
  .
  .
Initial dt verified
Driver_abortC called
Error: negative return from memspace H5Screate_simple
Calling MPI_Abort for immediate shutdown

and:

  .
  .
  .
HDF5-DIAG: Error detected in HDF5 library version 1.6.6-snap21 MPI-process 
0. Back trace follows.
#000: H5S.c line 1595 in H5Screate_simple(): zero sized dimension 
for non-unlimited dimension
major(01): Function arguments
minor(05): Bad value
MPI: On host columbia1, Program 
/u/.realmounts/home7/nieusma/Flash/Planet_IR, Rank 0, Process 8697
1 called MPI_Abort(<communicator>, 1)


A similar error was reported on the FLASH mailing list in the past and a 
suggestion was made that this could be due to (min refinement) > (max 
refinement) in flash.par. We checked the refinement levels and this does 
not appear to be a problem in the above run so the cause of this problem 
is probably different in our case. Did anybody encounter similar issues?


     thanks,
       Mateusz & Justin




From bastian.arnold at univie.ac.at  Wed Oct  7 03:46:00 2009
From: bastian.arnold at univie.ac.at (Bastian Arnold)
Date: Wed, 7 Oct 2009 10:46:00 +0200
Subject: [FLASH-USERS] unrecognized changes to eos routine
Message-ID: <f838c9eea33bfd9716b8fdd840e1613d.squirrel@webmail.univie.ac.at>

Dear Flash users,

I am Bastian, relatively new to Flash, and I am using version 3.1. There
is the following problem I cannot entangle so far:

I have set up a gas cloud (with several species, i.e. including
multigamma) in pressure equilibrium with the ambient medium. For that
purpose, I needed to change the Eos.F90 routine, such that a
temperature-dependent mean molecular weight is calculated by an external
routine. Everything is fine after Simulation_initBlock.F90 has finished
and BEFORE Eos_wrapped.F90 is applied to the initial block. But
afterwards, i.e. from the initial plot file on, the pressure is not
consistent with a temperature-dependent mean molec. weight, but instead as
it would be after applying the default Eos.F90 routine. It seems to me,
that Eos_wrapped.F90 (which calls Eos.F90) does not recognize my changes
to Eos.F90.

To overcome the problem I have commented the call to Eos_wrapped.F90 in
routine gr_expandDomain.F90 (since I have AMR enabled). Now the pressure
values in the first plot file are fine, but from the second plot file on
the problem is still present.

So can anyone provide me an idea how I can make a consistent change to
Eos.F90, such that Eos_wrapped.F90 realizes this change?

Thanks in advance for any helpful guess!

Regards,
Bastian





From klaus at flash.uchicago.edu  Wed Oct  7 08:34:26 2009
From: klaus at flash.uchicago.edu (Klaus Weide)
Date: Wed, 7 Oct 2009 08:34:26 -0500 (CDT)
Subject: [FLASH-USERS] unrecognized changes to eos routine
In-Reply-To: <f838c9eea33bfd9716b8fdd840e1613d.squirrel@webmail.univie.ac.at>
References: <f838c9eea33bfd9716b8fdd840e1613d.squirrel@webmail.univie.ac.at>
Message-ID: <Pine.LNX.4.64.0910070828350.9568@flash.uchicago.edu>

On Wed, 7 Oct 2009, Bastian Arnold wrote:

> Dear Flash users,
> 
> I am Bastian, relatively new to Flash, and I am using version 3.1. There
> is the following problem I cannot entangle so far:
> 
> I have set up a gas cloud (with several species, i.e. including
> multigamma) in pressure equilibrium with the ambient medium. For that
> purpose, I needed to change the Eos.F90 routine, such that a
> temperature-dependent mean molecular weight is calculated by an external
> routine. Everything is fine after Simulation_initBlock.F90 has finished
> and BEFORE Eos_wrapped.F90 is applied to the initial block. But
> afterwards, i.e. from the initial plot file on, the pressure is not
> consistent with a temperature-dependent mean molec. weight, but instead as
> it would be after applying the default Eos.F90 routine. It seems to me,
> that Eos_wrapped.F90 (which calls Eos.F90) does not recognize my changes
> to Eos.F90.

Dear Bastian,

This sounds like your Eos.F90 may not even be linked into the object 
directory, so that you are still compiling a different version of Eos into
the executable.

Please check that the Eos.F90 *in the object directory* (i.e., where you
'make' FLASH) is actually the one that you intended.  If not, you may have
to run setup again.  (For convenience, try "./resetup -noclobber" in the 
object directory.)


Klaus


From gaspmax at hotmail.com  Mon Oct 12 04:32:46 2009
From: gaspmax at hotmail.com (Massimo Gaspari)
Date: Mon, 12 Oct 2009 11:32:46 +0200
Subject: [FLASH-USERS] Slice interpolation and integration
Message-ID: <COL106-W225448ECBE39CD026E3A5DBEC80@phx.gbl>


Dear Flash users,

I am Massimo and I am using Flash since 3.1 (now 3.2). I'm working on AGN jet-outflows 3D hydro-simulations in galaxy clusters and groups. 

I have two questions:

 

1) I made some 2D slices of the outputs using Quickflash, but I noticed that the final map (which is at max "refine level") is not interpolated; so I'm seeing all the AMR level squares. Is there a method/flag/option to produce interpolated maps? Has someone implemented this feature yet?

I've written some code lines to interpolate the data, but cell-centered data cannot be "caught" on 2D slices, so I've to do 3D trilinear interpolation (takes 8 hours just to produce a map), searching through all 3D (available) points for the eight corners, for every level. Is there a smarter way?

 

2) I want also to generate surface brightness maps, so I have to integrate all the slices (dens and temp) along line of sight (let's say x-direction). In this case, if I take non interpolated slices the results are very bad, so I really (?) need to interpolate. Again, has someone done something in this direction, producing good SB maps (with QF or other codes)?

 

Thank you in advance.

Regards,

Massimo
 		 	   		  
_________________________________________________________________
Gossip, Sport, Notizie... Accendi la Messenger TV!
http://www.messenger.it/messengerTV.aspx
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://flash.uchicago.edu/pipermail/flash-users/attachments/20091012/b047f685/attachment.html 

From alexander.y.wagner at googlemail.com  Mon Oct 19 07:10:25 2009
From: alexander.y.wagner at googlemail.com (Alexander Wagner)
Date: Mon, 19 Oct 2009 23:10:25 +1100
Subject: [FLASH-USERS] SPECIES_FLUX_BEGIN and some RHD questions
Message-ID: <4ADC5731.1050509@gmail.com>

Dear FLASH users,

Certain lines in my Flash.h file generated in FLASH3.2 read

#define SPECIES_FLUX_BEGIN (PROP_FLUX_END + CONSTANT_ONE)
#define SPECIES_FLUX_END (FLUX_SPECIES_BEGIN + NSPECIES_FLUX - CONSTANT_ONE)

Isn't FLUX_SPECIES_BEGIN in the second line supposed to be 
SPECIES_FLUX_BEGIN? I couldn't find anything on this in the user mail 
archives.

Also, does anyone know whether the update for species in 
hy_rhd_sweep.F90 works? (That long expression that includes vint).

Is the RHD unit currently being developed?

Many thanks,
Alex




From alexander.y.wagner at googlemail.com  Tue Oct 20 19:38:04 2009
From: alexander.y.wagner at googlemail.com (Alexander Wagner)
Date: Wed, 21 Oct 2009 11:38:04 +1100
Subject: [FLASH-USERS] SPECIES_FLUX_END and some RHD questions
Message-ID: <4ADE57EC.9080307@gmail.com>

Dear FLASH users,

Certain lines in my Flash.h file generated in FLASH3.2 read

#define SPECIES_FLUX_BEGIN (PROP_FLUX_END + CONSTANT_ONE)
#define SPECIES_FLUX_END (FLUX_SPECIES_BEGIN + NSPECIES_FLUX - 
CONSTANT_ONE)

Isn't FLUX_SPECIES_BEGIN in the second line supposed to be 
SPECIES_FLUX_BEGIN? I changed it in the template file and it works.

Another thing I noticed in unitCfg.py, is that the line

return {}, 
r'MASS_SCALAR\s+(?P<msname>\w+)(\s+)?(?:(NORENORM|RENORM:(?P<gpname>\S*)|' 
+ strEosmapRE + ')\s*)*\s*$'

in def initparseMASS_SCALAR(self):, I think, should read

return {}, 
r'MASS_SCALAR\s+(?P<msname>\w+)(\s+)?(?:(NORENORM|RENORM:\s*(?P<gpname>\S*)|' 
+ strEosmapRE + ')\s*)*\s*$'

that is, there is a \s* missing after RENORM:


RHD questions:

Does anyone know whether the update for species in hy_rhd_sweep.F90 
works? (That long expression that includes vint).

Is the RHD unit currently being developed?

I also noticed in the manual that the runtime parameter for the 
reconstruction type in the RHD unit is listed as being called 
"rhd_reconType", but the one read in in RHD/Hydro_init.F90 is named 
"reconType".

Many thanks,
Alex




From flash-users at flash.uchicago.edu  Wed Oct 28 07:56:26 2009
From: flash-users at flash.uchicago.edu (VIAGRA Official Site)
Date: Wed, 28 Oct 2009 07:56:26 -0500 (CDT)
Subject: [FLASH-USERS] [SPAM] Dear flash-users@flash.uchicago.edu 56% 0FF on
	Pfizer !
Message-ID: <20091028125626.E8F52282C13D@flash.uchicago.edu>

An HTML attachment was scrubbed...
URL: http://flash.uchicago.edu/pipermail/flash-users/attachments/20091028/e3cf9455/attachment.html