From sgomez at cab.cnea.gov.ar  Tue Sep 11 01:13:18 2018
From: sgomez at cab.cnea.gov.ar (=?UTF-8?Q?Santiago_G=C3=B3mez?=)
Date: Mon, 10 Sep 2018 20:13:18 -0300
Subject: [mcstas-users] rssa source file
Message-ID: <426785d3030a25e2ca4f71b9882f7a20@cab.cnea.gov.ar>

Hello

I am trying to read a rssa source file (generated with SSW card in MCNP5 
v1.6) with the component Virtual_mcnp_ss_input.comp using the instrument 
file Test_SSW_SSR.instr, but i get this error message:

Instrument parameters for test_SSW_SSR (source.instr)
[test_SSW_SSR] Initialize
MCNP_input: testread: opening MCNP SSW file 'rssa' )
At line 63 of file neutronics-subs.f (unit = 123, file = 'rssa')
Fortran runtime error: I/O past end of record on unformatted file


I also get the same error with a rssa source file generate by MCNP6 
v1.0.


I did the compilation of the neutronics-subs.o file with gfortran 
version 4.8.3.


Can any one help?


I saw in the maillis a similar issue reported with 
Virtual_mcnp_ss_input.comp and fortran compiler, and a possible solution 
proposed by Peter Willendrup was to convert SSW to MCPL.

I tried this and it is working good, but I would like know if it is 
better to keep using converted MCPL than SSW MCNP source files?


Best regards


-- 
Santiago Miguel G?mez
Departamento de F?sica de Reactores y Radiaciones
Centro At?mico Bariloche
Int.: 5971

"Si queremos un mundo de paz y de justicia hay que poner decididamente 
la
inteligencia al servicio del amor." Antoine de Saint-Exupery

From pkwi at fysik.dtu.dk  Tue Sep 11 16:23:09 2018
From: pkwi at fysik.dtu.dk (=?utf-8?B?UGV0ZXIgS2rDpnIgV2lsbGVuZHJ1cA==?=)
Date: Tue, 11 Sep 2018 14:23:09 +0000
Subject: [mcstas-users] rssa source file
In-Reply-To: <426785d3030a25e2ca4f71b9882f7a20@cab.cnea.gov.ar>
References: <426785d3030a25e2ca4f71b9882f7a20@cab.cnea.gov.ar>
Message-ID: <2967651A-E764-4127-8CF7-35AB9EE2995F@fysik.dtu.dk>

Dear Santiago,


First of all good to hear that the MCPL solution seems to work for you.

As I remember it could be the case that the _mcnp_ss_ components implement SSR/SSW I/O wrt. MCNPX - which could easily explain that the record sizes do not match up. (Esben Klinkby (cc:) may know the full details.)

In any case, I would indeed suggest a complete move to MCPL as that code makes diagnostics on the file format to ensure correct import of MCNP5, MCNP6 and MCNPX.

(I will even consider obsoleting the _mcnp_ss_ component files for next McStas release)


Best,

Peter Willendrup



Peter Kj?r Willendrup
Forskningsingeni?r, Speciakonsulent
N?stformand for DTU Fysik LSU



DTU Physics

[cid:image001.gif at 01CCCAF1.5E6331F0]

Technical University of Denmark



[cid:image002.gif at 01CCCAF1.5E6331F0]


Department of Physics
Fysikvej

Building 307

DK-2800 Kongens Lyngby

Direct +45 2125 4612

Mobil +45 2125 4612
Fax +45 4593 2399


pkwi at fysik.dtu.dk<mailto:pkwi at fysik.dtu.dk>



On 11 Sep 2018, at 01.13, Santiago G?mez <sgomez at cab.cnea.gov.ar<mailto:sgomez at cab.cnea.gov.ar>> wrote:

I am trying to read a rssa source file (generated with SSW card in MCNP5 v1.6) with the component Virtual_mcnp_ss_input.comp using the instrument file Test_SSW_SSR.instr, but i get this error message:

Instrument parameters for test_SSW_SSR (source.instr)
[test_SSW_SSR] Initialize
MCNP_input: testread: opening MCNP SSW file 'rssa' )
At line 63 of file neutronics-subs.f (unit = 123, file = 'rssa')
Fortran runtime error: I/O past end of record on unformatted file


