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

[Tommasini, Riccardo]

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