checked drifts and mirror. Individually working

init_g.f90

@@ -255,8 +255,6 @@ subroutine ginit_default
             end do
          end if
       end if
-
-      phi(1,1,:) = 0.0 
 
       ! need better way to initialise for full flux surface cases
       ia = 1

mirror_terms.f90

@@ -214,21 +214,21 @@ subroutine init_invert_mirror_operator
 
       !> corresponds to sign of mirror term positive on RHS of equation
       a(2:, 1) = -0.5 * (1.0 - 2.0 * vpa_upwind) / dvpa
-      b(1:, 1) = -2.0 * vpa_upwind / dvpa
-      c(1:nvpa - 1, 1) = 0.5 * (1.0 + 2.0 * vpa_upwind) / dvpa
+      b(2:, 1) = -2.0 * vpa_upwind / dvpa
+      c(2:nvpa - 1, 1) = 0.5 * (1.0 + 2.0 * vpa_upwind) / dvpa
       !> must treat boundary carefully
       !> assumes fully upwinded at outgoing boundary
-      !b(1, 1) = -1.0 / dvpa
-      !c(1, 1) = 1.0 / dvpa
+      b(1, 1) = -1.0 / dvpa
+      c(1, 1) = 1.0 / dvpa
 
       !> corresponds to sign of mirror term negative on RHS of equation
       a(2:nvpa, -1) = -0.5 * (1.0 + 2.0 * vpa_upwind) / dvpa
       b(:nvpa, -1) = 2.0 * vpa_upwind / dvpa
       c(:nvpa - 1, -1) = 0.5 * (1.0 - 2.0 * vpa_upwind) / dvpa
       !> must treat boundary carefully
       !> assumes fully upwinded at outgoing boundary
-      !a(nvpa, -1) = -1.0 / dvpa
-      !b(nvpa, -1) = 1.0 / dvpa
+      a(nvpa, -1) = -1.0 / dvpa
+      b(nvpa, -1) = 1.0 / dvpa
 
       !> time_upwind = 0.0 corresponds to centered in time
       !> time_upwind = 1.0 corresponds to fully implicit (upwinded)

time_advance.f90

@@ -1351,6 +1351,7 @@ subroutine advance_wstar_explicit(phi, gout)
 
          !> calculate d<phi>/dy in k-space
          call get_dgdy(g_gyro, g0)
+
          !> transform d<phi>/dy from ky-space to y-space
          do ivmu = vmu_lo%llim_proc, vmu_lo%ulim_proc
             do iz = -nzgrid, nzgrid