// Jupiter Auroral Simulation
// GEANT4-10.6.2
// Gabriel Bridges Aug, 2020

// B-field definition


//includes
// user headers
#include "JovianField.hh"
#include "global.h"
// GEANT headers
#include "G4SystemOfUnits.hh"
#include "globals.hh"

JovianField::JovianField()
: G4MagneticField()
{
  //jovianMessenger = new JupiterMessenger(this);
}

JovianField::~JovianField(){
  //delete jovianMessenger;
}

void JovianField::GetFieldValue(const G4double point[4], double *bField) const
{
  bField[0] = 0.;
  bField[1] = 0.;
  bField[2] = 6.1e5*nT -(6.57e3*nT/km)*point[2];
}