Commit ddf0908f authored by Caneses Marin, Juan Francisco's avatar Caneses Marin, Juan Francisco
Browse files

Added finite spread on pitch angle

parent e64a32a3

src/MoveParticlePack.f90

+4 −4
Original line number Diff line number Diff line
@@ -374,11 +374,11 @@ SUBROUTINE loadParticles(zp,kep,xip,in0) 
      xip = uparArray/uArray

  elseif (in0%xip_InitType .EQ. 2) then

      ! Beam pitch angle

      xip = in0%xip_init

      !xip = in0%xip_init

      ! Beam pitch angle with a Gaussian distribution in pitch angle space:

      !call random_number(RmArray1)

      !call random_number(RmArray2)

      !xip = in0%xip_init_std*sqrt(-2.*log(RmArray1))*cos(2.*pi*RmArray2)  +  in0%xip_init

      call random_number(RmArray1)

      call random_number(RmArray2)

      xip = in0%xip_init_std*sqrt(-2.*log(RmArray1))*cos(2.*pi*RmArray2)  +  in0%xip_init

  end if

return

END SUBROUTINE loadParticles