[CIG-LONG] Segmentation Fault

[Walter Landry]

Hi Johnny,

Cc'ing the list.

Johnny Seales <johnny.seales at gmail.com> wrote:
> Walter,
> 
> I'm currenlty building a layered plateau model and cannot get my code to
> run except for very small nx and ny (2 is the max for each). I have checked
> the input file over and over to see if something is wrong. The input file
> is attached and the message I receive is below. Thanks.

The shapes of the crust and mantle do not cover the entire domain.
The easiest way to ensure that is to make the crust, for example,
cover everything that the mantle does not cover.  I am attaching a
fixed input file.

Cheers,
Walter Landry
-------------- next part --------------
{
    "EulerDeform":
    {
        "systems": [
            {
                "mesh": "v-mesh",
                "p-mesh": "p-mesh",
                "remesher": "velocityRemesher",
                "velocityField": "VelocityField",
                "wrapTop": "True"
            }
        ]
    },
    "components":
    {
        "buoyancyForceTerm":
        {
            "Type": "BuoyancyForceTerm",
            "ForceVector": "mom_force",
            "Swarm": "gaussSwarm",
            "gravity": "gravity"
        },
        "backgroundShape":
        {
            "Type": "EquationShape",
            "equation": "1"
        },
	// "crustLithosphereShape":
	// {
	//     "Type": "PolygonShape",
	//     "vertices": [
	// 	"asciidata",
	// 	["x", "y"],
	// 	0.0, 25.0,
        //         300.0, 25.0,
        //         300.0, 19.0,
        //         800.0, 19.0,
        //         800.0, 25.0,
        //         1100.0, 25.0,
	// 	1100.0, 68.0,
        //         800.0, 68.0,
        //         800.0, 70.0,
        //         300.0, 70.0,
        //         300.0, 68.0,
        //         0.0, 68.0
	//     ]
	// },
        "crustLithosphereShape":
        {
            "Type": "Intersection",
            "shapes":
            [
                "backgroundShape",
                "!mantleLithosphereShape"
            ]
        },
        "mantleLithosphereShape":
        {
            "Type": "PolygonShape",
            "vertices": [
                "asciidata",
                ["x", "y"],
                0.0, 0.0,
                1100.0, 0.0,
                1100.0, 25.0,
                800.0, 25.0,
                800.0, 19.0,
                300.0, 19.0, 
                300.0, 25.0, 
                0.0, 25.0
            ]
        },
        "crustLithosphereViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "1.0e4"
        },
        "crustLithosphereStrainWeakening":
        {
            "Type": "StrainWeakening",
            "TimeIntegrator": "timeIntegrator",
            "MaterialPointsSwarm": "materialSwarm",
            "softeningStrain": "0.1",
            "initialDamageFraction": "0.0",
            "initialDamageWavenumber": "0.5",
            "initialDamageFactor": "0.5",
            "healingRate": "0.0"
        },
        "crustLithosphereYielding":
        {
            "Type": "DruckerPrager",
            "PressureField": "PressureField",
            "VelocityGradientsField": "VelocityGradientsField",
            "MaterialPointsSwarm": "materialSwarm",
            "Context": "context",
            "StrainWeakening": "strainWeakening",
            "StrainRateField": "StrainRateField",
            "cohesion": "1.0",
            "cohesionAfterSoftening": "0.0001",
            "frictionCoefficient": "0.0",
            "frictionCoefficientAfterSoftening": "0.0",
            "minimumViscosity": "1.0e-4"
        },
        "crustLithosphereViscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "crustLithosphereShape",
            "density": "1.0",
            "Rheology": [
                "crustLithosphereViscosity",
                "crustLithosphereYielding",
                "storeViscosity",
                "storeStress"
            ]
        },        
        "mantleLithosphereViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "1.0"
        },
        "mantleLithosphereViscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "mantleLithosphereShape",
            "density": "1.222",
            "Rheology": [
                "mantleLithosphereViscosity",
                "storeViscosity",
                "storeStress"
            ]
        }       
    },
    "velocityBCs" :{
        "type": "CompositeVC",
        "vcList": [
            {
                "type": "WallVC",
                "wall": "left",
                "variables": [
                    {
                        "name": "vx",
                        "value": "0"
                    }
                ]
            },
            {
                "type": "WallVC",
                "wall": "right",
                "variables": [
                    {
                        "name": "vx",
                        "value": "0"
                    }
                ]
            },
            {
                "type": "WallVC",
                "wall": "bottom",
                "variables": [
                    {
                        "name": "vy",
                        "value": "0"
                    }
                ]
            }
        ]
    },

    "FieldVariablesToCheckpoint": [
        "StrainRateInvariantField",
        "VelocityField",
        "PressureField"
    ],
    "maxTimeSteps": "10",
    "outputPath": "./CPC/output",
    "dim": "2",
    "minX": "0",
    "minY": "0",
    "minZ": "0",
    "maxX": "1100",
    "maxY": "70",
    "maxZ": "0.3",
    "nx": "16",
    "ny": "4",
    "nz": "4",
    "particlesPerCell": "40",
    "seed": "13",
    "checkpointEvery": "1",
    "gravity": "1.0" 
}