anchor.diff:  				write(WRITE_UNIT_1,'(F24.15,A)') totalMassGHills(j), ' ; ...'
   			enddo
   			write(WRITE_UNIT_1,'(F24.15,A)') totalMassGHills(timesteps), ' ] ;'
  -
  +			
  +			write(WRITE_UNIT_1,'(A,F24.15,A)') 'tracerVar = [ ', tracerVar(0), ' ; ...'
  +			do j = 1, timesteps-1
  +				write(WRITE_UNIT_1,'(F24.15,A)') tracerVar(j), ' ; ...'
  +			enddo
  +			write(WRITE_UNIT_1,'(F24.15,A)') tracerVar(timesteps), ' ] ;'
  +			
   		endif
   		close(WRITE_UNIT_1)
   
  
positive.code: 	call AdvectionRK4Timestep(timekeeper, sphere, dt, t, procRank, numProcs, LauritzenEtAlNonDivergentWind)
  
  	totalMassGHills(timeJ+1) = ( TotalMass(sphere, tracerID) - mass0 ) / mass0
  
  	t = real( timeJ+1, kreal) * dt
  
  
negative.code: !
  ! test case variables
  !
  real(kreal), allocatable :: totalMassGHills(:), tracerVar(:)
  real(kreal) :: sphereL2, sphereLinf, panelsLinf, particlesLinf, phiMax, phiMin, deltaPhi, phimax0, phimin0
  real(kreal) :: mass0, var0
  
  !
  ! logging
  