advance_field_adapt.F90
a58f14c587be3D_parallel

advance_field_adapt.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 grid
			use prec_const
			IMPLICIT NONE

			COMPLEX(dp), DIMENSION ( ikxs:ikxe, ikys:ikye, ntimelevel ) :: f
			COMPLEX(dp), DIMENSION ( ikxs:ikxe, ikys:ikye, ntimelevel ) :: f_rhs
			REAL(dp) :: error
			INTEGER :: istage

			SELECT CASE (updatetlevel)

			CASE(1)
			SELECT CASE (numerical_scheme)

			CASE ('DOPRI5_ADAPT')
			error = 0._dp
			DO iky=ikys,ikye
			DO ikx=ikxs,ikxe
			fs = f(ikx,iky,1)
			DO istage=1,ntimelevel
			f(ikx,iky,1) = f(ikx,iky,1) + dtb_E(istage)f_rhs(ikx,iky,istage)
			fs = fs + dtb_Es(istage)f_rhs(ikx,iky,istage)
			END DO
			IF ( ABS(f(ikx,iky,1) - fs) .GT. error ) THEN
			error = ABS(f(ikx,iky,1) - fs)
			ENDIF
			END DO
			END DO
			IF (error > TOL)
			CASE DEFAULT
			DO iky=ikys,ikye
			DO ikx=ikxs,ikxe
			fs = f(ikx,iky,1)
			DO istage=1,ntimelevel
			f(ikx,iky,1) = f(ikx,iky,1) + dtb_E(istage)f_rhs(ikx,iky,istage)
			END DO
			END DO
			END DO
			END SELECT

			CASE DEFAULT
			DO iky=ikys,ikye
			DO ikx=ikxs,ikxe
			f(ikx,iky,updatetlevel) = f(ikx,iky,1);
			DO istage=1,updatetlevel-1
			f(ikx,iky,updatetlevel) = f(ikx,iky,updatetlevel) + &
			dtA_E(updatetlevel,istage)f_rhs(ikx,iky,istage)
			END DO
			END DO
			END DO
			END SELECT

			END SUBROUTINE advance_field

			END MODULE advance_field_routine