Constraints on melt and fluid in the continental crust from electrical conductivity anisotropy

Codes: Random: 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 secondy number represents the resisitance vaules of secondary mineral phases; the third number represents the number of secondary mineral phases in the network. lines 6 to 22 : the distribution functions of electrical anisotropy for first minerals (Qtz or Pl). Layering (Middle crust): the code generates compenents20x and components20y. line 4 : the first number represents the network size (n); the secondy number represents the resisitance vaules of secondary mineral phases; the third number is none; the 4-7th numbers represent the numbers secondary mineral phases in different layers. Layering (Lower crust): the code generates a resistor network consisting of 10 layers. Fluid/melt-bearing models is similar to layering calculations.