I also get the same error with a rssa source file generate by MCNP6 v1.0.


I did the compilation of the neutronics-subs.o file with gfortran version 4.8.3.


Can any one help?


I saw in the maillis a similar issue reported with Virtual_mcnp_ss_input.comp and fortran compiler, and a possible solution proposed by Peter Willendrup was to convert SSW to MCPL.

I tried this and it is working good, but I would like know if it is better to keep using converted MCPL than SSW MCNP source files?


Best regards


--
Santiago Miguel G?mez
Departamento de F?sica de Reactores y Radiaciones
Centro At?mico Bariloche
Int.: 5971

"Si queremos un mundo de paz y de justicia hay que poner decididamente la
inteligencia al servicio del amor." Antoine de Saint-Exupery
_______________________________________________
mcstas-users mailing list
mcstas-users at mcstas.org<mailto:mcstas-users at mcstas.org>
https://mailman2.mcstas.org/mailman/listinfo/mcstas-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180911/65165c71/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.gif
Type: image/gif
Size: 58 bytes
Desc: image001.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180911/65165c71/attachment-0002.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.gif
Type: image/gif
Size: 1055 bytes
Desc: image002.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180911/65165c71/attachment-0003.gif>

From ruvinskaia.e.e at mail.ru  Sun Sep 16 21:56:25 2018
From: ruvinskaia.e.e at mail.ru (=?UTF-8?B?0JXQutCw0YLQtdGA0LjQvdCwICDQoNGD0LLQuNC90YHQutCw0Y8=?=)
Date: Sun, 16 Sep 2018 22:56:25 +0300
Subject: [mcstas-users] =?utf-8?q?elliptic_mirror_problem?=
Message-ID: <1537127785.775917555@f405.i.mail.ru>

Dear all,

I have a problem. I need to neutron beam pass through elliptical mirrors system. I tried approach through component "Guide_tapering" and component "Mirror_Elliptic". But in both cases neutron beam does not move after the mirror. Could you say what the source of the problem and how I can fix this. Code with "Guide_tapering" for example are attached.

Thanking in advance.
Yours faithfully, Ekaterina .


DEFINE INSTRUMENT NSE_res(double slit_w = 0.1)


DECLARE
%{

double slit_h;
int iter;
double arm_lenght = 5;
%}


INITIALIZE
%{

%}


TRACE


COMPONENT Origin = Progress_bar()
? AT (0,0,0) ABSOLUTE

COMPONENT source =? Source_gen(
?? ?yheight=0.00001,xwidth=slit_w, lambda0=9, dlambda=0.000000000000001,I1=1e13,
?? ?verbose=1,focus_xw=slit_w,focus_yh=0.0001, dist = arm_lenght/4)
?? ?AT (0, 0, 0)? RELATIVE Origin
? EXTEND
? %{
???? t = rand0max(1e-20); // set time from 0 to 1 ms for TOF instruments.
? %}
?? ?

COMPONENT mon_TOF_source = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "time, limits=[0.000000001 0.3], bins=180, lambda, limits=[8.9 9.1], bins=180")
?? ?AT (0, 0, 0.00001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_source = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "lambda, limits = [8.99, 9.01], bins = 180")
?? ?AT (0, 0, 0.00001) RELATIVE PREVIOUS


COMPONENT first_arm = Arm()
?? ?AT (0, 0, 0.000001) RELATIVE PREVIOUS

?? ?COMPONENT slit_first = Slit(xwidth = slit_w, yheight = 0.00001)
?? ??? ?AT (0, 0, 1E-4) RELATIVE first_arm

?? ?COMPONENT elliptical_1 = Guide_tapering(
?? ??? ?option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
?? ??? ?linw = 0.5, loutw = 0.5,? linh = 0.5, louth = 0.5,
?? ??? ? segno = 100, R0 = 0.995, 
?? ??? ?Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
?? ??? ?alphay = 6.07, W = 0.003, mx = 3, my = 3)
?? ??? ?AT (-0.5+slit_w/2, 0,1E-3) RELATIVE first_arm
?? ??? ?EXTEND
?? ??? ?%{
? ??? ??? ?if(!SCATTERED) ABSORB;
?? ??? ?%}


