#!/usr/bin/env python
# coding: utf-8

# In[1]:


from em_examples.DC_cylinder import *
from IPython.display import display
get_ipython().run_line_magic('matplotlib', 'inline')


# # Purpose 
# 
# For a direct current resistivity (DCR) survey, currents are injected to the earth, and flow. 
# Depending upon the conductivity contrast current flow in the earth will be distorted, and these changes 
# can be measurable on the sufurface electrodes. 
# Here, we focus on a cylinder target embedded in a halfspace, and investigate what are happening in the earth when static currents are injected. Different from a sphere case, which is a finite target, "coupling" among Tx, target (conductor or resistor), and Rx will be significanlty different upon various scenarios. 
# By investigating changes in currents, electric fields, potential, and charges upon different geometry of cylinder and survey, Tx and Rx location, we understand geometric effects of the target for DCR survey. 

# # Setup

# <img src="https://github.com/geoscixyz/em_apps/blob/master/images/DCR_Setup Cylinder.png?raw=true"></img>

# # Question

# - Is the potential difference measured by a dipole over a conductive (/resisitive) target higher or lower compared to the half-space reference?
# - how do the field lines bend in presence of a conductive (/resistive) target?
# - Compared to the positive and negative sources (A and B), how are oriented the positive and negative accumulated charges around a conductive (/resistive) target?
# - How would you describe the secondary fields pattern? Does it remind you one fundamental object of electromagnetics?

# # Cylinder app
# 
# ## Parameters:

#  - **survey**: Type of survey
#  - **A**: (+) Current electrode  location
#  - **B**: (-) Current electrode  location
#  - **M**: (+) Potential electrode  location
#  - **N**: (-) Potential electrode  location
#  - **r**: radius of cylinder
#  - **xc**: x location of cylinder center
#  - **zc**: z location of cylinder center
#  - **$\rho_1$**: Resistivity of the halfspace
#  - **$\rho_2$**: Resistivity of the cylinder
#  - **Field**: Field to visualize
#  - **Type**: which part of the field
#  - **Scale**: Linear or Log Scale visualization

# In[2]:


app = cylinder_app()
display(app)


# In[ ]: