advance_field.F90
9a072e2a6205HeLaZ2

advance_field.F90
View Options

			MODULE advance_field_routine

			USE prec_const
			implicit none

			CONTAINS

			SUBROUTINE advance_time_level

			USE basic
			USE time_integration
			use prec_const
			IMPLICIT NONE

			call set_updatetlevel(mod(updatetlevel,ntimelevel)+1)

			END SUBROUTINE advance_time_level



			SUBROUTINE advance_field( f, f_rhs )

			USE basic
			USE time_integration
			USE array
			USE space_grid
			use prec_const
			IMPLICIT NONE

			real(dp), DIMENSION ( ikrs:ikre,ikzs:ikze,ntimelevel ) :: f
			real(dp), DIMENSION ( ikrs:ikre,ikzs:ikze,ntimelevel ) :: f_rhs
			INTEGER :: ikr, ikz, istage

			SELECT CASE (updatetlevel)
			CASE(1)
			DO ikz=ikzs,ikze
			DO ikr=ikrs,ikre
			DO istage=1,ntimelevel
			f(ikr,ikz,1) = f(ikr,ikz,1) + dtb_E(istage)f_rhs(ikr,ikz,istage)
			END DO
			END DO
			END DO
			CASE DEFAULT
			DO ikz=ikzs,ikze
			f(ikr,ikz,updatetlevel) = f(ikz,1);
			DO ikr=ikrs,ikre
			DO istage=1,updatetlevel-1
			f(ikr,ikz,updatetlevel) = f(ikr,ikz,updatetlevel) + &
			dtA_E(updatetlevel,istage)f_rhs(ikr,ikz,istage)
			END DO
			END DO
			END DO
			END SELECT
			END SUBROUTINE advance_field

			END MODULE advance_field_routine