memory.F90
56b9dafb082cdivided_output

memory.F90
View Options

			SUBROUTINE memory
			! Allocate arrays (done dynamically otherwise size is unknown)

			USE array
			USE basic
			USE fields
			USE grid
			USE time_integration
			USE model, ONLY: CO, NON_LIN, A0KH

			USE prec_const
			IMPLICIT NONE

			! Moments and moments rhs
			CALL allocate_array( moments_e, ips_e,ipe_e, ijs_e,ije_e, ikrs,ikre, ikzs,ikze, 1,ntimelevel )
			CALL allocate_array( moments_i, ips_i,ipe_i, ijs_i,ije_i, ikrs,ikre, ikzs,ikze, 1,ntimelevel )
			CALL allocate_array( moments_rhs_e, ips_e,ipe_e, ijs_e,ije_e, ikrs,ikre, ikzs,ikze, 1,ntimelevel )
			CALL allocate_array( moments_rhs_i, ips_i,ipe_i, ijs_i,ije_i, ikrs,ikre, ikzs,ikze, 1,ntimelevel )

			! Electrostatic potential
			CALL allocate_array(phi, ikrs,ikre, ikzs,ikze)

			! Collision matrix
			IF (CO .EQ. -1) THEN
			CALL allocate_array( Ceepj, 1,(pmaxe+1)(jmaxe+1), 1,(pmaxe+1)(jmaxe+1))
			CALL allocate_array( CeipjT, 1,(pmaxe+1)(jmaxe+1), 1,(pmaxe+1)(jmaxe+1))
			CALL allocate_array( CeipjF, 1,(pmaxe+1)(jmaxe+1), 1,(pmaxi+1)(jmaxi+1))

			CALL allocate_array( Ciipj, 1,(pmaxi+1)(jmaxi+1), 1,(pmaxi+1)(jmaxi+1))
			CALL allocate_array( CiepjT, 1,(pmaxi+1)(jmaxi+1), 1,(pmaxi+1)(jmaxi+1))
			CALL allocate_array( CiepjF, 1,(pmaxi+1)(jmaxi+1), 1,(pmaxe+1)(jmaxe+1))
			ENDIF

			! Non linear terms and dnjs table
			IF ( NON_LIN .OR. (A0KH .NE. 0) ) THEN
			CALL allocate_array( Sepj, ips_e,ipe_e, ijs_e,ije_e, ikrs,ikre, ikzs,ikze )
			CALL allocate_array( Sipj, ips_i,ipe_i, ijs_i,ije_i, ikrs,ikre, ikzs,ikze )
			CALL allocate_array( dnjs, 1,maxj+1, 1,maxj+1, 1,maxj+1)
			ENDIF

			END SUBROUTINE memory