Re: [fluka-discuss]: magnetic field
Ok, I was bitten by the 8 character region name feature!
What doesn't make sense is that when printing region numbers out for
"rRMVVac" and "rRMHVac", I saw the correct region number for "rRMHVac"
region but -1 for "rRMVVac". What doesn't make sense is that both
region names are 7 characters and I was not putting a space after the
"c" in the region name for either but was seeing the correct region
number for "rRMHVac"!
Mina
=============================
* Set dipole filed in tesla (B field is 250 Gauss)
DIPFLD = 0.0250D+00
* Some useful printout
WRITE (LUNOUT,*) ''
WRITE (LUNOUT,*) ' magfld.f called for EHDT Raster Magnets'
WRITE (LUNOUT,*) ' Dipole field: ',DIPFLD,' T'
WRITE (LUNOUT,*) ''
* Retrieve the indices of the magnetic regions starting from their
* names, and store them
CALL GEON2R ("rRMVVac ",NRMDV, IERR)
CALL GEON2R ("rRMHVac ",NRMDH, IERR)
WRITE (LUNOUT,*) ''
WRITE (LUNOUT,*) ' "rRMVVac" region number is: ', NRMDV
WRITE (LUNOUT,*) ' "rRMHVac" region number is: ', NRMDH
WRITE (LUNOUT,*) ''
LFIRST=.FALSE.
END IF
On 14-02-12 02:25 PM, Mina Nozar wrote:
> Hello everyone,
>
>
> How does one know whether a magnetic field has actually been applied
> to a region? I can't see anything in the output files that conveys
> this message.
>
> In my first attempt to set the magnetic field for a dipole magnet, I
> modified magfld.f, assigned the appropriate B field to a dipole
> sensitive region (vacuum) and in the ASSIGNMA card, and set magnetic
> field on. I compiled magfld and built the executable linking to
> magfld.o. The magnetic field is 250 G, so I put in 0.025 T for the
> intensity in magfld.f. The results didn't show any deflections in the
> backscattered e+/e- in the beam pipe from a 75 MeV beam dump. Most of
> these e+/e- have energies less than 10 MeV (with the peak around 1-2
> MeV).
>
> Reading through the manual, I realized I forgot to include the
> MGNFIELD card in the input file. So I included this card, recompiled
> and ran again. And I can't see any difference.
>
> MGNFIELD 0.0 0.0 0.0
>
>
> Thanks,
> Mina
>
Received on Thu Feb 13 2014 - 04:08:50 CET
This archive was generated by hypermail 2.3.0
: Thu Feb 13 2014 - 04:08:51 CET