COMPONENT second_arm = Arm()
?? ?AT (0, 0, arm_lenght) RELATIVE first_arm

?? ?COMPONENT elliptical_2 = Guide_tapering(
?? ??? ?option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
?? ??? ?linw = 0.5, loutw = 0.5,? linh = 0.5, louth = 0.5,
?? ??? ? segno = 10, R0 = 0.995, 
?? ??? ?Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
?? ??? ?alphay = 6.07, W = 0.003, mx = 3, my = 3)
?? ??? ?AT (0.5-slit_w/2, 0,1E-3) RELATIVE second_arm
?? ??? ?EXTEND
?? ??? ?%{
? ??? ??? ?if(!SCATTERED) ABSORB;
?? ??? ?%}


?? ?COMPONENT slit_f = Slit(xwidth = slit_w, yheight = 0.00001)
?? ??? ?AT (0, 0, arm_lenght) RELATIVE second_arm

COMPONENT mon_TOF_end = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "time bins=180, lambda, limits=[29 31], bins=180")
?? ?AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_end = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "lambda, limits = [8.99, 9.01], bins = 180")
?? ?AT (0, 0, 0.001) RELATIVE PREVIOUS


COMPONENT monitor_time_end_9 = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "time, limits = [0.022784, 0.022792] bins = 200")
?? ?AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_time_end_30 = Monitor_nD(
?? ?xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
?? ?options = "time, limits = [0.0455305, 0.0455307] bins = 200")
?? ?AT (0, 0, 0.001) RELATIVE PREVIOUS



FINALLY
%{
%}

END


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180916/3989071d/attachment.html>

From pkwi at fysik.dtu.dk  Mon Sep 17 13:15:08 2018
From: pkwi at fysik.dtu.dk (=?utf-8?B?UGV0ZXIgS2rDpnIgV2lsbGVuZHJ1cA==?=)
Date: Mon, 17 Sep 2018 11:15:08 +0000
Subject: [mcstas-users] elliptic mirror problem
In-Reply-To: <1537127785.775917555@f405.i.mail.ru>
References: <1537127785.775917555@f405.i.mail.ru>
Message-ID: <CC98D58E-BD02-4911-B835-48400212F2E5@fysik.dtu.dk>

Dear Ekatarina,


I have had a quick look at your instrument, which describes the attached geometry (top view) when visualised with mcdisplay.

Is the large horizontal offset between the two guides intentional?

Could you supply some more information about what you are trying to simulate?


Best,

Peter Willendrup

[cid:6BF1FF73-E488-4562-91BE-32509504C243 at win.dtu.dk]
On 16 Sep 2018, at 21.56, ????????? ????????? <ruvinskaia.e.e at mail.ru<mailto:ruvinskaia.e.e at mail.ru>> wrote:

Dear all,

I have a problem. I need to neutron beam pass through elliptical mirrors system. I tried approach through component "Guide_tapering" and component "Mirror_Elliptic". But in both cases neutron beam does not move after the mirror. Could you say what the source of the problem and how I can fix this. Code with "Guide_tapering" for example are attached.

Thanking in advance.
Yours faithfully, Ekaterina.


DEFINE INSTRUMENT NSE_res(double slit_w = 0.1)


DECLARE
%{

double slit_h;
int iter;
double arm_lenght = 5;
%}


INITIALIZE
%{

%}


TRACE


COMPONENT Origin = Progress_bar()
  AT (0,0,0) ABSOLUTE

COMPONENT source =  Source_gen(
    yheight=0.00001,xwidth=slit_w, lambda0=9, dlambda=0.000000000000001,I1=1e13,
    verbose=1,focus_xw=slit_w,focus_yh=0.0001, dist = arm_lenght/4)
    AT (0, 0, 0)  RELATIVE Origin
  EXTEND
  %{
     t = rand0max(1e-20); // set time from 0 to 1 ms for TOF instruments.
  %}


COMPONENT mon_TOF_source = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits=[0.000000001 0.3], bins=180, lambda, limits=[8.9 9.1], bins=180")
    AT (0, 0, 0.00001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_source = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "lambda, limits = [8.99, 9.01], bins = 180")
    AT (0, 0, 0.00001) RELATIVE PREVIOUS


