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

# In[ ]:


from geoscilabs.dcip.DCWidgetResLayer2_5D import ResLayer_app


# In[ ]:


from matplotlib import rcParams
rcParams['font.size'] = 16


# # Effects of a highly resisitive surface layer

# # 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 below a highly resistive surface layer, and investigate what are happening in the earth when static currents are injected. Different from a sphere case, which is a finite target, the resistive layer will also impact the illumination of the target (conductor or resistor).
# By investigating changes in currents, electric fields, potential, and charges upon different geometry, Tx and Rx location, we understand geometric effects of the resistive layer for DCR survey. 

# # Setup

# <img src="https://github.com/geoscixyz/geosci-labs/blob/main/images/em/DC_ResLayer_Setup.png?raw=true"> </img>

# # Question
# 
# - How does the cylinder affect the apparent resistivity without the resistive layer?
# - How does the resistive layer affect the apparent resistivity? Is there a difference if you add or remove the cylinder target?

# ## Plate model

#  - **survey**: Type of survey
#  - **A**: Electrode A (+) location
#  - **B**: Electrode B (-) location
#  - **M**: Electrode A (+) location
#  - **N**: Electrode B (-) location
#  - **$dz_{layer}$**: thickness of the resistive layer
#  - **$zc_{ayer}$**: z location of the resistive layer
#  - **xc**: x location of cylinder center
#  - **zc**: z location of cylinder center
#  - **$\rho_{1}$**: Resistivity of the half-space
#  - **$\rho_{2}$**: Resistivity of the layer
#  - **$\rho_{3}$**: Resistivity of the cylinder
#  - **Field**: Field to visualize
#  - **Type**: which part of the field
#  - **Scale**: Linear or Log Scale visualization
#  
# ###  **Do not forget to hit Run Interact to update the figure after you made modifications**

# In[ ]:


ResLayer_app()


# In[ ]: