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

[Tommasini, Riccardo]

I'll do. 
Thanks again,

Riccardo

-----Original Message-----
From: Norbert Flocke <flocke at flash.uchicago.edu>
Date: Thursday, October 24, 2013 8:58 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

>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