COMPONENT first_arm = Arm()
    AT (0, 0, 0.000001) RELATIVE PREVIOUS

    COMPONENT slit_first = Slit(xwidth = slit_w, yheight = 0.00001)
        AT (0, 0, 1E-4) RELATIVE first_arm

    COMPONENT elliptical_1 = Guide_tapering(
        option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
        linw = 0.5, loutw = 0.5,  linh = 0.5, louth = 0.5,
         segno = 100, R0 = 0.995,
        Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
        alphay = 6.07, W = 0.003, mx = 3, my = 3)
        AT (-0.5+slit_w/2, 0,1E-3) RELATIVE first_arm
        EXTEND
        %{
           if(!SCATTERED) ABSORB;
        %}


COMPONENT second_arm = Arm()
    AT (0, 0, arm_lenght) RELATIVE first_arm

    COMPONENT elliptical_2 = Guide_tapering(
        option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
        linw = 0.5, loutw = 0.5,  linh = 0.5, louth = 0.5,
         segno = 10, R0 = 0.995,
        Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
        alphay = 6.07, W = 0.003, mx = 3, my = 3)
        AT (0.5-slit_w/2, 0,1E-3) RELATIVE second_arm
        EXTEND
        %{
           if(!SCATTERED) ABSORB;
        %}


    COMPONENT slit_f = Slit(xwidth = slit_w, yheight = 0.00001)
        AT (0, 0, arm_lenght) RELATIVE second_arm

COMPONENT mon_TOF_end = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time bins=180, lambda, limits=[29 31], bins=180")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_end = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "lambda, limits = [8.99, 9.01], bins = 180")
    AT (0, 0, 0.001) RELATIVE PREVIOUS


COMPONENT monitor_time_end_9 = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits = [0.022784, 0.022792] bins = 200")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_time_end_30 = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits = [0.0455305, 0.0455307] bins = 200")
    AT (0, 0, 0.001) RELATIVE PREVIOUS



FINALLY
%{
%}

END


_______________________________________________
mcstas-users mailing list
mcstas-users at mcstas.org<mailto:mcstas-users at mcstas.org>
https://mailman2.mcstas.org/mailman/listinfo/mcstas-users


Peter Kj?r Willendrup
Forskningsingeni?r, Speciakonsulent
N?stformand for DTU Fysik LSU

DTU Physics

[cid:b6339554-7f28-441d-9f31-5944f811397c at win.dtu.dk]


Technical University of Denmark


[cid:e108e480-fcab-46e2-9531-b38165079572 at win.dtu.dk]



Department of Physics
Fysikvej
Building 307
DK-2800 Kongens Lyngby
Direct +45 2125 4612
Mobil +45 2125 4612
Fax +45 4593 2399
pkwi at fysik.dtu.dk<mailto:pkwi at fysik.dtu.dk>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2018-09-17 at 13.11.26.png
Type: image/png
Size: 66203 bytes
Desc: Screen Shot 2018-09-17 at 13.11.26.png
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.gif
Type: image/gif
Size: 58 bytes
Desc: image001.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0002.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.gif
Type: image/gif
Size: 1055 bytes
Desc: image002.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0003.gif>

From shapirod at mail.ru  Sun Sep 30 15:19:13 2018
From: shapirod at mail.ru (=?UTF-8?B?0JTQuNC80LAg0KjQsNC/0LjRgNC+?=)
Date: Sun, 30 Sep 2018 16:19:13 +0300
Subject: [mcstas-users] =?utf-8?q?Debye-Waller_factor?=
Message-ID: <1538313553.534361282@f177.i.mail.ru>

Hello everyone

I have a question concerning Debye-Waller factor in the component PowderN. What should I substitude here? 
In manual it's written that DW factor is exp(-2W). Is this W the same as thermal parameter B, which is connected with root-mean square diversion from equilibrium of atoms <u^2>?
How should I take into account dependence of DW on scattering vector q (i.e. on scattering angle)? 

Thanks in advance
Dmitrii Shapiro (PNPI, Gatchina)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180930/7e8a4586/attachment.html>