The abnormal storm fronts and epic lightning storms continue across the Earth but the eye of this activity has been the area surrounding the town of Lakeview. A few months have passed and your colleagues at HQ have designed and constructed a new age measurement and prediction device called the Mapping N$\vec{E}$twork Sensory Array (MNSA) to try and understand the source of these storms. This sophisticated measuring equipment works off of a large number of golf ball sized sensors buried in the ground in the surrounding area by HQ. These sensors send signals back to the main relay, located atop Stormchaser HQ. The MNSA was designed for predicting lightning strikes in the surrounding area as well as calibrating instrumentation used by field teams. This is a process you and your team are very familiar with, so your colleagues showed no hesitation calling you and asking for your inputs.

Using your existing knowledge of electric field computer modeling, your Lakeview colleagues have commissioned your team to help with the maiden voyage of the MNSA. In the sky on this strangely calm but partly cloudy day, are three clouds: a large cloud (Q = -150C), a small cloud (q = -50C), and an even smaller cloud (qq = 50 C). Your colleagues have insisted that this is the perfect day to launch their new network because the weather is calm. Some of their programmers have started to work on the model, but due to their inexperience with this type of work, they were unable to finish. So far, this is what they have:

##Scene Setup
scene = display(width=1000, length=1000, height = 500)

## Parameters and Initial Conditions
k = 9e9
Q = -150 ##C
q = -50 ##C
qq = 50 ##C
cloud_height = 600 ##m

## Objects
ground = box(pos = vec(0,0,0), width=5000, length=5000, height=0.1, color=color.green)
mountains = [cone(pos=vec(-1800,0,0), axis=vec(0,1000,0), radius=700, color=vec(1,0.7,0.2)),
            cone(pos=vec(-1800,0,-1600), axis=vec(0,1000,0), radius=700, color=vec(1,0.7,0.2)), 
            cone(pos=vec(-1800,0,1600), axis=vec(0,1000,0), radius=700, color=vec(1,0.7,0.2))]
HQ = box(pos = vec(-100,100,0), width = 15, length = 15, height= 200, color = color.blue)
large_cloud = sphere(pos = vec(-1000,cloud_height,0), radius = 100, color = color.white)
small_cloud = sphere(pos = vec(900,cloud_height,100), radius = 40, color = color.white)
smaller_cloud = sphere(pos = vec(800, cloud_height, 800), radius = 25, color = color.white)

##Calculation
E = vec(0,0,0)
point = vec(HQ.pos.x,HQ.pos.y,0)
field = arrow(pos = point, axis = E, color=color.yellow)

##Calculation Loop
x = -200
dx = 20
xmax = 200
while x<= xmax:
    z = -200
    dz = 20
    zmax = 200
    point = vec(x,0,z)
    E = vec(0,0,0)
    field = arrow(pos = point, axis = E, color=color.yellow)
    x = x + dx

Before you can start calibrating the field sensors, the main relay needs to be adjusted to the electric fields caused by the nearby clouds. The main sensor relay is located at the top of the HQ. Once the main relay is programmed properly, you and your team will need to complete this program to map out the electric field on the surface of the ground.

## Creating the scene for the code to run in
scene.range = 2

## Constants
TotalCharge = 15 #C
pointcharge = TotalCharge/7  
k = 9e9  
vscale = 1e-4

## Objects
charge1 = sphere(pos=vec(-3,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge2 = sphere(pos=vec(-2,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge3 = sphere(pos=vec(-1,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge4 = sphere(pos=vec(0,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge5 = sphere(pos=vec(1,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge6 = sphere(pos=vec(2,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charge7 = sphere(pos=vec(3,0,0), Q=pointcharge,  color=color.red,  size=5e-1*vec(1,1,1))
charges = [charge1, charge2, charge3, charge4, charge5, charge6, charge7]

## Calculation Loop 1
E = vec(0,0,0)
point = vec(0,0,0)
for c in charges:
    r = point - c.pos
field = arrow(pos=point, axis=vscale*E, color = color.cyan)

## Calculation Loop 2
x = -5
dx = 0.5
xmax = 5
while x<=xmax:
    theta = 0
    dtheta = 0.1
    R = 0
    while theta < 2*pi:
        E = vec(0,0,0)
        point = vec(x, R*sin(theta), R*cos(theta))
        field = arrow(pos=point, axis = E*vscale, color = color.green)
        theta += dtheta
    x+=dx

Complete the program above to first represent the total electric field just to the right and left of the line of charges. Then, calculate the total electric field surrounding the line of charges.

*/