anchor.diff:  public InitSolidBodyRotation, SetSolidBodyRotationOnMesh
   public SOLID_BODY_NINT, SOLID_BODY_NREAL
   public NullVorticity
  +public InitSingleGaussianVortex, SetSingleGaussianVortexOnMesh
  +public SINGLE_VORTEX_NINT, SINGLE_VORTEX_NREAL
   
   !
   !----------------
  
positive.code: 
  end subroutine
  
  subroutine NullVorticity(aMesh,nullVort)
  	type(SphereMesh), intent(inout) :: aMesh
  	type(BVESetup), intent(in) :: nullVort
  
negative.code: 	SolidBodyX = 2.0_kreal*rotationRate*xyz(3)/EARTH_RADIUS
  end function
  
  function GeneralGaussian(xyz,xyzCent, bb, maxVal)
  	real(kreal) :: GeneralGaussian
  	real(kreal), intent(in) :: xyz(3), xyzCent(3), bb, maxVal
  	GeneralGaussian = maxVal * exp( -2.0_kreal*bb*bb*( EARTH_RADIUS*EARTH_RADIUS - sum(xyz*xyzCent)))
  end function
  
  subroutine InitLogger(aLog,rank)
  	type(Logger), intent(inout) :: aLog
  