Thank you for providing this really useful program!

I have one question about the example code you've offered (split_model.m).
According to the code of generating elasticities for lower layer (L2), it contains 2-L2_faln term in MS_VRH function, which is different with that of upper layer (L1) 1-L1_faln.

However, since L1_faln and L2_faln are both defined as volume fractions (up to my understanding) given as 0.3 in examples, I was confused why 2-L2_faln was used for lower layer instead of 1-L2_faln.

I attach below the part of the code related to my inquiry:

`%  ** generate layer elasticities:
%        This is a Voigt-Reuss-Hill average of the appropriately rotated olivine
%        tensor and its isotropic equivalent.
      [L1_C,~] = MS_VRH([L1_faln 1-L1_faln],...
         MS_rot3(Cani,0,-L1_dip,L1_aaz,'order',[3 2 1]),rh, Ciso, rh) ;

      [L2_C,~] = MS_VRH([L2_faln 2-L2_faln],...
         MS_rot3(Cani,0,-L2_dip,L2_aaz,'order',[3 2 1]),rh, Ciso, rh) ;`