add_explicit_terms Subroutine

public subroutine add_explicit_terms(g1, g2, g3, phi, apar, bpar, istep, bd)

FIXME : Add documentation
! NEED TO DO SOMETHING HERE... BD GGH

Arguments

Type IntentOptional AttributesName
complex, intent(inout), dimension (-ntgrid:,:,g_lo%llim_proc:):: g1
complex, intent(inout), dimension (-ntgrid:,:,g_lo%llim_proc:):: g2
complex, intent(inout), dimension (-ntgrid:,:,g_lo%llim_proc:):: g3
complex, intent(in), dimension (-ntgrid:,:,:):: phi
complex, intent(in), dimension (-ntgrid:,:,:):: apar
complex, intent(in), dimension (-ntgrid:,:,:):: bpar
integer, intent(in) :: istep
real, intent(in) :: bd

Contents

Source Code


Source Code

  subroutine add_explicit_terms (g1, g2, g3, phi, apar, bpar, istep, bd)
    use theta_grid, only: ntgrid
    use gs2_layouts, only: g_lo
    use gs2_time, only: save_dt_cfl
    use job_manage, only: time_message
    implicit none
    complex, dimension (-ntgrid:,:,g_lo%llim_proc:), intent (in out) :: g1, g2, g3
    complex, dimension (-ntgrid:,:,:), intent (in) :: phi,    apar,    bpar
    integer, intent (in) :: istep
    real, intent (in) :: bd
    real :: dt_cfl
    logical, parameter :: nl = .true.
    call time_message(.false., time_add_explicit_terms, 'Explicit terms')
    select case (nonlinear_mode_switch)
    case (nonlinear_mode_none)
!!! NEED TO DO SOMETHING HERE...  BD GGH
       dt_cfl = 1.e8
       call save_dt_cfl (dt_cfl)
#ifdef LOWFLOW
       if (istep /=0) &
            call add_explicit (g1, g2, g3, phi, apar, bpar, istep, bd)
#endif
    case (nonlinear_mode_on)
       if (istep /= 0) then
         call add_explicit (g1, g2, g3, phi, apar, bpar, istep, bd, nl)
       endif
    end select
    call time_message(.false., time_add_explicit_terms, 'Explicit terms')
  end subroutine add_explicit_terms