[FLASH-USERS] T electron error from 3D laser slab example

[Norbert Flocke]

Ok. Great to know. If you have any more questions about the laser feel 
free to contact me.

On Thu, 24 Oct 2013, Tommasini, Riccardo wrote:

> Hi Norbert,
>
> Great tip and good catch! That did it.
> Thanks a lot.
>
> Best Regards,
>
> Riccardo
>
>
> -----Original Message-----
> From: Norbert Flocke <flocke at flash.uchicago.edu>
> Date: Thursday, October 24, 2013 6:37 AM
> To: Riccardo Tommasini <tommasini2 at llnl.gov>
> Cc: "flash-users at flash.uchicago.edu" <flash-users at flash.uchicago.edu>
> Subject: Re: [FLASH-USERS] T electron error from 3D laser slab example
>
>> Hi Riccardo,
>>
>> What you need to do is the following: in the flash.par file for the
>> LaserSlab add the lines
>>
>>   ed_enforcePositiveTele = .true.
>>   ed_enforcePositiveNele = .true.
>>
>> Probably the FLASH version you have has these runtime parameters set as
>> false as a default. Let me know if that solves the problem.
>>
>> Regards,
>> Norbert
>>
>>
>> On Wed, 23 Oct 2013, Tommasini, Riccardo wrote:
>>
>>> Hi All,
>>>
>>>
>>>
>>> I am new to FLASH4 therefore I apologize in advance if this issue has
>>> been already discussed ( I  was not able to find it in the FLASH user
>>> archives though).
>>>
>>> I had no problems running the Laser Slab example in 2D, however when I
>>> run the 3D case  the simulation ended up with the error reported below,
>>> indicating negative electron temperature was generated.
>>>
>>> "
>>> *** Wrote checkpoint file to lasslab_hdf5_chk_0007 ****
>>> Š.
>>> Š.
>>> Driver_abort called. See log file for details.
>>> Error message is [ed_traceBlockRays3DRec] ERROR: Tele <= 0 for a cell
>>> (target)
>>> Calling MPI_Abort() for shutdown in   2 seconds!
>>>
>>> application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
>>> "
>>>
>>> The setup command I used for the simulation is:
>>>
>>>
>>> ./setup.py -auto LaserSlab -noc -3d +pm4dev +hdf5typeio -nxb=16 -nyb=16
>>> -nzb=16 +mtmmmt +laser +mgd +uhd3t species=cham,targ mgd_meshgroups=6
>>> -parfile=example3d.par -objdir=obj_LaserSlab3D
>>>
>>> None of the files has been modified from the original distribution.
>>>
>>> Any ideas about what I could be missing?
>>>
>>> Thanks,
>>>
>>> Riccardo
>