Constraints on melts and fluids in the continental crust from electrical conductivity anisotropy simulations

Codes:

Random (3 phases): The code can be run once to generate two text files (componentsx and componentsy), i.e. a resistor networks in both directions. Importing these files into the circuit simulator allows a calculation to be performed.

line 4: the first number represents the network size (n); the second number represents the resistance values of second mineral phases (Plagioclase or Clinopyroxene); the third number represents the number of second mineral phases in the network.

line 5: the first number represents the resistance values third mineral phases (Amphibole or Quartz); the second number represents the number of third mineral phases in the network.

lines 7 to 28: the distribution functions of electrical anisotropy for first minerals (Qtz or Plag).

Layering (Middle crust)-three phases: the code generates compenents20x and components20y.

line 4: the first number represents the network size (n); the second number represents the resistance values of second mineral phases (Plagioclase); the third number is none; the 4-7th numbers represent the number of second mineral phases in different layers.

line 5: the first number represents the resistance values of third mineral phases (Amphibole); the 2-5th numbers represent the number of third mineral phases in different layers.

lines 7 to 33: the distribution functions of electrical anisotropy for Quartz.

Layering (Lower crust)-three phases: the code generates a resistor network consisting of 10 layers.

line 4: the parameters for second mineral phases (Clinopyroxene).

line 5: the parameters for third mineral phases (Quartz).

lines 7 to 23: the distribution functions of electrical anisotropy for Plagioclase.

Fluid/melt-bearing models is similar to layering calculations.