Underworld modeling of large-scale folds in ice sheets

Code for the revising manuscript, Zhang et al., Formation Mechanisms of Large-scale Folding in Greenland’s Ice Sheet. Note in the code: the z-axis is the horizontal inflow direction; the x-axis is the horizontal outflow direction; and the y-axis is the vertical direction.

Model 1: isotropic ice, free-slip basal boundary Model 2: anisotropic ice, free slip Model 3: anisotropy, flat bedrock surface Model 4: anisotropy, bumpy bedrock surface (max 400m) Model 5: isotropy, bumpy bedrock surface (max 400m) Model 6: anisotropy, bumpy bedrock surface (max 400m), softness in ice viscosities Model 7: isotropy, bumpy bedrock surface (max 400m), softness in ice viscosities model 8: larger along-flow model size based on Model 4 Model 9: smaller bumps based on Model 4 Model 10: larger bumps based on Model 4 Model 11: larger anisotropy based on Model 4 Model 12: with ice precipitation based on Model 11 Model 13: pre-existed folds (300 m), flat bedrock surface, anisotropy Model 14: pre-existed folds (300 m), flat bedrock surface, anisotropy Model 15: pre-existed folds (300 m), flat bedrock surface, isotropy