<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hi Norbert,<div><br></div><div>I am not sure. I found it quite confusing, so to check things I thought I'd look at it in visit. The code that doesn't run defines the lasers using:</div><div><br></div><div><div>--------------------------------------------------------------------------------------------------------</div><div>--------------------------------------------------------------------------------------------------------</div><div><div><div>useEnergyDeposition = .true.</div><div>ed_maxRayCount = 100000</div><div>ed_gradOrder = 2</div><div><br></div><div># Activate 3D-in-2D ray trace:</div><div>ed_laser3Din2D = .true.</div><div>ed_laser3Din2DwedgeAngle = 0.5</div><div><br></div><div>### LASER IO OPTIONS ###</div><div>ed_useLaserIO = .true.</div><div>ed_laserIOMaxNumberOfPositions = 10000</div><div>ed_laserIOMaxNumberOfRays = 128</div><div><br></div><div>ed_numberOfPulses = 3</div><div>ed_numberOfBeams = 3</div></div><div><br></div><div># Define Beam 1: 2/7 of the power 2.84e12*(2/7)</div><div># approximate lens location for OMEGA (R_cham=3.33m) define angles relative to target normal</div><div>#---------- CONE 1 at 32 degrees, y=R*cos(theta), x=R*sin(theta)</div><div>ed_numberOfSections_1 = 4</div><div>ed_time_1_1 = 0.0</div><div>ed_time_1_2 = 0.1e-09</div><div>ed_time_1_3 = 1.1e-09</div><div>ed_time_1_4 = 1.2e-09</div><div># power can be adjusted to test cone dependence in the future</div><div>ed_power_1_1 = 0.0</div><div>ed_power_1_2 = 0.81e+12</div><div>ed_power_1_3 = 0.81e+12</div><div>ed_power_1_4 = 0.0</div><div># set up the beam profile</div><div>ed_lensX_1 = 176.5 # radial position</div><div>ed_lensY_1 = 282.4 # axial position</div><div>ed_lensZ_1 = 0.0</div><div>ed_lensSemiAxisMajor_1 = 14.925</div><div>ed_targetX_1 = 0.0e-04</div><div>ed_targetZ_1 = 0.0e-04</div><div>ed_targetY_1 = 0.0e-04</div><div># define major and minor radii at the target</div><div>ed_targetSemiAxisMajor_1 = 500.0e-04</div><div>ed_targetSemiAxisMinor_1 = 500.0e-04</div><div>ed_pulseNumber_1 = 1</div><div>ed_wavelength_1 = 0.351</div><div>ed_crossSectionFunctionType_1 = "gaussian2D"</div><div>ed_gaussianExponent_1 = 4.3</div><div>ed_gaussianRadiusMajor_1 = 250.0e-04</div><div>ed_gaussianRadiusMinor_1 = 250.0e-04</div><div>ed_numberOfRays_1 = 5000</div><div>ed_gridType_1 = "radial2D"</div><div>ed_semiAxisMajorTorsionAngle_1= 0.0</div><div>ed_semiAxisMajorTorsionAxis_1 = "x"</div><div><br></div><div><br></div><div># Define Beam 2: 3/7 of the power 2.84e12*(3/7)</div><div># approximate lens location for OMEGA (R_cham=3.33m) define angles relative to target normal</div><div>#-------------- CONE 2 at 42 degrees, y=R*cos(theta), x=R*sin(theta)</div><div>ed_numberOfSections_2 = 4</div><div>ed_time_2_1 = 0.0</div><div>ed_time_2_2 = 0.1e-09</div><div>ed_time_2_3 = 1.1e-09</div><div>ed_time_2_4 = 1.2e-09</div><div># power can be adjusted to test cone dependence in the future</div><div>ed_power_2_1 = 0.0</div><div>ed_power_2_2 = 1.22e+12</div><div>ed_power_2_3 = 1.22e+12</div><div>ed_power_2_4 = 0.0</div><div># set up the beam profile</div><div>ed_lensX_2 = 222.8 # radial position</div><div>ed_lensY_2 = 247.5 # axial position</div><div>ed_lensZ_2 = 0.0</div><div>ed_lensSemiAxisMajor_2 = 14.925</div><div>ed_targetX_2 = 0.0e-04</div><div>ed_targetZ_2 = 0.0e-04</div><div>ed_targetY_2 = 0.0e-04</div><div># define major and minor radii at the target</div><div>ed_targetSemiAxisMajor_2 = 500.0e-04</div><div>ed_targetSemiAxisMinor_2 = 500.0e-04</div><div># laser profile from with SG4 profile</div><div>ed_pulseNumber_2 = 2</div><div>ed_wavelength_2 = 0.351</div><div>ed_crossSectionFunctionType_2 = "gaussian2D"</div><div>ed_gaussianExponent_2 = 4.3</div><div>ed_gaussianRadiusMajor_2 = 250.0e-04</div><div>ed_gaussianRadiusMinor_2 = 250.0e-04</div><div>ed_numberOfRays_2 = 5000</div><div>ed_gridType_2 = "radial2D"</div><div>ed_semiAxisMajorTorsionAngle_2= 0.0</div><div>ed_semiAxisMajorTorsionAxis_2 = "x"</div><div><br></div><div><br></div><div># Define Beam 3: 2/7 of the power 2.84e12*(2/7)</div><div># approximate lens location for OMEGA (R_cham=3.33m) define angles relative to target normal</div><div>#------------- CONE 3 at 50 degrees, y=R*cos(theta), x=R*sin(theta)</div><div>ed_numberOfSections_3 = 4</div><div>ed_time_3_1 = 0.0</div><div>ed_time_3_2 = 0.1e-09</div><div>ed_time_3_3 = 1.1e-09</div><div>ed_time_3_4 = 1.2e-09</div><div># power can be adjusted to test cone dependence in the future</div><div>ed_power_3_1 = 0.0</div><div>ed_power_3_2 = 0.81e+12</div><div>ed_power_3_3 = 0.81e+12</div><div>ed_power_3_4 = 0.0</div><div># set up the beam profile</div><div>ed_lensX_3 = 255.1 # radial position</div><div>ed_lensY_3 = 214.0 # axial position</div><div>ed_lensZ_3 = 0.0</div><div>ed_lensSemiAxisMajor_3 = 14.925</div><div>ed_targetX_3 = 0.0e-04</div><div>ed_targetZ_3 = 0.0e-04</div><div>ed_targetY_3 = 0.0e-04</div><div># define major and minor radii at the target</div><div>ed_targetSemiAxisMajor_3 = 500.0e-04</div><div>ed_targetSemiAxisMinor_3 = 500.0e-04</div><div>ed_pulseNumber_3 = 3</div><div>ed_wavelength_3 = 0.351</div><div>ed_crossSectionFunctionType_3 = "gaussian2D"</div><div>ed_gaussianExponent_3 = 4.3</div><div>ed_gaussianRadiusMajor_3 = 250.0e-04</div><div>ed_gaussianRadiusMinor_3 = 250.0e-04</div><div>ed_numberOfRays_3 = 5000</div><div>ed_gridType_3 = "radial2D"</div><div>ed_semiAxisMajorTorsionAngle_3= 0.0</div><div>ed_semiAxisMajorTorsionAxis_3 = "x"</div><div><br></div><div><div>--------------------------------------------------------------------------------------------------------</div><div>--------------------------------------------------------------------------------------------------------</div></div></div></div><div><br></div><div>however when I change the lines:</div><div><br></div><div>ed_laser3Din2D = .false.</div><div>ed_crossSectionFunctionType_X = "gaussian1D"</div><div>ed_gridType_X = "regular1D"</div><div><br></div><div>The code runs and the attached image shows a sample output early in the simulation. The geometry at least looks correct in that there are 3 beams coming in at different angles onto the target. I don't quite understand why they don't seem to be centered on 0,0, but perhaps I am not saving enough rays for this to look correct.</div><div><br></div><div>For what its worth, I have also compiled the LaserSlab simulation on the system I'm using with the standard setup script in the 4.2.2 manual. When I run LaserSlab with the example.par setup, it reproduces the behavior by hanging up right after the output line 'Initial dt verified'. That to me suggests something more fundamental is going wrong on my system. </div><div><br></div><div>Any input is very much welcomed.</div><div><br></div><div>Cheers,</div><div>Mario</div><div><br><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">______________________________________</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"> <span class="Apple-converted-space"> </span></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Dr. Mario Manuel<br>Einstein Postdoctoral Fellow, Smithsonian Astrophysical Observatory<br><a href="mailto:mmanuel@umich.edu">mmanuel@umich.edu</a><br>617-899-8381<br><br>University of Michigan<br>Space Research Building<br>2455 Hayward Street<br>Ann Arbor, MI 48109<br>Office: 734-764-4585<br></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">______________________________________</div></div>
</div>
<img apple-inline="yes" id="90871976-4062-425C-9158-73F412DB901C" height="742" width="1024" apple-width="yes" apple-height="yes" src="cid:2E11C16D-78EC-4DBD-8EA1-665A826AF772@adsroot.itcs.umich.edu"><div><br></div><div>\</div><div><br><div><div>On Apr 21, 2016, at 4:15 PM, Norbert Flocke <<a href="mailto:flocke@flash.uchicago.edu">flocke@flash.uchicago.edu</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hi Manuel,<br><br>Are you sure you are setting up the beams correctly for this 3D in 2D simulation? The reason I am asking is because it is easy to get confused<br>between the cylindrical FLASH coordinates (x,y) -> (R,z) and the 3D beam coordinates you need to use. In your application (cylindrical) the FLASH y coordinate corresponds to the beam 3D z coordinate.<br><br>Best,<br>Norbert<br><br><br>On Thu, 21 Apr 2016, Mario Manuel wrote:<br><br><blockquote type="cite">Hello FLASH users,<br><br>I am running a FLASH4.2.2 simulation that is a combination of the LaserSlab and FieldLoop<br>example simulations. It is a laser-matter interaction (multi-material, 3T, tabulated EOSs) with MHD as well. I recently implemented a more complex laser setup with multiple beams and wanted to use the ed_laser3Din2D = .true. for a more accurate ray tracing in the R-Z geometry that I am using. After fixing all the initialization errors, I got it to run, but it never gets in to the first step. The simulation hangs up, but doesn't throw an error. The output looks like:<br><br><br>-------------------------------------------------------------------------------------------------------------------------------<br>-------------------------------------------------------------------------------------------------------------------------------<br>Wed Apr 20 20:57:51 PDT 2016<br>NODES= 40<br>NTASKS= 1280<br>[Driver_initParallel]: Called MPI_Init_thread - requested level 2, given level 2<br>[Driver_initParallel]: Number of OpenMP threads in each parallel region 2<br>NOTE: Enabling curvilinear, cartesian_pm/cylindrical_pm/spherical_pm/polar_pm so far was F T F F<br>MaterialProperties initialized<br>Cosmology initialized<br>[eos_tabBrowseIonmix4Tables] IONMIX4 file found: he-imx-005.cn4<br>in eos_inittabulated, tableName = he-imx-005.cn4<br>in eos_inittabulated, tableName = he-imx-005.cn4<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, tableName = he-imx-005.cn4<br>in eos_inittabulated, groupName = -none-<br>[eos_tabBrowseIonmix4Tables] IONMIX4 file found: polystyrene-imx-002.cn4<br>in eos_inittabulated, tableName = he-imx-005.cn4<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, tableName = polystyrene-imx-002.cn4<br>in eos_inittabulated, tableName = polystyrene-imx-002.cn4<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, tableName = polystyrene-imx-002.cn4<br>in eos_inittabulated, groupName = -none-<br>in eos_inittabulated, tableName = polystyrene-imx-002.cn4<br>in eos_inittabulated, groupName = -none-<br>Source terms initialized<br> iteration, no. not moved = 0 0<br>refined: total leaf blocks = 3<br>refined: total blocks = 3<br>INFO: Grid_fillGuardCells is ignoring masking.<br> iteration, no. not moved = 0 2<br> iteration, no. not moved = 1 0<br>refined: total leaf blocks = 12<br>refined: total blocks = 15<br> iteration, no. not moved = 0 13<br> iteration, no. not moved = 1 2<br> iteration, no. not moved = 2 0<br>refined: total leaf blocks = 48<br>refined: total blocks = 63<br> iteration, no. not moved = 0 60<br> iteration, no. not moved = 1 13<br> iteration, no. not moved = 2 0<br>refined: total leaf blocks = 192<br>refined: total blocks = 255<br> iteration, no. not moved = 0 251<br> iteration, no. not moved = 1 60<br> iteration, no. not moved = 2 0<br>refined: total leaf blocks = 768<br>refined: total blocks = 1023<br> iteration, no. not moved = 0 283<br> iteration, no. not moved = 1 101<br> iteration, no. not moved = 2 0<br>refined: total leaf blocks = 216<br>refined: total blocks = 287<br> iteration, no. not moved = 0 50<br> iteration, no. not moved = 1 12<br> iteration, no. not moved = 2 0<br>refined: total leaf blocks = 96<br>refined: total blocks = 127<br> Finished with Grid_initDomain, no restart<br>Ready to call Hydro_init<br>[Hydro_init] The upwind electric field construction is NOT allowed for diffusion.<br>Hydro initialized<br>Gravity initialized<br>Initial dt verified<br>*** Wrote checkpoint file to /p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_chk_0000 ****<br>*** Wrote plotfile to /p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_plt_cnt_0000 ****<br>Initial plotfile written<br>Driver init all done<br>srun: got SIGCONT<br>srun: forcing job termination<br>-------------------------------------------------------------------------------------------------------------------------------<br>-------------------------------------------------------------------------------------------------------------------------------<br><br>and the last lines of the log file look like:<br><br>-------------------------------------------------------------------------------------------------------------------------------<br>-------------------------------------------------------------------------------------------------------------------------------<br>[ 04-20-2016 21:18:56.632 ] [GRID amr_refine_derefine]: redist. phase. tot blks requested: 287<br>[GRID amr_refine_derefine] min blks 0 max blks 1 tot blks 287<br>[GRID amr_refine_derefine] min leaf blks 0 max leaf blks 1 tot leaf blks 216<br>[ 04-20-2016 21:18:56.654 ] [GRID amr_refine_derefine]: refinement complete<br>[ 04-20-2016 21:18:56.847 ] [GRID gr_expandDomain]: iteration=6, create level=5<br>[ 04-20-2016 21:18:56.891 ] [GRID amr_refine_derefine]: initiating refinement<br>[ 04-20-2016 21:18:56.896 ] [GRID amr_refine_derefine]: redist. phase. tot blks requested: 127<br>[GRID amr_refine_derefine] min blks 0 max blks 1 tot blks 127<br>[GRID amr_refine_derefine] min leaf blks 0 max leaf blks 1 tot leaf blks 96<br>[ 04-20-2016 21:18:56.919 ] [GRID amr_refine_derefine]: refinement complete<br>[ 04-20-2016 21:18:57.110 ] [GRID gr_expandDomain]: iteration=7, create level=5<br>[ 04-20-2016 21:18:57.116 ] [gr_ensureValidNeighborInfo] found mpi_pattern_id: -10<br>[ 04-20-2016 21:18:57.124 ] [mpi_amr_comm_setup]: buffer_dim_send=157, buffer_dim_recv=1<br>[ 04-20-2016 21:18:57.130 ] [gr_hypreInit]: HYPRE_RELEASE_NAME=hypre, HYPRE_RELEASE_VERSION=2.8.0b, HYPRE_RELEASE_DATE=2011/11/14<br>[ 04-20-2016 21:18:57.136 ] memory: bg heap use (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:57.141 ] memory: bg heap peak (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:57.145 ] memory: bg heap free (MB): 324.20 (min) 332.95 (max) 327.56 (avg)<br>[ 04-20-2016 21:18:57.155 ] memory: bg heap use (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:57.169 ] memory: bg heap peak (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:57.263 ] memory: bg heap free (MB): 324.20 (min) 332.95 (max) 327.56 (avg)<br>[ 04-20-2016 21:18:58.199 ] memory: bg heap use (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:58.205 ] memory: bg heap peak (MB): 126.46 (min) 135.21 (max) 131.85 (avg)<br>[ 04-20-2016 21:18:58.209 ] memory: bg heap free (MB): 324.26 (min) 333.01 (max) 327.63 (avg)<br>[ 04-20-2016 21:18:58.407 ] [IO_writeCheckpoint] open: type=checkpoint name=/p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_chk_0000<br>[ 04-20-2016 21:18:59.943 ] [io_writeData]: wrote 127 blocks<br>[ 04-20-2016 21:19:00.596 ] [IO_writeCheckpoint] close: type=checkpoint name=/p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_chk_0000<br>[ 04-20-2016 21:19:00.682 ] [IO_writePlotfile] open: type=plotfile name=/p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_plt_cnt_0000<br>[ 04-20-2016 21:19:01.269 ] [io_writeData]: wrote 127 blocks<br>[ 04-20-2016 21:19:01.695 ] [IO_writePlotfile] close: type=plotfile name=/p/lscratchv/manuel3/ACSEL15B_CH/ACSEL15B_hdf5_plt_cnt_0000<br>[ 04-20-2016 21:19:01.706 ] memory: bg heap use (MB): 134.27 (min) 135.21 (max) 134.74 (avg)<br>[ 04-20-2016 21:19:01.711 ] memory: bg heap peak (MB): 134.27 (min) 151.21 (max) 134.77 (avg)<br>[ 04-20-2016 21:19:01.716 ] memory: bg heap free (MB): 324.26 (min) 325.20 (max) 324.73 (avg)<br>[ 04-20-2016 21:19:01.722 ] [Driver_evolveFlash]: Entering evolution loop<br>[ 04-20-2016 21:19:01.727 ] step: n=1 t=0.000000E+00 dt=1.000000E-15<br>[ 04-20-2016 21:19:01.746 ] [mpi_amr_comm_setup]: buffer_dim_send=27933, buffer_dim_recv=1<br>[ 04-20-2016 21:19:01.783 ] [hy_uhd_unsplit]: gcNeed(MAGI_FACE_VAR,MAG_FACE_VAR) - FACES<br>[ 04-20-2016 21:19:01.792 ] [mpi_amr_comm_setup]: buffer_dim_send=2333, buffer_dim_recv=1<br>[ 04-20-2016 21:19:01.886 ] [mpi_amr_comm_setup]: buffer_dim_send=27933, buffer_dim_recv=1<br>[ 04-20-2016 21:19:01.913 ] [mpi_amr_comm_setup]: buffer_dim_send=25757, buffer_dim_recv=1<br>-------------------------------------------------------------------------------------------------------------------------------<br>-------------------------------------------------------------------------------------------------------------------------------<br><br>Everything runs fine when I switch the flag to ed_laser2Din3D = .false.<br><br>Has anyone come across this type of thing before? Are there issues with the 2Din3D package and some other packages I may have implemented? Any thoughts are much welcomed.<br><br>Cheers,<br>Mario<br><br>______________________________________<br><br>Dr. Mario Manuel<br>Einstein Postdoctoral Fellow, Smithsonian Astrophysical Observatory<br><a href="mailto:mmanuel@umich.edu">mmanuel@umich.edu</a><br>617-899-8381<br><br>University of Michigan<br>Space Research Building<br>2455 Hayward Street<br>Ann Arbor, MI 48109<br>Office: 734-764-4585<br>______________________________________<br><br><br></blockquote></blockquote></div><br></div></div></body></html>