inital.F90
815393243914divided_output

inital.F90
View Options

			SUBROUTINE inital

			use basic

			use prec_const
			implicit none

			WRITE(*,'(a/)') '=== Set initial conditions ==='

			CALL init_profiles

			END SUBROUTINE inital


			SUBROUTINE init_profiles
			! Set initial conditions

			USE basic
			USE fourier_grid
			USE fields
			use initial_par
			USE time_integration

			use prec_const
			IMPLICIT NONE

			INTEGER :: ip,ij, ikr,ikz
			INTEGER, DIMENSION(12) :: iseedarr
			REAL(dp) :: noise

			! Seed random number generator
			iseedarr(:)=iseed
			CALL RANDOM_SEED(PUT=iseedarr)

			CALL set_updatetlevel(1)

			DO ikr=ikrs,ikre
			DO ikz=ikzs,ikze

			DO ip=ips_e,ipe_e
			DO ij=ijs_e,ije_e
			CALL RANDOM_NUMBER(noise)
			moments_e( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp) ! Re
			CALL RANDOM_NUMBER(noise)
			moments_e( ip,ij, ikr,ikz, :) = imagu * (initback_moments + initnoise_moments*(noise-0.5_dp)) ! Im
			END DO
			END DO

			DO ip=ips_i,ipe_i
			DO ij=ijs_i,ije_i
			CALL RANDOM_NUMBER(noise)
			moments_i( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp) ! Re
			CALL RANDOM_NUMBER(noise)
			moments_i( ip,ij, ikr,ikz, :) = imagu * (initback_moments + initnoise_moments*(noise-0.5_dp)) ! Im
			END DO
			END DO

			END DO
			END DO

			CALL poisson ! To set phi

			END SUBROUTINE init_profiles